《实验3利用matlab分析连续系统及离散系统的复频域特性》由会员分享，可在线阅读，更多相关《实验3利用matlab分析连续系统及离散系统的复频域特性（11页珍藏版）》请在金锄头文库上搜索。
1、1广东技术师范学院实验报告学院： 自动化学院 专业： 班级： 成绩：姓名： 学号： 计算机编号：实验地点： 微机 6 室 实验日期： 指导教师签名：实验 （三） 项目名称：利用 MATLAB 分析连续系统及离散系统的复频域特性 一实验目的1.掌握 Laplace 变换的意义、基本性质及应用。2.掌握拉普拉斯变换的三维可视化表示。3.理解系统函数的零、极点分布(极、零图) 决定系统时间原函数的特性。4.掌握系统冲激响应。5. H(z)部分分式展开的 MATLAB 实现6. H(z)的零极点与系统特性的 MATLAB 计算二实验原理1Laplace 变换和逆变换定义为( 4 1 )jststdeF
2、tffs)(21)(0在 Matlab 中实现 Laplace 变换有两个途径：直接调用指令 laplace 和 ilaplace 进行；根据定义式 ( 4 1 )，利用积分指令 int 实现。相较而言，直接利用 laplace 和 ilaplace 指令实现机器变换要简洁一些。调用格式：L=laplace(F)F=ilaplace(L)2实现拉普拉斯曲面图及其可视化的步骤如下：a定义两个向量 x 和 y 来确定绘制曲面图的复平面横座标和纵座标的范围。b调用 meshgrid 函数产生包含绘制曲面图的 s 平面区域所有等间隔取样点的复矩阵。c计算复矩阵 s 定义的各样点处信号拉氏变换 F(s)
3、的函数值，并调用 abs 函数求其模。d调用 mesh 函数绘出其幅度曲面图。3在连续系统的复频域分析中，系统函数起着十分重要的作用，它包含了连续系统的固有特性。通过系统函数可以对系统的稳定性、时域特性、系统频率响应等系统特性进行分析。若连续系统的系统函数的零极点已知，系统函数便可确定下来，即系统函数预习情况 良好 操作情况 良好 考勤情况 全勤 数据处理情况 良好2H（s）的零极点分布完全决定了系统的特性。系统函数的零点和极点位置可以用matlab 的多项式求根函数 roots()来求得。用 roots()函数求得系统函数 H(s)的零极点后，就可以用 plot 命令在复平面上绘制出系统函数
4、的零极点图。4系统冲激响应 h（t）的时域特性完全由系统函数 H（s）的极点位置决定，H（s）的每一个极点将决定 h（t）的一项时间函数。显然， H（s）的极点位置不同，h（t）的时域特性也完全不同。用函数 residue（）求出 H（s）部分分式展开的系数后，便可根据其极点位置分布情况直接求出 H（s）的拉普拉斯反变换 h（t ） 。且利用绘制连续时间系统冲激响应曲线的 matlab 函数 impulse（） ，将系统冲激响应 h（t）的时域波形绘制出来。5利用 tf()函数、pole()函数、zero()函数和 pzmap()函数，能方便地求出系统函数的零极点，并绘出其零极点分布图。调用格
5、式：sys=tf(b,a); %b 为系统函数分子多项式系数构成的行向量；a 为分母多项式系数构成的行向量；sys 为系统函数对象。p=pole(sys); %输出参量 p 为返回包含系统函数所有极点位置的列向量。z=zero(sys);pzmap(sys);%用于绘制系统函数零极点分布图和计算系统函数的零极点位置6部分分式展开的 MATLAB 实现r,p,k=residuez(num,den)num,den 分别为 X(z)分子多项式和分母多项式的系数向量。r 为部分分式的系数，p 为极点， k 为多项式的系数。若为真分式，则 k 为零。7H(z) 的零极点与系统特性的 MATLAB 计算利
6、用 tf2zp 函数计算 H(z)的零极点，调用形式为z,p,k=tf2zp(b,a)b 和 a 分别为 H(z)分子多项式和分母多项式的系数向量。返回值 z 为零点、p 为极点、 k 为增益常数。H(z)零极点分布图可用 zplane 函数画出，调用形式为zplane(b,a)三实验内容1 试用 MATLAB 求函数 的拉普拉斯变换，绘出其零极)t(ue)t()t(h2点分布图。第 1 步：求出 h(t)的拉普拉斯变换syms t;h=exp(-t)+exp(-2*t);H=laplace(h)结果为：1/(s + 1) + 1/(s + 2)= 32s3-2 -1.5 -1 -0.5 0
7、0.5 1 1.5 2-2-1.5-1-0.500.511.52位位位位位位位位位位位位位（第 1 题）第 2 步：绘出其零极点分布图A=1 3 2;B=2 3;p=roots(A);q=roots(B);x=max(abs(p;q);x=x+0.1;y=x;hold on;axis(-x x -y y);axis(square);plot(-x x,0 0,m); plot(-0 0,-y y,m); plot(real(p),imag(p),x);plot(real(q),imag(q),o);title(连续时间系统的零极点图 );text(0,3.9,虚轴);text(3.7,0.2,
8、实轴);42 使用 Matlab 绘出下列信号拉普拉斯变换的三维曲面图。a. ()cos()2tfetu第一步：求出 f(t)的拉普拉斯变换syms t;h=exp(-t)*cos(t*pi/2);H=laplace(h)结果为：(s + 1)/(s + 1)2 + pi2/4)第二步：绘制拉普拉斯变换的三维曲面图x=-1:0.1:0.5;y=-5:0.1:5;x,y=meshgrid(x,y);s=x+i*y;H=(s + 1)./(s + 1).2 + pi2/4);F=abs(H);surf(x,y,F)colormap(spring);title(拉普拉斯变换曲面图 );xlabel(
9、实轴 );ylabel(虚轴 );-1-0.500.5-50505101520位位位位位位位位位位位位位第 2 题（a）5b. ()2sin()(4fttut第一步：求出 f(t)的拉普拉斯变换syms t;h=2*sin(2*t-pi/4);H=laplace(h)结果为：-(2*(2(1/2)*s)/2 - 2(1/2)/(s2 + 4)第二步：绘制拉普拉斯变换的三维曲面图x=-1:0.1:0.5;y=-5:0.1:5;x,y=meshgrid(x,y);s=x+i*y;H=-(2*(2(1/2)*s)/2 - 2(1/2)./(s.2 + 4);F=abs(H);surf(x,y,F)c
10、olormap(lines);title(拉普拉斯变换曲面图 );xlabel(实轴 );ylabel(虚轴 );-1-0.500.5-505024681012位位位位位位位位位位位位位第 2 题（b）63已知系统函数如下，试用 Matlab 绘出其零极点分布图，求出冲激响应，并判断系统是否稳定。 25431() 2sHss第一步：求零极点分布图A=1 2 -3 3 3 2;B=1 0 1;p=roots(A)q=roots(B)x=max(abs(p;q)x=x+0.1y=xhold onaxis(-x x -y y)axis(square)plot(-x x,0 0,r) plot(-0
11、0,-y y,r) plot(real(p),imag(p),x)plot(real(q),imag(q),o)title(连续时间系统的零极点图 )text(0,3,虚 轴)text(3,0.1,实轴)-3 -2 -1 0 1 2 3-3-2-10123位位位位位位位位位位位位位位位第 3 题（1）由于系统函数的极点有位于 s 右半平面部分，故系统不稳定74利用 Matlab 的 residuez 函数求下式的部分分式展开及对应的 hk。432216563()81zzH第 2 步：求出冲激响应A=1 2 -3 3 3 2;B=1 0 1;impulse(B,A)0 10 20 30 40 5
12、0 60 70-2.5-2-1.5-1-0.500.51x 1028 Impulse ResponseTime (seconds)Amplitude第 3 题（2）方法 2：A=1 2 -3 3 3 2;B=1 0 1;figure(1);zplane(B,A);xlabel(实轴 );ylabel(虚轴 );title(零极点分布图);-3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1-1.5-1-0.500.511.53位位位位位位位位位位85试画出系统的零极点分布图，求其单位冲激响应 hk和频率响应 H(ej) 。第一步：求出系统的零极点分布图A=3 3 -15 18 -12
13、;B=2 16 44 56 32;p=roots(A);q=roots(B);x=max(abs(p;q);x=x+0.1;y=x;hold on;axis(-x x -y y);axis(square);plot(-x x,0 0,g); plot(-0 0,-y y,g); plot(real(p),imag(p),x);plot(real(q),imag(q),o);title(连续时间系统的零极点图 );text(0,3.9,虚轴);text(3.7,0.2,实轴);grid on;A=3 3 -15 18 -12;B=2 16 44 56 32;r,p,k=residue(B,A)a
14、ngle,mag=cart2pol(real(r),imag(r)运行结果为：k=0.6667p =-3.2361 + 0.0000i 1.2361 + 0.0000i 0.5000 - 0.8660i 0.5000 + 0.8660ir =0.0574 + 0.0000i 11.7321 + 0.0000i -3.5614 + 1.4890i -3.5614 - 1.4890iangle =0 0 -2.7456 2.7456mag =0.0574 11.7321 3.8601 3.8601由此可得： 86.051.386.051.326.173.5467.)(H 74.274.2jzejz
15、ezzz jj 所以 )7456.2.cos(72.)(72.)(0.)(.0)(h 5.0361.2361. keuekuekk kk9-4 -3 -2 -1 0 1 2 3 4-4-3-2-101234位 位 位 位 位 位 位 位 位 位 位位 位位 位法 2：求出系统的零极点分布图A=3 3 -15 18 -12;B=2 16 44 56 32;figure(1);zplane(B,A);xlabel(实轴 );ylabel(虚轴 );title(零极点分布图);-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1-2-1.5-1-0.500.511.52位位位位位位位位位位第 5 题（1）10第 2 步：求其单位冲激响应 hk和频率响应 H(ej)A=3 3 -15 18 -12;B=2 16 44 56 32;h=impz(B,A);figure(2);stem(h)title(单位冲击响应)axis(1 12 -2500 1000);H,w=freqz(B,A);figure(3);plot
