sifouri.m

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

function [t,tdata]=sifouri(f,fdata)
%f=[0 1/3 2/3]; fdata=[6	-1.5 + 0.86603i	-1.5 - 0.86603i];

N=length(fdata);
ts=1/(f(2)-f(1));
dt=ts/N;
t=[0:dt:dt*(N-1)];
tdata=ifft(fdata);
tdata=real(tdata);
subplot(2,1,1);
plot(f,abs(fdata));

xlabel('频率');
ylabel('幅度');

subplot(2,1,2);
plot(t,tdata);
axis([0 max(t) -max(tdata) max(tdata)]);
xlabel('时间');
ylabel('幅度');



pause;