matlab教程第四章數(shù)值計(jì)算(參考版)
2024-08-25 13:33本頁(yè)面
【正文】 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)。D=0。C=[ 0。 B=[0 0 4。 。 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)圖（圖 ）。 P1=minreal(P1)。H3=1。 H1=tf([1 1],[1 2])。 G3=tf([1 1],[1 4 4])。 圖 多環(huán)系統(tǒng) G1=tf(1,[1 10])。 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ù)。0 ]。1 0 1]39。 8]。 A=[ 。,0) 16 0 21012Input Output 圖 10階 Butterworth濾波器的濾波效果 系統(tǒng)分析 線性時(shí)不變對(duì)象 LTI 【例 】已知一個(gè)兩輸入兩輸出系統(tǒng)的狀態(tài)方程四對(duì)組，據(jù)此建立 LTI 對(duì)象“狀態(tài)方程子類”模型。,39。,10) legend(39。,39。,t,y,39。 plot(t,x,39。 [B,A]=butter(10,30/wn)。 x=sin(2*pi*10*t)+cos(2*pi*100*t)+*randn(size(t))。,1) ws=1000。 clear,randn(39。) 0 5 10 15 200200400600800100012001400Frequency Rad/s 圖 受噪聲污染信號(hào)幅頻譜的局部放大 數(shù) 字濾波 【例 】設(shè)計(jì)一個(gè)低通濾波器，從受噪聲干擾的多頻率混合信號(hào) )(tx 中獲取 10Hz的信號(hào)，見圖 。 plot(w(ii),Ya(ii)) grid,xlabel(39。 Ya=abs(Y(1:length(t)/2))。 Ts=t(2)t(1) Ws=2*pi/Ts。 y=3*sin(5*t)6*cos(9*t)+5*randn(size(t))。state39。在此， )5,0(~)( NtN 。) 14 0 5 10 15 20 25 3002468nncc0 5 10 15 20 25 30101error 圖 變 換法和直接法求卷積結(jié)果比較 【例 】 ff
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1