MATLAB練習(xí)題和答案

限制系統(tǒng)仿真試驗(yàn)

Matlab部分試驗(yàn)結(jié)果

目錄

試驗(yàn)一MATLAB基本操作......................................................1

試驗(yàn)二Matlab編程.........................................................5

試驗(yàn)三Matlab底層圖形限制.....................................................6

試驗(yàn)四限制系統(tǒng)古典分析.......................................................12

試驗(yàn)五限制系統(tǒng)現(xiàn)代分析........................................................15

試驗(yàn)六PID限制器的設(shè)計(jì)......................................................19

試驗(yàn)七系統(tǒng)狀態(tài)空間設(shè)計(jì).......................................................23

試驗(yàn)九直流雙閉環(huán)調(diào)速系統(tǒng)仿真.................................................25

試驗(yàn)一MATLAB基本操作

1用MATLAB可以識(shí)別的格式輸入下面兩個(gè)矩陣

1233

-1+4/43678

2357

233554+2z

A=1357B=

26+7z5342

3239

189543

1894-

再求出它們的乘積矩陣C,并將C矩陣的右下角2X3子矩陣賦給D矩陣。賦值完成后，調(diào)

用相應(yīng)的叮囑查看MATLAB工作空間的占用狀況。

A=[l,2,3,3;2,3,5,7;l,3,5,7;3,2,3,9;l,8,9,4];

B=[l+4iz4z3z6,7,8;2,3,3,5,5,4+2i;2,6+7i,5,3,4,2;l,8,9,5,4,31;

C=A*B;

D=C(4:5Z4:6);

whos

NameSizeBytesClassAttributes

A5x4160double

B4x6384doublecomplex

C5x6480doublecomplex

D2x396doublecomplex

2選擇合適的步距繪制出下面的圖形

sin(l〃)，其中re(-1,1)

t=[-l：0.1：l]；

y=sin(l./t);

plot(t,y)

3對(duì)下面給出的各個(gè)矩陣求取矩陣的行列式、秩、特征多項(xiàng)式、范數(shù)、特征根、特征向量和

逆矩陣。

-7.53.500--5765-

8334.1071087

A=,B=

09103--1.568109

_003.719.3_57910

」23丁-3-3-24-

56785-518

c=,D=

91011121185-7

131415165-1-3-1

A=[7.5,3.5,0,0;8,33,4.1,0;0,9,103,-1.5;0,0,3.7,19.3];

B=[5,7,6,5;7,10,8,7;6,8,10,9;5,7,9,10];

C=[l:4;5:8;9:12;13:lrtf6];

det(A);det(B);det(C);det(D);

rank(A);

rank(B);

rank(C);

rank(D);

a=poly(A);

b=poly(B);

c=poly(C);

d=poly(D);

norm(A);

norm(B);

norm(C);

norm(D);

[v,d]=eig(Aznobalance');

[vzd]=eig(B,nobalance*);

1

[vzd]=eig(C,nobalance);

[v,d]=eig(D,nobalance1);

m=inv(A);

n=inv(B);

p=inv(C);

q=inv(D);

4求解下面的線性代數(shù)方程，并驗(yàn)證得出的解真正滿足原方程。

721-241321390

9153-27721-264

(a)X=，(b)X=

115-19153-2117

132130-2-2115-2-1

(a)

A=[7,2,1,-2;9,15,3,-2;-2,-2,11,5;1,3,2,13];

B=[4;7;-1;0];

X=A\B;

C=A*X;

(b)

A=[l,3,2,13;7,2,1,-2;9,15,3,-2;-2,-2,11,5];

B=[9,0;6,4;ll,7;-2,-l];

X=A\B;

C=A*X;

5.

(1)初始化一10*10矩陣，其元素均為1

ones(10,10);

(2)初始化一10*10矩陣，其元素均為0

zeros(10,10);

(3)初始化-10*10對(duì)角矩陣

v=[l:10];

diag(v);

(4)輸入A=/7/5；256；315],B=[lI1;222;333],執(zhí)行下列叮囑，理解其含義

A(2,3)表示取A矩陣第2行、第3列的元素；

A(:,2)表示取A矩陣的第2列全部元素；

A(3,:)表示取A矩陣第3行的全部元素；

A(:,1:2:3)表示取A矩陣第1、3列的全部元素；

A(:,3).*B(:,2)表示A矩陣第3列的元素點(diǎn)乘B矩陣第2列的元素

A(:,3)*B(2,:)表示A矩陣第3列的元素乘以B矩陣第2行

A*B矩陣AB相乘

A.*B矩陣A點(diǎn)乘矩陣B

AA2矩陣A的平方

A32矩陣表示求矩陣A的每一個(gè)元素的平方值

B/A表示方程AX=B的解X

B./A表示矩陣B的每一個(gè)元素點(diǎn)除矩陣A的元素

6在同一坐標(biāo)系中繪制余弦曲線y=cos(t-0.25)和正弦曲線y=sin(t-O.5),te[0,2Jt],

用不同顏色，不同線的類型予以表示，留意坐標(biāo)軸的比例限制。

t=[0:0.01:2*pi];

yl=cos(t-0.25);

plot(t,yl,1r-1)

holdon

y2=sin(t-0.5);

1r

plot(try2rk)

試驗(yàn)二Matlab編程

1分別用for和while循環(huán)結(jié)構(gòu)編寫程序，求出

63

==1+2+22+23++262+263

/?=0

并考慮一種避開循環(huán)的簡潔方法來進(jìn)行求和。

(a)j=l;n=0;sum=l;

forn=n+l：63

fori=l:n

j=j*2;

end

sum=sum+j;

j=l;

end

sum

(b)j=l;n=l;sum=l;

whilen?=64

i=l；

whilei<n+l

j=j*2；

i=i+l;

end

n=n+l;

sum=sum+j;

j=l;

end

Sum

(c)i=0:63;k=sum(2.Ai);

2計(jì)算l+2+-+n<2000時(shí)的最大n值

s=0;m=0;while(s<=2000),m=m+l;s=s+m;end,m

3用MATLAB語言實(shí)現(xiàn)下面的分段函數(shù)

h,x>D

y=f(x)=<hlDx,|x|<D

—h,x<—D

存放于文件ff.m中，令D=3,h=l求出，如1.5),40.5),45).

