南郵Matlab實(shí)驗(yàn)報(bào)告三.doc

Matlab上機(jī)實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)名稱：MATLAB的數(shù)值運(yùn)算班級： 自動(dòng)化二班 學(xué)號： B11050216 姓名： 李鵬飛 南京郵電大學(xué) 2013年4月9日一、實(shí)驗(yàn)?zāi)康?1學(xué)習(xí)MATLAB的基本矩陣運(yùn)算； 2學(xué)習(xí)MATLAB的點(diǎn)運(yùn)算；3學(xué)習(xí)復(fù)雜運(yùn)算。二、實(shí)驗(yàn)基本知識 1.基本矩陣運(yùn)算； 2.多項(xiàng)式運(yùn)算； 3.代數(shù)方程組；4.數(shù)值微積分。 三、上機(jī)練習(xí)1. 給a,b,c賦如下數(shù)據(jù):1)求a+b,a*b,a.*b,a/b,a./b,a2,a.2,length(c)的結(jié)果.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.求解方程，并對所得結(jié)果作出解釋。程序代碼： 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.已知有理分式，其中，。求該分式的商多項(xiàng)式和余多項(xiàng)式。程序代碼： 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=商多項(xiàng)式為 ;cr=余多項(xiàng)式為 ;disp(cq,poly2str(q,s),disp(cr,poly2str(r,s)商多項(xiàng)式為 0.6 s - 1.44余多項(xiàng)式為 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)四、實(shí)驗(yàn)小結(jié)：在實(shí)驗(yàn)三：MATLAB的符號運(yùn)算中，我們學(xué)會(huì)了用MATLAB解決實(shí)際的數(shù)學(xué)和符號運(yùn)算問題，但是在使用軟件的過程中，我還是遇到了不少的問題，有很多的代碼不熟悉，這需要我在日后的學(xué)習(xí)過程中不斷加強(qiáng)深入學(xué)習(xí)，多增加動(dòng)手操作，爭取早日熟悉MATLAB的基本操作要求。在實(shí)驗(yàn)第四題中，用到了M文件的使用，同