北科大Matlab_数学实验报告1~6次(全)

《数学实验》报告

实验名称 Matlab 基础知识 学 院 专业班级 姓 名 学 号

2014年 6月

一、 【实验目的】

1.认识熟悉Matlab这一软件，并在此基础上学会基本操作。 2.掌握Matlab基本操作和常用命令。

3.了解Matlab常用函数，运算符和表达式。 4.掌握Matlab工作方式和M文件的相关知识。 5.学会Matlab中矩阵和数组的运算。

二、 【实验任务】

P16 第4题

编写函数文件，计算 P27第2题

n

k!，并求出当k=20时表达式的值。

k 1

123 468 矩阵A=456，B=556，计算A*B,A.*B,并比较两者的区别。 789 322

P27第3题 已知矩阵A=

52

，B=91 12

92 ，做简单的关系运算A>B,A==B,A<B,并做逻辑运算

(A==B)&(A<B),(A==B)&(A>B)。

P34 第1题 用

111

1 ……公式求 的近似值，直到某一项的绝对值小于10-6为止。 4357

三、 【实验程序】

P16 第4题

function sum=jiecheng(n) sum=0; y=1;

for k=1:n for i=1:k y=y*i; end

sum=sum+y; end sum P27第2题

>>A=[1 2 3;4 5 6;7 8 9]

>>B=[4 6 8;5 5 6;3 2 2] >>A*B

P27第3题

>> A=[5 2;9 1];B=[1 2;9 2]; >>A>B >>A==B >>A<B

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

P34 第1题 t=1; pi=0;

n=1; s=1;

while abs(t)>=1e-6 pi=pi+t; n=n+2; s=-s; t=s/n; end

pi=4*pi;

四、 【实验结果】

P16 第4题

P27第2题

两者的区别：A*B是按正规算法进行矩阵的计算, A.*B是对应元素相乘。

P27第3题

P34 第1题 >> pi pi=

3.141590653589692e+000

五、 【实验总结】

这次实验是第一次接触Matlab这个软件，所以有些生疏，花的时间也比较多，但功夫不怕有心人，而且当一个程序弄出来后感觉也特别开心，以后再继续努力学习。

了解并掌握matlab的基本绘图

二、 【实验任务】

P79页 1,3,5题

三、 【实验程序】

1.

clf;

x=0:pi/50:4*pi;

y1=exp(x/3).*sin(3*x);y2=exp(x/3);y3=-exp(x/3); plot(x,y1,'b*',x,y2,'r-.',x,y3,'r-.'),grid on legend('y1=exp(x/3).*sin(3*x)','y2=+-exp(x/3)')

3.

clf;

x1=-pi:pi/50:pi; x2=pi:pi/50:4*pi; x3=1:0.1:8; y1=x1.*cos(x1);

y2=x2.*tan(x2.^(-1)).*sin(x2.^3); y3=exp(x3.^(-1)).*sin(x3);

subplot(2,2,1),plot(x1,y1,'m.'),grid on,title('y=x*cosx') xlabel('xÖá'),ylabel('yÖá')

gtext('y=x*cosx'),legend('y=x*cosx')

subplot(2,2,2),plot(x2,y2,'r*'),grid on,title('y=x*tan(1/x)*sin(x^3)') xlabel('xÖá'),ylabel('yÖá')

gtext('y=x*tan(1/x)*sin(x^3)'),legend('y=x*tan(1/x)*sin(x^3)') subplot(2,2,3),plot(x3,y3,'bp'),grid on,title('y=e(1/x3)*sinx') xlabel('xÖá'),ylabel('yÖá')

gtext('y=e(1/x3)*sinx'),legend('y=e(1/x3)*sinx')

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.

五、 【实验总结】

一、 【实验目的】

1. 学会用Matlab进行三维的曲线绘图； 2. 掌握绘图的基本指令和参数设置 二、 【实验任务】 P79 习题5

绘制圆锥螺线的图像并加标注，圆锥螺线的参数方程为；

x tcot 6

y tsit (0 t 20 )。

6

z 2t

P79 习题9

22

画三维曲线z 5-x-y(-2 x,y 2)与平面z=3的交线。

三、 【实验程序】 习题5：

clf;

t=0:pi/100:20*pi; x=t.*cos(t.*pi/6); y=t.*sin(t.*pi/6); z=2*t;

plot3(x,y,z)

title('圆锥螺线')

xlabel('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;

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进行数值方法计算定积分 二、 【实验任务】 P114 习题12

随机输入一个六阶方阵，并求其转置、行列式、秩，以及行最简式。 P114 习题14

2 1 1

1 2 1求矩阵A 的特征多项式、特征值和特征向量。 1 1 2

P115 习题20

求下列线性方程组的通解：

x1-x2-x3 x4 0 x1 x2 2x3-x4 0

(1) -x1 x2 3x3 0 (2) x1-x2 x3-3x4 1 2x-3x 4x-x 0 1234 x1-x2-2x3 3x4 -1

2

P167 习题17

用三种方法求下列积分的数值解： (2) 0

xsinx

1 cosx

2

P167 习题18

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

1

dx，并与精确值2进行比较，观察

1-sinx

习题12

>> A=[1 9 5 3 6 5;2 4 6 8 1 0;3 4 6 9 7 2;4 6 7 8 10 4;5 0 7 3 2 1;3 8 6 3 1 9] >> A'

>> det(A) >> rank(A) >> rref(A)

习题14：

>> B=[2 1 1;1 2 1;1 1 2] >> p=poly(B) >> [V D]=eig(B)

习题20： (1)

>> A=[1 1 2 -4;-1 1 3 0;2 -3 4 -1] >> rref(A) (2)

>> B=[1 -1 -1 1;1 -1 1 -3;1 -1 -2 3] >> rref(B)

>> C=[1 -1 -1 1 0;1 -1 1 -3 1;1 -1 -2 3 -1/2] >> rref(C)

习题17：

（2）

function y=jifen(x)

y=x.*sin(x)./(1+cos(x).^2); h=0.01; x=0:h:pi;

y0=1+cos(x).^2; y1=x.*sin(x)./y0; t=length(x);

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-sq …… 此处隐藏：2128字，全部文档内容请下载后查看。喜欢就下载吧 ……