D=3;h=l;

x=-2*D:l/2:2*D;

y=-h*(x<-D)+h/D./x.*((x>=-D)&(x<=D))+h*(x>D);

plot(x,y);

gridon

fl=y(find(x==-1.5))

f2=y(find(x==0.5))

f3=y(find(x==5))

試驗(yàn)三Matlab底層圖形限制

1在MATLAB叮囑行中編程得到y(tǒng)=sin(t)和yl=cos⑴函數(shù),plot(t,y);自gure(10);plot(t,yl);

?t=[-pi:O.O5:pi];

?y=sin(t);

?yl=cos(t);

?plot(t,y)

?figure(lO);

?plot(t,yl)

2在MATLAB叮囑行中鍵入h=get(O),查看根屏幕的屬性，h此時(shí)為根屏幕句柄的符號(hào)表

示，0為根屏幕對(duì)應(yīng)的標(biāo)號(hào)。

?h=get(O)

BeingDeleted:'off

BusyAction:'queue'

ButtonDownFcn:"

Callbackobject:[]

Children:[2x1double]

Clipping:'on'

CommandWindowSize:[8927]

CreateFcn:"

CurrentFigure:1

DeleteFcn:"

Diary:'off

DiaryFile:'diary*

Echo:'off

FixedWidthFontName:'CourierNew'

Format:'short'

FormatSpacing:loose1

HandleVisibility:'on'

HitTest:'on'

Interruptible:'on'

Language:'zh_cn.gbk'

MonitorPositions:[1114409001

More:'off

Parent:[]

PointerLocation:[1048463]

PointerWindow:0

RecursionLimit:50()

ScreenDepth:32

ScreenPixelsPerlnch:96

ScreenSize:[111440900]

Selected:'off

SelectionHighlight:'on'

ShowHiddenHandles:'off

Tag:"

Type:‘root'

UIContextMenu:[]

Units:'pixels'

UserData:[]

Visible:'on'

3hl=get(l);h2=get(I0),1,10分別為兩圖形窗口對(duì)應(yīng)標(biāo)號(hào)，其中1為Matlab自動(dòng)支配，標(biāo)號(hào)

10已在figure。0)中指定。查看hl和h2屬性，留意CurrentAxes和CurrenObject屬性。

?hl=get(l)

hl=

Alphamap:[1x64double]

BeingDeleted:'off

BusyAction:'queue*

ButtonDownFcn:"

Children:170.0012

Clipping:'on'

CloseRequestFcn:'closereq'

