资源描述：
《信号系统实验代码》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、四、实验内容1、求出下图中周期方波信号的频谱，并参照例3-1,并画出频谱图。(A=l,t=0.5,Ti=l)是在把周期信号在时域上右移了0.25个单位，频域不变MATLAB实现傅里叶级数计算的程序如下：dt=0.01;Tl=1;wl=2*pi/Tl;t=-Tl:dt:Tl;tau=0.5;A=1;f=A*(heaviside(t+0)-heaviside(t-tau));subplot(2，l，l)plot(t，f)axis([-Tl，Tl，-0.1,1.1])title(’f(t)时域波形i)N=10;n=-N:N;Fn=f^exp(-
2、j*t'*w1*n)*dt/T1;subplot(2，l，2)stem(n，Fn)holdondw=0.01;w=-N*wl:dw:N*wl;F=A*tau/Tl*sinc(w*tau/pi);plot(w/wl,F，'--’)title(’傅里叶级数Fn*)f(t)时域波形10.80.60.40.20•1-0.8-0.6-0.4-0.200.20.40.60.811)脉冲宽度t=0.5保持不变，分别取TOMt和T0=8t:,分别绘制相应的频谱图，并讨论周期乃与频谱的关系。T,越大，频谱图的幅度越小，衰减速度增快。f(t)时域波形10.8
3、0.60.40.20•4-3-2-101234T,=2tTl=ldt=0.01;T1=1;wl=2*pi/Tl;t=-Tl:dt:Tl;tau=0.5;A=1;f=A*(heaviside(t+0)-heaviside(t-tau));subplot(2,l,l)plot(t,f)axis([-Tl,Tl,-0.1,1.1])titleCf⑴时域波形•)N=10;n=-N:N;Fn=f^expGj^t’^wPnpdt/Tl;subplot(2,1，2)stem(n，Fn)holdondw=0.01;w=-N*wl:dw:N*wl;F=A*
4、tau/Tl*sinc(w*tau/pi);plot(w/wl,F，'--’)title(’傅里叶级数F_n')f(t)吋域波形1I.-J0.80.60.40.20-1-0.8-0.6-0.4-0.200.20.40.60.81T0=4tT0=2T!=2tTl=ldt=0.01;T1=2;wl=2*pi/Tl;t=-Tl:dt:Tl;tau=0.5;A=1;f=A*(heaviside(t+0)-heaviside(t-tau));subplot(2，l，l)plot(t，f)axis([-Tl,Tl，-0.1，1.1])titlc(T(
5、t)时域波形胃)N=10;n=-N:N;Fn=f*exp(-jw1*n)*dt/T1;subplot(2,l,2)stem(n，Fn)holdondw=0.01;w=-N*wl:dw:N*wl;F=A*tau/Tl*sinc(w*tau/pi);plot(w/wl,F，'--’)titleC傅里叶级数Fn*)f(t)时域波形10.80.60.40.20-2-1.5-1-0.500.511.52T0=8tt=0.5t0=4dt=0.01;T1=4;wl=2*pi/Tl;t=-Tl:dt:Tl;tau=0.5;A=1;f=A*(heavisi
6、de(t+0)-heaviside(t-tau));subplot(2，l，l)plot(t，f)axis([-Tl,Tl,-0.1,1.1])titleff(t)时域波形’)N=10;n=-N:N;Fn=Pexp(-j*t'*w1*n)*dt/T1;subplot(2,1,2)stem(n,Fn)holdondw=0.01;w=-N*wl:dw:N*wl;F=A*tau/Tl*sinc(w*tau/pi);plot(w/wl，F，’--’)titlef傅里叶级数F_n’)2)脉冲周期T1=l保持不变，分别取t=0.75、t=0.5和t=
7、0.25，分别绘制相应的频谱图，并讨论脉冲宽度t与频谱的关系。W吋域波形10.80.60.40.201-0.8-0.6-0.4-0.200.20.40.60.81T越大。幅度越大dt=0.01;T1=1;wl=2*pi/Tl;t=-Tl:dt:Tl;tau=0.75;A=1;f=A*(heaviside(t+0)-heaviside(t-tau));subplot(2，l，l)plot(t，f)axis([-Tl，Tl，-0.1，1.1])titlc(’f(t)时域波形’)N=10;n=-N:N;Fn=f*exp(-j*t'*wl*n)*
8、dt/Tl;subplot(2，l，2)stem(n,Fn)holdondw=0.01;w=-N*wl:dw:N*wl;F=A*tau/Tl*sinc(w*tau/pi);plot(w/wl,
