概率求积分面积.docx

Xxxx实践6实践名称：用概率算法计算定积分问题的实验时间安排：3学时一、实验目的通过上机实验，要求用概率算法计算定积分问题的描述、算法设计思想、程序设计和算法复杂性分析等。二、实验环境装有Visual C＋＋6.0的计算机。本次实验共计3学时。三、实验内容1、熟悉概率算法思想掌握如何进行问题的算法表示和描述。掌握计算定积分问题描述。2、掌握如何编译程序理解编译过程中的错误信息，并掌握如何排错。3、掌握如何对算法进行时间复杂度的分析掌握时间复杂度的概念和用渐进符号进行时间复杂度的评估。4、实验题：求曲线 和直线x=0、x= 及y=0所围成图形面积A。要求：实现该实验结果。通过该实验题，解决计算定积分问题。一、算法分析产生一个随即点(即两个随机数，一个作为x，另一个作为y)，我们统计其落在指定范围内的频率来估算面积。以估算一个圆的面积为例： Figure 1如上图一个矩形中心在原点，边长为2,且边于坐标轴平行。我们产生一个随机点：生成一个随机数范围在 [-1,1]作为这个点的x坐标，同样的方法生成另一个随机数。这个点必在这个正方形内。重复此步，生成足够多的点。判断这个点是否在圆内：若X^2+Y^2=1,则在园内，统计所有在园内的点个数记为s1.若X^2+Y^2 1,则不在圆内，但仍在正方形内，统计不在圆内的点个数记为s2.3.通过统计估算圆的面积S_arer,正方形的面积记为S_square。 (1)4.同理可求sinx函数面积Figure 2x随机数的范围在[0,π]y随机数的范围在[0, 1]ysin(x),则统计为有效点，否则无效。二、代码清单#includeiostream#include stdlib.h#includemath.h#include windows.h#include time.husing namespace std;int main(){srand( (unsigned)time(NULL)); //以时间为种子产生随即数int maxN //随机数试验次数int top=0; //记录在sin 上方不记为面积int bot=0; //在sin函数下面，记为面积double PI = 2*acos(0); //求根据反三角函数求PI的值//printf(%f,PI);int MIN=0;int MAX=1000; double s=0; for (int i = 0; i maxN; i++){ double dy=(double)(rand() % (MAX + 1 - MIN) + MIN )/MAX;//积分区间y的范围 0-1 产生随机数 double dx=(double)(rand() % (MAX + 1 - MIN) + MIN )/MAX*PI;//积分区间x取值范围0-pi内，产生随机数 double sinx=sin(dx); if(dy=sinx)//在sin函数值下面，算为面积 { bot++; } else//在上面，不算面积 { top++; } s=PI*bot/(bot+top);//根据面积求出面积 }printf(随机%d次： 面积： %f\n\n,maxN, s);}三、实验结果Figure 3Figure 4Figure 5Figure 6四、结果分析我们发现随着试验次数增多，实验结果逼近真实结果2，但是任然不是2。猜测这是因为计算机生成的是伪随机数列而不是真随机。我们生成0-1000的随机数实验结果如下：Figure 7所以本次试验结果随着实验次数增多之后不会再逼近真实结果，原因是rand() 函数为伪随机数列。