Matlab数学实验报告03

《数学实验》报告

实验名称 Matlab绘图及线性代数相关运算 学 院 机械工程学院 专业班级 工业101 姓 名 于洋 学 号 41040318

2011年 10月18

一、 【实验目的】

进一步练习Matlab的三维绘图，并掌握线性代数相关运算。 二、 【实验任务】

1、 用mesh与surf命令绘制三维曲面z?x2、 画三维曲面z?5-x222?3y的图像，并使用不同的着色效果及光照效果。

2?y(?2?x,y?2)与平面z=3的交线。

3、 生成如下数列：

生成一个从-10到10的步长是2 的等差数列。 生成一个从0到-20 的步长是-2的等差数列。 4、 已知矩阵A???-2?1144??，实现下列操作： 7?(1) 提取矩阵A的第一行元素并生成以此为主对角线元素的对角阵 （提示：用diag命令生成对角矩阵。）

(2) 在矩阵A后添加第三行元素(4 7 10)，构成矩阵C. (3) 生成矩阵D=(B C)， 。 (4) 删除矩阵C的第一列。

三、 【实验程序】 1、

t=-5:0.1:3;

[x,y]=meshgrid(t); z=x^2+3*y^3;

subplot(1,3,2);colordef none;

mesh(x,y,z),colormap(hsv),light('position',[-5,6,1]); subplot(1,3,3);colordef black; mesh(x,y,z),colormap(jet);

subplot(1,3,1);colordef white; mesh(x,y,z),colormap(jet)

t=-2:0.1:2;

[x,y]=meshgrid(t); z1=-(x.^2+y.^2)+5;

2

subplot(1,3,1),mesh(x,y,z1),title('oˉêy1') z2=3*ones(size(x));

subplot(1,3,2),mesh(x,y,z2),title('oˉêy2') r0=abs(z2-z1)<=1;

zz=r0.*z2;yy=r0.*y;xx=r0.*x;

subplot(1,3,3),plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0)),title('????') 2、

t=-2:0.1:2;

[x,y]=meshgrid(t); z1=-(x.^2+y.^2)+5;

subplot(1,3,1),mesh(x,y,z1),title('oˉêy1') z2=3*ones(size(x));

subplot(1,3,2),mesh(x,y,z2),title('oˉêy2') r0=abs(z2-z1)<=1;

zz=r0.*z2;yy=r0.*y;xx=r0.*x;

subplot(1,3,3),plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0)),title('????')

3、-10:2:10

0:-2:-20

4、A=[-2 1 4;1 4 7]; B=diag(A(1,:))

A=[-2 1 4;1 4 7]; C=[A;[4 7 10]]

A=[-2 1 4;1 4 7]; B=diag(A(1,:)); C=[A;[4 7 10]]; D=[B,C] F=[B;C]

A=[-2 1 4;1 4 7]; A(:,1)=[]

3

四、 【实验结果】 1、

4

2、

3、 ans =

-10 -8 -6 -4 -2 0 2 4 6 8 10 ans =

0 -2 -4 -6 -8 -10 -12 -14 -16 -18 -20 4、 B =

-2 0 0 0 1 0 0 0 4

5

C =

-2 1 4 1 4 7 4 7 10 D =

-2 0 0 -2 1 4 0 1 0 1 4 7 0 0 4 4 7 10 F=

-2 0 0 0 1 0 0 0 4 -2 1 4 1 4 7 4 7 10 A =

1 4 4 7

五、 【实验总结】

同过实验，我掌握了Matlab的三维绘图及简单的线性代数的运算。收获不少。

6

C =

-2 1 4 1 4 7 4 7 10 D =

-2 0 0 -2 1 4 0 1 0 1 4 7 0 0 4 4 7 10 F=

-2 0 0 0 1 0 0 0 4 -2 1 4 1 4 7 4 7 10 A =

1 4 4 7

五、 【实验总结】

同过实验，我掌握了Matlab的三维绘图及简单的线性代数的运算。收获不少。

6

本文来源：https://www.bwwdw.com/article/0hmr.html