蒙特卡罗方法概率拓展1.ppt

%模拟投币 format long %只决定显示的位数 clear %清除 N=100; %试验次数 n=0; %正面的次数 p=0.5; %正面的概率 for i=1:N %循环 x=rand; %产生随机数 if xp, n=n+1; %满足件计一次数 end end y=n/N %统计平均值，频率 抛硬币试验 用蒲丰投针法在计算机上计算π值，取a=4、l=3。 a=4; % 设置两条平行线之间的距离 l=3; % 投针的长度 n=0; % 针与平行线相交的次数 N % 投掷次数 x=unifrnd(0,a/2,1,N); %产生n个(0,a/2)之间均匀分布的随机数 f=unifrnd(0,pi,1,N); % 产生n个(0,pi)之间均匀分布的随机数， for i=1:N if x(i)l*sin(f(i))/2 % 只要x小于l*sin(phi(i))/2,则相交 n=n+1; %计数 end end Pi=2*l*N/(a*n) % 从相交的频率总求的pi %运行结果 *注意函数的调用：unifrnd(A,B,M,N) 产生数值区间从A到B的随机数组，M行N列 作业 1、用蒙特卡罗方法模拟射击打靶问题 2、用计算机模拟计算，求连续掷两颗骰子，点数之和大于6且第一次掷出的点数大于第二次掷出点数的概率。 * MATLAB启动界面： MATLAB入门 常量、变量 数据类型 运算符与操作符 工作区操作命令 常用函数 常量与变量 常量 MATLAB使用下列名字表示常量： exp 指数e i或j 虚数单位 pi π inf 正无穷大 eps 无穷小的正数或机器的浮点运算误差限 nan 不定式 lasterr 存放必威体育精装版的错误信息 lastwarn 存放必威体育精装版的警告信息 应避免取与常量同名的变量。 常量与变量 变量 MATLAB变量命名的规则： （1）第一个字母必须是英文字母 （2）字母间不可留空格 （3）最多只能有31个字符，MATLAB会忽略多馀字母 （4）大小写不同 MATLAB变量可以通过赋值方式来定义。 方式1：赋值变量=赋值表达式 例如：b=0.2222,则定义一个变量b，它的值为0.2222。 变量可以通过重新赋值方式来覆盖。 方式2：[返回变量列表]=函数名(输入变量列表) 例如：[u,v] =max(X),返回数组X的最大值u和位置v。 数据类型 数值型数据、字符型数据、符号型数据和类与对象等格式。 数值量和数据格式显示 在默认情况下，当结果是整数时，MATLAB将它作为整数显示；当结果是实数，MATLAB以小数点后4位的精度近似显示。如果结果中的有效数字超出了这一范围，MATLAB以科学计数法来显示结果。 数据类型 常用命令： 说明 例：a=1/3 format short 显示小数点后4位（默认值） 0.3333 format long 显示14位 0.33333333333333 format bank 显示小数点后2位 0.33 format + 显示+，-，0 + format short e 5位科学记数法 3.3333e-001 format long e 16位科学记数法 3.333333333333333e-001 format rat 最接近的有理数显示 1/3 数据类型 字符型 在MATLA