一维信号小波去噪,傅里叶幅值和相位角计算
Spectral Subtraction:E=SIG; % 输入的一维信号
n=size(E);
fft_E = fft(E);
MEMSNoise=load('data_2015051206Noise.mat');
MEMSNoiseUD=MEMSNoise.MEMSUD;
MEMSNoiseUD=(MEMSNoiseUD-mean(MEMSNoiseUD))*100;
fft_n = fft(MEMSNoiseUD(1:n));
E_noise = sum(abs(fft_n)) / N;
mag_E = abs(fft_E); % 幅值
phase_E = angle(fft_E); % 相位角
mag_E = mag_E(1:2) - E_noise(1:2);
mag_E(mag_E<0)=0;
fft_E = mag_E .* exp(1i.*phase_E);
E = ifft(fft_E); % 相位角不变,改变幅值参考:去随机飘逸噪声
页:
[1]