南邮数学实验答案全.pdf

第一次练习

教学要求：熟练掌握Matlab软件的基本命令和操作，会作二维、三维几何

图形，能够用Matlab软件解决微积分、线性代数与解析几何中的计算问题。

补充命令

vpa(x,n)显示x的n位有效数字，教材102页

fplot(‘f(x)’,[a,b])函数作图命令，画出f(x)在区间[a,b]上的图形

在下面的题目中m为你的学号的后3位（1-9班）或4位（10班以上）

mxsinmxmxsinmx

1.1计算lim与lim

x0x3xx3

程序：

symsx

limit((1001*x-sin(1001*x))/x^3,x,0)

结果：

1003003001/6

程序：

symsx

limit((1001*x-sin(1001*x))/x^3,x,inf)

结果：

0

mx

x

1.2yecos，求y

1000

程序：

symsx

diff(exp(x)*cos(1001*x/1000),2)

结果：

-2001/1000000*exp(x)*cos(1001/1000*x)-1001/500*exp(x)*sin(1001/10

00*x)

word文档可自由编辑

11

22

xy

1.3计算edxdy

00

程序：

dblquad(@(x,y)exp(x.^2+y.^2),0,1,0,1)

结果：

2.13935019514228

x4

1.4计算dx

22

m4x

程序：

symsx

int(x^4/(1000^2+4*x^2))

结果：

1/12*x^3-1002001/16*x+1003003001/32*atan(2/1001*x)

x(10)

1.5yecosmx,求y

程序：

symsx

diff(exp(x)*cos(1000*x),10)

结果：

-1009999759158992000960720160000*exp(x)*cos(1001*x)-1009023999899

0319040000160032*exp(x)*sin(1001*x)

word文档可自由编辑

m

1.6给出x在x0的泰勒展式(最高次幂为4).

1000.0

程序：

symsx

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*1001

0^(1/2)*x^4

1.7Fibonacci数列{x}的定义是x1,x1，

n12

,xxx(n3,4,)用循环语句编程给出该数列的前20项（要求

nn1n2

将结果用向量的形式给出）。

程序：

x=[1,1];

forn=3:20

x(n)=x(n-1)+x(n-2);

end

x

结果：

Columns1through10

11235

813213455

Columns11through20

89