北科大Matlab数学实验报告16次全文档格式.docx
- 文档编号:14347225
- 上传时间:2022-10-22
- 格式:DOCX
- 页数:32
- 大小:518.08KB
北科大Matlab数学实验报告16次全文档格式.docx
《北科大Matlab数学实验报告16次全文档格式.docx》由会员分享,可在线阅读,更多相关《北科大Matlab数学实验报告16次全文档格式.docx(32页珍藏版)》请在冰豆网上搜索。
B)。
P34第1题
用公式求的近似值,直到某一项的绝对值小于为止。
三、【实验程序】
P16第4题
functionsum=jiecheng(n)
sum=0;
y=1;
fork=1:
n
fori=1:
k
y=y*i;
end
sum=sum+y;
end
sum
>
>
A=[123;
456;
789]
B=[468;
556;
322]
A*B
A、*B
A=[52;
91];
B=[12;
92];
A>
B
A==B
A<
(A==B)&
B)
t=1;
pi=0;
n=1;
s=1;
whileabs(t)>
=1e-6
pi=pi+t;
n=n+2;
s=-s;
t=s/n;
pi=4*pi;
四、【实验结果】
P16第4题
两者的区别:
A*B就是按正规算法进行矩阵的计算,A、*B就是对应元素相乘。
P34第1题
pi
pi=
3、1492e+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,'
),gridon
legend('
y1=exp(x/3)、*sin(3*x)'
'
y2=+-exp(x/3)'
)
3、
x1=-pi:
pi;
x2=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、'
),gridon,title('
y=x*cosx'
xlabel('
xÖ
á
'
),ylabel('
yÖ
gtext('
),legend('
subplot(2,2,2),plot(x2,y2,'
r*'
y=x*tan(1/x)*sin(x^3)'
subplot(2,2,3),plot(x3,y3,'
bp'
y=e(1/x3)*sinx'
5、
t=0:
20*pi;
x=t、*cos(t*pi/6);
y=t、*sin(t*pi/6);
z=2*t;
plot3(x,y,z)
1.
通过本次课程与作业,我初步了解了matlab在绘图方面的优势与重要性。
1、学会用Matlab进行三维的曲线绘图;
2、掌握绘图的基本指令与参数设置
P79习题5
绘制圆锥螺线的图像并加标注,圆锥螺线的参数方程为;
。
P79习题9
画三维曲线与平面z=3的交线。
习题5:
pi/100:
x=t、*cos(t、*pi/6);
y=t、*sin(t、*pi/6);
title('
圆锥螺线'
x轴'
y轴'
),zlabel('
z轴'
习题9:
t=-2:
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),'
、'
交线'
这次三维曲线(曲面)的绘制虽然不算复杂,但还就是要注意一些细节,而且要注意弄懂其中的原因,不能硬套书上的,否则很容易不明道理的出错。
1、学会用Matlab练习使用矩阵的基本运算;
2、掌握用Matlab运用矩阵的特征值、特征向量、特征多项式;
3、学会用Matlab解线性方程组;
4、掌握用Matlab进行数值方法计算定积分
P114习题12
随机输入一个六阶方阵,并求其转置、行列式、秩,以及行最简式。
P114习题14
求矩阵的特征多项式、特征值与特征向量。
P115习题20
求下列线性方程组的通解:
(1)
(2)
P167习题17
用三种方法求下列积分的数值解:
(2)
P167习题18
用多种数值方法计算定积分,并与精确值进行比较,观察不同方法相应的误差。
习题12
A=[195365;
246810;
346972;
4678104;
507321;
386319]
A'
det(A)
rank(A)
rref(A)
习题14:
B=[211;
121;
112]
p=poly(B)
[VD]=eig(B)
习题20:
(1)
A=[112-4;
-1130;
2-34-1]
B=[1-1-11;
1-11-3;
1-1-23]
rref(B)
C=[1-1-110;
1-11-31;
1-1-23-1/2]
rref(C)
习题17:
functiony=jifen(x)
y=x、*sin(x)、/(1+cos(x)、^2);
h=0、01;
h:
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:
y=1、/(1-sin(x));
pi/4;
formatlong
s3=trapz(x,y)
s4=quad('
0,pi/4)
formatshort
u1=s1-sqrt
(2)
u2=s2-sqrt
(2)
u3=s3-sqrt
(2)
u4=s4-sqrt
(2)
习题14
习题20
原方程对应的同解方程组为:
解得方程基础解系为:
所以方程组的通解为:
=
(2)
解对应的齐次方程组,可得一个基础解系:
原方程组对应的同解方程组为:
可找到一个特解为:
因此,此方程组的通解为:
在掌握线性代数相关运算与数值积分的理论基础上进行操作,学会了用Matlab相关指令与编程,并进行计算与误差分析,感觉原来很繁琐的计算用Matlab很方便!
实验目的】
1、学会用Matlab进行曲线拟合与使用插值函数;
2、掌握曲线拟合与插值处理的基本指令与参数设置
P130习题9
已知在某实验中测得某质点的位移s与速度v随时间t变化如下:
、
t
0、5
1、0
1、5
2、0
2、5
3、0
v
0、4794
0、8415
0、9975
0、9093
0、5985
0、1411
s
1
2
3
3、5
4
求质点的速度与位移随时间的变化曲线以及位移随速度变化曲线。
P130习题10
在某种添加剂的不同浓度之下对铝合金进行抗拉强实验,得到数据如下,现分别使用不同的插值方法,对其中没有测量的浓度进行推测,并估算出浓度X=18及26时的抗压强度Y的值。
浓度X
10
15
20
25
30
抗压强度Y
25、2
29、8
31、2
31、7
29、4
P130习题12
利用不同的方法对在(-3,3)上的二维插值效果进行比较。
0、5:
3;
v=[00、47940、84150、99750、90930、59850、1411];
s=[11、522、533、54];
p1=polyfit(t,v,2);
p2=polyfit(t,s,2);
p3=polyfit(s,v,2);
disp('
速度与时间函数'
),f1=poly2str(p1,'
t'
位移与时间的函数'
),f2=poly2str(p2,'
位移与速度的函数'
),
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 北科大 Matlab 数学 实验 报告 16