《南京邮电大学数学实验练习题参考答案.docx》由会员分享,可在线阅读,更多相关《南京邮电大学数学实验练习题参考答案.docx(9页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、南京邮电大学数学实验练习题参考答案第一次练习教学要求:熟练把握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-
2、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+10
3、03003001/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
4、)*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
5、?-?=?-?,求该矩阵的逆矩阵,特征值,特征向量,行列式,计算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?=?-此页面能否是列表页或首页?未找到适宜正文内容。
限制150内