南通大学罗开元
t=[t1,t2,t3];
f=[zeros(1,length(t1)),ones(1,length(t2)),zeros(1,length(t3))]; m0=5;
f1=[zeros(1,m0),f];f2=f1(1:length(f)) w=-2*pi:0.1:2*pi; F=f*exp(-j*t'*w)*dt;
subplot(2,2,1);plot(w,abs(F));grid; subplot(2,2,2);plot(w,angle(F));grid; F2=f2*exp(-j*t'*w)*dt;
subplot(2,2,3);plot(w,abs(F2));grid; subplot(2,2,4);plot(w,angle(F2));grid; 仿真结果:
19
南通大学罗开元
频移性质：
Command Window
dt=0.1; N=500; door__width=10; tao=door__width/2; t1=-(N-1):dt:-tao+1; t2=-tao:dt:tao; t3=tao-1:dt:N-1; t=[t1,t2,t3];
f=[zeros(1,length(t1)),ones(1,length(t2)),zeros(1,length(t3))]; w1=0.5*pi;
f2=f.*exp(-j*w1*t); w=-2*pi:0.1:2*pi; F=f*exp(-j*t'*w)*dt;
subplot(2,2,1);plot(w,abs(F));grid; subplot(2,2,2);plot(w,angle(F));grid; F2=f2*exp(-j*t'*w)*dt;
subplot(2,2,3);plot(w,abs(F2));grid; subplot(2,2,4);plot(w,angle(F2));grid; 仿真结果：
20
南通大学罗开元
(3)画出g（*g（的频谱图，验证时域卷积定理。 ?t）?t）Command Window
dt=0.1;N=50; door__width=8; tao=door__width/2; t1=-(N-1):dt:-tao+dt; t2=-tao:dt:tao; t3=tao-dt:dt:N-1; t=[t1,t2,t3];
f=[zeros(1,length(t1)),ones(1,length(t2)),zeros(1,length(t3))]; [k,f2]=myconv(f,f,t,t,dt) subplot(2,2,1);stairs(t,f); subplot(2,2,2);plot(k,f2); w=-2*pi:0.1:2*pi;
21
南通大学罗开元
F=f*exp(-j*t'*w)*dt;
subplot(2,2,3);plot(w,abs(F));grid; F2=f2*exp(-j*t'*w)*dt;
subplot(2,2,4);plot(w,abs(F2));grid; 得到波形如下图所示：
（4）画出g（和?Sa（? t/2）的频谱图，比较两者的联系，验证傅里叶变换的对?t）称性；
Command Window
dt=0.03;N=10; door__width=10; tao=door__width/2; t1=-(N-1):dt:-tao-dt; t2=-tao:dt:tao; t3=tao-dt:dt:N-1; t=[t1,t2,t3];
22
南通大学罗开元
f=[zeros(1,length(t1)),ones(1,length(t2)),zeros(1,length(t3))]; subplot(2,2,1);stairs(t,f);grid;title('门函数波形'); f2=10*sin(5*t)./(5*t);
subplot(2,2,2); plot(t,f2);grid;title('Sa函数波形'); w=-2*pi:0.1:2*pi; F=f*exp(-j*t'*w)*dt;
subplot(2,2,3);plot(w,real(F));grid;title('门函数频谱'); F2=f2*exp(-j*t'*w)*dt;
subplot(2,2,4);plot(w,real(F2));grid;title('Sa函数频谱'); 仿真结果：
(5)画sin(2?t)和g?(t)sin(2?t)的振幅频谱图，验证频域卷积定理。
Command Window
dt=0.01;N=20; door_width=5; tao=door_width/2;
23