matlab数学实验复习题有答案Word格式文档下载.docx

1、（做不出来）二、写出运行结果：1、eye（3，4）=2、size（1,2,3）=1；33、设b=round（unifrnd（-5，5，1，4），则=3 5 2 -5x,m=min（b）；x=-5；m=4,x,n=sort（b）-5 2 3 54 3 1 2mean（b）=1.25，median（b）=2.5，range（b）=104、向量b如上题，则any（b）,all（b2）,all（b5 6;7 87 8;5 6=6、若，则7、diag（diag（B）=8、4:-2:1.*-1,6=-4 129、acos（0.5）,atan（1）ans=.0471*0.78539816339744810、

2、norm（1,2,3）Ans=3.74165738677394111、length（1，3，-1）=312、x=0:0.4:2;plot（x,2*x,k*）13、zeros（3,1）;14、ones（3）=，vander（2,3,5）=16、floor（1：0.3：3）=1 1 1 1 2 2 218、subplot（2,2,1）; fplot（sin,0,2*pi）;subplot（2,2,2）;plot（1,2,-1）;x=linspace（0,6*pi）;subplot（2,2,3）;plot3（cos（x）,sin（x）,x）;subplot（2,2,4）;polar（x,5*sin（

3、4*x/3）;19、t=linespace（0，2，11）0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.020、a,b=binostat（15,0.2）a=3 b=2.4y1=binopdf（5,10,0.7）=0.1029,y2=binocdf（5,10,0.7）=0.150321、log10（1,10,100）=0 1 222、p=1;for k=2:3:9 p=p*k;end;p p=8023、s=0;9 s=s+k;s s=1524、Ans=3.864725、a1=norminv（0.6,3,4）a1=4.013426、unifinv（0.4,1,

4、5）,unifpdf（0.4,1,5）,unifpdf（2,1,5）Ans=2.6 0 0.2527、A=0 1-1;2 1 0;1-1 1;0 1 -11 -1 1A（1,3,:）A（3,1,:）=1 -1 1 0 1 -1A（2,:）=2 1 0-2*A（1,:）= 0 -2 228、quad（sin（x）,0,pi/2）=1.000029、trapz（3,4,6,1,2,3）=6.500030、 int（x-sin（x）,0,1）Ans=cos（1） - 1/231、round（3:5）,ceil（3:5）;floor（3:5）3 3 4 4 5 53 3 3 4 4 5limit（1+

5、1/（3*x）x,inf）=1diff（sin（3*x）+x3,2）=6*x-9*sin（3*x） taylor（exp（3*x）,5,1）:命令输入： y=taylor（exp（3*x）,x,1,Order,5） exp（3） + 3*exp（3）*（x - 1） + （9*exp（3）*（x - 1）2）/2 + （9*exp（3）*（x - 1）3）/2 + （27*exp（3）*（x - 1）4）/8a1=mod（15,4）,b1=rem（15,4）=3,3a2=mod（-15,-4）,b2=rem（-15,-4）=-3,-3a3=mod（15,-4）,b3=rem（15,-4）=-1

6、,-3a4=mod（-15,4）,b4=rem（-15,4）=1,-334、x=binornd（20,0.4,2,4）8 7 10 810 7 9 12sign（x）,1 1 1 1y=-poissrnd（8,2,4）-16 -10 8 -7-7 -8 -6 -9sign（y）-1 -1 -1 -135、a1,b1=binostat（20,0.4） a1=8 b1=4.8a2,b2=poisstat（8）ans=8,8a3,b3=chi2stat（15）ans=15 3036、运行M文件：chi2fign=5;a=0.9;xa=chi2inv（a,n）;0.1:15;y=chi2pdf（x,n

7、）;plot（x,y,b）;hold on;xf=0:xa;yf=chi2pdf（xf,n）;fill（xf,xa,yf,0,gtext（xa*1.01,0.005,num2str（xa）;text（2.5,0.05,alpha=0.9,fontsize,20）;text（9,0.09,Xchi2（4）,16）;37、t=linspace（0,2*pi）;polar（t,3*t,g*）38、quadl（exp（2*x）.*log（3*x）,1,3）ans = 398.635239、x0=0:2*pi/6:2*pi;y0=sin（x0）.*cos（x0）;x=linspace（0,2*pi,10

