sfouri.m

来自「语音信号通过滤波器分析」· M 代码 · 共 21 行

function [f,fdata]=sfour(t,tdata)
%t=[0 1 2];tdata=[1 -2 3];
N=length(tdata);
fs=1/(t(2)-t(1));
df=fs/N;
f=[0:df:df*(N-1)]%-fs/2;
fdata=fft(tdata);
%fdata=fftshift(fdata/fs);

subplot(2,1,1);
plot(t,tdata);
xlabel('时间');
ylabel('信号');

subplot(2,1,2);
plot(f,abs(fdata));
xlabel('频率');
ylabel('幅度');
axis([0 fs 0 max(abs(fdata))])
pause;