实验二连续时间的频域分析,浙江工业大学.docxVIP

Matlabs实验 浙江工业大学 信号与系统上机作业 光信1301 三、实验内容和要求 实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。 实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。 Q2-1 编写程序Q2_1，绘制下面的信号的波形图： 其中，?0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos(?0t)、cos(3?0t)、cos(5?0t) 和x(t) 的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。 抄写程序Q2_1如下： (1)M—文件： function f=y(t); n=input(请输入要累加的次数n的值); w=0.5*pi; f=0; for i=1:n; f=f+1/i*sin(i*pi/2)*cos(i*w*t); end; (2)命令窗口文件： w=0.5*pi;n=100;t=-pi:0.00001:pi; x1=cos(w*t);x2=cos(3*w*t);x3=cos(5*w*t); x4=y(t); n=100; subplot(221);plot(t,x1); title(函数cos（wt）的波形图);grid on;xlabel(时域t); xlabel(时域t); subplot(222);plot(t,x2); title(函数cos（3wt）的波形图);grid on;xlabel(时域t); subplot(223);plot(t,x3); title(函数cos（5wt）的波形图);grid on;xlabel(时域t); subplot(224);plot(t,x4); 请输入要累加的次数n的值100 执行程序Q2_1所得到的图形如下： Q2-2 反复执行程序Program2_2，每次执行该程序时，输入不同的N值，并观察所合成的周期方波信号。通过观察，你了解的吉伯斯现象的特点是： 了解到的吉布斯现象：在不连续点附近会有9%的过冲，且这个过冲大小不会随N的增加而变化，只是N越大，过冲会越靠近不连续点。 1、周期信号的傅里叶级数与GIBBS现象 给定如下两个周期信号： Q2-3 仿照程序Program2_2，编写程序Q2_3，计算并绘制出原始信号x1(t) 的波形图，用有限项级数合成的y1(t) 的波形图，以及x1(t) 的幅度频谱和相位频谱的谱线图。 编写程序Q2_3如下： T=2;dt=0.00001;t=-4:dt:4; x11=(t+1).*(u(-t)-u(-t-1))-(t-1).*(u(t-dt)-u(t-dt-1));%用单位阶跃和斜坡产生一个周期内的信号% x1=0;for m=-1:1; x1=x1+(t-m*T+1).*(u(-t+m*T)-u(-t+m*T-1))-(t-m*T-1).*(u(t-m*T-dt)-u(t-m*T-dt-1)); %通过时移，相加产生三个周期的信号x1（t）% end; w0=2*pi/T;%信号基波频率% N=input(Type in the number of the harmonic components N=:); Type in the number of the harmonic components N=:20 L=2*N+1;%总的谐波次数% for k=-N:1:N; ak(N+1+k)=(1/T)*x11*exp(-j*k*w0*t)*dt;%利用矩阵的乘法计算傅里叶积分的值% end phi=angle(ak);%一个angle就可以计算出ak的元素（傅里叶系数）的相位角% y=0; for q=1:L; y=y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T);%计算原信号的傅里叶级数表示结果% end; subplot(221);plot(t,x1);title(the original signal x(t));axis([-3,3,-1.5,1.5]);xlabel(time t);grid on; subplot(223);plot(t,y);title(the original signal y(t));axis([-3,3,-1.5,1.5]);xlabel(time t);grid on; subplot(222);k=-N:N;stem(k,abs(ak),k.);grid on; subplot(224);s