Solved

Random Distribution Number Generator Algorithm

Posted on 2001-06-07
9
557 Views
Last Modified: 2008-02-01
I need to build a random number generator. The number generated should follow some distribution like noemal distribution, Poisson distribution, exponential distribution and so on....
I'm using Visual C++ programming, but I will be appreciate if anybody can just let me know the algorithm.
I have gone through the exact generator I needed in Excel, but I can't get the source code.
Thanx for your help.
0
Comment
Question by:ronaldyoong
[X]
Welcome to Experts Exchange

Add your voice to the tech community where 5M+ people just like you are talking about what matters.

  • Help others & share knowledge
  • Earn cash & points
  • Learn & ask questions
9 Comments
 

Expert Comment

by:nir_arad
ID: 6162935
There are specific algorithm tips in
"Simulation Modeling and Analysis", 2nd Edition by Law and Kelton.
If you need something specific - I have the book at hand.

Example:
Bernoulli (p.472)
-----------------
1. Generate U~U(0,1)
2. If U <= p return X = 1. Otherwise return X = 0.


Uniform distribution can be obtained by the C library rand() function.
0
 
LVL 18

Accepted Solution

by:
deighton earned 50 total points
ID: 6163202
Here is exponential for starters, include math.h

double exponential_distribution(double dlamda)
{
     double dprop;

     do
     {          
          dprop = (double) rand() / (double) RAND_MAX;
     }
     while (dprop == 0);

     return -log(dprop) / dlamda;
}


void main()
{


     printf("\nexponential lamda = 3 %lf",exponential_distribution(3));
         
         
}
0
 
LVL 18

Expert Comment

by:deighton
ID: 6163206
include stdlib, stdio as well

and do a

// include time.h

srand( (unsigned)time( NULL ) );  

at program start to seed the random number if you want different numbers every time you run the program

0
Independent Software Vendors: We Want Your Opinion

We value your feedback.

Take our survey and automatically be enter to win anyone of the following:
Yeti Cooler, Amazon eGift Card, and Movie eGift Card!

 
LVL 18

Expert Comment

by:deighton
ID: 6163264
poisson distribution

long poisson_distribution(double lamda)
{
     double total=0;
     double dprop,dterm;
     long k=0;
     long c;
     
     do
     {
          dprop = (double) rand() / (double) RAND_MAX;
     } while (dprop == 1);

     do
     {
          dterm = (double) exp(-lamda);
          for(c = 1;c <= k ; c++)
          {
               dterm *= lamda / c;                    
          }

          total += dterm;
          k++;
     } while(total < dprop);
     k--;
     
     return k;
}
0
 
LVL 2

Expert Comment

by:yairy
ID: 6163720
A 'Real' random number generator is impossible to build with a deterministic machine.




0
 
LVL 18

Expert Comment

by:deighton
ID: 6167760
yairy,

I agree that the rand function produces pseudo random numbers, a mathematical sequence that appears to be random, but obviously isn't as such random.  However they are adequate for most tasks,  and it is possible to replicate in c++ anything that excel can do.
0
 
LVL 7

Expert Comment

by:KangaRoo
ID: 6170488
0
 
LVL 11

Expert Comment

by:griessh
ID: 6811672
I think you forgot this question. I will ask Community Support to close it unless you finalize it within 7 days.
Unless there is objection or further activity,  I will suggest to accept "deighton" comment(s) as an answer.

If you think your question was not answered at all, you can post a request in Community support (please include this link) to refund your points.
The link to the Community Support area is: http://www.experts-exchange.com/jsp/qList.jsp?ta=commspt

Please do not accept this comment as an answer!
======
Werner
0
 
LVL 5

Expert Comment

by:Netminder
ID: 6828409
Per recommendation comment force/accepted by

Netminder
Community Support Moderator
Experts Exchange
0

Featured Post

Independent Software Vendors: We Want Your Opinion

We value your feedback.

Take our survey and automatically be enter to win anyone of the following:
Yeti Cooler, Amazon eGift Card, and Movie eGift Card!

Question has a verified solution.

If you are experiencing a similar issue, please ask a related question

Unlike C#, C++ doesn't have native support for sealing classes (so they cannot be sub-classed). At the cost of a virtual base class pointer it is possible to implement a pseudo sealing mechanism The trick is to virtually inherit from a base class…
Templates For Beginners Or How To Encourage The Compiler To Work For You Introduction This tutorial is targeted at the reader who is, perhaps, familiar with the basics of C++ but would prefer a little slower introduction to the more ad…
The viewer will learn how to pass data into a function in C++. This is one step further in using functions. Instead of only printing text onto the console, the function will be able to perform calculations with argumentents given by the user.
The viewer will learn how to user default arguments when defining functions. This method of defining functions will be contrasted with the non-default-argument of defining functions.

739 members asked questions and received personalized solutions in the past 7 days.

Join the community of 500,000 technology professionals and ask your questions.

Join & Ask a Question