Color:[0.800()0.80000.8000J

Colormap:[64x3double]

CreateFcn:"

CurrentAxes:170.0012

CunentCharacter:"

Currentobject:[]

CunentPoint:[00]

DeleteFcn:"

DockControls:'on,

FileName:"

FixedColors:[10x3double]

Handlevisibility:'on,

HitTest:'on'

IntegerHandle:'on'

Interruptible:*on*

InvertHardcopy:'on'

KeyPressFcn:"

KeyReleaseFcn:"

MenuBar:'figure'

MinColormap:64

Name:"

NextPlot:'add'

NumberTitle:'on'

PaperOrientation:'portrait'

PaperPosition:[0.63456.345220.304615.2284]

PaperPositionMode:'manual1

PaperSize:[20.984029.6774]

PaperType:'A4'

PaperUnits:'centimeters'

Parent:0

Pointer:'arrow*

PointerShapeCData:[16x16double]

PointerShapeHotSpot:[11]

Position:[440378560420]

Renderer:'painters'

RendererMode:'auto'

Resize:'on'

ResizeFcn:"

Selected:'off

SelectionHighlight:'on'

SelectionType:'normal1

Tag:"

ToolBar:'auto,

Type:'figure'

UIContextMenu:[]

Units:'pixels'

UserData:f]

Visible:'on'

WindowButtonDownFcn:"

WindowButtonMotionFcn:"

WindowButtonUpFcn:"

WindowKeyPressFcn:"

WindowKeyReleaseFcn:"

WindowScrollWheelFcn:"

WindowStyle:'normal'

WVisual:*00(RGB32GDI,Bitmap,Window),

WVisualMode:'auto'

?h2=get(10)

h2

Alphamap:[1x64double]

BeingDeleted:*off

BusyAction:'queue1

ButtonDownFcn:"

Children:342,0011

Clipping:'on'

CloseRequestFcn:'closereq'

Color:[0.80000.80000.8000]

Colormap:[64x3double]

CreateFcn:"

CurrentAxes:342.0011

Currentcharacter:"

Currentobject:[1

CurrentPoint:[00]

DeleteFcn:"

DockControls:'on,

FileName:"

FixedColors:[10x3double]

Handlevisibility:'on'

HitTest:'on,

IntegerHandle:'on'

Interruptible:'on'

InvertHardcopy:'on'

KeyPressFcn:"

KeyReleaseFcn:"

MenuBar:'figure'

MinColormap:64

Name:"

NextPlot:'add'

NumberTitle:'on'

PaperOrientation:'portrait'

PaperPosition:[0.63456.345220.304615.2284]

PaperPositionMode:'manual*

PaperSize:[20.984029.6774]

PaperType:'A4'

PaperUnits:'centimeters'

Parent:0

Pointe亡'arrow'

PointerShapeCData:[16x16double]

PointerShapeHotSpot:[11]

Position:[440378560420]

Renderer:'painters*

RendererMode:'auto'

Resize:‘on'

ResizeFcn:

Selected:'off

SelectionHighlight:'on'

SelectionType:'normar

Tag:"

ToolBar:'auto'

Type:'figure'

UIContextMenu:[]

Units:'pixels'

UserData:f]

Visible:'on'

WindowButtonDownFcn:"

WindowButtonMotionFcn:"

WindowButtonUpFcn:"

WindowKeyPressFcn:"

WindowKeyReleaseFcn:"

WindowScrollWheelFcn:"

WindowStyle:'normal'

WVisual:'00(RGB32GDI,Bitmap,Window)1

WVisualMode:'auto'

4輸入h.Children,視察結(jié)果。

?h.Children

ans=

1

10

5鍵入gcf,得到當(dāng)前圖像句柄的值，分析其結(jié)果與h,hl,h2中哪個(gè)一樣，為什么？

ans=

I

結(jié)果與h的一樣

6鼠標(biāo)點(diǎn)擊Figure1窗口，讓其位于前端，在叮囑行中鍵入gcf,視察此時(shí)的值，和上一步

中有何不同，為什么？

ans=

1

7視察hl.Children和h2.Children,gca的值。

?hl.Children

ans=

170.0012

?h2.Children

ans=

342.0011

?gca

ans=

170.0012

8視察以下程序結(jié)果h3=hl.Children;set(h3,'Color*/green');h3_1=get(h3,'children,);set(h3_l,

,Color1,Yed*);其中h3_l為Figure1中線對(duì)象句柄，不能干脆接受h3_l=h3.Children叼囑獲

