南邮 数学试验答案（全）

本文格式为Word版，下载可任意编辑——南邮数学试验答案（全）专业姓名学号成绩

第一次练习

教学要求：熟练把握Matlab软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab软件解决微积分、线性代数与解析几何中的计算问题。

补充命令

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

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

在下面的题目中m为你的学号的后3位（1-9班）或4位（10班以上）1.1计算limmx?sinmxmx?sinmxlim与33x?0x??xx程序：

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

1.2y?ecosxmx，求y''1000程序：symsx

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

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

数学试验试验报告

1

专业姓名学号成绩

1.3计算

??1100ex2?y2dxdy

程序：

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

2.13935019514228

x4dx1.4计算?22m?4x程序：symsx

int(x^4/(1000^2+4*x^2))结果：

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

1.5y?excosmx,求y(10)

程序：symsx

diff(exp(x)*cos(1000*x),10)结果：

-1009999759158992000960720230000*exp(x)*cos(1001*x)-10090239998990319040000160032*exp(x)*sin(1001*x)

数学试验试验报告

2

专业姓名学号成绩

1.6给出m?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/1002023*10010^(1/2)*x^2+625000/1003003001*10010^(1/2)*x^3-390625000/1004006004001*10010^(1/2)*x^4

1.7Fibonacci数列{xn}的定义是x1?1,x2?1，

,xn?xn?1?xn?2(n?3,4,)用循环语句编程给出该数列的前20项（要求

将结果用向量的形式给出）。程序：x=[1,1];forn=3:20

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

结果：

Columns1through10

11235813213455

Columns11through20

891442333776109871597258441816765

数学试验试验报告

3

专业姓名学号成绩

????211???A??020?1.8对矩阵，求该矩阵的逆矩阵，特征值，特

??m??41?1000??6征向量，行列式，计算A，并求矩阵P,D（D是对角矩阵），使得

A?PDP?1。

程序与结果：

a=[-2,1,1;0,2,0;-4,1,1001/1000];inv(a)

0.50100100100100-0.00025025025025-0.5005005005005000.5000000000000002.00202300202300-0.50050050050050-1.00100100100100eig(a)

-0.49950000000000+1.32230849275046i-0.49950000000000-1.32230849275046i2.00000000000000[p,d]=eig(a)p=

0.3355-0.2957i0.3355+0.2957i0.2425000.97010.89440.89440.0000注：p的列向量为特征向量d=

-0.4995+1.3223i000-0.4995-1.3223i0002.0000a^6

11.968013.0080-4.9910064.0000019.9640-4.9910-3.0100

数学试验试验报告

4

专业姓名学号成绩

1.9作出如下函数的图形（注：先用M文件定义函数，再用fplot进行函数作图）：

1?2x0?x???2f(x)??

1?2(1?x)?x?11??20.9函数文件f.m:0.80.7functiony=f(x)

0.6if0<=x

0.4else1/2<x

0.2end

0.1

000.10.20.30.40.50.60.70.80.91程序：fplot(@f,[0,1])

1.10在同一坐标系下作出下面两条空间曲线（要求两条曲线用不同的颜色表示）

?x?cost?x?2cost??（1）?y?sint（2）?y?2sint

?z?t?z?t??程序：

t=-10:0.01:10;x1=cos(t);