定积分近似计算方法实验报告.pdfVIP

实验七定积分的近似计算

班级：统计0502学号：1303050206：欧钊锋

实验目的：

⚫学习、了解并掌握定积分近似计算的MeCarlo法（包括（1）随机投点法（2）均值

估计法）、梯形法、法

⚫对定积分的计算有进一步的认识和了解

⚫掌握上述几种方法的思想

实验理论：

MeCarlo法（包括（1）随机投点法（2）均值估计法）：

根据大数定律，“（x(i),y(i)）落在四分之一单位圆内“发生的频率k/n依概率收敛于（n

—∞）该发生的概率p。

梯形法、法：

在小区间上用矩形面积近似曲边梯形的面积，可取左端点函数值为小矩形的高，或取右端点

函数值为小矩形的高，于是在整个区间（a,b）构成台阶形，容易知道两个台阶的面积，将

两者平均，则每个小区间上的小矩形变为小梯形，整个区间上的结果即为梯形求积。为提高

精度可用分段二次插值函数作为f(x)的近似，由于每段要用到相邻两个小区间端点的三个函

数值，所以小区间的数目必须是偶数，在第k段的两个小区间上用三个节点做二次插值，然

后积分，求m段之和即为所求积分。

实验内容：

一正方形里有四分之一单位圆，如果向边长为1的正方形里随机投n块小石头，当n很大时

小石头会均匀的分布在正方形中，数一下落在四分之一圆里的小石头，假定有k个，那么

k/n就能看作是四分之一单位圆面积pi/4的近似值，于是有pi=4*k/n显然这可以看做近似计

算pi的法。试分别用MeCarlo法（包括（1）随机投点法（2）均值估计法）、梯

形法、法作其近似计算。

程序设计：

用（MeCarlo法）随机投点法作其近似计算。

n=0;

fori=1:1000000

s=rand(1,2);

x=s(1);

y=s(2);

ify=sqrt(1-x^2)

n=n+1;

end

end

s=n/1000000

pi=4*s

实验结果：

s=

0.7853

pi

3.1414

用（MeCarlo法）均值估计法作其近似计算。

程序设计：

m=0;

forn=1:1000000

x=rand;

f=sqrt(1-x^2);

m=m+f;

end

s=m/1000000

pi=4*s

实验结果：

s

0.7852

pi

3.1407

用梯形法作其近似计算。

程序设计：

functions=ozftixing(f,a,b,M)%在命令窗口输入ozftixing(ozf111,0,1,10000)即得结果

h=(b-a)/M;%feval函数的作用是调用定义的f函数在x初的值

s=0;

fork=1:(M-1)

x=a+h*k;

s=s+feval(f,x);

end

s=h*(feval(f,a)+feval(f,b))/2+h*s

pi=4*s

functiony=ozf111(x)

y=sqrt(1-x^2);

实验结果：

s

0.7854

pi

3.1416

用法作其近似计算。

程序设计：

functions=ozfxinpusen(f,a,b,M)