定积分Monte Carlo计算方法实现课件.ppt

《计算机高级语言》认知实习报告 ; ;一、设计任务与要求;;二、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