南邮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模型的状态方程，并图示.初值在(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文件的使用，同