一维信号小波去噪，傅里叶幅值和相位角计算

Halcom

Spectral Subtraction：

1.     E=SIG;   % 输入的一维信号
   
2.     n=size(E);
   
3.     fft_E = fft(E);
   
4.    
   
5.     MEMSNoise=load('data_2015051206Noise.mat');
   
6.     MEMSNoiseUD=MEMSNoise.MEMSUD;
   
7.     MEMSNoiseUD=(MEMSNoiseUD-mean(MEMSNoiseUD))*100;
   
8.     fft_n = fft(MEMSNoiseUD(1:n));
   
9.     E_noise = sum(abs(fft_n)) / N;
   
10.     mag_E = abs(fft_E);        % 幅值
    
11.     phase_E = angle(fft_E);    % 相位角
    
12.     mag_E = mag_E(1:2) - E_noise(1:2);
    
13.     mag_E(mag_E<0)=0;   
    
14.     fft_E = mag_E .* exp(1i.*phase_E);
    
15.     E = ifft(fft_E);     % 相位角不变，改变幅值

复制代码

参考：去随机飘逸噪声