信号与系统实验报告-wmy.docVIP

信号与系统实验报告 学 院: 核技术与自动化工程学院 专 业: 电气工程及其自动化 指导老师: 姓 名: 学 号: 日 期: 实验二 线性非时变系统的时域分析 实验目的 掌握在时域中对连续和离散时间线性非时变系统响应进行分析的方法。 实验内容与步骤 已知系统的微分方程如下，用MATLAB画出该系统的冲激响应及该系统在输入信号的零状态响应的波形。（改变取样的时间间隔p观察仿真的效果） 方程： 冲激响应 a=[1 3 2]; b=[3]; impulse(b,a) 在输入信号的零状态响应 a=[1 3 2]; b=[3]; p=0.01; t=0:p:10; x=exp(-2*t); lsim(b,a,x,t); 上述情况当p=1时，情况如图： 已知离散系统的差分方程为：y(n)+y(n-1)+0.25y(n-2)=x(n)用MATLAB画出该系统的单位函数响应。 单位函数响应： a=[1 1 0.25]; b=[1]; impz(b,a) 实验分析 实验总结 初步学习并且使用了MATLAB软件，对LTI系统进行了简单的时域分析，其中应当注意的问题如下： 首先，输入正确的函数调用格式，函数调用格式正确与否决定了连续系统的范围，以及冲击响应的数值解； 其次，注意MATLAB中一些常用函数的调用，例如：sqrt exp等等，能够用这些表达式正确的写出所编程序； 最后，根据不同的方程列写表达式时，注意系数的选择，缺项用0补齐。 实验四 系统的零极点分析 实验目的 掌握系统函数及零极点的概念； 掌握对连续和离散系统的稳定性进行分析方法。 实验原理和方法 连续系统的零极点分析 系统函数H(s) 在一个CTLTI系统可用一个线性常系 数微分方程的描述： 其中an…a0,bm…b0…均为常数 定义H（s） H(s)= =（计算时通常使an=1） 冲激函数h(t) 系统函数H(s) 系统零极点 H(s)=,假定为真分式 零点：使H(s)=0的根（即分子为0），z1z2……zm; 极点：H(s)=∞的根（即分母为零），p1p2……pn; H(s)= 可见:若系统的零极点已知，则 H(s)可知 因此，对系统的零极点进行分析不仅可确定系统的H(s),还可以判断系统的稳定性 零极点分布于系统稳定性的关系 判断系统稳定性的方法： 若输入有界，则输出有界； h(n)绝对可和； H(z)的所有极点在单位圆内； 因此，考察H(z)的极点分布，就可判断离散系统的稳定性。 实验内容 用MATLAB来绘制连续系统的零极点图 连续系统的零极点的位置可以用MATLAB中的多项式求根函数roots（）来求得，其调用格式为：p=roots(D):D是由多项式的系数构成的行向量。 用roots（）函数求出系统函数H（s）的零极点后，即可用plot命令在复平面上绘出系统函数的零极点。下面的函数可以求出系统函数的零极点并绘制零极点图。 Function[p,z]=ljdt(D,N) P=roots(D) %求系统的极点 Z=roots(N) %求系统的零点 p=p’; %将极点列向量转置为行向量 z=z’; %求将零点列向量转置为行向量 x=max(abs([p z])); %用来确定坐标轴的范围 x=x+0.1; y=x; hold on %重叠绘图 axis([-x x –y y]); %确定坐标轴显示范围 plot([-x x],[0 0]) %画出横坐标 plot([0 0],[-y y]) %画出纵坐标 plot(real(p),imag(p),’x’) %画出极点 plot(real(z),imag(z),’O’) %画出零点 title(‘’) %加标题 text(0.2,x-0.2,’虚轴’) %加文本标注 text(y-0.2,0.2,’实轴’) （1）已知以连续时间线性非系统的系统函数为 画出系统的零极点图并判断系统的稳定性。 程序： a=[1 2 -3 2 1]; b=[1 0 -4]; ljdt(a,b); p = -3.1300 0.7247 + 0.6890i 0.7247 - 0.6890i -0.3195 z = 2.0000 -2.0000 2）已知一离散时间线性非系统函数为画出系统的零极点图并判断系统的稳定性。 程序： a=[2