《数学实验与数学软件》课程实验报告.docVIP

广西工学院信计系 《数学实验与数学软件》课程实验报告告姓名 专业班级 学号 实验一 MATLAB的基本使用 1、设有分块矩阵，其中E、R、O、S分别为单位阵、随机阵、零矩阵和对角阵。试通过数值计算验证。 解：在在命令窗口输入： E=eye(3);R=rand (3,2);O=zeros (2,3);S=eye(2,2); A=[E R;O S];X=A^2,Y=[E R+R*S;O S^2]; X = 2.1813 1.0928 1.6534 1.0000 1.0000 1.1813 2.0928 1.6534 1.0000 1.0000 1.1813 1.0928 2.6534 1.0000 1.0000 0.9501 0.6068 0.8913 2.4483 2.4483 0.2311 0.4860 0.7621 1.4792 1.4792 Ｙ= 2.1813 1.0928 1.6534 1.0000 1.0000 1.1813 2.0928 1.6534 1.0000 1.0000 1.1813 1.0928 2.6534 1.0000 1.0000 0.9501 0.6068 0.8913 2.4483 2.4483 0.2311 0.4860 0.7621 1.4792 1.4792 X=Y，即得证 2、求下列矩阵方程的解：⑴；⑵。 其中，， 解：（1）在命令窗口输入： A=[2 1 -1;2 1 0;1 -1 1];B=[1;2;3];C=[1 -1 3;4 3 2]; X=A\B X = 1.3333 -0.6667 1.0000 所以矩阵方程AX=B的解为 （2）接着（1）在命令窗口继续输入： X=C/A X = 所以矩阵方程XA=C的解为 3、用两种方法在同一个坐标系下作出以下四条曲线的图形，并在图上添加各种标注（如函数名、图例等）：。 解：x=-1:0.01:1; y1=x.^2;y2=x.^3;y3=x.^4;y4=x.^5; plot(x,y1) hold on plot(x,y2) plot(x,y3) plot(x,y4) legend(x.^2,x.^3,x.^4,x.^5) title(xyw) x=linspace(-1,1,100); plot(x,x.^2,x,x.^3,x,x.^4,x,x.^5) grid on legend(x.^2,x.^3,x.^4,x.^5) title(xyw) 4、用subplot风别在不同的坐标系下作出以下曲线，为每幅图像加上标题： ⑴概率曲线；⑵四叶玫瑰线（极坐标） 解：x=-5:0.01:5; y=exp(-x.^2); subplot(1,2,1) plot(x,y) x=0:pi/180:2*pi; y=sin(2*x); subplot(1,2,2) polar(x,y) 5、作出以下曲面的图形：。 解：x=-2:0.01:2; y=-2:0.01:2; [x,y]=meshgrid(x,y); z=sin(pi*sqrt(x.^2+y.^2)); mesh(x,y,z) 6、建立一个M文件，求所有的“水仙花数”。所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为。 解：A=[]; for x=100:999 a=fix(x/100);b=fix((x-a*100)/10);c=x-a*100-b*10; if x==a^3+b^3+c^3 A=[A,x]; end end disp(A) 153 370 371 407 7、编写函数M文件SQRT.m：用迭代法求的值。其迭代公式为：。迭代的终止条件为。 解：建立M文件 %函数 mysqrt.m function a=mysqrt(x) If x==0;a=x; else y=abs(x) x1=y;x2=(x1+y/x1)/2; while sbs(x2-x1)=