东大matlab实践课三.docx

东大matlab实践课三.docx

《东大matlab实践课三.docx》由会员分享，可在线阅读，更多相关《东大matlab实践课三.docx（18页珍藏版）》请在冰豆网上搜索。

东大matlab实践课三

实验六

一

2

b=[1010];

a=[10];

[z,p,k]=tf2zp（b,a）

z=

-1

p=

0

k=

10

3

b=[10];

a=[1-5];

[z,p,k]=tf2zp（b,a）

z=

Emptymatrix:

0-by-1

p=

5

k=

10

二

1

>>z=0;

>>p=[-1-4];

>>k=1;

>>sys1=zpk（z,p,k）

Zero/pole/gain:

s

-----------

（s+1）（s+4）

2

>>Ns=[1,1];

>>Ds=[1,0,-1];

>>sys2=tf（Ns,Ds）

Transferfunction:

s+1

-------

s^2–1

3

>>Ns=[1,6,6,0];

>>Ds=[1,6,8];

>>sys3=tf（Ns,Ds）

Transferfunction:

s^3+6s^2+6s

-----------------

s^2+6s+8

实验七

三

1

>>z=[1,2];

>>p=[-1,-2];

>>zplane（z,p）

2

num=1;

den=[1,0];

[z,p,k]=tf2zp（num,den）

zplane（z,p）

z=

Emptymatrix:

0-by-1

p=

0

k=

1

3

num=[1,0,1];

den=[1,2,5];

[z,p,k]=tf2zp（num,den）

zplane（z,p）

z=

0+1.0000i

0-1.0000i

p=

-1.0000+2.0000i

-1.0000-2.0000i

k=

1

4

num=[1.8,1.2,1.2,3];

den=[1,3,2,1];

[z,p,k]=tf2zp（num,den）

zplane（z,p）

z=

-1.2284

0.2809+1.1304i

0.2809-1.1304i

p=

-2.3247

-0.3376+0.5623i

-0.3376-0.5623i

k=

1.8000

五

5

>>b=100*conv（[1,0],conv（[1,2],conv（[1,2],conv（[1,3,2],[1,3,2]））））;

>>a=conv（[1,-1],conv（[1,1],conv（[1,3,5,2],conv（[1,0,2,0,4],[1,0,2,0,4]））））;

>>[z,p,k]=tf2zp（b,a）

zplane（z,p）

z=

0

-2.0004+0.0004i

-2.0004-0.0004i

-1.9996+0.0004i

-1.9996-0.0004i

-1.0000+0.0000i

-1.0000-0.0000i

p=

1.0000

0.7071+1.2247i

0.7071-1.2247i

0.7071+1.2247i

0.7071-1.2247i

-1.2267+1.4677i

-1.2267-1.4677i

-0.7071+1.2247i

-0.7071-1.2247i

-0.7071+1.2247i

-0.7071-1.2247i

-1.0000

-0.5466

k=

100

不稳定

6

>>b=conv（[1,-1.414,1],[1,1]）;

>>a=conv（[1,0.9,0.81],[1,-0.3]）;

>>[z,p,k]=tf2zp（b,a）

zplane（z,p）

z=

-1.0000

0.7070+0.7072i

0.7070-0.7072i

p=

-0.4500+0.7794i

-0.4500-0.7794i

0.3000

k=

1

稳定

七

1

>>b=[1];

>>a=[1,3,2];

>>sys=tf（b,a）;

>>nyquist（sys）;

稳定

2

>>b=[1];

>>a=[1,4,4,0];

>>sys=tf（b,a）;

>>nyquist（sys）;

3

>>b=[1];

>>a=[1,2,2];

>>sys=tf（b,a）;

>>nyquist（sys）;

实验八

一

1

>>b=[1,3];

>>a=[1,3,2];

>>sys=tf（b,a）;

>>impulse（sys）;

>>B=[1,3];

>>A=[1,3,2];

>>[a,b,c,d]=tf2ss（B,A）

a=

-3-2

10

b=

1

0

c=

13

d=

0

>>sys=ss（a,b,c,d）;

>>x0=[4,-1];

>>initial（sys,x0）;

>>grid;

2

>>b=[1,3];

>>a=[1,2,2];

>>sys=tf（b,a）;

>>impulse（sys）;

>>B=[1,3];

>>A=[1,2,2];

>>[a,b,c,d]=tf2ss（B,A）;

>>sys=ss（a,b,c,d）;

>>x0=[1.6,-0.2];

>>initial（sys,x0）;

>>grid;

3

>>b=[1,3];

>>a=[1,2,1];

>>sys=tf（b,a）;

>>impulse（sys）;

>>B=[1,3];

>>A=[1,2,1];

>>[a,b,c,d]=tf2ss（B,A）;

>>sys=ss（a,b,c,d）;

>>x0=[7/4,-1/4];

>>initial（sys,x0）;

>>grid;

二

b=1;

a=[1,1,1,0];

sys=tf（b,a）;

subplot（1,2,1）;

impulse（sys）;

subplot（1,2,2）;

step（sys）;

subplot（1,2,1）;

impulse（sys）;

subplot（1,2,2）;

step（sys）;t=0:

0.08:

10

e=sin（t）;

sys=tf（b,a）;

lsim（sys,e,t）;

lsim（sys,e,t）;

三

1

>>clear;

>>b=[1,3];

>>a=[1,3,2];

>>t=0:

0.08:

10;

>>e=[exp（-3.*t）];

>>sys=tf（b,a）;

>>lsim（sys,e,t）;

2

b=[1,3];

a=[1,2,2];

t=0:

0.08:

10;

e=[exp（-0.*t）];

sys=tf（b,a）;

lsim（sys,e,t）;

3

b=[1,3];

a=[1,2,1];

t=0:

0.08:

10;

e=[exp（-2.*t）];

sys=tf（b,a）;

lsim（sys,e,t）;