《南邮 数学实验》由会员分享，可在线阅读，更多相关《南邮 数学实验（11页珍藏版）》请在金锄头文库上搜索。
1、第一次练习题1.求解下列各题：1) limit(1834*x-sin(1834*x)/(x3)ans =3084380852/3 2) diff(exp(x)*cos(1834*x/1000.0),10)ans =-426276065704904857519435573449/976562500000000000000000000*exp(x)*cos(917/500*x)+296986972130617301707171403/195312500000000000000000*exp(x)*sin(917/500*x) 3）int(x4/(18342+4*x2),x)ans =1/12*x3
2、-840889/4*x+771095213/4*atan(1/917*x)4）将在展开(最高次幂为8). taylor(sqrt(1834/1000.0+x),9,x)ans =1/50*4585(1/2)+5/917*4585(1/2)*x-625/840889*4585(1/2)*x2+156250/771095213*4585(1/2)*x3-48828125/707094310321*4585(1/2)*x4+2441406250/92629354652051*4585(1/2)*x5-915527343750/84941118215930767*4585(1/2)*x6+251770
3、0195312500/545237037828059593373*4585(1/2)*x7-1022815704345703125/499982363688330647123041*4585(1/2)*x82.求矩阵 的逆矩阵 及特征值和特征向量。逆矩阵： A=-2,1,1;0,2,0;-4,1,1834;inv(A)ans = -0.5005 0.2501 0.0003 0 0.5000 0 -0.0011 0.0003 0.0005特征值： A=-2,1,1;0,2,0;-4,1,1834;eig(A)ans = 1.0e+003 * -0.0020 1.8340 0.0020特征向量：
4、A=-2,1,1;0,2,0;-4,1,1834;P,D=eig(A)P = -1.0000 -0.0005 0.2425 0 0 0.9701 -0.0022 -1.0000 0.0000D = 1.0e+003 * -0.0020 0 0 0 1.8340 0 0 0 0.00203.已知分别在下列条件下画出的图形：,分别为(在同一坐标系上作图); x=-2:1/50:2;y1=1/(sqrt(2*pi)*1834/600)*exp(-x.2/(2*(1834/600)2);y2=1/sqrt(2*pi)*1834/600*exp(-(x+1).2/(2*(1834/600)2);y3=1
5、/sqrt(2*pi)*1834/600*exp(-(x-1).2/(2*(1834/600)2);plot(x,y1,x,y2,x,y3)4.画 （1） t=0:pi/1000:20;u=0:pi/10000:2;x=u.*sin(t);y=u.*cos(t);z=100.*t/1834;plot3(x,y,z) (2) ezmesh(sin(1834*x*y),0,3,0,3)（3） t,u=meshgrid(0:.01*pi:2*pi,0:.01*pi:2*pi);x=sin(t).*(1834./100+cos(u);y=cos(t).*(1834./100+cos(u);z=sin(
6、u);surf(x,y,z)5对于方程，先画出左边的函数在合适的区间上的图形，借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。最后写出你做此题的体会. subplot(2,1,1);ezplot(x5-1834/200*x-.1);subplot(2,1,2);ezplot(x5-1834/200*x-.1);axis(-4 4 -.1 .1);solve(x5-1834/200*x-.1=0) ans = -1.7379128266700894611391209979802 -.1089
7、3246204072582693798581490422e-1 .27232448901993127689441860269764e-2-1.7406575085305563272476413157894*i .27232448901993127689441860269764e-2+1.7406575085305563272476413157894*i 1.7433595830937634182950312074166第二次练习题1. 统计1到以内可以写为两个素数之和的偶数与奇数的个数.function vector=judge(x); vector=; i=2; while isempty(
8、vector)&(i format long;x=1;stopc=1;eps=1e-17;n=1;p=7+1834/1000;while stopceps x=x+1/np;n=n+1;stopc=1/np;end nxn = 85x = 2.002257920470133.设，数列是否收敛？若收敛，其值为多少？精确到6位有效数字。 format long;x=3;eps=1e-6;n=1;stopc=1;while x-stopceps stopc=x; x=(x+1834/x)/2; n=n+1;endnxn = 3x = 1.565686833062037e+0024.能否找到分式函数以
9、及分式函数,使它产生的迭代序列收敛到(对于为整数的学号，请改为求。收敛时要求精确到17位有效数字。有一个要求：必须全部是整数)?并研究如果迭代收敛,那么迭代的初值与收敛的速度有什么关系.写出你做此题的体会.function y=fex2_4_1(x,a,b,c,d,e) y=(a*x2+b*x+c)/(d*x+e);x=-10000; answer=; vpa(x,17) for a=-10000:10000 for b=-10000:10000 for c=-10000:10000 for d=-10000:10000 for e=-10000:100000 if vpa(fex2_4_1(
10、x,a,b,c,d,e),17)=x answer=x,a,b,c,d,e end end end end end end第三次练习题练习4A=4,2,;1,3; a=rand; b=rand; x=a;b t=; for i=1:40 x=A*x; t(i,1:2)=i,x(1,1)/x(2,1); end x = 0.95012928514718 0.23113851357429t = 1.00000000000000 2.59365859734352 2.00000000000000 2.21226129089303 3.00000000000000 2.08144691106100 4
11、.00000000000000 2.03205658249965 5.00000000000000 2.01274094675769 6.00000000000000 2.00508342517318 7.00000000000000 2.00203130487201 8.00000000000000 2.00081219198689 9.00000000000000 2.00032482403086 10.00000000000000 2.00012992117204 11.00000000000000 2.00005196711849 12.00000000000000 2.0000207
12、8663135 13.00000000000000 2.00000831461797 14.00000000000000 2.00000332584166 15.00000000000000 2.00000133033578 16.00000000000000 2.00000053213417 17.00000000000000 2.00000021285365 18.00000000000000 2.00000008514145 19.00000000000000 2.00000003405658 20.00000000000000 2.00000001362263 21.00000000000000 2.00000000544905 22.00000000000000 2.00000000217962 23.00000000000000 2.00000000087185 24.00000000000000 2.00000000034874 25.00000000000000 2.00000000013950 26.00000
