MATLAB全实验报告

  • 格式:doc
  • 大小:1.32 MB
  • 文档页数:32

下载文档原格式

  / 32
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

《数学实验》报告

实验名称 Matlab 基础知识

学院

专业班级

姓名

学号

2014年 6月

一、【实验目的】

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

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

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

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

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

二、【实验任务】

P16 第4题

编写函数文件,计算

1!

n

k

k =

∑,并求出当k=20时表达式的值。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)&(AB)。

P34 第1题

111

1

4357

π

=-+-+……公式求π的近似值,直到某一项的绝对值小于-6

10为止。

三、【实验程序】

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

>> (A==B)&(A

>> (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)