原创力文档- 1、本文档共53页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
数学实验答案
第一次练习
教学要求：熟练掌握Matlab软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab软件解决微积分、线性代数与解析几何中的计算问题。
补充命令
vpa(x,n) 显示x的n位有效数字，教材102页
fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形
在下面的题目中为你的学号的后3位（1-9班）或4位（10班以上）
1.1 计算与
程序：
syms x
limit((1001*x-sin(1001*x))/x^3,x,0)
结果：
1003003001/6
程序：
syms x
limit((1001*x-sin(1001*x))/x^3,x,inf)
结果：
0
1.2 ，求
程序：
syms x
diff(exp(x)*cos(1001*x/1000),2)
结果：
-2001/1000000*exp(x)*cos(1001/1000*x)-1001/500*exp(x)*sin(1001/1000*x)
1.3 计算
程序：
dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)
结果：
2.13935019514228
1.4 计算
程序：
syms x
int(x^4/(1000^2+4*x^2))
结果：
1/12*x^3-1002001/16*x+1003003001/32*atan(2/1001*x)
1.5
程序：
syms x
diff(exp(x)*cos(1000*x),10)
结果：
-1009999759158992000960720160000*exp(x)*cos(1001*x)-10090239998990319040000160032*exp(x)*sin(1001*x)
1.6 给出在的泰勒展式(最高次幂为4).
程序：
syms x
taylor(sqrt(1001/1000+x),5)
结果：
1/100*10010^(1/2)+5/1001*10010^(1/2)*x-1250/1002001*10010^(1/2)*x^2+625000/1003003001*10010^(1/2)*x^3-390625000/1004006004001*10010^(1/2)*x^4
1.7 Fibonacci数列的定义是，
用循环语句编程给出该数列的前20项（要求将结果用向量的形式给出）。
程序：
x=[1,1];
for n=3:20
x(n)=x(n-1)+x(n-2);
end
x
结果：
Columns 1 through 10
1 1 2 3 5 8 13 21 34 55
Columns 11 through 20
89 144 233 377 610 987 1597 2584 4181 6765
1.8 对矩阵，求该矩阵的逆矩阵，特征值，特征向量，行列式，计算，并求矩阵（是对角矩阵），使得。
程序与结果：
a=[-2,1,1;0,2,0;-4,1,1001/1000];
inv(a)
0.50100100100100 -0.00025025025025 -0.50050050050050
0 0.50000000000000 0
2.00200200200200 -0.50050050050050 -1.00100100100100
eig(a)
-0.49950000000000 + 1.32230849275046i
-0.49950000000000 - 1.32230849275046i
2.00000000000000
[p,d]=eig(a)
p =
0.3355 - 0.2957i 0.3355 + 0.2957i 0.2425
0 0 0.9701
0.8944 0.8944 0.0000
注：p的列向量为特征向量
d =
-0.
文档评论（0）