t=0:0.1:30;y=1-2*exp(-t).*sin(t);plot(t,y);('Time');('Amplitude');title('Decaying-oscillatingExponential');FOR語句：S=0;form=1:100ifS>=1000break;endS=S-m^2m=m-1

每個(gè)環(huán)節(jié)都對(duì)應(yīng)一個(gè)

WHILE語句：S=0;m=0;whileS<1000endS=S-m^2m=m-1K=6;Z=[-3];P=[-1;-2;-5];零點(diǎn)模傳函數(shù)模狀態(tài)空間模型tf(num,den)[A,B,C,D]=zp2ss(Z,P,K)ss(A,B,C,D)num1=[10];den1=[11];num2=[1];den2=[20.5];num3=[540];den3=[1];num4=[0.1];den4=[1];[na,da]=series(num1,den1,num2,den2);[nb,db]=feedback(na,da,num4,den4);cloop:單位饋系統(tǒng)閉函num2={11;10};den2={[13],[11];[11],1};G2=tf(num2,den2)

陣分母個(gè)矩陣遞數(shù)陣聯(lián)聯(lián)num=7*[23];den=conv([10],conv([10],conv([31],conv([12],conv([12],[508])))));tf(num,den)numo=[20];deno=[1836400];單位反饋系統(tǒng)閉環(huán)函數(shù)plot(t,y);

單位階躍響應(yīng)曲線，為仿真時(shí)disp(['超調(diào)量p=',num2str(Mp),出函數(shù)

數(shù)值轉(zhuǎn)換字符串wn=6;zeta=[0.2:0.2:1.0,2.0];figure(1);holdfori=zeta

hold在當(dāng)前圖的軸（坐標(biāo)）中畫了一幅圖，再畫一幅圖時(shí)，原來圖還在圖片注釋，標(biāo)題endnumo=[1];deno=conv([10],conv([0.11],[0.23]));單位反系統(tǒng)閉函數(shù)t=0:0.1:30;r=t;[y,x]=lsim(num,den,u,t)u為系統(tǒng)輸入信號(hào)t為仿時(shí)間針對(duì)線性時(shí)不變模型給定任意輸入得到輸出響應(yīng)num=[15];den=conv([10],conv([12],[13]));K=rlocfind(num,den)%移動(dòng)鼠標(biāo)到根軌跡與實(shí)軸的交點(diǎn)，得到增益Rlocus格式：[r,k]=rlocus(n,d)[r,k]=rlocus(g)[r,k]=rlocus(n,d,k)函數(shù)可算出或畫出系統(tǒng)根軌，中g(shù)(或?yàn)閷?duì)象模型，輸入變量為用戶已選擇增益向，當(dāng)k缺省時(shí)則為系統(tǒng)動(dòng)生成增向量,返回變量根軌跡各個(gè)點(diǎn)構(gòu)成的復(fù)矩陣.如果在函數(shù)調(diào)用中不返回任何參數(shù)，則rlocus函數(shù)在當(dāng)前口中畫出系統(tǒng)的根軌圖。numo=[12];deno=conv([143],[143]);figure(1);rlocus(numo,deno);K0=rlocfind(numo,deno)%移動(dòng)鼠標(biāo)到根軌跡與虛軸的交點(diǎn)，得到增益，使環(huán)系統(tǒng)定取值范圍為：0<K<K0。K1=10;[numc1,denc1]=cloop(K1*numo,deno);impulse是脈沖信號(hào)輸入這個(gè)遞函數(shù)后出來的脈沖應(yīng)函數(shù)圖（，）可使傳遞，也可是pk也可以狀態(tài)函數(shù)t為響應(yīng)時(shí)間gtext('K=10'

根軌圖上點(diǎn)了虛軸上的交之后，出現(xiàn)=10的脈沖響應(yīng)點(diǎn)擊圖中會(huì)出現(xiàn)的標(biāo)題K2=100;[numc2,denc2]=cloop(K2*numo,deno);impulse(numc2,denc2);gtext('K=100'

圖二點(diǎn)擊之會(huì)出現(xiàn)=100脈沖響應(yīng)，點(diǎn)擊圖中會(huì)現(xiàn)=100的標(biāo)題當(dāng)題目需要析穩(wěn)定性的候，點(diǎn)虛軸交點(diǎn)，脈沖響應(yīng)線發(fā)散和收斂的情況各不同。numo=[-12];deno=conv([10],[13]);rlocus(numo,deno);K1=rlocfind(numo,deno)K2=rlocfind(numo,deno)%動(dòng)鼠標(biāo)到根跡與虛軸上交點(diǎn)，得到益，再移動(dòng)鼠標(biāo)到根軌跡實(shí)軸上的分離點(diǎn)，得到增益%閉環(huán)系統(tǒng)響應(yīng)衰減振蕩的取值范圍：K2<K<K1。wn=6;w=logspace(-1,1,100);num=wn^2;figure(1);forζ用來代替按照為橫坐標(biāo)，生成bodem為頻，p為相頻;grid

依次點(diǎn)擊會(huì)出現(xiàn)的曲線endlogspace(a,b,n)，其中a、、分別表示開值（10^a）結(jié)束值(10^b)、元個(gè)數(shù)n)生成從0的次方的次方之間按對(duì)數(shù)等的n元素的行向量。如果省略則默認(rèn)值50。系統(tǒng)穩(wěn)定的充要條件為：線按逆針包圍臨界點(diǎn)的圈N，等于開環(huán)傳函數(shù)位于s右半平面極點(diǎn)數(shù)；否則閉環(huán)系統(tǒng)不穩(wěn)定，閉環(huán)正實(shí)特征根個(gè)數(shù)Z。若剛好過臨界點(diǎn)則系臨界穩(wěn)定num=[0.5];den=[1210.5];[z,p]=tf2zp(num,den)nyquist(num,den);N=0Z=P-NifZ==0'Stable'););end

