MATLAB全实验报告
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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)