Chinaunix首页 | 论坛 | 博客
  • 博客访问: 838900
  • 博文数量: 190
  • 博客积分: 2991
  • 博客等级: 少校
  • 技术积分: 2400
  • 用 户 组: 普通用户
  • 注册时间: 2012-09-24 18:11
文章分类

全部博文(190)

文章存档

2015年(3)

2014年(1)

2013年(65)

2012年(121)

我的朋友

分类: C/C++

2012-11-08 16:12:24

在用计算机的一些智能算法(GA,PSO,ANN etc.)仿真时经常需要随机生成初始种群(初始样本),看看中的这两个函数的伪随机数生成吧~~~

1. 生成[a,b]之间的一个实数和一个整数
  1. /* 
  2. 定义函数     int   rand(void); 
  3. 函数说明     rand()会返回一随机数值,范围在0至RAND_MAX   间。 
  4. 在调用此函数产生随机数前,必须先利用srand()设好随机数种子,如果未设随机数种子,rand()在调用时会自动设随机数种子为1。 
  5. 关于随机数种子请参考srand()。  
  6. 返回值     返回0至RAND_MAX之间的随机数值,RAND_MAX定义在stdlib.h,其值为2147483647。  
  7. 范例 : 
  8. */  
  9. #include    
  10. #include   
  11.   
  12. double  doubleRand(double a,double b); 
  13. int         intRand(int a,int b);  
  14.   
  15. int main(void)  
  16. {  
  17.     double i=doubleRand(2.0,9.0);  
  18.     int     j=intRand(2,9);  
  19.     printf("%f \n",i);  
  20.     printf("%d \n",j );  
  21.   
  22.     return 0;  
  23. }   
  24. double doubleRand(double a,double b)  
  25. {  
  26.     double r;  
  27.     r=(double)rand()/RAND_MAX;  
  28.     return a+r*(b-a);  
  29. }  
  30. int intRand(int a,int b)  
  31. {  
  32.     return (int)doubleRand(a,b);  
  33. }  
以上代码中的每个函数只能生成一个随机数,至于为什么呢?
2.生成多个随机数

       之所以rand()每次的随机数都一样是因为rand()函数使用不正确。各种编程语言返回的随机数(确切地说是伪随机数)实际上都是根据递推公式计算的一组数值,当序列足够长,这组数值近似满足均匀分布。如果计算伪随机序列的初始数值(称为种子)相同,则计算出来的伪随机序列就是完全相同的。这个特性被有的软件利用于加密和解密。加密时,可以用某个种子数生成一个伪随机序列并对数据进行处理;解密时,再利用种子数生成一个伪随机序列并对加密数据进行还原。这样,对于不知道种子数的人要想解密就需要多费些事了。当然,这种完全相同的序列对于你来说是非常糟糕的。要解决这个问题,需要在每次产生随机序列前,先指定不同的种子,这样计算出来的随机序列就不会完全相同了。你可以在调用rand()函数之前调用srand(   (unsigned)time(NULL)),这样以time函数值(即当前时间)作为种子数,因为两次调用rand函数的时间通常是不同的,这样就可以保证随机性了。你也可以使用srand函数来人为指定种子数。

好,那按照这样,我就这样写~~~

  1. #include        
  2. #include        
  3. #include        
  4.   
  5. int main()     
  6. {     
  7.     for(int   i=0;i <100000;i++)     
  8.     {     
  9.         srand(   (unsigned)time(   NULL   )   );     
  10.         printf("%d\n",rand() );  
  11.     }    
  12.     return 0;  
  13. }    
答:你的程序是有问题的,你每产生一个随机数之前,都调用一次srand,而由于计算机运行很快,所以你每次用time得到的时间都是一样的(time的时间精度较低,只有55ms)。这样相当于使用同一个种子产生随机序列,所以产生的随机数总是相同的。
 
你应该把srand放在循环外:
  1. /* 
  2. #include  
  3. void srand(unsigned seed);  
  4. #include  
  5. time_t  time(time_t *time);
  6. */  
  7. #include        
  8. #include        
  9. #include        
  10.   
  11. int main()     
  12. {     
  13.     int i;  
  14.      //以当前系统时间作为种子  
  15.     srand(   (unsigned)time(   NULL   )   );     
  16.     for(i=0;i <10;i++)     
  17.     {     
  18.         printf("%d\n",rand() );  
  19.     }    
  20.     return 0;  
  21. }    
3. 若要不重复呢?即种群中的粒子都是不同的~~~
先来个最笨的办法:就是我拿一个数组来存你生成的随机数,一个一个放进来,边放边检查,这样的复杂度随着个数成阶层增长~~~且时间是不可预测的,这对RTOS是不好的消息~~~
但是简单好实现,走一个先~~~
  1. #include        
  2. #include        
  3. #include        
  4. #define MAX_NUM 10000  
  5. /*when insert a Rand_Num then check it*/  
  6. int check(int a[],int i)  
  7. {  
  8.     int j;  
  9.     for(j=0;j
  10.         if(*(a+j)==*(a+i))  
  11.             return 0;  
  12.     return 1;  
  13. }  
  14. int main()     
  15. {     
  16.     int i;  
  17.     int a[MAX_NUM];  
  18.      //以当前系统时间作为种子  
  19.     srand(   (unsigned)time(   NULL   )   );     
  20.     for(i=0;i <10;i++)     
  21.     {     
  22.         a[i]=rand();  
  23.         if (check(a,i)==0)  
  24.         {  
  25.             i--;    
  26.             continue;      //the number is the same with of one of the array number,so once again  
  27.         }  
  28.         printf("%d\n",a[i] );  
  29.     }    
  30.     return 0;  
  31. }    
这个和洗牌算法很类似,但没有那么分牌规则那么严格,应用的地方不同~~~
 
话说其实,种群产生的粒子很多情况下可以可以重复(与具体问题模型有关)~~~
那么不笨的方法呢?这个不像洗牌算法那么多规则,再想想~~~
好了,至少PSO的粒子可以产生了~~~:)
源文地址:
阅读(813) | 评论(0) | 转发(0) |
给主人留下些什么吧!~~