m0708.m

来自「MATLAB仿真技术与实例应用教程」· M 代码 · 共 13 行

df=0.01;
fs=5;
ts=1/fs;
t=[-5:ts:5];
x=zeros(1,length(t));
x(16:26)=t(16:26)+2;
[X,x1,df1]=fft_seq(x,ts,df);
f=[0:df1:df1*(length(x1)-1)]-fs/2;
X1=X/fs;
plot(f,fftshift(abs(X1)));