MATLAB——非周期序列的傅里叶变换

题目：
有一非周期方波信号x(t)的脉冲宽度为1ms，信号持续时间为2ms，在0~2ms区间外信号为0。试求其含有20次谐波的信号的频谱特性。求其逆变换并与原时间信号的波形进行比较。
这个题目不难，基本上根据傅里叶变换和反变换的公式就可以了。如果信号比较好表达的话，就是直接利用ifourier和fourier。但之前要定义变量 syms t。
代码：

T=2;
f1=1/T;
N=256;

t=linspace(0,T,N);
x=[ones(1,N/2),zeros(1,N/2)];
dt=T/(N-1);
f=linspace((-20*f1),(20*f1),N);
w=2*pi*f;
dw=40*pi*f1/(N-1);



X=x*exp(-j*t'*w)*dt; 
subplot(1,2,1);
plot(f,abs(X));

x2=X*exp(j*w'*t)/(pi)*dw
subplot(1,2,2);
plot(t,x,'r',t,x2);

运行结果：