得。

9叮囑行中鍵入plot(t,sin(t-pi/3)),視察曲線出現(xiàn)在哪個(gè)窗口。h4=h2.Children;axes(h4);

plot(t,sin(t-pi/3)),看看此時(shí)曲線顯示在何窗口。

pk>t(t,sin(t-pi/3))后,曲線出現(xiàn)在figure1窗口;h4=h2.Children;axes(h4);plot(t,sin(t-pi/3))后,

曲線出現(xiàn)在figurelO

試驗(yàn)四限制系統(tǒng)古典分析

10

3.已知二階系統(tǒng)G(s)=

52+25+10

(1)編寫程序求解系統(tǒng)的階躍響應(yīng);

a=sqrt(10);

zeta=(1/a);

num=[10];

den=[12*zeta*a10];

sys=tf(num,den);

t=0:0.01:3;

figure(1)

step(sysft);grid

修改參數(shù)，實(shí)現(xiàn)J=1和J=2的階躍響應(yīng);

"1時(shí)：

a=sqrt(10);

zeta=l;

num=[10];

den=[12*zeta*a10];

sys=tf(num,den);

t=0:0.01:3;

figure(1)

step(sysrt);grid

4=2時(shí)：

a=sqrt(10);

zeta=2;

num=[10];

den=[12*zeta*a10];

sys=tf(num,den);

t=0:0.01:3;

figure(1)

step(sys,t);grid

修改參數(shù)，實(shí)現(xiàn)/〃1=和口〃2=的階躍響應(yīng)(G“=JI6)

時(shí)：

a=sqrt(10);

zeta=(1/a);

num=[0.25];

den=[12*zeta*0.5*a0.25];

sys=tf(num,den);

t=0:0.01:3;

figure(1)

step(sys,t);grid

%2=2%時(shí)：

a=sqrt(10);

zeta=(1/a);

num=[40];

den=[12*zeta*2*a40];

sys=tf(num,den);

t=0:0.01:3;

figure(1)

step(sys,t);grid

(2)試做出以下系統(tǒng)的階躍響應(yīng)，并比較與原系統(tǒng)響應(yīng)曲線的差別與特點(diǎn)，作出相應(yīng)的試驗(yàn)

分析結(jié)果。

2s+10/+0.5s+10s2+0.5s

G(s)=G(5)=；G.3(S)=

+2s+102s2+25+1052+2s+10

G(5)=

252+25+10

要求：分析系統(tǒng)的阻尼比和無阻尼振蕩頻率對(duì)系統(tǒng)階躍響應(yīng)的影響;

分析響應(yīng)曲線的零初值、非零初值與系統(tǒng)模型的關(guān)系：

分析響應(yīng)曲線的穩(wěn)態(tài)值與系統(tǒng)模型的關(guān)系；

分析系統(tǒng)零點(diǎn)對(duì)階躍響應(yīng)的影響；

a=sqrt(10);

zeta=(1/a);

num=[10];

den=[12*zeta*a10];

sys=tf(num,den);

t=0:0.01:3;

step(sys,t);

holdon

numl=[0210];

sysl=tf(numl,den);

step(sysl,t);

num2=[10.510];

sys2=tf(num2zden);

step(sys2zt);

num3=[10.50J;

sys3=tf(num3,den);

step(sys3,t);

num4=[010];

sys4=tf(num4,den);

step(sys4zt);grid

?ZjFigure1□00

FileEditVie?InsertToolsDeskiopWindovHelp

口目金)玲，，門@*房?③口國?國

StepResponse

System:sys1

Time(sec):0.646

Amplitude:1.35

System:sys

Time(sec):0.827

■Afriplifude:T26--

3

.8

pg6

snsSystem:sys2

dTime(sec):0.731

與4AmnOhMle*A80S

sSystem:sys4

2Time(sec):0.676

Amplitude:0.152

o.o

System:sys3

-0.2-------Time(sec):0.904

Amplitude:-0.396

11.5

Time(sec)

5

k(s+1)

已知G(s)=

52(0.15+1)

令k=l作Bode圖，應(yīng)用頻域穩(wěn)定判據(jù)確定系統(tǒng)的穩(wěn)定性，并確定使系統(tǒng)獲得最大相位

裕度的增益k值。

G=tf([11],[0.1100]);

figure(1)

margin(G);grid

Figure1

FileEditViewInsertToolsDesktopWindowHelp

