北理工信号与系统实验5

《北理工信号与系统实验5》由会员分享，可在线阅读，更多相关《北理工信号与系统实验5（13页珍藏版）》请在装配图网上搜索。

1、优质文本实验5 连续时间系统的复频域分析一、 实验目的1 掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB实现方法。2 学习和掌握连续时间系统函数的定义及复频域分析方法。3 掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。二、 实验原理与方法1. 拉普拉斯变换连续时间信号x(t)的拉普拉斯变换、反变换定义为MATLAB中，可以采用符号数学工具箱的laplace函数和ilaplace函数进行拉氏变换和拉氏反变换。L=laplace(F)符号表达式F的拉氏变换，F中时间变量为t，返回变量为s的结果表达式。L=laplace(F,t)用t替换结果中的变量t。F=ilaplace(L

2、)以s为变量的符号表达式L的拉氏反变换，返回时间变量为t的结果表达式。F=ilaplace(L,x)用x替换结果中的变量t。2. 连续时间系统的系统函数连续时间系统的系统函数是系统单位冲击响应的拉氏变换连续时间系统的系统函数还可以由输入和输出信号的拉氏变换之比得到单位冲击响应h(t)反映了系统的固有性质，而H(s)从复频域反映了系统的固有性质。由上式描述的连续时间系统，其系统系数为s的有理函数3. 连续时间系统的零极点分析系统的零点指使上式的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s平面上，零点用表示，极点用表示，得

3、到零极点分布图。由零点定义可知，零点和极点分别指上式分子和分母多项式的根，利用MATLAB求多项式的根可以通过roots来实现，该函数调用格式为：r=roots(c) c为多项式的系数向量，返回值r为多项式的根向量。分别对上式的分子、分母多项式求根即可求得零极点。此外，在MATLAB中还提供了pzmap函数来求取零极点和绘制系统函数的零极点分布图，该函数的调用格式为：Pzmap(sys) 绘出系统模型sys描述的系统的零极点分布图。p,z=pzmap(sys) 这种调用方法返回极点和零点，而不绘出零极点分布图。其中sys为系统传函模型，由t命令sys=tf(b,a)实现，b、a为传递函数的分子

4、和分母多项式的系数向量。MATLAB还提供了两个专用函数tf2zp和zp2tf来实现系统传递函数模型和灵机电网增益模型的转换，其调用格式为z,p,k=th2zp(b,a) b,a=th2zp(z,p,k)其中b、a为传递函数的分子和分母多项式的系数向量，返回值z为零点列向量，p为极点列向量，k为系统函数零极点形式的增益。三、 实验内容1.系统的冲激响应，输入信号，试采用复频域的方法求解系统的响应，编写MATLAB程序实现。解：MATLAB程序如下：h=sym(heaviside(t)-heaviside(t-2);H=laplace(f)x=sym(heaviside(t);X=laplace

5、(x)Y=H*X结果如下：H =1/s - exp(-2*s)/sX =1/sY =-(exp(-2*s)/s - 1/s)/s2.因果连续时间系统的系统函数分别如下： 试采用MATLAB画出零极点分布图，求解系统的冲激响应和频率响应，并判断系统是否稳定。1 零极点分布图b=1;a=1 2 2 1;sys=tf(b,a);p,z=pzmap(sys)pzmap(sys)结果如下：p = -1.0000 + 0.0000i -0.5000 + 0.8660i -0.5000 - 0.8660iz = Empty matrix: 0-by-1冲激响应b=1;a=1 2 2 1;impulse(b,

6、a)频率响应b=1;a=1 2 2 1;H,w=freqs(b,a);subplot(211);plot(w/pi,abs(H);xlabel(omega(rad/s);ylabel(Magnitude);title(|H(jomega)|);grid on;subplot(212);plot(w/pi,angle(H)/pi);xlabel(omega(rad/s);ylabel(Phase);title(|phi(omega)|);grid on;2零极点分布图b=1 0 1;a=1 2 -3 3 3 2;sys=tf(b,a);p,z=pzmap(sys)pzmap(sys)结果如下：p

7、 = -3.1704 + 0.0000i 0.9669 + 0.9540i 0.9669 - 0.9540i -0.3817 + 0.4430i -0.3817 - 0.4430iz = 0.0000 + 1.0000i 0.0000 - 1.0000i冲激响应b=1 0 1;a=1 2 -3 3 3 2;impulse(b,a)频率响应b=1 0 1;a=1 2 -3 3 3 2;H,w=freqs(b,a);subplot(211);plot(w/pi,abs(H);xlabel(omega(rad/s);ylabel(Magnitude);title(|H(jomega)|);grid

8、on;subplot(212);plot(w/pi,angle(H)/pi);xlabel(omega(rad/s);ylabel(Phase);title(|phi(omega)|);grid on;3.连续时间系统函数的极点位置分别如下所示 试用MATLAB绘制上述6种不同情况下，系统函数的零极点分布图，并绘制相应冲激响应的时域波形，观察并分析系统函数极点位置对冲激响应时域特性的影响。1b=1;a=1 0;sys=tf(b,a);subplot(211)pzmap(sys);subplot(212)impulse(sys);2b=1;a=1 2;sys=tf(b,a);subplot(21

9、1)pzmap(sys);subplot(212)impulse(sys);3b=1;a=1 -2;sys=tf(b,a);subplot(211)pzmap(sys);subplot(212)impulse(sys);4b=1;a=1 0 4;sys=tf(b,a);subplot(211)pzmap(sys);subplot(212)impulse(sys,10);axis(0 10 -1 1)5b=1;a=1 2 17;sys=tf(b,a);subplot(211)pzmap(sys);subplot(212)impulse(sys);6b=1;a=1 -2 17;sys=tf(b,a

10、);subplot(211)pzmap(sys);subplot(212)impulse(sys);4. 连续时间系统的系统函数分别如下： 上述3个系统具有相同的极点，只是零点不同，试用MATLAB分别绘出系统的零极点分布图及相应冲激响应的时域波形，观察并分析系统函数零点位置对冲激响应时域特性的影响。1b=1;a=1 2 17;sys=tf(b,a);subplot(211)pzmap(sys);subplot(212)impulse(sys);2b=1 8;a=1 2 17;sys=tf(b,a);subplot(211)pzmap(sys);subplot(212)impulse(sys)

11、; 3b=1 -8;a=1 2 17;sys=tf(b,a);subplot(211)pzmap(sys);subplot(212)impulse(sys);四、 实验总结通过本次实验掌握了拉普拉斯变换及其反变换的定义，并掌握了MATLAB实现方法；学习和掌握了连续时间系统函数的定义及复频域分析方法；掌握了系统零极点的定义，加深理解了系统零极点分布与系统特性的关系，包括零极点分布与冲激响应的关系、零极点分布与系统频率响应的关系、零极点分布与系统稳定性的关系等。通过一系列特定函数完成了实验中要求的各类变换，顺利画出频率响应图像。实验中的一些内容，以类比的形式将几道题并列给出，对于程序的编写有一定的训练作用，通过相似的程序编写，对函数的使用加深了印象。