南郵Matlab實驗報告三.doc

Matlab上機實驗報告實驗名稱：MATLAB的數值運算班級： 自動化二班 學號： B11050216 姓名： 李鵬飛 南京郵電大學 2013年4月9日一、實驗目的 1學習MATLAB的基本矩陣運算； 2學習MATLAB的點運算；3學習復雜運算。二、實驗基本知識 1.基本矩陣運算； 2.多項式運算； 3.代數方程組；4.數值微積分。 三、上機練習1. 給a,b,c賦如下數據:1)求a+b,a*b,a.*b,a/b,a./b,a2,a.2,length(c)的結果.2)求c中所有元素的平均值、最大值.3)求d=b(2:3,1,3).程序代碼：（1） a=1,3,4;5,12,44;7,8,27; b=-7,8,4;12,24,38;68,-5,3; c=1,0,8,3,6,2,-4,23,46,6; d=a+bd = -6 11 8 17 36 82 75 3 30 e=a*be = 301 60 130 3101 108 608 1883 113 413 f=a.*bf = -7 24 16 60 288 1672 476 -40 81 g=a/bg = 0.0966 0.0945 0.0080 -3.6125 1.5838 -0.5778 -1.9917 0.9414 -0.2682 h=a./bh = -0.1429 0.3750 1.0000 0.4167 0.5000 1.1579 0.1029 -1.6000 9.0000 i=a2i = 44 71 244 373 511 1736 236 333 1109 j=a.2j = 1 9 16 25 144 1936 49 64 729 k=length(c)k =10（2） l=max(c)l =46 m=mean(c)m =9.1000（3） d=b(2:3,1,3)d = 12 3868 32.求解方程，并對所得結果作出解釋。程序代碼： A=1,2,3;4,5,6;7,8,9;b=1,0,-1;ra=rank(A)ra = 2 rab=rank(A,b)rab = 2 xs=Ab;Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.541976e-018. xg=null(A);c=rand(1);ba=A*(xs+c*xg)ba = 1.0000 -0.0000 -1.0000 norm(ba-b)ans = 1.7764e-0153.已知有理分式，其中，。求該分式的商多項式和余多項式。程序代碼： format ratp1=conv(3,0,1,0,1,0,0,0.5);p2=5 12 -6 -3 2 -2;q,r=deconv(p1,p2);cq=商多項式為 ;cr=余多項式為 ;disp(cq,poly2str(q,s),disp(cr,poly2str(r,s)商多項式為 0.6 s - 1.44余多項式為 21.88 s4 - 5.34 s3 - 5.52 s2 + 4.58 s - 2.88 qp2=conv(q,p2);pp1=qp2+r;pp1=p1ans = 1 1 0 0 1 1 14.求解Lorenz模型的狀態(tài)方程，并圖示.初值在(0,1)中取值程序代碼：M-文件：function dx = myfun(t,x) dx = zeros(3,1); dx(1)=16*(x(2)-x(1); dx(2)=45*x(1)-x(2)-x(1)*x(3); dx(3)=x(1)*x(2)-4*x(3) ;comand window中：T,Y = ode45(myfun,0 1,1 1 1);plot(T,Y)tspan=0,1;x0=0.5;0.4;0.6;tt,xx=ode45(myfun,tspan,x0);subplot(3,1,1),plot(tt,xx(:,1),-b)hold onxlabel(t)ylabel(x1(t)title(Lorenz)legend(fontsize14it x1(t)hold onsubplot(3,1,2),plot(tt,xx(:,2),r*)xlabel(t)ylabel(x2(t)legend(fontsize14it x2(t)subplot(3,1,3),plot(tt,xx(:,3),go)xlabel(t)ylabel(x3(t)legend(fontsize14it x3(t)四、實驗小結：在實驗三：MATLAB的符號運算中，我們學會了用MATLAB解決實際的數學和符號運算問題，但是在使用軟件的過程中，我還是遇到了不少的問題，有很多的代碼不熟悉，這需要我在日后的學習過程中不斷加強深入學習，多增加動手操作，爭取早日熟悉MATLAB的基本操作要求。在實驗第四題中，用到了M文件的使用，同