在线等：怎么把输出的x和y的一一对应的数保存为txt呢？

xiezhh

回复 14 # ChaChing 的帖子

ChaChing兄说的没错，dlmwrite可能需要用户设置精度。

mingmingtree

回复 11 # xiezhh 的帖子

1. clear all；
   
2. clc;
   
3. a = 5.8;
   
4. b1 = 1528;b2=1535;b3=1560;b4=1583;
   
5. c1 = 0.3;c2=0.2;c3=0.25;c4=0.3;
   
6. g1 = @(t)a*exp(-4*log(2)*(((t-b1)/c1).^2));
   
7. g2 = @(t)a*exp(-4*log(2)*(((t-b2)/c2).^2));
   
8. g3 = @(t)a*exp(-4*log(2)*(((t-b3)/c3).^2));
   
9. g4 = @(t)a*exp(-4*log(2)*(((t-b4)/c4).^2));
   
10. x = 1510:0.005:1592;
    
11. y1 = g1(x);
    
12. y2 = g2(x);
    
13. y3 = g3(x);
    
14. y4 = g4(x);
    
15. s1=y1+y2+y3+y4;
    
16. size(s1);%纯信号
    
17. z=awgn(s1,-5,'measured');%加入噪声
    
18. s=z(1:16002);%去前16002个点作为程序的原始含噪信号
    
19. ls=length(s);
    
20. subplot(2,2,1) ,plot(x,z);title('含噪声的信号');grid;
    
21. 
    
22. %用db1小波对原始信号进行3层分解并提取系数
    
23. [c,l]=wavedec(s,3,'db1');
    
24. ca3=appcoef(c,l,'db1',3);
    
25. cd3=detcoef(c,l,3);
    
26. cd2=detcoef(c,l,2);
    
27. cd1=detcoef(c,l,1);
    
28. 
    
29. %用三种阈值分别进行去噪处理
    
30. 
    
31. %对信号进行强制性消噪处理并图示结果
    
32. cdd3=zeros(1,length(cd3));
    
33. cdd2=zeros(1,length(cd2));
    
34. cdd1=zeros(1,length(cd1));
    
35. c1=[ca3 cdd3 cdd2 cdd1];
    
36. s1=waverec(c1,l,'db1');
    
37. subplot(2,2,2);
    
38. plot(s1);
    
39. title('强制消噪后的信号');grid;
    
40. %用默认阈值对信号进行消噪处理并图示结果
    
41. %用ddencmp函数获得信号的默认阈值,使用wdencmp命令函数实现消噪过程
    
42. [thr,sorh,keepapp]=ddencmp('den','wv',s);
    
43. s2=wdencmp('gbl',c,l,'db1',3,thr,sorh,keepapp);
    
44. subplot(2,2,3);
    
45. plot(s2);
    
46. title('默认阈值消噪后的信号');grid;
    
47. %用给定的软阈值进行消噪处理
    
48. cd1soft=wthresh(cd1,'s',1.465);
    
49. cd2soft=wthresh(cd2,'s',1.823);
    
50. cd3soft=wthresh(cd3,'s',2.768);
    
51. c2=[ca3 cd3soft cd2soft cd1soft];
    
52. s3=waverec(c2,l,'db1');
    
53. subplot(2,2,4);
    
54. plot(s3);
    
55. title('给定软阈值消噪后的信号');grid
    
56. 
    
57. %function y=snr(s,s1);%s是原始信号，s1，s2，s3分别是三种阈值去噪后的信号
    
58. y1=sum(s1.^2); y2=sum(s2.^2);y3=sum(s3.^2);
    
59. y12=sum((s-s1).^2);y22=sum((s-s2).^2);y32=sum((s-s3).^2);%y2=sum((abs(x1)-abs(x2)).^2)
    
60. y13=10*log10((y1/y12))
    
61. y23=10*log10((y2/y22))
    
62. y33=10*log10((y3/y32))

复制代码

mingmingtree

回复 16 # xiezhh 的帖子

按照您上次教我的，我对自己生成的含噪声信号用小波的三种阈值进行去噪，可是为何信噪比都是负数呢？是程序错了？改怎么纠正呢？当您有空读到我的请教时，能帮忙指点一下不呢？万分感激哦！！又来麻烦您了。

mingmingtree

回复 16 # xiezhh 的帖子

如果说信噪比可以为负数，那么怎么区别哪种去噪方法更好呢？
y13 =

   -3.0545


y23 =

   -3.0545


y33 =

   -1.9478

若凡

matlab就是边练习边提高，阅读帖子也能积累知识！:handshake