q通|巧?要/，|曰|口目Ind

BodeDiagram

Gm=-InfdB(at0rad/sec),Pm-44.5deg(at1.26rad/sec)

m

p

)

①

p

n

v

u

o

e

w

-150

Frequency(rad/sec)

試驗(yàn)五限制系統(tǒng)現(xiàn)代分析

(2)Bode圖法推斷系統(tǒng)穩(wěn)定性：

已知兩個(gè)單位負(fù)反饋系統(tǒng)的開環(huán)傳遞函數(shù)分別為：

2727

G(s)=3…“；G,(s)=3一，

53+5s-+4s's+5s--4s

用Bode圖法推斷系統(tǒng)閉環(huán)的穩(wěn)定性。

Gl=tf([2.7],[1540]);

figure(1)

margin(G);grid

G2=tf([2.7]z[15-40]);

figure(2)

margin(G2);grid

igure1

FileEditViewInsertToolsDesktopWindowHelp

□ddQ+、、爭⑨0/□SI■O

BodeDiagram

Gm=-InfdB(at0rad/sec),Pm=44.5deg(at1.26rad/sec)

100

m50

p

)

g

p

n

u0

g

e

w

-100

-120

(

aB

)p

8

s

e

q

d

Frequency(rad/sec)

Figure2［二］叵|區(qū)

FileEditVie*InsertToolsDesktopWindowHelp

□t3J4+、，C?編/?學(xué)□IDl?國

BodeDiagram

Gm=Inf,Pm=-58.1deg(at0.535rad/sec)

(

B

p

)

0p>

n

w

罩

se

ie

di

Frequency(rad/sec)

2

系統(tǒng)能控性、能觀性分析

已知連續(xù)系統(tǒng)的傳遞函數(shù)模型:

s+a

G(s)=

?+10J2+275+18

當(dāng)a分別取一1,0,+1時(shí)，判別系統(tǒng)的能控性與能觀性。

當(dāng)Q取-1時(shí)：

num=[1-1];

den=[l102718];

G=tf(num,den);

Gl=ss(G)

a=[-10-3.375-2.25;800;010];

b=[0.5;0;0];

A

Uc=[bza*b,a2*b];

rank(Uc)

a=

xlx2x3

X1-10-3.375-2.25

x2800

x3010

b=

ul

xl0.5

x20

x30

c=

xlx2x3

yl00.25-0.25

d=

ul

yl0

Continuous-timemodel.

rankUc=

3

rankUo=

3

由此，可以得到系統(tǒng)能控性矩陣Uc的秩是3,等于系統(tǒng)的維數(shù)，故系統(tǒng)是能控的。能觀性矩

陣U。的秩是3,等于系統(tǒng)的維數(shù)，故系統(tǒng)能觀測(cè)的。

當(dāng)a取0時(shí)：

a

xlx2x3

xl-10-3.375-2.25

x2800

x3010

ul

xl0.25

x20

x30

x1x2x3

yl00.50

ul

yl0

Continuous-timemodel.

rankUc=

3

rankUo=

3

由此，可以得到系統(tǒng)能控性矩陣Uc的秩是3,等于系統(tǒng)的維數(shù)，故系統(tǒng)是能控的。能觀性矩

陣U。的秩是3,等于系統(tǒng)的維數(shù)，故系統(tǒng)能觀測(cè)的。

當(dāng)a取1時(shí):

a=

xlx2x3

xl-10-3.375-2.25

x2800

x3010

b=

ul

xl0.5

x20

x30

c

xlx2x3

yl00.250.25

d=

ul

yl0

rankUc=

3

rankUo=

2

由此，可以得到系統(tǒng)能控性矩陣Uc的秩是3,等于系統(tǒng)的維數(shù)，故系統(tǒng)是能控的。能觀性矩

陣U。的秩是2,小于系統(tǒng)的維數(shù)，故系統(tǒng)不能觀測(cè)的。

試驗(yàn)六PID限制器的設(shè)計(jì)

1.已知三階對(duì)象模型G(s)=l/(s+l)3,利用MATLAB編寫程序，探討閉環(huán)系統(tǒng)在不同限

制狀況下的階躍響應(yīng)，并分析結(jié)果。

(1)8,%->0時(shí)，在不同KP值下，閉環(huán)系統(tǒng)的階躍響應(yīng)；

S=tf(*S1);

G=l/(s+1)A3;

forK=[0:0.5:2];holdon

