小波滤波下的信号频谱分析

Halcom

1. clc,clear,close all
   
2. warning off
   
3. %% 随机分析一个信号
   
4. load('data.mat')  % 加载数据
   
5. x = RC304(:, 2);  % 信号
   
6. Ts = 7.8125e-5;   % 采样时间
   
7. Fs = 1/Ts;        % 采样频率
   
8. %% 小波滤波
   
9. [thr,sorh,keepapp,crit] = ddencmp('cmp','wp',x);
   
10. [xc,wpt,perf0,perfl2] = wpdencmp(x,sorh,3,'db2',crit,thr,keepapp);
    
11. [c,l] = wavedec(xc,3,'db1');   % 多尺度小波分解
    
12. %% FFT分析
    
13. L = length(c);
    
14. NFFT = 2^nextpow2(L); % Next power of 2 from length of y
    
15. c_fft = fft(c, NFFT)/L;
    
16. figure(1)
    
17. plot( Fs/2*linspace(0, 1, NFFT/2+1), abs( c_fft(1:NFFT/2+1) ) )

复制代码