2012-07-08 2 views
-1

평균 제로 및 분산 2.25의 가우시안 분포로 이산 난수 (1000 개의 앙상블에 대해)를 생성하고 싶습니다. NUMERICAL RECIPES in Fortran 90, 7 장, 7.2 절에서이 기능을 사용했는데이 기능을 사용하는 방법을 모르겠습니다.fortran 90의 이산 Guassian 모델

앙상블 수와 평균 및 분산을 어떻게 사용할 수 있습니까?

답변

1

IMO, NR에서 함수를 사용하는 방법을 가장 NR 자체로 작성 감사합니다.

정상적으로 배포 된 실수를 생성해야하는 경우 "Box-Muller transform"에 대해 google을 사용하거나 Wiki page 및/또는 수십 개의 질문과 답변을 여기에서 바로 찾으십시오.

귀하의 요구 사항이 정수를 생성하는 경우

, 나는 가장 쉬운 방법 is to compute the cdf 말 것 (당신이 약 이산 난수를? 물어 것을 의미했다).

+0

내가 이산 난수를 생성하려면 (이 숙제 인 경우 BTW, 제대로 같은 태그되어야한다) 제로 평균과 분산 2.25로 가우시안 분포 (1000 앙상블). 어떻게해야합니까? – user1272138

0

gasdev 수익률은 일반적으로 자주 N(0, 1)로 작성 0.0의 평균과 1.0의 분산이 임의의 변수를 배포했습니다. 원하는 분포와 일치 시키려면 분산을 사용하여 결과의 ​​크기를 조정해야합니다. xN(0, 1) (gasdev에 의해 반환 됨) 인 경우 y = a + b*xN(a, b^2), 즉 평균이 a이고 분산이 b^2 인 정규 분포의 확률 변수입니다.

+0

그건 그렇고, 이것은 숙제가 아닙니다. 나는 생물학 박사 과정 학생이며, 세포 집단을 모델링 할 때 이것을 필요로합니다. 1000 앙상블에 대해서, 나는이 숫자를 gasdev 함수에 넣을 곳이 궁금합니다. – user1272138

+0

gasdev 함수에서 idum이란 무엇입니까? – user1272138

+0

'gasdev'가 사용하는 의사 난수 생성 프로그램의 시드입니다. 'gasdev'에 대한 첫 번째 호출에서 음의 값을 전달한 다음'0'을 전달해야합니다. –