原创力文档- 1、本文档共26页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
定积分Monte Carlo计算方法实现课件
《计算机高级语言》 认知实习报告 ; ;一、设计任务与要求;;二、Monte Carlo算法引导 ;;三、Monte-Carlo思想计算定积分的实现;1.定积分的定义;;2、Monte Carlo积分法的原理;四、详细设计及实现;五 测试;n;107;六 总结;;附录：源程序;;double ran2(long *idum) /*产生随机数*/
{
int j;
long k;
static long idum2=123456789;
static long iy=0;
static long iv[NTAB];
double temp;
if(*idum=0); {
if(-(*idum)1)
*idum=1;
else
*idum=-(*idum);
for(j=NTAB+7;j=0;j--)
{
k=(*idum)/IQ1;
*idum=IA1*(*idum-k*IQ1)-k*IR1;
if(*idum0)
*idum+=IM1;
if(jNTAB)
iv[j]=*idum;
}
iy=iv[0];
};k=(*idum)/IQ1;
*idum=IA1*(*idum-k*IQ1)-k*IR1;
if(*idum0)
*idum+=IM1;
k=idum2/IQ2;
idum2=IA2*(idum2-k*IQ2)-k*IR2;
if(idum20)
idum2+=IM2;
j=iy/NDIV;
iy=iv[j]-idum2;
iv[j]=*idum;
if(iy1)
iy+=IMM1;
if((temp=AM*iy)RNMX)
return RNMX;
else
return temp;
};double fun1(double x) /*被积函数*/
{
return(sin(x));
};ko1=ko2=1;
for(i=0;in;i++)
{
t=ran2(ko1);
s=ran2(ko2);
x=a+(b-a)*t;
y=(-1)*fMAX+2*fMAX*s;
if(y=0)
{
if(y==0)
sy++;
if(y=(*p)(x))
sx++;
};else
{
if(y=(*p)(x))
vx++;
}
}
return ((sx-vx)/n);
}
main()
{
double n,fMAX,a,b,Mont,A,B,h,i=0,N;
long time0,time1;
printf(Enter a=:);scanf(%lf,a);;printf(\nEnter b=:);scanf(%lf,b);
printf(\nEnter N=:);scanf(%lf,N); /*N表示区间的个数*/
printf(\nEnter n=:);scanf(%lf,n);
printf(\nEnter fMAX(此函数的界)=:);scanf(%lf,fMAX);
h=(b-a)/N;
time0=clock();
; for(i=0;iN;i++)
{
A=a+i*h;
B=A+h;
Mont+=Monte(n,A,B,fun1,fMAX)*(B-A)*2*fMAX;
}
printf(\nintegration=%lf,Mont);
time1=clock();
printf(\n耗时：%ldms,time1-time0);
get
文档评论（0）