LTI系统复频域分析的MATLAB实现.docVIP

实验项目名称：LTI系统复频域分析的MATLAB实现 上机实验题目： 拉氏变换与Z变换的基本性质在系统分析中的应用 实验项目的目的和任务： 掌握拉氏变换、Z变换的基本性质及其在系统分析中的典型应用 实验题目： 第9章 第10章 实验过程 9a_1 b=[1,5 0]; a=[1,2,3]; zs=roots(b) ps=roots(a) plot(real(zs),imag(zs),o); hold on plot(real(ps),imag(ps),x); grid axis([-5 2 -2 2]); 运行截图： 9a_2 b=[2,5,12]; a=[1,2,10]; zs=roots(b) ps=roots(a) plot(real(zs),imag(zs),o); hold on plot(real(ps),imag(ps),x); grid axis([-10 10 -10 10]); 运行截图： 9a_3 b=[2,5,12]; a=[1,4,14,20]; zs=roots(b) ps=roots(a) plot(real(zs),imag(zs),o); hold on plot(real(ps),imag(ps),x); grid axis([-10 10 -10 10]); 运行截图： 9c、经过拉氏变换之后得到的系统函数为： b=[1,2,5]; a=[1,-3]; zs=roots(b) ps=roots(a) plot(real(zs),imag(zs),o); hold on plot(real(ps),imag(ps),x); grid axis([-10 10 -10 10]); 运行截图： 10、dpzplot函数代码： function dpzplot(b,a) la=length(a); lb=length(b); if (lalb) b=[b zeros(1,la-lb)]; elseif (lbla) a=[a zeros(1,lb-la)]; end ps=roots(a); zs=roots(b); mx=max(abs([ps zs .95]))+0.05; clg axis([-mx mx -mx mx]); axis(equal); hold on w=[0:0.01:2*pi]; plot(cos(w),sin(w),.); plot([-mx mx],[0 0]); plot([0 0],[-mx mx]); text(0.1,1.1,Im,sc); text(1.1,0.1,Re,sc); plot(real(ps),imag(ps),rx); plot(real(zs),imag(zs),ro); numz=sum(abs(zs)==0); nump=sum(abs(ps)==0); if numz1 text(-.1,-.1,num2str(numz)); elseif nump1 text(-.1,-.1,num2str(nump)); end hold off end 10a b=[1 -1 0]; a=[1 3 2]; dpzplot(b,a); 运行截图： 10b经过Z变换之后可以得到： b=[1 0 0]; a=[1 1 0.5]; dpzplot(b,a); 运行截图： 10c经过Z变换之后可以得到 b=[1 0.5 0 0]; a=[1 -1.25 0.75 -0.125]; dpzplot(b,a); 运行截图： 实验总结 这次实验较以往的实验都容易，基本上就是按照课本上的程序来重新输入一遍就几乎能得到结果了，所以能很快完成。但是还需要大家对拉普拉斯变换和z变换有一定的了解，而且需要知道roots函数的用法。 致谢 转眼间，三年时|’日J已过去，看着眼前这厚厚的一叠文稿，想想这一路走来，实 在接受了太多人的关心与帮助，心中充满感激。希望拙文的的完成，能为我这三年 的研究生学习生涯画上一个圆满的句号。在此，向多年来在求学路上关心我、支持 我的老师、亲人、同窗，表达我最诚致的谢意。 感谢我尊敬的导师给予我在学习上、生活上无微不至的指导与关心。我的每一点进步、每一份 收获，无不包含着黄老师的谆谆教诲和殷殷期望。感谢黄老师，再一次表达我最真挚 的谢意和敬意! 感谢各位老师的指导和帮助。在三年的成长过程中，每一位老师都给了 我无价的知识和真诚的关爱，让我从懵懂无知的青年逐渐成熟睿智，你们的教诲将 成为我一生宝贵的精神财富财富，受益终尘。 感谢与我共同成长的同学们，是你们的情谊与帮助，令我的生活充满温暖和快乐。 感谢生我养我父母。我所迈出的第一步，都凝聚着你们的心血和汗水。谁言寸 草心，报得三春晖，我惟有尽我的努力让你们能为我骄