信号频谱图实验报告.doc

信号频谱图实验报告.doc:..大连理工大学实验报告学院(系):专业:班级:姓 名:学号:组:实验时间:实验室:实验台:指导教师签字:成绩:实验一信号的频谱图一、、实验内容:已知周期三角信号如图所示,试求出该信号的傅里叶级数,利用MATLAB编程实现其各次谐波的叠加,并验证其收敛性。MATLAB源程序为:»t=-3::3;»nO=-3;n1=-1;t0=2;»fori=0:2tl=nO::nO+tO/2;xl=tl-nO;t2=nl-t0/2::nl;x2=-t2+nl;plot(tlzxlz'r',t2,x2/r');holdon;nO=nO+tO;nl=nl+tO;end»n_max二[1371531];»N=length(n_max);»fork=l:Nn=1;sum=O;while(n<(n_max(k)+l))b=4./pi/pi/n/n;y=b*cos(n*pi*t);sum=sum+y;n=n+2;endfigure;nO=-3;nl=-l;tO=2;fori=0:2tl=n0::n0+t0/2;xl=tl-n0;t2=nl-t0/2::nl;x2=-t2+nl;holdon;nO=nO+tO;nl=nl+tO;endy=sum+;plot(t,y;b');xlabeft'),ylabel('wove‘);holdoff;axis([--]);gridon;title(['themax=',num2str(n_max(k))])end运行结果:ihemax"3-hemax=15ihemax"=31-3 -2-10 12 3t8765432■■•••••。当周期三角信号的周期和三角信号的宽度变化时,试观察其频谱的变化。MATLAB源程序为:»fork=l:3;n=-30:30;tao=k;T=2*k; w=2*pi/T;x=n*tao*=sinc(x/pi);fn=tao*fnl.*fnl;subplot(3/l,k),stem(n*w/fn);gridontitle(['T=',num2str(2*k)]);axis([-30300k]);end运行结果:X=Columns1through9- - - - - - - --・ ・ - - - - - --- --------- -- & —Columns1through16-30 -29 -28 -27 -26-25 -24-23-22-21・20-18 -17-16 -15Columns17through32-14 -13 -12 -11 -10・9 -8・7-6・5-4・2-101Columns33through482 3 45 678910111213 1415 16 17Columns49through6118 19 202********** 272829 30X-Columns1through9・ ・--・-・・--31.