我希望对信号进行幅度和相位重建。我的信号带宽为40KHz,起始频率为70Khz至110KHz。假设我的信号为x
。
nfft = length(x);
res = fft(x,nfft)/ nfft; % normalizing the fft
f = fs/2*linspace(0,1,nfft/2+1); % choosing correct frequency
res = res(1:nfft/2+1); % amplitude of fft
res2 = fft(res);
我想用
绘制频率与幅度的关系figure
plot(f,abs(res2))
其中振幅应位于70kHz至110kHz之间,频率与相位相似
figure
plot(f,angle(res2))
相位应分布在70kHz至110kHz之间。那么如何选择正确的频率轴呢。
答案 0 :(得分:0)
如果你想绘制信号x
的频率内容而不是已经转换的信号res
的频率内容,那么你应该摆脱这条线
res2 = fft(res);
然后plot(f,abs(res))
和plot(f,angle(res))
会正确显示从0到fs/2
的信号。
然后您可以使用以下方式更精确地放大70kHz-110kHz的感兴趣频率范围(假设为fs/2>110000
):
axis([70000 110000]);