我正在将所谓的维纳滤波器应用于心跳心电信号。我的问题是傅里叶和傅里叶逆变换,理论上应该很容易。 获得Pxx和Pxxn后。用于构造的期望信号的PSD和我们的信号中的噪声在中心移位零频率。我很确定我做对了。
% frequency response function
for z= 1:Nfft
Hf(z)=Pxx(z)/(Pxx(z)+(Pxxn(z)));
end
% fourier transform on one ecg cycle
ecgf=fftshift(abs(fft(ecg(100:800),Nfft)));
% getting the output signal
Yf=ecgf.*Hf;
figure;
plot(f,Yf,'green');
% inverse fourier transform
Y=real(ifft(fftshift((Yf)),Nfft));
figure;
plot(Y);
我做错了什么?并且可以将ifft的输出长度更改为700以匹配输入信号。