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

实验项目名称：LTI系统复频域分析的MATLAB实现上机实验题目：拉氏变换与Z变换的基本性质在系统分析中的应用实验项目的目的和任务：掌握拉氏变换、Z变换的基本性质及其在系统分析中的典型应用实验题目：第9章第10章实验过程9a_1b=1,5 0;a=1,2,3;zs=roots(b)ps=roots(a)plot(real(zs),imag(zs),o);hold onplot(real(ps),imag(ps),x);gridaxis(-5 2 -2 2);运行截图：9a_2b=2,5,12;a=1,2,10;zs=roots(b)ps=roots(a)plot(real(zs),imag(zs),o);hold onplot(real(ps),imag(ps),x);gridaxis(-10 10 -10 10);运行截图：9a_3b=2,5,12;a=1,4,14,20;zs=roots(b)ps=roots(a)plot(real(zs),imag(zs),o);hold onplot(real(ps),imag(ps),x);gridaxis(-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 onplot(real(ps),imag(ps),x);gridaxis(-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);endps=roots(a);zs=roots(b);mx=max(abs(ps zs .95)+0.05;clgaxis(-mx mx -mx mx);axis(equal);hold onw=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);endhold offend10ab=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函数的用法。