Matlab第二次作业.doc

Matlab 第二次作业 班级：通信11001班 姓名：张安娜学号：201006091题目：用Matlab说明频谱“栅栏效应”的影响。 已知三角波信号如图所示，采样间隔T=0.1s，选取NO=64、128和256。 所谓“栅栏效应”就是在一组离散的频率上观察加窗采样信号的频谱。从效果上看，就像通过一个栅栏来观看频谱一样。如果在某些频率上存在较大的频谱分量或者谱分量是0，那么这些谱分量不会如同脉冲的采样信号的频谱上显示出来。减小栅栏效应产生的频谱失真的方法是在保持采样频率不变的情况下，提高DFT的点数，这样可使“栅栏”更密。即频率分辨率更高。如果原始信号为有限区间信号，则可以通过补零来增加NO。但补零并不能提高准确度，因为它并未增加任何新的信号信息，它只是减小了频谱间距。解：Matlab的程序： %栅栏效应的影响T=0.1;NO=64;TO=NO*T;t=0:T:4;f=1-abs(0.5*(t-2),zeros(1,NO-length(t);%f=8*(u(t)-u(t-1);%f(1)=4;f(33)=4;Fr=T*fft(f);F_r=fftshift(abs(Fr);r=-NO/2:NO/2-1;w=linspace(-3*pi,3*pi,500);F=2*sinc(w/pi).2;figure(1)plot(w/pi,abs(F),r:),hold onw=r*2*pi/TO;stem(w/pi,F_r,.);hold offaxis(-3 3 0 inf);xlabel(omega xpi(rad/s)text(0.7,1.6,T=,num2str(T),s NO=,num2str(NO)text(0.7,1.4,fO=,num2str(1/TO),Hz)运行结果：Matlab的程序：%栅栏效应的影响T=0.1;NO=128;TO=NO*T;t=0:T:4;f=1-abs(0.5*(t-2),zeros(1,NO-length(t);%f=8*(u(t)-u(t-1);%f(1)=4;f(33)=4;Fr=T*fft(f);F_r=fftshift(abs(Fr);r=-NO/2:NO/2-1;w=linspace(-3*pi,3*pi,500);F=2*sinc(w/pi).2;figure(1)plot(w/pi,abs(F),r:),hold onw=r*2*pi/TO;stem(w/pi,F_r,.);hold offaxis(-3 3 0 inf);xlabel(omega xpi(rad/s)text(0.7,1.6,T=,num2str(T),s NO=,num2str(NO)text(0.7,1.4,fO=,num2str(1/TO),Hz)运行结果：Matlab的程序： %栅栏效应的影响T=0.1;NO=256;TO=NO*T;t=0:T:4;f=1-abs(0.5*(t-2),zeros(1,NO-length(t);%f=8*(u(t)-u(t-1);%f(1)=4;f(33)=4;Fr=T*fft(f);F_r=fftshift(abs(Fr);r=-NO/2:NO/2-1;w=linspace(-3*pi,3*pi,500);F=2*sinc(w/pi).2;figure(1)plot(w/pi,abs(F),r:),hold onw=r*2*pi/TO;stem(w/pi,F_r,.);hold offaxis(-3 3 0 inf);xlabel(omega xpi(rad/s)text(