【文章內(nèi)容簡介】 1 ? fprintf(fid,39。%39。,T(i,:))。 ? fprintf(fid,39。\r\n39。)。 ? end ? fclose(fid)。 ? disp(39。系統(tǒng)的狀態(tài)空間表達(dá)式 (狀態(tài)空間模型 )： 39。)。 ? G=ss(A1,B1,C1,D1) %求狀態(tài)空間表達(dá)式 ? disp(39。用矩陣符號(hào)表示的狀態(tài)空間模型如下： 39。)。 ? disp(39。dx/dt=Ax+Bu (狀態(tài)方程 )39。)。 ? disp(39。y=Cx+Du (輸出方程 )39。)。 ? disp(39?？煽乇鎰e矩陣 T為： 39。) ? T=ctrb(A1,B1) %計(jì)算可控判別矩陣 ? disp(39?？煽嘏袆e矩陣 T的秩為： 39。) ? RT=rank(T) %計(jì)算可控判別矩陣的秩 ? if RT==3 ? disp(39。系統(tǒng)可控 !39。)。 ? else ? disp(39。系統(tǒng)不可控！ 39。)。 ? end ? disp(39?？捎^辨別矩陣 N為： 39。) ? N=obsv(A1,C1) %計(jì)算可觀判別矩陣 ? disp(39?？捎^判別矩陣 N的秩為： 39。) ? RN=rank(N) %計(jì)算可觀判別矩陣的秩 ? HS2=size(N,1)。 ? fprintf(fid,39。%s39。,39?？捎^判別矩陣 :39。)。 ? fprintf(fid,39。\r\n39。)。 ? for i=1:HS2 ? fprintf(fid,39。%39。,N(i,:))。 ? fprintf(fid,39。\r\n39。)。 ? end ? fclose(fid)。 ? figure(1)。 ? t=0::10。 ? subplot(1,2,1)。 ? impulse(G,t)。 %計(jì)算脈沖響應(yīng)并添加?xùn)鸥? ? grid。 ? title(39。脈沖響應(yīng) 39。)。 ? subplot(1,2,2)。 ? step(G,t)。 ? grid %計(jì)算階躍響應(yīng)并添加?xùn)鸥? ? title(39。階躍響應(yīng) 39。)。 ? y=1。 ? Y=2。 ? n=3。 ? N=4。 ? b=input(39。是否用 simulink仿真運(yùn)行脈沖響應(yīng)和階躍響應(yīng) [y/n]: 39。)。 ? if b==1|b==2 ? open_system(39。untitled839。)。 ? set_param(39。untitled8/Transfer F39。,39。Numerator39。,39。num39。,39。Denominator39。,39。den39。)。 ? save_system ? sim(39。untitled839。)。 ? open_system(39。untitled939。)。 ? set_param(39。untitled9/StateSpace39。,39。A39。,39。A139。,39。B39。,39。B139。,39。C39。,39。C139。,39。D39。,39。D139。)。 ? save_system ? sim(39。untitled939。)。 ? end ? disp(39。你是否還要繼續(xù)使用該傳遞函數(shù)模型， 39。) ? c=input(39。判別可控性和可觀性 [y/n]: 39。)。 ? end ? x=input(39。你是否還想返回菜單，繼續(xù)進(jìn)行別的題的運(yùn)行 [y/n]: 39。)。 ? if x==1|x==2 ? clc ? disp(39。： G(s)=s+a/s^3+10s^2+27s+18,39。)。 ? disp(39。當(dāng) a分別取 1， 0， 1時(shí)，判別系統(tǒng)的能控性和能觀測性。 39。)。 ? disp(39。： 39。)。 ? disp(39。G(s)=3s^4+2s^3+s^2+4s+2/3s^5+5s^4+s^3+2s^2+2s+139。)。 ? disp(39。判斷系統(tǒng)的穩(wěn)定性，并給出不穩(wěn)定極點(diǎn)。 39。)。 ? disp(39。， 39。)。 ? disp(39。G(z)=5z^5+4z^4+z^3+^23z+^539。)。 ? disp(39。判斷單位負(fù)反饋系統(tǒng)的穩(wěn)定性。 39。)。 ? disp(39。： 39。)。 ? disp(39。dx/dt=Ax+Bu39。)。 ? disp(39。y=Cx39。)。 ? disp(39。其中， A=[2 2 1。0 2 0。1 4 0]。B=[0。0。1]。C=[1 0 0]。39。)。 ? disp(39。(1) 利用 MATLAB建立控制系統(tǒng)的數(shù)學(xué)模型； 39。)。 ? disp(39。(2) 分析系統(tǒng)的可控性和可觀性； 39。)。 ? disp(39。(3) 計(jì)算系統(tǒng)的傳遞函數(shù)； 39。)。 ? disp(39。(4) 根據(jù)傳遞函數(shù)，計(jì)算系統(tǒng)的零極點(diǎn)，并判斷系統(tǒng)的穩(wěn)定性； 39。)。 ? disp(39。(5) 利用 LTI Viewer工具繪制系統(tǒng)的階躍響應(yīng)曲線，驗(yàn)證系統(tǒng)的穩(wěn)定性。 39。)。 ? x=input(39。請問您是否還要繼續(xù)做出選擇 ?[y/n] 39。)。 ? end ? end ? %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%第一題%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ? %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%第二題%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ? if cho==2 ? clc ? disp(39。已知閉環(huán)系統(tǒng)的傳遞函數(shù)為： 39。)。 ? disp(39。3s^4+2s^3+s^2+4s+2/3s^5+5s^4+s^3+2s^2+2s+139。)。 ? disp(39。判斷系統(tǒng)的穩(wěn)定性，并給出不穩(wěn)定極點(diǎn)。 39。)。 ? num=[3 2 1 4 2]。 ? den=[3 5 1 2 2 1]。 ? disp(39。系統(tǒng)的傳遞函數(shù)為： 39。)。 ? G=tf(num,den) ? %get(G) ? disp(39。系統(tǒng)的零極點(diǎn)以及增益為： 39。)。 ? [z,p,k]=tf2zp(num,den) %傳遞函數(shù)模型轉(zhuǎn)換為零極點(diǎn)模型 ? hs1=size(p,1)。 ? hs2=size(z,1)。 ? fid=fopen(39。第二題的 txt文檔和截圖 \零極點(diǎn)及其增益 .txt39。,39。w39。)。 ? fprintf(fid,39。%s39。,39。以下是該系統(tǒng)的零極點(diǎn)和增益： 39。)。 ? fprintf(fid,39。\r\n39。)。 ? fprintf(fid,39。%s39。,39。極點(diǎn) p為 :39。)。 ? fprintf(fid,39。\r\n39。)。 ? for i=1:hs1 ? b=p(i,:)。 ? aa=num2str(b)。 ? fprintf(fid,39。%s39。,aa)。 ? fprintf(fid,39。\r\n39。)。 ? end ? fprintf(fid,39。%s39。,39。零點(diǎn) z為： 39。)。 ? fprintf(fid,39。\r\n39。)。 ? for i=1:hs2 ? bb=num2str(z(i,:))。 ? fprintf(fid,39。%s39。,bb)。 ? fprintf(fid,39。\r\n39。)。 ? end ? fprintf(fid,39。%s39。,39。系統(tǒng)增益 K為： 39。)。 ? fprintf(fid,39。%39。,k)。 ? fprintf(fid,39。\r\n39。)。 ? fclose(fid)。 ? disp(39。由以上的零極點(diǎn)得出該系統(tǒng)不穩(wěn)定 39。)。 ? disp(39。系統(tǒng)不穩(wěn)定極點(diǎn)為： 39。)。 ? for n=1:hs1 ? if real(p(n,:))0 ? p(n,:) ? fprintf(39。\r\n39。)。 ? end ? end ? disp(39。系統(tǒng)的零極點(diǎn)傳遞函數(shù)模型為： 39。)。 ? G1=zpk(G) ? figure(1)。 %繪制零極圖與根軌跡圖 ? subplot(2,1,1)。 ? pzmap(G) ? grid on。 ? title(39。零極點(diǎn)分布圖 39。)。 ? subplot(2,1,2)。 ? rlocus(num,den)。 ? title(39。根軌跡分布圖 39。)。 ? saveas(figure(1),39。第二題的 txt文檔和截圖 \零極點(diǎn)與根軌跡分布圖 .bmp39。)。 ? figure(2)。 %繪制階躍響應(yīng)和脈沖響應(yīng) ? subplot(1,2,1)。 ? step(G)。 ? axis([0 10 25 10])。 ? title(39。系統(tǒng)階躍響應(yīng) 39。)。 ? grid on ? subplot(1,2,2)。 ? impulse(G)。 ? axis([0 10 5 5])。 ? title(39。系統(tǒng)脈沖響應(yīng) 39。)。 ? grid on ? saveas(figure(2),39。第二題的 txt文檔和截圖 \系統(tǒng)的階躍響應(yīng)和脈沖響應(yīng) .bmp39。)。 ? figure(3)。 %繪制 Bode圖和 Nyquist圖 ? subplot(1,2,1)。 ? bode(G)。 %Bode圖的繪制 ? grid on ? subplot(1,2,2)。 ? nyquist(G)。 %Nyquist圖的繪制 ? hold on ? grid on ? %line([0,0],[1,1])。 %畫一條起點(diǎn)為（ 0， 1），終點(diǎn)為（ 0， 1）的直線 ? saveas(figure(3),39。第二題的 txt文檔和截圖 \頻域特性曲線 .bmp39。)。 ? %disp(39。顯示裕度及對應(yīng)頻率 : 39。)。 ? %[Gm,Pm,Wcg,Wcp]=margin(G) %顯示裕度及對應(yīng)頻率 ? %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ? y=1。 ? Y=2。 ? n=3。 ? N=4。 ? b=input(39。是否用 simulink仿真運(yùn)行系統(tǒng)的階躍響應(yīng)和脈沖響應(yīng) [y/n]: 39。)。 ? if b==1|b==2 ? open_system(39。untitled139。)。 ? sim(39。untitled139。)。 ? open_system(39。untitled239。)。 ? sim(39。untitled239。)。 ? end ? cb=input(39。是否啟動(dòng) matlab自帶的根軌跡分析與設(shè)計(jì)工具 rltool[y/n]: 39。)。 ? if cb==1|cb==2 ? rltool(G)。 %根軌跡分析與設(shè)計(jì)工具 rltool ? end ? cc=input(39。是否啟動(dòng) matlab自帶的線性時(shí)不變系統(tǒng)的分析工具 LTI Viewer[y/n]: 39。)。 ? if cc==1|cc==2 ? ltiview %matlab線性時(shí)不變系統(tǒng)圖形用戶界面 LTI Viewer ? end ? x=input(39。你是否還想返回菜單，繼續(xù)進(jìn)行別的題的運(yùn)行 [y/n]: 39。)。 ? if x==1|x==2 ? clc ? disp(39。： G(s)=s+a/s^3+10s^2+27s+18,39。)。 ? disp(39。當(dāng) a分別取 1， 0， 1時(shí)，判別系統(tǒng)的能控性和能觀測性。 39。)。 ? disp(39。： 39。)。 ? disp(39。G(