matlab教程第四章數(shù)值計算(完整版)
【正文】 rei=norm(A*xib)/norm(b) ti = eri = rei = （ 3） tic。state39。 x=ones(100,1)。A(:,1)=[]。 [VR,DR]=cdf2rdf(V,D) VR = 0 DR = 3 【例 】對虧損矩陣進行 Jordan分解。 A_Ap=A.^ A_Mp=A^ A_Ap = A_Mp = + + + + + 4 【例 】標量的數(shù)組乘方和矩陣乘方的比較。 cq=39。)]) 商多項式為 s + 5 余多項式為 5 s^2 + 4 s + 3 【例 】求 3階方陣 A的特征多項式。 P=poly(R) PR=real(P) PPR=poly2str(PR,39。 b=[3,2,5,4,6]。b39。 x0=0::1。,x0,y0,39。（ B） 假設(shè) A 序列中最后 4 個非零值未知，而成為截尾序列，求卷積并圖示。nc1=n1+n3。 nx=nn2+n4。filled39。,0) A=randn(1000,4)。 CCX=Xmx39。 [X,Y]=meshgrid(x,y)。x39。 ss10=dt*sum(y)。cumsum39。r39。cumtrapz39。 clf,plot(x,Y,39。)。 ff=inline(39。 數(shù)值積分 【例 】求 dxeI x? ?? 10 2，其精確值為 ? 。 y0=[1。a([1,2,3])=0。 cc=real(ifft(CF))。),ylabel(39。在此， )5,0(~)( NtN 。 Ya=abs(Y(1:length(t)/2))。,1) ws=1000。,t,y,39。,0) 16 0 21012Input Output 圖 10階 Butterworth濾波器的濾波效果 系統(tǒng)分析 線性時不變對象 LTI 【例 】已知一個兩輸入兩輸出系統(tǒng)的狀態(tài)方程四對組，據(jù)此建立 LTI 對象“狀態(tài)方程子類”模型。0 ]。 H1=tf([1 1],[1 2])。 。 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)。 B=[0 0 4。H3=1。 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)的傳
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1