step(feedback(G*K,1))

end

I3回區(qū)

E.xleEditView工nsertToolsDesktopWindowHelp

□aj\-xo?c口園-回

StepResponse

1.4

1.2

8

o

p

n

M

d

lo.

u

vo.6

比例環(huán)節(jié):成比例地反映限制系統(tǒng)的偏差信號(hào)e⑴，偏差一旦產(chǎn)生限制器立刻產(chǎn)生限制作用,

以減小偏差。比例系數(shù)越大，誤差越小。

(2)K〃=l,7；->0時(shí)，在不同7；值下，閉環(huán)系統(tǒng)的階躍響應(yīng)；

s=tf(*s1);

G=l/(s+1)A3;

forTi=[1:0.5:3];holdon

Gl=l+tf(1,[Tiz0]);

step(feedback(G*G1Z1))

end

13間國

FileEditViewInsertToolsDesktoplindowH?lp

□aJ口國■□

StepResponse

1.6

1.4

1.2

o

p

n

s8

d

uso.

/

v

6

51015202530354045

Time(sec)

積分環(huán)節(jié).?主要用于消退靜塞提高系統(tǒng)的無差度。積分作用的強(qiáng)弱取決于積分時(shí)間常數(shù)"

力越大，積分作用越弱,反之則越強(qiáng)。

(3)Kp=7；=l時(shí)，在不同7；值下，閉環(huán)系統(tǒng)的階躍響應(yīng)；

s=tf(*s*);

G=l/(s+1)A3;

forTd=[1:0.5:3];holdon

Gl=l+tf(1,[1,0])+tf([Td0],[Td/101]);

step(feedback(G*G1,1))

end

微分環(huán)節(jié):反映偏差信號(hào)的變更趨勢(shì)（變更速率），并能在偏差信號(hào)變得太大之前，在系統(tǒng)中引

入一個(gè)有效的早期修正信號(hào)，從而加快系統(tǒng)的動(dòng)作速度，削減調(diào)整時(shí)間。

2.被控對(duì)象同上，選擇合適的參數(shù)進(jìn)行模擬PID限制(PID參數(shù)整定)

S=tf(*S1);

G=l/(s+1)A3;

[KIzLIrT1rG1]=getfolpd(1FG);G1

N=10;K=1.014;T=2.101;L=l.l;a=K*L/T;

Kp=l.2/arTi=2*L,Td=0.5*LrG2=Kp*(1+tf(1,[Ti,0])+tf([Td0],[Td/N

1]))；

step(sysl=feedback(G*G2,1));grid

holdon

step(feedback(G,1))

Figure1□00

FileEditViewInsertToolsDesktopWindowHelp

QT百Hfe+、4c⑨組/，曼□0：TQ

S8

o

np

dw

/u

vo.6

0.2

1012

Time(sec)

試驗(yàn)七系統(tǒng)狀態(tài)空間設(shè)計(jì)

3.已知對(duì)象模型

--0.30.1-0.05112]

x(t)-10.10x(Z)+0

—1.5—8.9-0.054

y=[l23]x

⑴假如我們想將閉環(huán)系統(tǒng)的極點(diǎn)配置到-1,2-3,利用MATLAB設(shè)計(jì)限制器，并繪出閉環(huán)

系統(tǒng)的階躍響應(yīng)曲線。(說明：用兩種方法配置極點(diǎn))

接受Ackermann公式計(jì)算：

A=[-0.30.1-0.05;l0.10;-1.5-8.9-0.05];

B=[2;0;4];

C=[l23];

D=0;

P=[-l-2-3];

k=acker(A,B,P)

Ac=A-B*k

eig(Ac)

G=ss(Ac,B,C,D)

Gl=tf(G)

t=0:0.01:3;

step(G1tt);grid

Figure1;□iX

FileEditViewInsertToolsDesktopWindowHelp

□后UJ玲+、入踮黝0/，顯□國■Q

Time(sec)

接受魯棒極點(diǎn)配置算法：

A=[-0.30.1-0.05;l0.10;-l.5-8.9-0.05];

B=[2;0;4];

C=[l23];

D=0;

P=[-l-2-3];

k=place(A,B,P)

Ac=A-B*k

eig(Ac)1

G=ss(Ac,B,C,D)

Gl=tf(G)

t=0:0.01:3;

step(G1ft);grid

?JFigure113回區(qū)

FileEditVie?InsertToolsDesktopWindowHelp

StepResponse

2