信号和系统实验报告实验3周期信号的频谱分析

 end phi=angle(ak); plot(t,x1); axis([-4 4 0 1.2]); grid on; title('The signal x1(t)'); xlabel('Time t (sec)'); ylabel('signal x1(t)');(5)仿照程序3_1,编写程序Q3_5,以计算x2(t) 的傅里叶级数的系数(不绘图)。程序如下:clc,clear,closeallT=2;dt=0.00001;t=-3:dt:3;x=ut(t+0.2)-ut(t-0.2-dt);x2=0;form=-1:1x2=x2+ut(t+0.2-m*T)-ut(t-0.2-m*T)-ut(t-0.2-m*t-dt);endw0=2*pi/T;N=10;L=2*N+1fork=-N:N;ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t')*dt;endphi=angle(ak);plot(t,x2);axis([-2.52.501.2]);gridon;title('Thesignalx2(t)');xlabel('Timet(sec)');ylabel('signalx2(t)'); (6)仿照程序3_2,编写程序Q3_6,计算并绘制出原始信号x1(t) 的波形图,用有限项级数合成的y1(t) 的波形图,以及x1(t) 的幅度频谱和相位频谱的谱线图。程序如下:clc,clear,closeallT=2;dt=0.00001;t=-3:dt:3;x=(t+1).*(ut(t+1)-ut(t))-(t-1).*(ut(t)-ut(t-1));x1=0;form=-2:2x1=x1+(t+1-m*T).*(ut(t+1-m*T)-ut(t-m*T))-(t-1-m*T).*(ut(t-m*t)-ut(t-1-m*t));endw0=2*pi/T;N=10;