南京邮电大学数学实验练习题参考答案.docx
南京邮电大学数学实验练习题参考答案第一次练习教学要求：熟练把握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计算30sinlimxmxmxx-与3sinlimxmxmxx-程序：symsxlimit(1001*x-sin(1001*x)/x3,x,0)结果：1003003001/6程序：symsxlimit(1001*x-sin(1001*x)/x3,x,inf)结果：01.2cos1000xmxye=，求''y程序：symsxdiff(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计算221100xyedxdy+?程序：dblquad(x,y)exp(x.2+y.2),0,1,0,1)结果：2.139350195142281.4计算4224xdxmx+?程序：symsxint(x4/(10002+4*x2)结果：1/12*x3-1002001/16*x+1003003001/32*atan(2/1001*x)1.5(10)cos,xyemxy=求程序：symsxdiff(exp(x)*cos(1000*x),10)结果：-1009999759158992000960720210000*exp(x)*cos(1001*x)-10090239998990319040000160032*exp(x)*sin(1001*x)1.60x=的泰勒展式(最高次幂为4).程序：symsxtaylor(sqrt(1001/1000+x),5)结果：1/100*10010(1/2)+5/1001*10010(1/2)*x-1250/1002001*10010(1/2)*x2+625000/1003003001*10010(1/2)*x3-390625000/1004006004001*10010(1/2)*x41.7Fibonacci数列nx的定义是121,1xx=，12,(3,4,)nnnxxxn-=+=用循环语句编程给出该数列的前20项要求将结果用向量的形式给出。程序：x=1,1;forn=3:20x(n)=x(n-1)+x(n-2);endx结果：Columns1through1011235813213455Columns11through208914423337761098715972584418167651.8对矩阵211020411000Am?-?=?-?，求该矩阵的逆矩阵，特征值，特征向量，行列式，计算6A，并求矩阵,PDD是对角矩阵，使得1APDP-=。程序与结果：a=-2,1,1;0,2,0;-4,1,1001/1000;inv(a)0.50100100100100-0.00025025025025-0.5005005005005000.5000000000000002.00200200200200-0.50050050050050-1.00100100100100eig(a)-0.49950000000000+1.32230849275046i-0.49950000000000-1.32230849275046i2.00000000000000p,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.0000a611.968013.0080-4.9910064.0000019.9640-4.9910-3.01001.9作出如下函数的图形注：先用M文件定义函数，再用fplot进行函数作图：1202()12(1)12xxfxxx?=?-此页面能否是列表页或首页？未找到适宜正文内容。
