数学实验报告4.docVIP

数学实验报告 实验序号： 日期：2013年5月30日 学号 姓名：陈楠楠 班级：11数本 实验名称：概率与频率 实验目的：通过对本实验的学习，加深频率和概率等概念的 理解和认识，并帮助掌握一些概率统计的 原理. 实验内容：1、通过实验，填写完成表格2~6的数据. 2、用Monte Carlo方法求两平面曲线与及轴所围成的区域面积. 下面是与之相关的两个Matlab程序： 程序甲 程序乙 n=100;a=sqrt(2)/2;b=1;S=a*b; m=0; for i=1:n if rand(1)*b=(rand(1)*a)^2rand(1)*b1-(rand(1)*a)^2 m=m+1; end end m/n*S quad(1-2*x.^2,0,a) n=100;a=sqrt(2)/2;b=1;S=a*b; m=0; for i=1:n x=rand(1)*a;y=rand(1)*b; if y=1-x^2yx^2 m=m+1; end end m/n*S quad(1-2*x.^2,0,a) 实验值（=0.4911）与实际值（=0.4714） 相距甚远（偏大） 实验值（=0.4719）与实际值（=0.4714） 十分接近 试分析【程序甲】和【程序乙】的不同之处.试问：哪一个程序是对的？为什么？ 实验程序： 第1题（1）、 n=1000;m=0; for i=1:n; y=unidrnd(2)-1; if y==0; m=m+1; end end m/n 第1题（2）、 n=100;m1=0;m2=0;m3=0;m4=0;m5=0;m6=0; for i=1:n y=unidrnd(6); switch y case 1 m1=m1+1; case 2 m2=m2+1; case 3 m3=m3+1; case 4 m4=m4+1; case 5 m5=m5+1; otherwise m6=m6+1; end end disp([num2str(m1/n),,,num2str(m2/n),,,num2str(m3/n),,,num2str(m4/n),,,num2str(m5/n),,,num2str(m6/n)]) 第1题（3）、 n=100;a=2;m=0; for i=1:n x=rand(1)*a/2;y=rand(1)*a/2; if x^2+y^2=(a/2)^2; m=m+1; end end 4*m/n 第1题（4）、 n=100;l=0.5;m=0;d=1; for i=1:n alpha=rand(1)*pi;y=rand(1)*d/2; if y=l/2*sin(alpha) m=m+1; end end m/n 2*n*l/(m*d) 第1题（5）、 %%实验值计算 n=1000;p=0;m=50; for t=1:n a=[];q=0; for k=1:m b=randperm(365); a=[a,b(1)]; end c=unique(a); if length(a)~=length(c); p=p+1; end end p/n %%理论值计算 p1=1:365; p2=1:(365-m); p2=[p2,ones(1,m)*365]; pp=p1./p2; pp=1-prod(pp) 第2题、 顺序问题（2）( 显示的结果 n=100000;a=sqrt(2)/2;b=1;S=a*b; m=0; for i=1:n if rand(1)*b=1-(rand(1)*a)^2rand(1)*b(rand(1)*a)^2