MATLAB数学实验报告.doc

《MATLAB数学实验报告.doc》由会员分享，可在线阅读，更多相关《MATLAB数学实验报告.doc（19页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、MATLAB数学实验报告MATLAB数学实验报告数学实验报告题目：matlab软件基本操作实验一.题目下面给出一些自然数的特点，请读者编程完成相应的任务。(1)水仙花数若一个三位自然数的各位数字的立方和等于该数本身，则称该自然数为花数，例如：153=13+53+33,所以153就是水仙花数,编程计算出所有的水仙花数.。(2)回文数对于一个自然数,若将各位数字到序排出,加到原来的数字上,反复这样多次，若能得到一个从左到右与从右到左读完全一样的数,则称该自然数能产生回文数或对称数,例如:195就可以产生一个回文数9339,因为:195+951=786786+687=14731473+3741=52

2、145214+4125=9339通过编程计算,你能找到多少个能产生回文数的数,又能找出多少个不能产生回文数的数?二者的最小数是多少?(3)完美数对于一个自然数而言,若它的真因子之和小于数本身,则称该数为盈数,若它的真因子之和大于数本身,则称该数为亏数,若它的真因子之和等于数本身,则称该数为完美数(完全数).例如:12是亏数,14是盈数,6是完全数.通过编程计算,你能找出多少个完美数?(完美数的两个特点是:它可以写成若干连续自然数之和;它的全部因子的倒数和都是2)(4)亲和数亲和数又叫友好数,是指两个自然数,其中每个自然数的真因子之和等于另一个数.例如:220与284就是一对亲和数.令人吃惊的事

3、实是:第二对最小的亲和数(1184,1210)是意大利的一个16岁是男孩在19世纪后期发现的.通过编程计算,你能算出多少对亲和数.二.问题分析（1）分别求出每一个三位数的百位，十位和个位上的数字，然后算出它们的立方和与原数进行比较。依照这种方法穷举出所有的三位数，就能求出所有的水仙花数。（2）利用循环将一个数回文数求解出来，然后将该数与原数进行比较，若不相同，则将该数与原数相加而得到新数，再重复上述操作，将新数的回文数求解出来直至得到的数与原数相同或数字超过范围（当得到的新数末尾数字为零时，该数的回文数无法表示，这时也认为初始数据无法产生回文数），进而根据结果判断原始数据是否为能产生回文数的数

4、。（3）对于一个自然数，可以通过检测所有比他小的正整数来确定出它所有的真因数，因此，利用穷举法求出某一个自然数的真因数后，再对其进行求和操作，将所得的结果与原数进行比较就可以确定出该数是否是完美数。（4）利用（3）中的结论，首先选取两个数，依次算出两个数的所有真因数的和，分别与另一个数进行比较，然后判断出所选取的两个数是否为亲和数。三.假设约定:（1）第一题三位数的范围为100到999。（2）第二题范围在10000以内。（3）第三题，约定查找范围为2到10000内的所有完美数。（4）第四题，约定范围为2到1500内的所有亲和数。四.模型分析：（1）第一题应用for语句，将100到999内的所有

5、三位数依次判断是否是水仙花数。（2）第二题运用多种系统函数结合分步实现要求范围内的回文数的查找。(3)第三题运用for语句嵌套及if语句判断，分别实现真因子的查找和完美数条件的判断。（4）第四题应用和第三题相同的方法找到真因子，并通过多重for循环的嵌套和if语句的配合实现目的。五.模型求解（1）问题一所用的代码如下：i=100:999j=mod(i,10);k=mod(i-j,100)/10;l=(i-j-k*10)/100;ifi=j3+k3+l3iend运行结果：i=153i=370i=371i=407（2）问题二所用的代码如下：a=input(enteranumber:);b=floo

