信号与线性系统第二次试验选读.docVIP

河南农业大学 信号与线性系统分析 上 机 实 验 报 告 学院：机电工程学院 科目：信号与系统 班级：13级电信一班 姓名：王驿钊 学 号：1304101029 【实验二】实验基本知识与操作 （1）掌握使用 MATLAB 实现信号的卷积运算、卷积的可视化。 （2）用 MATLAB 分析周期矩形脉冲、三角波脉冲的频谱。 ①设计要求：某LTI系统的冲激响应h(t)=e-0.1t,输入u(t)如图2-1（a初始条件为零，求系统的响应y(t)。 图一 卷积的计算 设计要求：设二阶连续系统，其特性可用常微分方程表示：d2y/dt2+2(dy/dt)+8y=u 求其冲激响应。若输入为u=3t+cos(0.1t),求其零状态响应。 图二 LTI系统的零状态响应设计要求： 周期电流、电压（统称其为信号）f(t)可展开为直流与各次谐波之和，即 式中(=2(/T是基波角频率，T为周期。 周期信号的有效值定义为 若用各谐波有效值 则表示为 全波整流电压Us(t)的波形，用傅立叶级数可求得 可写出其展开式为（它只含直流和偶次谐波，令k=2n） 若Um=100V,频率f=50Hz,(相应的T=0.02S,(1=100( rad/s),分别计算其有效值Us1和Us2（取至六次谐波），并求Us2的误差。 图三 周期信号的频谱 设计要求：如图2-5(a)的矩形脉冲信号，求其在(=-40rad/s~40rad/s区间的频谱。 图四 非周期信号的频谱分析 【实验三】实验基本知识与操作 （1）傅里叶性质 ①方波分解为多次正弦波之和 设计要求：如图3-1所示的周期性方波，其傅立叶级数为： f(t)=4/([sint+1/3(sin3t)+…+(1/(2k-1))sin(2k-1)t+…] k=1,2,… 用MATLAB演示谐波合成情况。 t=0:0.01:2*pi; y=sin(t);plot(t,y),figure(gcf),pause y=sin(t)+sin(3*t)/3;plot(t,y),pause y=sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/9; plot(t,y) y=zeros(10,max(size(t)));x=zeros(size(t)); for k=1:2:19 x=x+sin(k*t)/k; y((k+1)/2, : )=x; end pause, figure(1),plot(t,y(1:9, : )),grid line([0,pi+0.5],[pi/4,pi/4]) text(pi+0.5,pi/4,pi/4) halft=ceil(length(t)/2);pause, figure(2),mesh(t(1:halft),[1:10],y(: ,1:halft)) ②周期信号的滤波 设计要求：如图3-3滤波电路，已知L=400mH，C=10(F,R1=200(.如激励电压us(t)为全波整流信号，Um=100V,(1=100(rad/s,求负载R两端的直流和各次谐波（它只含偶次谐波）分量。 图3-3滤波电路 ①利用MATLAB计算H(Z)的零极点与系统稳定性。 设计要求：已知一离散因果LTI系统的系统函数为： 利用MATLAB画出该系统的零极点分布图，求系统的单位冲激响应h[k]和幅频响应|H（ej(）|，并判断系统是否稳定。 分析：根据已知的H（z）,利用MATLAB中的zplane函数即可画出系统的零极点分布图。而利用impz函数和freqz函数求系统的单位冲激响应和频率响应时，一般需要将H（z）改写为下列形式： 参考程序： ②求离散系统在各种输入下的响应 设计要求：二阶巴特沃斯低通数字滤波器的频率响应。 二阶巴特沃斯低通数字滤波器的系统函数（传递函数）为 求其频率响应并做图（0~2(）。 参考程序： b=[1,2 ,1];a=[2+sqrt(2),0,2-sqrt(2)]; N=input(N=); w=[0:N-1]*pi/N; H=polyval(b,exp(i*w))./polyval(a,exp(i*w)); figure(1) subplot(211),plot(w,abs(H)),grid subplot(212),plot(w,unwrap(angle(H))),grid figure(2) subplot(211),semilogx(w,20*log10(abs(H))),grid subplot(212),sem