哈工大_数学实验报告

《哈工大_数学实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、作者：王卓学号：1013600414班级：01360041数学实验报告第17页共17页作者：王卓学号：1013600414班级：0136004实验一　Matlab的使用1.上机实验各种数据输入方法：第17页共17页作者：王卓学号：1013600414班级：0136004程序语句：a=[123;456;7,8,9]程序语句：linspace(1,10,5)等等…………计算结果：a=计算结果：ans=1234561.00003.25005.50007.750010.0000789第17页共17页作者：王卓学号：1013600414班级：01360042.(1)(a

2、)方法：(b)方法：程序语句：程序语句：a=[-3508;1-82-1;0-593;-70-45];a=[-3508;1-82-1;0-593;-70-45];b=[0;2;-1;6];b=[0;2;-1;6];ab计算结果：计算结果：ans=ans=-0.6386-0.6386-0.4210-0.4210-0.3529-0.35290.02370.0237(2)4个矩阵的生成语句：矩阵a的生成语句：e=eye(3,3);a=[er;os]r=rand(3,2);验证语句：o=zeros(2,3);a^2s=diag([1,2]);%此为一个任取的2X2矩阵

3、b=[er+r*s;os^2]计算结果相同：ans=1.0000001.90031.457901.000000.46232.6739001.00001.21372.28630001.0000000004.00003.生成多项式的语句：poly([2,-3,1+2i,1-2i,0,-6])计算结果：ans=15-9-172-1800计算x＝0.8,-x=-1.2之值的指令与结果：指令：polyval([1,5,-9,-1,72,-180,0],0.8)结果：ans=-100.2179指令：polyval([1,5,-9,-1,72,-180,0],-1.2)结果

4、：ans=293.29004.求a的指令与结果：指令：a=compan([1,0,-6,3,-8])结果：a=06-38第17页共17页作者：王卓学号：1013600414班级：0136004100001000010求a的特征值的指令与结果：roots(p)的指令与结果为：指令：eig(a)指令：roots([1,0,-6,3,-8])结果：结果：ans=ans=-2.8374-2.83742.46922.46920.1841+1.0526i0.1841+1.0526i0.1841-1.0526i0.1841-1.0526i结论：利用友元阵函数a=compan

5、y(p)和eig(a)可以与roots(p)有相同的作用，结果相同。1.(1)图形：作图指令：P3.5x=0:0.01:1.5;y=[x.^2;x.^3;x.^4;x.^5];plot(x,y)(2)图形：P3.5(2)作图指令：x=0:0.01:10;y1=x.^2;y2=x.^3;y3=x.^4;y4=x.^5;subplot(2,2,1),plot(x,y1),title('x^2')subplot(2,2,2),plot(x,y2),title('x^3')subplot(2,2,3),plot(x,y3),title('x^4')subplot(2,

6、2,4),plot(x,y4),title('x^5')第17页共17页作者：王卓学号：1013600414班级：01360041.fplot命令：建立函数：funfplot(x)functiony=funfplot(x)y=sin(1./tan(pi*x));%此文件以funfplot.m存盘绘图语句：fplot('funfplot',[-0.1,0.1],2e-4)2.用作图法求4sinx–x-2=0的根的近似值P4.7绘图语句为：x=-3*pi:0.01:pi;y1=4*sin(x);y2=x+2;plot(x,y1,x,y2)对x，y轴加标记并加网格：

7、gridxlabel('x方向')ylabel('y方向')P4.83.作曲面z=x2–y2的三维图形绘图语句(a)：x=-10:0.5:10;y=x;[x,y]=meshgrid(x,y);z=x.^2-y.^2;mesh(x,y,z)第17页共17页作者：王卓学号：1013600414班级：0136004P4.8绘图语句(b)：x=-10:0.2:10;y=x;[x,y]=meshgrid(x,y);z=x.^2-y.^2;plot3(x,y,z)1.（a）自然数的阶乘：求阶乘的主程序：建立阶乘函数：jie(x)n=input('你想求几的阶乘？');fu

8、nctiony=jie(n);y=1;