希望能帮到你

1. N=500;fs=1000;  %设置采样点及采样频率
   
2. n=0:N-1;
   
3. t=n/fs;
   
4. s=sin(50*2*pi*n/fs)+sin(300*2*pi*n/fs);  %信号生成
   
5. s1=awgn(s,2);   %信噪比2dB
   
6. s_ft1=fftshift(abs(fft(s1,N)))/fs;
   
7. s2=awgn(s,0);   %信噪比0dB
   
8. s_ft2=fftshift(abs(fft(s2,N)))/fs;
   
9. s3=awgn(s,-2);  %信噪比-2dB
   
10. s_ft3=fftshift(abs(fft(s3,N)))/fs;
    
11. w=fs*(n-N/2)/N;   %频域坐标转换
    
12. subplot(311);plot(w,s_ft1);title('信噪比为2dB频谱图');
    
13. subplot(312);plot(w,s_ft1);title('信噪比为0dB频谱图');
    
14. subplot(313);plot(w,s_ft1);title('信噪比为-2dB频谱图');
    
15. 
    
16. 
    
17. 
    
18. 第二个
    
19. 
    
20. N=500;fs=1000;  %设置采样点及采样频率
    
21. n=0:N-1;
    
22. t=n/fs;
    
23. s_i=3*sin(100*pi*t)+2*sin(68*pi*t)+5*cos(72*pi*t);  %原始信号
    
24. s_n=wgn(1,N,3);          %噪声
    
25. s=s_i+s_n;               %信号声称
    
26. scale=[1 1.2 1.4 3];     %设置尺度
    
27. 
    
28. %% 四种小波变换
    
29. subplot(221);s_c_h1=cwt(s,scale,'haar','plot');title('Haar wavelet');
    
30. subplot(222);s_c_h2=cwt(s,scale,'mexh','plot');title('Mexican hat wavelet');
    
31. subplot(223);s_c_h3=cwt(s,scale,'morl','plot');title('Morlet wavelet');
    
32. subplot(224);s_c_h4=cwt(s,scale,'meyr','plot');title('Meyer wavelet');
    
33. 
    
34. figure;   %以Mexican hat小波为例
    
35. subplot(221);plot(s_c_h2(1,:));title('尺度为1');
    
36. subplot(222);plot(s_c_h2(2,:));title('尺度为1.2');
    
37. subplot(223);plot(s_c_h2(3,:));title('尺度为1.4');
    
38. subplot(224);plot(s_c_h2(4,:));title('尺度为3');

复制代码

评分

参与人数 1威望 +5 积分 +5 +5 收起 理由 A670521546 + 5 + 5 + 5 帮助他人

查看全部评分

学习小波变换ing

正学习小波，多谢楼主分享。