sdma.rar_SDMA_visual c

#include <stdio.h> #include <math.h> #include <time.h> #include <stdlib.h> #define N 5 //可用服务信道的数目，则复用全部后可用的为2N #define P 1000000 //假设1000000个用户来接受服务 #define ps 0.6 //信道复用成功的概率，用以产生0-1分布的随机数 double exp(double e) //产生指数分布随机数的函数 { double E=(rand()+1.0)/(RAND_MAX+1.0); return -e*log(E); } double fun() //产生0到1之间的随机数的函数,且小于ps的概率为0.6,大于ps表示创建复用信道失败，小于ps表示创建复用信道成功 { double x=(rand()+1.0)/(RAND_MAX+1.0); return x; } main() { FILE *fp1, //保存阻塞率和A值的文件指针 *fp2, *fp3; //保存用于仿真数据的文件指针 double A, //每个通话时间到达的用户数 H=100, //平均通话时间，假设为100s，用以产生指数分布的随机数 servetime, //用户的通话时间 LM, //用户到达的平均时间间隔，用以产生泊松分布的随机数 intervaltime, //到达时间间隔 primary[N]={0.0}, //为原始信道开辟数组 second[N]={0.0}, //为复用信道开辟数组 FDblock; //阻塞率 int sign[N], //标志复用的成功有否，0表示没有复用信道，1表示有复用信道 sum, //阻塞的用户数 i,j; srand((unsigned)time(NULL)); //初始化伪随机时间变量种子 if((fp1=fopen("FDblockA.txt","w"))==NULL) //打开保存阻塞率和A值的文件 { printf("file can not open!\n"); exit(0); } if((fp2=fopen("FDblock.txt","w"))==NULL) //打开用于仿真的文件 { printf("file can not open!\n"); exit(0); } if((fp3=fopen("A.txt","w"))==NULL) { printf("file can not open!\n"); exit(0); } for(LM=0.005;LM<=0.08;LM+=0.005) //由不同的的用户到达的平均时间间隔率产生不同的用户到达的状况,假设每次来P个用户 { A=LM*H; sum=0; for(i=0;i<N;i++) //初始化 { primary[i]=0.0; second[i]=0.0; sign[i]=0; } for(i=0;i<P;i++) { intervaltime=exp(1/LM); servetime=exp(H); for(j=0;j<N;j++) //初始化原始信道，释放已使用的复用信道 { primary[j]=primary[j]-intervaltime; if(primary[j]<=0) //原始信道空闲 { primary[j]=0.0; if(sign[j]==1) //已经使用的复用信道 { second[j]=second[j]-intervaltime; if(second[j]<=0) //已经使用完的复用信道释放 { second[j]=0.0; sign[j]=0; } else //未使用完的复用信道赋予原始信道并释放 { primary[j]=second[j]; sign[j]=0; } } } else //原始信道繁忙 { if(sign[j]==1) { second[j]=second[j]-intervaltime; if(second[j]<=0) //已经使用完的复用信道释放 { second[j]=0.0; sign[j]=0; } } } } for(j=0;j<N;j++) { if(primary[j]>0) //原始信道繁忙 { if(sign[j]==1) continue; //有复用信道 else //没有复用信道 { if(fun()<=ps) //创建复用信道成功，服务 { second[j]=servetime; sign[j]=1; break; } else continue; //创建复用信道失败 } } else //原始信道空闲,服务 { primary[j]=servetime; break; } } if(j==N) sum++; } FDblock=(double)sum/(P); printf("FDblock：%lf\n",FDblock); fprintf(fp1,"FDblock：%lf\n",FDblock); fprintf(fp2,"%lf\n",FDblock); printf("A：%lf\n",A); fprintf(fp1,"A：%lf\n",A); fprintf(fp3,"%lf\n",A); } fclose(fp3); fclose(fp2); fclose(fp1); }