南邮数学实验报告word版本 18页Word格式.docx
- 文档编号:14905071
- 上传时间:2022-10-25
- 格式:DOCX
- 页数:14
- 大小:18.39KB
南邮数学实验报告word版本 18页Word格式.docx
《南邮数学实验报告word版本 18页Word格式.docx》由会员分享,可在线阅读,更多相关《南邮数学实验报告word版本 18页Word格式.docx(14页珍藏版)》请在冰豆网上搜索。
1003003001/6
limit((627*x-sin(627*x))/x^3,x,inf)
1.2y?
ecosxmx,求y'
'
1000
diff(exp(x)*cos(627*x/1000),2)
-201X/1000000*exp(x)*cos(1001/1000*x)-1001/500*exp(x)*sin(1001/10
00*x)
1.3计算?
1?
y2
0?
10ex2dxdy
dblquad(@(x,y)exp(x.^2+y.^2),0,1,0,1)
2.139********228
1.4计算?
x4
m2?
4x2dx
int(x^4/(627^2+4*x^2))
1/12*x^3-100201X/16*x+1003003001/32*atan(2/627*x)
1.5y?
excosmx,求y(10)
diff(exp(x)*cos(627*x),10)
-9389137388146839380380277888*cos(627*x)*exp(x)
149759579095532896918284384*sin(627*x)*exp(x)
1.6
x?
0的泰勒展式(最高次幂为4).
-
taylor(sqrt(627/1000+x),4)
(62500*627^(1/2)*1000^(1/2)*x^3)/246491883-
(125*627^(1/2)*1000^(1/2)*x^2)/393129+
(627^(1/2)*1000^(1/2)*x)/1254+(627^(1/2)*1000^(1/2))/1000
1.7Fibonacci数列{xn}的定义是x1?
1,x2?
1,
xn?
xn?
2(n?
3,4,)用循环语句编程给出该数列的前20项(要求
将结果用向量的形式给出)。
x=[1,1];
forn=3:
20
x(n)=x(n-1)+x(n-2);
end
x
Columns1through10
11235813213455
Columns11through20
891442333776109871597258441816765
211?
A?
020?
1.8对矩阵,求该矩阵的逆矩阵,特征值,特?
m?
41?
1000?
6征向量,行列式,计算A,并求矩阵P,D(D是对角矩阵),使得
PDP?
1。
程序与结果:
a=[-2,1,1;
0,2,0;
-4,1,627/1000];
inv(a)
0.22830.0679-0.3642
00.50000
1.4567-0.3642-0.7283
eig(a)
-0.6865+1.5082i
-0.6865-1.5082i
2.0000
[p,d]=eig(a)
p=
0.2937-0.3372i0.2937+0.3372i0.2425000.97010.89440.89440.0000注:
p的列向量为特征向量
d=
-0.6865+1.5082i000-0.6865-1.5082i0002.0000
a^6
11.968013.0080-4.9910
064.00000
19.9640-4.9910-3.0100
1.9作出如下函数的图形(注:
先用M文件定义函数,再用fplot进行函数作图):
2x0?
2f(x)?
1?
2(1?
x)?
1
2
函数文件f.m:
functiony=f(x)
if0<
=x&
x<
=1/2
y=2.0*x;
else1/2<
x&
=1
y=2.0*(1-x);
fplot(@f,[0,1])
1.10在同一坐标系下作出下面两条空间曲线(要求两条曲线用不同的颜色表示)
cost?
2cost?
(1)?
y?
sint
(2)?
2sint
z?
t?
t=-10:
0.01:
10;
x1=cos(t);
y1=sin(t);
z1=t;
plot3(x1,y1,z1,'
k'
);
holdon
x2=cos(2*t);
y2=sin(2*t);
z2=t;
plot3(x2,y2,z2,'
r'
holdoff
篇二:
南邮201XMATLAB数学实验答案(全)
熟练掌握Matlab软件的基本命令和操作,会作二维、三维几何图形,能够用Matlab软件解决微积分、线性代数与解析几何中的计算问题。
sinmx与limx?
x3x3
limit((902*x-sin(902*x))/x^3)
ans=
366935404/3
limit((902*x-sin(902*x))/x^3,inf)
diff(exp(x)*cos(902*x/1000),2)
(46599*cos((451*x)/500)*exp(x))/250000-(451*sin((451*x)/500)*exp(x))/250
e0011x2?
y2dxdy
2.1394
dx1.4计算?
2m?
4x2
int(x^4/(902^2+4*x^2))
(91733851*atan(x/451))/4-(203401*x)/4+x^3/12
diff(exp(x)*cos(902*x),10)
-356485076957717053044344387763*cos(902*x)*exp(x)-395232302427764249482201X8
84*sin(902*x)*exp(x)
taylor(sqrt(902/1000+x),5,x)
-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602
+(15625*451^(1/2)*500^(1/2)*x^3)/91733851
-(125*451^(1/2)*500^(1/2)*x^2)/406802+(451^(1/2)*500^(1/2)*x)/902
+(451^(1/2)*500^(1/2))/500
1,xn?
3,4,?
)用循环语句编
程给出该数列的前20项(要求将结果用向量的形式给出)。
x=
1.8对矩阵,求该矩阵的逆矩阵,特征值,特征向量,行列式,?
计算A,并求矩阵P,D(D是对角矩阵),使得A?
PDP。
A=[-2,1,1;
-4,1,902/1000];
inv(A)
0.41070.0223-0.4554
1.8215-0.4554-0.9107
eig(A)
-0.5490+1.3764i
-0.5490-1.3764i
det(A)
4.39206?
[P,D]=eig(A)
P=%特征向量
0.3245-0.3078i0.3245+0.3078i0.2425
000.9701
0.89440.89440.0000
D=
-0.5490+1.3764i00
0-0.5490-1.3764i0
002.0000
P*D^6*inv(P)%A^6的值
15.366112.1585+0.0000i-5.8531
23.4124-5.8531+0.0000i-1.6196
m文件:
functiony=fenduan(x)
ifx<
y=2*x
elsex<
y=2-2*x
执行函数:
fplot('
fenduan'
[0,1]);
gridon
title('
第1.9题图'
)
得下图:
第1.9题图
0.0.0.0.0.0.0.0.0.
plot3(x1,y1,z1);
m'
第1.10题图'
第1.10题图
4?
22?
134?
1.11已知A?
305?
B?
20?
3?
,在MATLAB命令窗口中建立A、B矩阵并?
15m3?
2?
11?
对其进行以下操作:
(1)计算矩阵A的行列式的值det(A)
(2)分别计算下列各式:
2A?
B,A*B,A.*B,AB?
1,A?
1B,A2,AT解:
A=[4,-
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 南邮数学实验报告word版本 18页 数学 实验 报告 word 版本 18