北科大Matlab数学实验分析报告次全

精心整理《数学实验》报告

实验名称Matlab基础知识

学院

一、【实验目的】

1.认识熟悉Matlab这一软件，并在此基础上学会基本操作。

2.掌握Matlab基本操作和常用命令。

3.了解Matlab常用函数，运算符和表达式。

4.掌握Matlab工作方式和M文件的相关知识。

5.学会Matlab中矩阵和数组的运算。

二、【实验任务】

P16第4题

P27

矩阵

P27

已运算P34

π

用

4

P16

for

for

end

sum=sum+y;

end

sum

P27第2题

>>A=[123;456;789]

>>B=[468;556;322]

>>A*B

>>A.*B

P27第3题

>>A=[52;91];B=[12;92];

>>A>B

>>A==B

>>A

>>(A==B)&(A>(A==B)&(A>B)

P34第1题

t=1;

pi=0;

n=1;

s=1;

end

P27

P27

P34

>>pi

pi=

了解并掌握matlab的基本绘图二、【实验任务】

P79页1,3,5题

三、【实验程序】

1.

clf;

3.

clf;

5.

t=0:pi/50:20*pi;

x=t.*cos(t*pi/6);

y=t.*sin(t*pi/6);

z=2*t;

plot3(x,y,z)

四、【实验结果】

1．

3.

5.

通过本次课程和作业，我初步了解了matlab在绘图方面的优势和重要性。

1.学会用Matlab 进行三维的曲线绘图；

2.掌握绘图的基本指令和参数设置

二、 【实验任务】

P79习题5

⎪⎪⎪⎩

⎪⎪⎪⎨⎧===z y x P79xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')

习题9：

clf;

t=-2:0.1:2;

[x,y]=meshgrid(t);

z1=5-x.^2-y.^2;

subplot(1,2,1),mesh(x,y,z1),title('曲面z1=5-x.^2-y.^2')

z2=3*ones(size(x));

r0=abs(z1-z2)<=0.05;

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

subplot(1,2,2),plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'.') title('交线')

四、【实验结果】

习题5：

习题9：

1.学会用Matlab 练习使用矩阵的基本运算；

2.掌握用Matlab 运用矩阵的特征值、特征向量、特征多项式；

3.学会用Matlab 解线性方程组；

4.掌握用Matlab 进行数值方法计算定积分

用三种方法求下列积分的数值解：

(2)dx x x x ⎰+π02cos 1sin

P167习题18

用多种数值方法计算定积分⎰

40sin -11π

dx x

，并与精确值2进行比较，观察不同方法相应的误差。

习题12

>>A=[195365;246810;346972;4678104;507321;386319] >>A'

>>det(A)

>>rank(A)

>>rref(A)

习题14：

（2

s1=sum(y1(1:(t-1)))*h

s2=sum(y1(2:t))*h

s3=trapz(x,y1)

s4=quad('jifen',0,pi)

习题18：

function y=jifen(x)

y=1./(1-sin(x));

h=0.01;

x=0:h:pi/4;

y=1./(1-sin(x));

t=length(x);

format long

s1=sum(y1(1:(t-1)))*h s2=sum(y1(2:t))*h s3=trapz(x,y)

s4=quad('jifen',0,pi/4) format short

u1=s1-sqrt(2)

u2=s2-sqrt(2)

u3=s3-sqrt(2)

（1⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡4321x x x x ⎥⎥⎦

⎢⎢⎣125(2)

解对应的齐次方程组⎩⎨⎧=+=434212x x x x x ，可得一个基础解系：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1212ε

原方程组对应的同解方程组为：⎪⎪⎩⎪⎪⎨

⎧+=++=2122143421x x x x x ，可找到一个特解为：⎥⎥⎥⎥⎥⎥⎥⎦

⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡=125125*

η 因此，此方程组的通解为：

习题17： （2） 习题很