6、r(log10(a)+1;v=mod(floor(a./10.(b-1:-1:0),10);w=fliplr(v);ifall(v=w)error(thisisahuiwennumber);endp=sort(v);q=-sort(-v);whilenot(all(v-w=0)n=polyval(q,10)+polyval(p,10);c=polyval(q,10);d=polyval(p,10);fprintf(%d+%d=%dn,c,d,n)b=floor(log10(n)+1;v=mod(floor(n./10.(b-1:-1:0),10);w=fliplr(v);p=sort(v);q

7、=-sort(-v);end（3）问题三所用的代码如下：sum=0;fori=2:10000forj=2:i0.5ifmod(i,j)=0sum=sum+j+i/j;endendifi=sum+1iendsum=0;end运行结果：i=6i=28i=96i=8128（4）问题四所用的代码如下：sum=0;fori=2:1500forj=i+1:1500forh=2:j0.5ifmod(j,h)=0sum=sum+h+j/h;endendifi=sum+1sum=0;forh=2:i0.5ifmod(i,h)=0sum=sum+h+i/h;endendifj=sum+1ijendendsum=0

8、;endend运行结果：i=220j=284i=1184j=1210六.实验总结体会这四个题目考察了matlab软件中for语句的灵活运用，以及各种语句之间的配合，体现了软件灵活丰富的编程功能。实验锻炼了我们分析实际问题、转化问题、并用所学解决问题的能力。各成员相互分工合作，协调配合，增进了彼此的了解与沟通。并且，通过实践我们发现了自己的许多问题，例如将实际问题转化为数学问题的能力，这也将在我们今后通过学习实践努力改进。扩展阅读：Matlab数学实验报告-1-MATLAB作业MATLAB与科学计算实验报告指导老师：易昆南实验时间：大二年级上学期学院：数学科学与计算技术学院专业班级：统计0902

9、姓名：丁克明学号：1304090112-2-MATLAB作业目录1.学生成绩管理32.函数极限问题53.蛛网模型74.水塔模型105.混沌问题126.koch曲线147.拉格朗日与三次样条插值以及曲线拟合178.解线性方程组209.矩阵对角化2310.矩阵标准化.2411.摸球实验2612.钉板问题2613.火车问题28-3-MATLAB作业一、学生成绩管理学号1304090112实验题目班级统计0902姓名丁克明指导教师评分易昆南学生成绩管理程序1、设计（实习）目的：1结合实际问题展现MATLAB在生活和学习方面的广泛应用2学会利用MATLAB编程并求解实际问题3学会并运用for循环和swi

10、tch结构，以及MATLAB中已有函数如sum4了解单元数组cell和结构数组struct的作用，学会创建这些数组.2、实验内容：编写一个学生成绩管理的函数程序：将学生成绩按五个等级划分：优秀（90以上），良好（80-90），中等（70-79），及格（60-69），不及格（60以下）对输入的学生成绩按五个等级划分后，打印学生姓名，得分，等级，并计算全班最高分、最低分、平均分。3、详细设计clearfori=1:10ai=89+i;bi=79+i;ci=69+i;di=59+i;endc=d,c;Name=input(pleaseinputname:);Score=input(pleaseinp

11、utscore:);n=length(Score);Rank=cell(1,n);S=struct(Name,Name,Score,Score,Rank,Rank);fori=1:nswitchS(i).Scorecase100S(i).Rank=满分;caseaS(i).Rank=优秀;casebS(i).Rank=良好;casecS(i).Rank=及格;casedS(i).Rank=不及格;endenddisp(学生姓名,得分,等级)；-4-MATLAB作业fori=1:ndisp(S(i).Name,blanks(6),num2str(S(i).Score),blanks(6),S(i

12、).Rank);ends=0;fori=1:ns=S(i).Score+s;endaverscore=s/n;t=S(1).Score;fori=1:(n-1)if(S(i).ScoreS(i+1).Score)m=S(i+1).Score;endenddisp(平均成绩);disp(averscore);disp(最高分);disp(t);disp(最低分);disp(m);4、实验结果：-5-MATLAB作业5、实验总结：这是我第一次接触matlab并运用matlab解决实际问题，在以前学习C+的时候，解决一个问题的代码会很繁琐。但相同的问题，在matlab中却变得很简单。这就是matla

13、b的强大之处。运用它可以方便的解决许多实际问题。知道了这一点，我决心以后会认真的学习这样一门课程。同时，它也给我们枯燥的数学公式的学习带来了动手解决实际问题的机会。是将理论与实际相结合的方法。让我们更加体会到知识的力量是强大的，我们应该更好的掌握科学技术和相关理论知识，并能够将其运用于实际生活当中。在解决这个问题的时候也确实遇到了一些问题。比如如何输入符号数组和数值数组的问题，当时一直不理解。当通过查询资料，和同学讨论，最后明白了。输入符号时，符号要用单引号，而数据就不需要了。通过这样一次实践，我更加明白了动手的必要性。只有理论知识是远远不够的。所以以后一定要加强自己的动手能力，勤动手。二、函

14、数极限问题学号1304090112实验题目班级统计0902姓名丁克明指导教师评分易昆南函数极限问题1、设计（实习）目的：1结合实际问题展现MATLAB在生活和学习方面的广泛应用2学会利用MATLAB编程并求解实际问题3学会并运用movie和moviein函数，以及plot结构。4理解getframe以及学会运用axis调节坐标抽2、实验内容：自选函数，运用plot进行绘图。并使用movie以及moviein，制作函数极限动画。3、详细设计clearm=moviein(16);holdonfori=1:50forj=1:it1=(j-1)*5;t2=j*5;t=t1:0.01:t2;axis(0

15、,40,-1,1)%限制动画的坐标显示大小plot(t,0,.r)plot(t,1./t,.)endm(:,i)=getframe;endmovie(m,50)clearm=moviein(16)-6-MATLAB作业holdonfori=1:10forj=1:it1=(j-1)*pi/10;t2=j*pi/10;t=t1:0.01:t2;axis(0,5*pi/2,-1,1)%限制动画的坐标显示大小plot(pi-0.1,t,.k)plot(pi-0.1,-t,.k)plot(pi+0.1,-t,.k)plot(pi+0.1,t,.k)plot(t,0,.r)plot(t,sin(t),.)

16、t1=(20-j)*pi/10;t2=(21-j)*pi/10;t=t1:0.01:t2;axis(0,5*pi/2,-1,1)%限制动画的坐标显示大小plot(t,0,.r)plot(t,sin(t),.)endm(:,i)=getframe;endmovie(m,50)4、实验结果：-7-MATLAB作业5、实验总结：使用matlab进行绘图觉得很有意思，简单的代码，却能绘出各种颜色的图形动画。在学习的过程中还是有很多的问题。还有很多方面的东西不够了解，只能边尝试边查询资料，让自己能够更了解，更明白。以便绘出更好的图形。三、蛛网模型学号1304090112实验题目班级统计0902姓名丁克明

17、指导教师评分易昆南蛛网模型1、设计（实习）目的：1.学会运用Matlab解决实际问题2.进一步掌握Matlab的一些基本操作3.通过范例体会Matlab的初步建模过程4学会用做动画2、实验内容：已知前两年的猪肉的产量和猪肉的价格分别为：39吨，28吨，12元/公斤，17元/公斤，根据前一年的猪肉价格影响后一年猪肉产量当年猪肉产量影响当年猪肉价格的线性关系，编写程序，利用动画原理画出前十年猪肉的产量价格的动态图形。（参数设置为，-8-MATLAB作业c1=39,c2=28,c3=36,r1=12,r2=17,k=16)。3、详细设计clear%c1为产量1,c2为产量2,c3为产量3,r1为%肉

18、价1,r2为肉价2,k为K年后产量与肉价%是否稳定holdoffc1=39;c2=28;c3=36;r1=12;r2=17;k=16;a1=c11;c21;%系数矩阵b1=r1,r2;%列向量a2=r11;r21;b2=c2,c3;a=a1b1;b=a2b2;%x0(1)=c1;forn=1:30y0(n)=a(1)*x0(n)+a(2);%a(1)a(2)为矩阵元素的引用，下行类似x0(n+1)=b(1)*y0(n)+b(2);x(n)=x0(n);y(n)=x0(n+1);endplot(x,y0,-g,y,y0,-b)holdonm=moviein(100);forn=1:kforj=1

19、:30t1=x0(n)+(j-1)*(x0(n+1)-x0(n)/30;t2=x0(n)+j*(x0(n+1)-x0(n)/30;ift2-9-MATLAB作业t2=t;elseendt=t1:0.01:t2;plot(x(n+1),t,.r)%划竖线endm(:,n)=getframe;endmovie(m,20)4、实验结果：ans=Columns1through639.000028.000036.000030.181834.413231.3358Columns7through1233.573931.946233.130032.269132.895232.4398Columns13thro

20、ugh1832.771032.530232.705332.577932.670632.6032Columns19through2432.652232.616632.642532.623632.637432.6274Columns25through3032.634632.629432.633232.630432.632432.6310Column3132.63205、实验图像：87.87.67.47.276.86.66.46.263535.53636.53737.53838.53939.5-10-MATLAB作业6、实验总结：通过做此实验，让我对MATLAB有更进一步的了解，学会怎样才能正确运用

21、MATLAB求解实际问题，了解如何利用数学模型去解释和分析社会经济问题，特别是这个典型经济问题的求解。我对MATLAB还不是特别熟悉，过程中遇到了很多问题，经过与同学讨论得到了良好的解决，希望以后会有进步。四、水塔模型学号1304090112实验题目班级统计0902姓名丁克明指导教师评分易昆南水塔问题1、设计（实习）目的：1掌握用MATLAB来求最值2学会在一个闭区间求最值3学会解决实际问题。2、实验内容：水塔问题：在地面上建有一座圆柱形水塔,水塔内部的直径为d,并且在地面处开了一个高为H的小门.现在要对水塔内部进行维修施工,施工方案要求把一根长为l(ld)的水管运到水塔内部.请问水塔的门高H

22、多高时,才有可能成功地把水管搬进水塔内?3、详细设计一：Clearfunctionx=lt523(l,d,h)k1=0;a=0l=7;d=3ymax=0;k2=0;b=pi/2;h=3;while(b-a)10(-2)k2=k2+1;m=0;a=0;ifk1=0n=ceil(b/h)-1elsen=ceil(b/h);endfori=1:n+1;-11-MATLAB作业x(i)=a+(i-1)*(b-a)/n;y(i)=l*sin(x(i)-d*tan(x(i);endfori=1:nify(i)ymaxymax=y(i);a=x(i-1);elsey(i)-12-MATLAB作业plot(x

23、,8*sin(x)-2*tan(x),-b)4、实验图像：5、实验总结：通过做此实验，让我对MATLAB有更进一步的了解，学会怎样才能正确运用MATLAB求解实际问题，了解如何利用数学模型去解释和分析社会经济问题，特别是这个典型经济问题的求解。我对MATLAB还不是特别熟悉，过程中遇到了很多问题，经过与同学讨论得到了良好的解决，希望以后会有进步。五、混沌问题学号1304090112实验题目1、设计（实习）目的：1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：利用迭代方法绘制混沌图形3、详细设计：x0=0.3;holdonfork=1:100班级统计

24、0902姓名丁克明指导教师评分易昆南混沌问题-13-MATLAB作业a=3+(k-1)/100;fori=1:300x0=a*x0*(1-x0);ifi100plot(a,x0,r.)endendend4、实验图像：5、实验总结：实际的绘图进一步让我体会到了matlab功能的强大。他能在实际生活中解决很多人工无法解决的复杂问题。通过做此实验，让我对MATLAB有更进一步的了解，学会怎样才能正确运用MATLAB求解实际问题，了解如何利用数学模型去解释和分析社会经济问题，特别是这个典型经济问题的求解。我对MATLAB还不是特别熟悉，过程中遇到了很多问题，经过与同学讨论得到了良好的解决，希望以后会有

25、进步。-14-MATLAB作业六、koch曲线学号1304090112实验题目1、设计（实习）目的：1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：绘制雪花曲线和其他一个美丽图形3、详细设计：1、雪花曲线p=00;100;n=2;A=cos(pi/3)-sin(pi/3);sin(pi/3)cos(pi/3);fork=1:4d=diff(p)/3;m=4*n-3;q=p(1:n-1,:);p(5:4:m,:)=p(2:n,:);p(2:4:m,:)=q+d;p(3:4:m,:)=q+d+d*A;p(4:4:m,:)=q+2*d;n=m;endpl

26、ot(p(:,1),p(:,2)axis(010010)班级统计0902姓名丁克明指导教师评分易昆南Koch曲线2、p=010;100;0-10;-100;010;n=5;A=0-1;10;fork=1:5d=diff(p)/3;m=4*n-3;q=p(1:n-1,:);p(5:4:m,:)=p(2:n,:);p(2:4:m,:)=q+d;p(3:4:m,:)=q+2*d+d*A;p(4:4:m,:)=q+2*d;n=m;endplot(p(:,1),p(:,2)axis(-1010-1010)3、花草树木-15-MATLAB作业p=50;510;n=2;A=cos(pi/3)-sin(pi/

27、3);sin(pi/3)cos(pi/3);B=cos(-pi/3)-sin(-pi/3);sin(-pi/3)cos(-pi/3);fork=1:4d=diff(p)/3;d1=d(1:2:n,:);m=5*n;q1=p(1:2:n-1,:);p(10:10:m,:)=p(2:2:n,:);p(1:10:m,:)=p(1:2:n,:);p(2:10:m,:)=q1+d1;p(3:10:m,:)=p(2:10:m,:);p(4:10:m,:)=q1+d1+d1*A;p(5:10:m,:)=p(2:10:m,:);p(6:10:m,:)=q1+2*d1;p(7:10:m,:)=p(6:10:m,

28、:);p(8:10:m,:)=q1+2*d1+d1*B;p(9:10:m,:)=p(6:10:m,:);n=m;endplot(p(:,1),p(:,2)axis(010010)4、实验图像：1、雪花曲线-16-MATLAB作业2、3、-17-MATLAB作业5、实验总结：通过做此实验，让我对MATLAB有更进一步的了解，学会怎样才能正确运用MATLAB求解实际问题，了解如何利用数学模型去解释和分析社会经济问题，特别是这个典型经济问题的求解。我对MATLAB还不是特别熟悉，过程中遇到了很多问题，经过与同学讨论得到了良好的解决，希望以后会有进步。七、拉格朗日插值与三次样条插值以及曲线拟合学号13

29、04090112实验题目班级统计0902姓名丁克明指导教师易昆南评分拉格朗日插值与三次样条插值以及曲线拟合1、设计（实习）目的：1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：拉格朗日插值和曲线拟合（拟合）用下面一组数据拟合c(t)abe0.0.2kt中的参数a，b，ktj1002021004005006007008009001000ci4.544.995.355.655.906.106.266.396.506.593、详细设计：1.拉格朗日插值x1=1:0.1:5;x0=1,3,5;y0=1,9,25;x=2,2.5,3,3.5,4,4.5;fo

30、ri=1:6s=0;fork=1:3p=1;forj=1:3ifj=kp=p*(x(i)-x0(j)/(x0(k)-x0(j);endend-18-MATLAB作业s=p*y0(k)+s;endy(i)=s;endplot(x1,x1.2,b,x,y,*r)2.三次样条插值x0=0:0.1:5;y0=x0.2;x=0:0.5:5;y=interp1(x0,y0,x,spline);plot(x0,y0,b,x,y,*r)3.曲线拟合f=inline(x(1)+x(2)*exp(0.02*x(3)*tdata),x,tdata);tdata=100:100:1000;cdata=4.54,4.9

31、9,5.35,5.65,5.90,6.10,6.26,6.39,6.50,6.59;x,JM=lsqcurvefit(f,111,tdata,cdata);xJMplot(tdata,cdata)4、实验图像：1.拉格朗日插值-19-MATLAB作业2.三次样条插值3.拟合-20-MATLAB作业5、实验总结：这次的试验个人觉得就拉格郎日的振荡比较难做，虽然写出来了一眼看的就很清楚，不过其思想还是找了挺多资料后在同学的帮助下才得以完成；至于其他的试验让我学习到了更多的知识，开阔了眼界，明白了多种方法，对于以后的解题思路有了更多的思考余地，这次试验中用到了拉格郎日插值法，分段线性插值，三次样条插

32、值，拟合的基本原理八、解线性方程组学号1304090112实验题目班级统计0902姓名丁克明指导教师评分易昆南解线性方程组1、设计（实习）目的：1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：1/比较用逆矩阵法、除法、克拉默法则解方程Ax=b的用时和误差2/求解方程组最多零解3、详细设计：1、A=rand(100)*1.e2;x=ones(100,1);b=A*x;ticy=inv(A)*b;tocerr=norm(y-x)res=norm(A*y-b)ticy=Ab;tocerr=norm(y-x)res=norm(A*y-b)tica=det(

33、A)fori=1:1-21-MATLAB作业B=A;B(1:100,i)=b;y(i)=det(B)/a;endtocerr=norm(y-x)res=norm(A*y-b)2、functionC=solution(A)l,u=lu(A);M1=;M2=;e=;f=;C=;r=rank(A);rf=size(A);n=rf(1,2);k=1;fori=1:rwhilek-22-MATLAB作业Elapsedtimeis0.004000seconds.err=5.9797e-012res=8.1628e-010Elapsedtimeis0.003000seconds.err=5.0572e-01

34、2res=2.1790e-011a=-1.3898e+225Elapsedtimeis0.083000seconds.err=5.0860e-012res=1.4895e-0092、-23-MATLAB作业5、实验总结：通过本次实验，我基本上掌握了MATLAB求矩阵的秩、行列式和逆矩阵的命令，并通过编程比较可知用逆矩阵法、除法、克拉默法则求解方程Ax=B时除法最省时，最精确；而且利用了以上方法求解给定的恰定，超定和欠定线性方程组，可以很快的求出结果；可见，利用MATLAB求解线性方程组是非常方便的，很有利于我们各方面涉及到这类问题的学习和工作。而由于本来对线性代数的知识就掌握得不是很好，所以做

35、起这样的题目来确实很吃力，不过使用matlab软件解决了计算烦琐的问题，所以相对来说，题目也基本上能得到解决。在做实验过程中，我对线性代数的知识又有了进一步的认识，而且对于它结合matlab在现实问题中的应用也有了初步的了解，这也是一个不小的收获吧。九、矩阵对角化学号1304090112实验题目1、设计（实习）目的：1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：矩阵对角化3、详细设计：functiony=trigle(A)y=1;c=size(A);ifc(1)=c(2)y=0;return;ende=eig(A);n=length(A);whi

36、le1ifisempty(e)return;endd=e(1);f=sum(abs(e-d)-24-MATLAB作业y=0;return;ende(find(abs(e-d)-25-MATLAB作业1、设计（实习）目的：1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：对矩阵进行标准化3、详细设计：A=1-11;-1-3-3;1-30P,D=eig(A)symsy1y2y3;y=y1;y2;y3f=y*D*y4、实验结果：A=1-11-1-3-31-30P=-0.0580-0.8018-0.5948-0.8554-0.26730.4437-0.514

37、70.5345-0.6703D=-4.87300000.00000002.8730y=y1y2y3f=-342905718471287/70368744177664*conj(y1)*y1+3196254119/19342813113834066795298816*conj(y2)*y2+3234691681855341/1125899906842624*conj(y3)*y35、实验总结：这个实验进一步让我感受到了matlab的实用性，以及其在数学无论是代数还是几何整个领域里的极大用处。个人觉得，学习数学专业的我们大家应该很好的学学这样一个多功能的软件。以便在以后的生活工作中游刃有余。-26

38、-MATLAB作业11.摸球实验学号1304090112实验题目1、设计（实习）目的：1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：不同情况下，计算条件摸球的概率。一共有是个球：三个黑球，七个白球。不放回摸球。1、第三次摸到黑球2、第三次才摸到黑球3、三次都摸到黑球3、详细设计：a=rand(1000000,3);a(:,1)=round(a(:,1)-0.2);a(:,2)=round(a(:,2)*0.9-0.2-0.1*(a(:,1)-1);a(:,3)=round(a(:,3)*0.8-0.2-0.1*(a(:,1)-1)-0.1*(a(

39、:,2)-1);fori=1:6b(i)=sum(a(1:10i,3)/(10i);m=a(1:10i,1)&a(1:10i,2)&a(1:10i,3);c(i)=sum(m)/10i;d=a(1:10i,1)&a(1:10i,2)&a(1:10i,3);e(i)=sum(d)/(10i);endbce4、实验结果：b=0.40000.22021.26400.30070.30130.3004c=0.10000.14000.14500.17560.17580.1752e=000.00300.00800.00790.00835、实验总结:第一次将软件和正在学习的概率论理论知识联系在一起。将理论在计

40、算机中模拟实现。觉得非常有意思。更加加强了我学习matlab的信心。虽然这个学期快要结束，虽然matlab课程快要结束，但我相信自己会继续学习这个软件。班级统计0902姓名丁克明指导教师评分易昆南摸球实验12.钉板问题学号1304090112实验题目1、设计（实习）目的：班级统计0902姓名丁克明指导教师评分易昆南钉板问题-27-MATLAB作业1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：钉板问题3、详细设计：clearclfm=100;n=5;y0=2;ballnum=zeros(1,n+1);p=0.5;q=1-p;fori=n+1:-1:

41、1x(i,1)=0.5*(n-i+1);y(i,1)=(n-i+1)+y0;forj=2:ix(i,j)=x(i,1)+(j-1)*1;y(i,j)=y(i,1);endendmm=moviein(m);fori=1:ms=rand(1,n);xi=x(1,1);yi=y(1,1);k=1;l=1;forj=1:nplot(x(1:n,:),y(1:n,:),o,x(n+1,:),y(n+1,:),.-);axis(-2n+20y0+n+1),holdonk=k+1;ifs(j)pl=l+0;elsel=l+1;endxt=x(k,l);yt=y(k,l);h=plot(xi,xt,yi,yt

42、);axis(-2n+20y0+n+1)xi=xt;yi=yt;endballnum(l)=ballnum(l)+1;ballnum1=3*ballnum./m;bar(0:n,ballnum1),axis(-2n+20y0+n+1)mm(i)=getframe;holdoffendmovie(mm,3)-28-MATLAB作业4、实验结果：5、实验总结：这次试验比较困难，但是同时也是让我们自己去探索，去找打解决办法，还是很有意义的。也感觉到模拟这个问题很实用，和现实联系的很紧密。13.火车问题学号1304090112实验题目班级统计0902姓名丁克明指导教师评分易昆南火车问题1、设计（实习）

43、目的：1.了解MATLAB在实际问题中的应用2.学会利用MATLAB做图并求解实际问题2、实验内容：一列火车从A站并往B站，某人每天赶往B站上这趟火车.他已经了解到火车从A站到B站的运行时间是均值为30分钟；标准差为2分钟的随机变量；火车在下午大约1点离开A站，离开时刻的频率分布如下；出发时刻频率午后1:000.7午后1:050.2午后1:100.-29-MATLAB作业此人到达B站的时刻的频率分布为时刻频率午后1l:280.3午后1:300.4午后1:320.2午后1:340.1问他能赶上火车的概率是多少?模拟求解答案。3、详细设计：functionhuocheabk=0;n=10000;fori=1:nr1=rand;r2=rand;x=-2*log(r1)(1/2)*cos(2*pi*r2);t2=2*x+30;r=rand;d=rand;ifr-30-MATLAB作业5、实验总结：这次的实验很难，有很多新颖的思路，我感觉系统模拟这个方法非常的有用，可以用来模拟很多数学上难于直接求出来的东西，而且是接近我们的生活，故有很大的用途,这个方法应该要学习好，对我们以后的工作有很大的用途。希望以后有时间能够继续对这门理论有更深入的研究。第 19 页 共 19 页