mengtekaluo.txt.txt

蒙特卡洛方法的一些应用

#include <stdio.h> 
#include <stdlib.h> 
#include <time.h> 
void main() 
{ 
int i,k,n=0; 
int N=10000;
float a,b;
float an[10000],bn[10000];
srand((unsigned)time(NULL)); 
for(k=0;k<N;k++)
{i = rand() % 1000; 
an[k]=i/1000.0;} 
for(k=0;k<N;k++)
{i = rand() % 1000; 
bn[k]=(i/1000.0);} 

for(k=0;k<=N;k++)
{ a=an[k];
  b=bn[k];
	if(a*a+b*b<=1) n++;}
 printf("pi=%f\n",4.0*n/N);
}

#include <math.h>
#include <stdio.h>
#define f(x) x*sin(x)/(1+cos(x)*cos(x))
static double xn=1;
double rand()
{double lambda=1687,p=2147483647;
 xn=fmod(lambda*xn,p);
 return xn/p;}
void srand(double seed)
{xn=seed;}
main()
{int i,n=10000;
 double x,z,v,s=0.0,pi=3.141592;
 srand(4);
 for(i=0;i<=n;i++)
 {x=rand();
  x=x*pi;
  s=s+f(x);}
  v=pi*s/n;
printf("V=%10.6f\n",v);
}


#include<stdio.h>
#include<math.h>
#define f(x) x*exp(x)/(1+x)/(1+x)
static double xn=1;
double rand()
{double lambda=16807,p=2147483647;
 xn=fmod(lambda*xn,p);
 return xn/p;
}
void srand(double seed)
{xn=seed;}
main()
{int i,n=10000;
 double x,v,s=0;
 srand(4);
 for(i=0;i<n;i++)
{x=rand();s=s+f(x);}
v=s/n;
printf("V=%10.6f\n",v);
}

#include<math.h>
#include<stdio.h>
#define f(x,y,z) x*y*z
static double xn=1;
double rand()
{double lambda=16807,p=2147483647;
 xn=fmod(lambda*xn,p);
 return xn/p;
}
void srand(double seed)
{ xn=seed;}
main()
{long i,n=1000,m=0;
 double x,y,z,ita,v;
 srand(256);
 for(i=0;i<n;i++)
{ x=rand();
  y=rand();
  z=rand();
  ita=rand();
  if(ita<=f(x,y,z)) m=m+1;
}
v=(double)m/n;
printf("N=%ld           v=%10.6f\n",n,v);


#include<math.h>
#include<stdio.h>
#define f(x,y,z) x*x*y+z*log(y)+exp(x)
#define lambda 16807
#define p 2147483647
static double xn=1;
double rand()
{xn=fmod(lambda*xn,p);
 return xn/p;}
void srand(double seed)
{ xn=seed;}
main()
{int i,n=1000;
 double x,y,z,v,s=0;
 srand(256);
 for(i=0;i<n;i++)
{ x=2*rand()-1;
  y=3*rand()+3;
  z=2*rand();
  s=s+f(x,y,z);}
  v=2*3*2*s/n;
printf("N=%8d           v=%10.6f\n",n,v);
}

#include<math.h>
#include<stdio.h>
void rand(int m,double x[])
{ int i;
  double q[9];
  static double r=0;
  q[1]=0.4142135623;q[2]=0.7320508057;
  q[3]=0.2360679774;q[4]=0.6457513110;
  q[5]=0.3166247904;q[6]=0.6055512755;
  q[7]=0.1231056256;q[8]=0.3588989435;
  for(r++,i=1;i<=m;i++)
  x[i]=fmod(r*q[i],1.0);
}
double f(int m,double x[])
{int i;
 double z=0.0;
 for(i=1;i<=m;i++)
 z=z+x[i]*x[i];
 return z;
}
main()
{int i,k,m=8,n=10000;
 double x[9],v,s;
 for(k=1;k<=m;k++)
 {for(s=0,i=0;i<n;i++)
 {rand(k,x);
  s=s+f(k,x);
  }
v=s/n;
printf("m=%4d       n=%8d       s=%12.6f\n",k,i,v);
}
}


#include"stdio.h"
#include"stdlib.h"
#include"math.h"
#include"time.h"
#include"conio.h"
 float random(float lanmda,int i)
{float y,t;
 srand(i);
y=(float)(rand()%100)/100;
t=(-lanmda)*log(1-y);
return t;
}
void main()
{float t1,t2,time,jgt[300],lanmda,sum=0,a[300][4],at[300],st[300],lt[300];
int h,min,i,N,M,head=0,tail=0,ql;

printf("请输入时间间隔平均值t1=");
scanf("%f",&t1);
lanmda=t1;
for(i=1;i<=3000;i++)
{jgt[i]=random(lanmda,i);
sum=sum+jgt[i];
if(sum>=480)
{N=i;
break;}
}

printf("请输入服务时间平均值t2=");
scanf("%f",&t2);
lanmda=t2;
sum=0;
for(i=1;i<=3000;i++)
{st[i]=random(lanmda,i);
sum=sum+st[i];
if(sum>=480)
{M=i;
break;}
}
if(M<N) N=M;

at[1]=jgt[1];
for(i=2;i<=N;i++)
{at[i]=at[i-1]+jgt[i];}

lt[1]=at[1]+st[1];
for(i=2;i<=N;i++)
{lt[i]=lt[i-1]+st[i];
   if(lt[i-1]<at[i])
    lt[i]=at[i]+st[i];}

for(i=1;i<=N;i++)
{a[i][1]=at[i];
a[i][2]=st[i];
a[i][3]=lt[i];}
printf("求任意时刻队列长度");
printf("请输入8：00-16：00之间任意时间：\n");
printf("输入整点是o'clock=");
scanf("%d",&h);
printf("请输入分minutes=");
scanf("%d",&min);
time=(h-8)*60+min;
for(i=1;i<=N;i++)
{if(a[i][1]>=time) {tail=i;break;}
tail=N;}
for(i=1;i<=N;i++)
{if(a[i][3]>=time) {head=i;break;}
head=N;}
ql=tail-head;
printf("%d点%d分的队长是%d\n",h,min,ql);
printf("今天窗口共服务%d个顾客\n",N);
for(i=1;i<=N;i++)
{printf("第%d个客户到来时间是%d:%d\n",i,8+(int)(a[i][1]/60),(int)(a[i][1])%60);
printf("服务时间是%f minutes\n",a[i][2]);
printf("第%d个客户离开时间是%d:%d\n\n",i,8+(int)(a[i][3]/60),(int)(a[i][3])%60);
}
}