8、0）;y=sin（x）.*cos（x）;y1=spline（x0,y0,x）;x;y;y1k,x,y1,b-注：此处省略100组数据40、A=round（unifrnd（0,100,3,3）;L,U=lu（A）L = 0.9897 0.4699 1.0000 0.1649 1.0000 0 1.0000 0 0U = 97.0000 80.0000 92.0000 0 35.8041 26.8247 0 0 -89.656841、a=sparse（1 3 3,2 3 5,1 2 3,4,5）;s=full（a）s = 0 1 0 0 0 0 0 0 0 0 0 0 2 0 3三、编程1、分别用

9、矩形公式、梯形公式、辛普森公式、Gauss-Lobatto公式及随机模拟方法计算数值积分，并与符号运算计算的结果进行比较。format long0.01:pi/2;y=exp（3*x）.*sin（2*x）;s1=sum（y）*0.01;s2=trapz（x,y）;s3=quad（exp（3*x）.*sin（2*x）,0,pi/2）;s4=quadl（n=10000;x=unifrnd（0,pi/2,1,n）;y=unifrnd（0,exp（5.5）,1,n）;k=0;for i=1:n if y（i）=1disp（迭代不收敛for k=1:x=B*x+f;x;if norm（A*x-b）tol

10、breakm=k;高斯-赛德尔迭代M文件；function x,m=ga（A,b,x0,tol,n）B=（D-L）U;f=（D-L）b;end 3、用欧拉方法和龙格库塔方法求下列微分方程初值问题的数值解:向前欧拉M文件：function z=foeula（f,a,b,y0,h）m=floor（b-a）/h）;x（1）=a;y（1）=y0;for n=1:m x（n+1）=x（1）+n*h; y（n+1）=y（n）+h*feval（f,x（n）,y（n）;z=y;改进欧拉M文件：function z=adveula（f,a,b,y0,h）;x=a:h:b; k1=feval（f,x（n）,y（n

11、）; k2=feval（f,x（n+1）,y（n）+h*k1）; y（n+1）=y（n）+h*（k1+k2）/2;函数调用M文件：function dy=ode121（x,y）dy=x2+y2;Returnz1=foeula（ode121,0,1,1,0.1）z2=adeveula（x,y=ode45（,0:1,1）或者用直接用“inline”符号函数 z1=foeula（inline（x2+y2）,0,1,1,0.1）z2=adeveula（inline（x,y=ode45（inline（）,0:4、用牛顿切线法求的根，要求相对误差不超过，并输出解和迭代次数。function x,m=new

12、ton（f,df,x0,n,tol）x（1）=x0;x（k+1）=x（k）-feval（f,x（k）/feval（df,x（k）;if abs（x（k+1）-x（k）/x（k） a（k+1）=（1+r）*a（k）-1000;plot（a）;grid;n=fix（k/12）;m=mod（k,12）;a=adisp（每岁末养老金的余额为：（1:n）,a（12:12:n*12）disp（老人养老金用完时的年龄为：）;time=num2str（n+60）, 岁 ,num2str（m）, 月B=0;p=1000;240 p=p/（1+r）;B=B+p;disp（若老人养老金想用到80岁，每月取1000元

13、，则60岁时应存入银行的钱为:）;money=num2str（B）,元在命令窗口键入文件名，计算机运行结果为： ex2xt2每岁末养老金余额为：ans = 1.00 92639.47 2.00 84917.75 3.00 76817.13 4.00 68319.02 5.00 59403.89 6.00 50051.30 7.00 40239.78 8.00 29946.80 9.00 19148.74 10.00 7820.82基金用完时老人的年龄为：time =70 岁 8 个月若老人想用到80岁，每月取1000元，则老人60岁时应存入的钱为：money =154093.3029 元三、结果分析：老人60岁时存入10万元，每月取1千元，那么70岁零8个月时钱用完；若想用到80岁，则60岁时应存入15409.30元。10、由某商店过去的销售记录知道，某种商品每月的销售数可以用参数的泊松分布来描述，为了有95%以上的把握不使商品脱销，问商店在每月月底应进该种商品多少件？y=poissinv（0.95,25）y=33P85实验练习：第三题