求右半平面軸的點(diǎn)數(shù)z=[];p=[0-1-10];k=10;[numo,deno]=zp2tf(z,p,k);[Gm,Pm]=margin(numo,deno)

margin函可從率應(yīng)據(jù)計(jì)出：幅裕度，角裕度及對(duì)應(yīng)頻率，值在頻率為時(shí)，頻率響的實(shí)裕度和相角裕度針對(duì)開環(huán)系統(tǒng)而言，指示斯特在頻率為時(shí)，頻響應(yīng)出系統(tǒng)閉環(huán)的相穩(wěn)定性，可在當(dāng)前圖形窗應(yīng)頻率（單位：）口中繪制出帶有度以及相應(yīng)頻率顯示的，頻率響應(yīng)的幅波特圖，其中幅裕度已分貝為單位

值波圖在頻率為時(shí)，頻率響應(yīng)的相位w：響頻率（單位）：值裕度；：角裕度；：角界率；截止頻率；對(duì)于線性時(shí)變系統(tǒng)如果對(duì)取定初始時(shí)刻的一個(gè)非零初始狀態(tài)，存一個(gè)時(shí)t1,(t1>t0)，和一個(gè)無約束的容許控u(t)，∈[t0,t1]，狀態(tài)由轉(zhuǎn)移到時(shí)，則此0t0時(shí)刻為能控的A=[10;-10];B=[1;0];

陣的最大長度能控性if函數(shù)是求矩陣秩個(gè)向量的極大線性無關(guān)組所含向量的個(gè)數(shù)稱為這向量組的秩狀態(tài)能控else系統(tǒng)狀態(tài)能控);end對(duì)于線性時(shí)變系統(tǒng)如果對(duì)取定初始時(shí)刻的一個(gè)非零初始狀態(tài)，存在一個(gè)有限刻，對(duì)所有t∈[t0,t1]，有，則稱此0在時(shí)刻為不能觀測的。A=[010;001;-2-4-3];C=[00-1;11];n=length(A);Vo=obsv(A,C);ifrank(Vo)==n'系統(tǒng)狀態(tài)能觀測);系統(tǒng)狀態(tài)不能觀測);endz=[];p=[0-1-2];k=10;n=length(A);ifrc==n

極點(diǎn)型轉(zhuǎn)狀態(tài)間先判斷能控性系狀能控函為系統(tǒng)的狀態(tài)空間矩陣量中是望的閉環(huán)極位置回值是增益向。end

這邊一定注意上頭引號(hào)不要掉！要轉(zhuǎn)置成一個(gè)列向量??！Tf=20;表左列是，右列是是10，數(shù)為陣列長的矩陣sim(‘剛才建立的模型名稱要和文件放在一個(gè)文件夾里面,行時(shí)間Tf,選項(xiàng)可以不寫ut);plot(t,y)sxG(

s

間den=[14100];kesi=den(2)/2/wn;kesi<1此系二阻統(tǒng)（分）N=length(t);ss=y

（分）ss)/yss<=0.05i=i-1;

uu

xux

yy

u

G()

yy

Ts=t(i)

（分）

y

G(

u

以、為入，、、。其中

,

,(s

s

,G(s

ss1=[1910;3127911;5121314];6;24;22;1021;802eros(2,2);（分）num2=[032;1025];5[A,B,C,D]=append(A1,B1,C1,D1,A2,B2,C2,D2);

（分）Z3=-1;K3=2;[A3,B3,C3,D3]=zp2ss(Z3,P3,K3);[Aa,Ba,Ca,Da]=append(A,B,C,D,A3,B3,C3,D3);%將個(gè)輸入輸出統(tǒng)接組合一多入輸出系，….sysn)Q=[3-5;444];[Ac,Bc,Cc,Dc]=connect(Aa,Ba,Ca,Da,Q,ins,outs);

（分）已狀：穩(wěn)最。A=[12-12;2630;4-8-5;721D=7;)>0);

x

udisp('Sy

y

xuelsedisp('Syisisaphaseone.');elsedisp('Syisphaseone.');10032x10032x1別利用for和循：ii

2

①循：S=0;i=0:100S已：(s

（分1)(0.25

②循：i=0;i<=100i=i+1;（分S10G()G)0.5求整個(gè)系統(tǒng)的傳num1=[16.71]);num2=[10];den2=[1

G()1

2num3=[0.5];den3=[1];

（）[nb,db]=series(num1,den1,na,da);

(s)3

（分已脈()

3z

z

z

z

要性2-20[numc,denc]=cloop(numo,deno);ii=find(abs(r)>1);disp('Sy(r(ii));

（）（分（4）已狀極j2;。1,2A=[00-4;0116];0;12;0rc=rank(ctrb(A,B))P=[-2+2*jK=place(A,B,P)end

0（分1161（分隨16、4、2整。用-else結(jié)x=input('pleaseinputanumber');ifbedividedbyelseifrem(x,4)==0be

elseifcandividedbeby16,4andend()H()：

(ss(ss)定并說明使系統(tǒng)產(chǎn)生臨界的值繪制此時(shí)的單位脈沖響應(yīng)。num=[1den=poly([