matlab教程第四章數(shù)值計(jì)算-資料下載頁(yè)
2024-08-21 13:33本頁(yè)面
【導(dǎo)讀】，進(jìn)行三種解法比較。其中A取MATLAB庫(kù)中的特殊。簡(jiǎn)單實(shí)陣的特征值問(wèn)題。把例征值對(duì)角陣D轉(zhuǎn)換成實(shí)數(shù)塊對(duì)角陣，使VR*DR/VR=A。對(duì)虧損矩陣進(jìn)行Jordan分解。sssss的“商”及“余”多項(xiàng)式。cq='商多項(xiàng)式為';cr='余多項(xiàng)式為';兩種多項(xiàng)式求值指令的差別。對(duì)于給定數(shù)據(jù)對(duì)x0,y0，求擬合三階多項(xiàng)式，并圖示擬合情況。以上例所給數(shù)據(jù)，研究一維插值，并觀察插值與擬合的區(qū)別。
【正文】 lot(t,x,39。b39。,t,y,39。r.39。,39。MarkerSize39。,10) legend(39。Input39。,39。Output39。,0) 16 0 21012Input Output 圖 10階 Butterworth濾波器的濾波效果 系統(tǒng)分析 線性時(shí)不變對(duì)象 LTI 【例 】已知一個(gè)兩輸入兩輸出系統(tǒng)的狀態(tài)方程四對(duì)組，據(jù)此建立 LTI 對(duì)象“狀態(tài)方程子類”模型。然后再?gòu)拇四Ｐ瞳@取傳遞函數(shù)二對(duì)組。 A=[ 。 。 8]。 B=[0 0 4。1 0 1]39。 C=[ 0。0 ]。 D=0。 S_ss=ss(A,B,C,D) a = x1 x2 x3 x1 x2 x3 8 b = u1 u2 x1 0 1 x2 0 0 x3 4 1 c = x1 x2 x3 y1 0 y2 0 d = u1 u2 y1 0 0 y2 0 0 Continuoustime model. 17 S_tf=tf(minreal(S_ss)) [Num,Den]=tfdata(S_ss) Transfer function from input 1 to output... 4 1: s^3 + 8 s^2 + 5 s + 4 2 s^2 + s 2: s^3 + 8 s^2 + 5 s + 4 Transfer function from input 2 to output... s^2 + s + 1: s^3 + 8 s^2 + 5 s + 4 s^2 + s 2: s^3 + 8 s^2 + 5 s + 4 Num = [1x4 double] [1x4 double] [1x4 double] [1x4 double] Den = [1x4 double] [1x4 double] [1x4 double] [1x4 double] Num{1,1},Den{1,1} ans = 0 0 0 ans = 【例 】求方框圖（圖 ）所描述的那多環(huán)系統(tǒng)的傳遞函數(shù)。已知此方框圖中 G s G s G ss s G ss H ss H H1 2 3 2 4 1 2 31 10 1 1 14 4 16 12 2 1? ? ? ? ? ?? ? ? ?? ? ?? ? ?, , , , , ,。 圖 多環(huán)系統(tǒng) G1=tf(1,[1 10])。G2=tf(1,[1 1])。 G3=tf([1 1],[1 4 4])。G4=tf([1 1],[1 6])。 H1=tf([1 1],[1 2])。H2=2。H3=1。 P1=G3*G4/(1H1*G3*G4)。 P1=minreal(P1)。 18 P2=minreal(G2*P1/(1+G2*P1*H2/G4))。 P3=minreal(G1*P2/(1+G1*P2*H3)) Transfer function: s^2 + 3 s + 2 s^5 + 21 s^4 + 157 s^3 + 564 s^2 + 1004 s + 712 系統(tǒng)的時(shí)域和頻域分析 【例 】繪制一個(gè)零初始狀態(tài)的二輸入二輸出系統(tǒng)各“通道”的單位階躍響應(yīng)（圖 ）和頻率響應(yīng)圖（圖 ）。 clear,A=[ 。 。 8]。 B=[0 0 4。1 0 1]39。C=[ 0。0 ]。D=0。 S_ss=ss(A,B,C,D)。 step(S_ss) Time (sec.)AmplitudeStep Response01234From: U(1)To: Y(1)From: U(2)0 5 10 15 2010To: Y(2)0 5 10 15 20 圖 二輸入二輸出系統(tǒng)的階躍響應(yīng) bode(S_ss) 19 Frequency (rad/sec)Phase (deg)。 Magnitude (dB)Bode Diagrams2000200From: U(1) From: U(2)4002000To: Y(1)500501011001011025000500To: Y(2)101100101102 圖 二輸入二輸出系統(tǒng)的頻率響應(yīng) 習(xí)題 4
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1