关于 Matlab 中 FFT 的补零
About zero-padding for the FFT in the Matlab
我对 fft
的零填充有疑问。我 运行 fft
使用零填充和不使用零填充进行比较。
sf = 100; %sampling frequency
dt=1/sf; %time sampling interval
L = 10; %Length of signal
t = linspace(0,L,L/dt+1);
%zero-padding
nfft = 2^nextpow2(length(t)); %expansion of the data length for fft
t(length(t)+1:nfft) = 0; L = dt*nfft; t = linspace(0,L,L/dt+1);
t(end)=[];
fr = 4; %frequency
data = cos(2*pi*fr1*t);
df = sf/length(data); %frequency increment
f = (0:length(data)/2)*df;
fft_result =fft(data)/length(data);
spec_fft = 2*abs(fft_result); %spectrum
pha_fft = angle(fft_result); %phase
pha_fft = rad2deg(pha_fft);
subplot(2,1,1);
stem(f,spec_fft(1:length(f)));
subplot(2,1,2);
stem(f,pha_fft(1:length(f)));
而且我可以看到两个结果图像之间的差异。
当我fft
不补零时,振幅显示清晰,相位也清晰(我认为其他频率的相位是由于非常小的振幅值而不是零)。但是当我使用零填充时,我可以看到我输入的附近频率 (4Hz) 的振幅显示出不同的方面,并且在我看来相位的结果是 st运行ge 。当我使用零填充时我的代码有什么问题吗?
*Cris Luengo 评论的附加问题
我试图通过扩展长度向数据填充零。
nfft = 2^nextpow2(length(t)); %expansion of the data length for fft
data(length(t)+1:nfft) = 0;
当我绘制数据时,我得到了
如您所见,超过 10 的值为零。
我得到了这个结果。
不知道我的成绩好不好
如果 FFT 幅度结果接近于零,则可能只是数值噪声(随机量化和舍入误差)。这种数字噪音的相位是无稽之谈。我通常将相位设置为零,除非相应的幅度高于实际非零频谱的某个阈值。
我对 fft
的零填充有疑问。我 运行 fft
使用零填充和不使用零填充进行比较。
sf = 100; %sampling frequency
dt=1/sf; %time sampling interval
L = 10; %Length of signal
t = linspace(0,L,L/dt+1);
%zero-padding
nfft = 2^nextpow2(length(t)); %expansion of the data length for fft
t(length(t)+1:nfft) = 0; L = dt*nfft; t = linspace(0,L,L/dt+1);
t(end)=[];
fr = 4; %frequency
data = cos(2*pi*fr1*t);
df = sf/length(data); %frequency increment
f = (0:length(data)/2)*df;
fft_result =fft(data)/length(data);
spec_fft = 2*abs(fft_result); %spectrum
pha_fft = angle(fft_result); %phase
pha_fft = rad2deg(pha_fft);
subplot(2,1,1);
stem(f,spec_fft(1:length(f)));
subplot(2,1,2);
stem(f,pha_fft(1:length(f)));
而且我可以看到两个结果图像之间的差异。
当我fft
不补零时,振幅显示清晰,相位也清晰(我认为其他频率的相位是由于非常小的振幅值而不是零)。但是当我使用零填充时,我可以看到我输入的附近频率 (4Hz) 的振幅显示出不同的方面,并且在我看来相位的结果是 st运行ge 。当我使用零填充时我的代码有什么问题吗?
*Cris Luengo 评论的附加问题
我试图通过扩展长度向数据填充零。
nfft = 2^nextpow2(length(t)); %expansion of the data length for fft
data(length(t)+1:nfft) = 0;
当我绘制数据时,我得到了
我得到了这个结果。
不知道我的成绩好不好
如果 FFT 幅度结果接近于零,则可能只是数值噪声(随机量化和舍入误差)。这种数字噪音的相位是无稽之谈。我通常将相位设置为零,除非相应的幅度高于实际非零频谱的某个阈值。