通信建模作业程序集锦.docxVIP

1.正态分布:#include math.h #include stdlib.h #include time.h #include fstream #include iostream using namespace std; #define pi 3.1415926 //Generate a randon number followed by uniform distribution between the inteval [min,max] templatetypename T T rand(T min, T max) { return min+(max-min)*rand()/(RAND_MAX+1.0); } //calculate the normal density distribution values double normal(double x, double miu,double sigma) { return 1.0/sqrt(2*pi)/sigma*exp(-1*(x-miu)*(x-miu)/(2*sigma*sigma)); } //Generate a distributed random output numberdouble randn(double miu,double sigma, double min ,double max) { double x,y,dScope; do{ x=rand(min,max); y=normal(x,miu,sigma); dScope=rand(0.0,normal(miu,miu,sigma)); }while(dScopey); return x; } //Generate 10000 distributed random output numbers and put into a text called123int main(int argc,char* argv[]) { double min,max,miu,sigma; srand((unsigned)time( NULL )); cinminmaxmiusigma; ofstream outfile(123.txt); for (int i=0;i10000;i++) { outfile randn(miu,sigma,min,max) endl; } return 0; } 2.瑞利分布:#include math.h #include stdio.h #include stdlib.h #include time.h #define pi 3.1415926// generate the random numbers between the interval [0,1]double randn(){double f;f = (double)(rand()%1001); //[0,1000]return f/1000.0; //[0,1]}//calculate the Rayleigh density function values by using the following generated random numbers double ruili(double sigma) { Return randn()/(sigma*sigma)*exp(-randn()*randn()/(2*sigma*sigma)); }//The main function can generate 10000 distributed random output numbers which are put into a text called”321”int main(void) { double sigma; srand((unsigned)time( NULL )); printf(input sigma:\n); scanf(%lf,sigma); for (int i=0;i10000;i++) { FILE *fp;fp=fopen(321.txt,a+);fprintf(fp,%lf,,ruili(sigma));printf(%lf\n,randn());fclose(fp); } return 0;]3.泊