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

《《数学实验与数学软件》课程实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、广西工学院信计系《数学实验与数学软件》课程实验报告告姓名专业班级学号实验一MATLAB的基本使用1、设有分块矩阵，其中E、R、O、S分别为单位阵、随机阵、零矩阵和对角阵。试通过数值计算验证。解：在在命令窗口输入：>>E=eye(3);R=rand(3,2);O=zeros(2,3);S=eye(2,2);>>A=[ER;OS];X=A^2,Y=[ER+R*S;OS^2];X=2.18131.09281.65341.00001.00001.18132.09281.65341.00001.00001.18131.09282.65341.00001.00000.

2、95010.60680.89132.44832.44830.23110.48600.76211.47921.4792Ｙ=2.18131.09281.65341.00001.00001.18132.09281.65341.00001.00001.18131.09282.65341.00001.00000.95010.60680.89132.44832.44830.23110.48600.76211.47921.4792X=Y，即得证2、求下列矩阵方程的解：⑴；⑵。其中，，解：（1）在命令窗口输入：>>A=[21-1;210;1-11];B=[1;2;3];C

3、=[1-13;432];>>X=ABX=1.3333广西工学院信计系-0.66671.0000所以矩阵方程AX=B的解为（2）接着（1）在命令窗口继续输入：>>X=C/AX=所以矩阵方程XA=C的解为3、用两种方法在同一个坐标系下作出以下四条曲线的图形，并在图上添加各种标注（如函数名、图例等）：。解：x=-1:0.01:1;y1=x.^2;y2=x.^3;y3=x.^4;y4=x.^5;>>plot(x,y1)>>holdon>>plot(x,y2)>>plot(x,y3)>>plot(x,y4)>>legend('x.^2','x.^3','x.^4'

4、,'x.^5')title('xyw')>x=linspace(-1,1,100);>>plot(x,x.^2,x,x.^3,x,x.^4,x,x.^5)>>gridon>>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)

5、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=[];forx=100:999a=fix(x/100);b=fix((x-a*100)/10);c=x-a*100-b*10;ifx==a^3+b^3+c^3A=[A,x];endend

6、disp(A)1533703714077、编写函数M文件SQRT.m：用迭代法求的值。其迭代公式为：。迭代的终止条件为。解：建立M文件%函数mysqrt.mfunctiona=mysqrt(x)Ifx==0;a=x;elsey=abs(x)x1=y;x2=(x1+y/x1)/2;whilesbs(x2-x1)>=1e-5x1=x2;x2=(x1+y/x1)/2;endifx>0;a=x2;elsea=x2*i;endend