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

《LTI系统复频域分析的MATLAB实现》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验项目名称：LTI系统复频域分析的MATLAB实现上机实验题目：拉氏变换与Z变换的基本性质在系统分析中的应用实验项目的目的和任务：掌握拉氏变换、Z变换的基本性质及其在系统分析中的典型应用实验题目：第9章第10章实验过程9a_1b=[1,50];a=[1,2,3];zs=roots(b)ps=roots(a)plot(real(zs),imag(zs),'o');holdonplot(real(ps),imag(ps),'x');gridaxis([-52-22]);运行截图：9a_2b=[2,5,12];a=[1,2,10];zs=roots(b)ps

2、=roots(a)plot(real(zs),imag(zs),'o');holdonplot(real(ps),imag(ps),'x');gridaxis([-1010-1010]);运行截图：9a_3b=[2,5,12];a=[1,4,14,20];zs=roots(b)ps=roots(a)plot(real(zs),imag(zs),'o');holdonplot(real(ps),imag(ps),'x');gridaxis([-1010-1010]);运行截图：9c、经过拉氏变换之后得到的系统函数为：b=[1,2,5];a=[1,-3];z

3、s=roots(b)ps=roots(a)plot(real(zs),imag(zs),'o');holdonplot(real(ps),imag(ps),'x');gridaxis([-1010-1010]);运行截图：10、dpzplot函数代码：functiondpzplot(b,a)la=length(a);lb=length(b);if(la>lb)b=[bzeros(1,la-lb)];elseif(lb>la)a=[azeros(1,lb-la)];endps=roots(a);zs=roots(b);mx=max(abs([ps'zs'.

4、95]))+0.05;clgaxis([-mxmx-mxmx]);axis('equal');holdonw=[0:0.01:2*pi];plot(cos(w),sin(w),'.');plot([-mxmx],[00]);plot([00],[-mxmx]);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);i

5、fnumz>1text(-.1,-.1,num2str(numz));elseifnump>1text(-.1,-.1,num2str(nump));endholdoffend10ab=[1-10];a=[132];dpzplot(b,a);运行截图：10b经过Z变换之后可以得到：b=[100];a=[110.5];dpzplot(b,a);运行截图：10c经过Z变换之后可以得到b=[10.500];a=[1-1.250.75-0.125];dpzplot(b,a);运行截图：实验总结这次实验较以往的实验都容易，基本上就是按照课本上的程序来重新输入一遍就

6、几乎能得到结果了，所以能很快完成。但是还需要大家对拉普拉斯变换和z变换有一定的了解，而且需要知道roots函数的用法。