南邮matlab实验报告三

《南邮matlab实验报告三》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、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,a^2,a.^2,length(c)的结果.2)求c中所有元素的平均

2、值、最大值.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=-611817368275330>>e=a*be=3016013031011086081883113413>>f=a.*bf=-72416602881672476-4081>>g=a/bg=0.09660.09450.0080-3.61251.5838-0.5778-

3、1.99170.9414-0.2682>>h=a./bh=-0.14290.37501.00000.41670.50001.15790.1029-1.60009.0000>>i=a^2i=447124437351117362363331109>>j=a.^2j=19162514419364964729>>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=12386832.求解方程，并对所得结果作出解释。程序代码：>>A=[1

4、,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:Matrixisclosetosingularorbadlyscaled.Resultsmaybeinaccurate.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.已知有理分式，其中

5、，。求该分式的商多项式和余多项式。程序代码：>>formatratp1=conv([3,0,1,0],[1,0,0,0.5]);p2=[512-6-32-2];[q,r]=deconv(p1,p2);cq='商多项式为';cr='余多项式为';disp([cq,poly2str(q,'s')]),disp([cr,poly2str(r,'s')])商多项式为0.6s-1.44余多项式为21.88s^4-5.34s^3-5.52s^2+4.58s-2.88>>qp2=conv(q,p2);pp1=qp2+r;pp1==

6、p1ans=11001114.求解Lorenz模型的状态方程，并图示.初值在(0,1)中取值程序代码：M-文件：functiondx=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);comandwindow中：[T,Y]=ode45(@myfun,[01],[111]);plot(T,Y)tspan=[0,1];x0=[0.5;0.4;0.6];[tt,xx]=ode45(

7、@myfun,tspan,x0);subplot(3,1,1),plot(tt,xx(:,1),'-b')holdonxlabel('t')ylabel('x1(t)')title('LorenzÄ£ÐÍͼ')legend('fontsize{14}itx1(t)')holdonsubplot(3,1,2),plot(tt,xx(:,2),'r*')xlabel('t')ylabel('x2(t)')legend('fontsize{14}itx2(t)')subplot(3,1,3),plot(tt,xx

8、(:,3),'go')xlabel('t')ylabel('x3(t)')legend('fontsize{14}itx3(t)')四、实验小结：在实验三：MATLAB的符号运算中，我们学会了用MATLAB解决实际的数学和符号运算问题，但是在使用软件的过程中，我还是遇到了不少的问题，有很多的代码不熟悉，这需要我在日后的学习过程中不断加强深入