<rt id="bn8ez"></rt>
<label id="bn8ez"></label>
  • 

    • <span id="bn8ez"></span>
    

    <label id="bn8ez"><meter id="bn8ez"></meter></label>
    

    
    

    




    







												


        


    

    
	
	



	

    

    


    
	
		
    
			


    
posts - 134,comments - 22,trackbacks - 0


    
	

    


    
	
					


    
	
	
    使用了 G. E. P. Box、M. E. Muller 和 G. Marsaglia 的極坐標法 (polar method)生成符合高斯分布的隨機數
    

     1#include <time.h>
    
     2#include <stdio.h>
    
     3#include <stdlib.h>
    
     4#include <math.h>
    
     5
    
     6/*生成符合0-1均勻分布的隨機數*/
    
     7double randf()
    
     8{
    
     9    return (double)rand()/RAND_MAX;
    
    10}
    
    11
    
    12
    
    13/* 高斯分布隨機數生成器 */
    
    14/* 均值 m, 標準差 s */
    
    15double randomGaussian(double m, double s)    
    
    16{                        
    
    17    double x1, x2, w, y1;
    
    18    static double y2;
    
    19    static bool haveNext= false;
    
    20
    
    21    if (haveNext)                
    
    22    {
    
    23        y1 = y2;
    
    24        haveNext = false;
    
    25    }
    
    26    else
    
    27    {
    
    28        do 
    
    29        {
    
    30            x1 = 2.0 * randf() - 1.0;
    
    31            x2 = 2.0 * randf() - 1.0;
    
    32            w = x1 * x1 + x2 * x2;
    
    33        }
    
    34        while ( w >= 1.0 || w==0);
    
    35        
    
    36        w = sqrt( (-2.0 * log( w ) ) / w );
    
    37        y1 = x1 * w;
    
    38        y2 = x2 * w;
    
    39        haveNext = true;
    
    40    }
    
    41    
    
    42    return( m + y1 * s );
    
    43}
    
    44
    
    45void main()
    
    46{
    
    47    srand((unsigned)time( NULL )); //初始化隨機種子
    
    48
    
    49    //生成10個服從均值為0  標準差為1的高斯分布的隨機數
    
    50    double tmp;
    
    51
    
    52    for(int i=0;i<10;i++)
    
    53    {
    
    54        tmp=randomGaussian(0,1);
    
    55        printf("%f\n",tmp);
    
    56    }
    
    57}
    
結果圖：
    

    
參考：http://www.taygeta.com/random/gaussian.html 

	
    posted on 2009-02-26 15:22 何克勤 閱讀(1352) 評論(0)  編輯  收藏  
    

    







    
	    
    




    





    


    

	
    
		
    
				
		
	
		
    
		
		
    只有注冊用戶登錄后才能發表評論。
    
		
    
			
		
    
			
		
    
		
		
    
		


    

		    		
		
    
		
    			
		
		
		
		
    
		
    
		
		
    
		
    
			網站導航:
		
		
    
		
    
			
		
    	
		
    
			 
		
    
		
    
			
		
    
		
    
			
				
		
    		
	
    	

    



    

				

    



    







    
        
	



主站蜘蛛池模板：
国内精品免费在线观看|
久久久久亚洲AV成人网人人软件|
亚洲国产精品日韩在线观看|
宅男666在线永久免费观看|
黄网站免费在线观看|
日韩毛片在线免费观看|
亚洲一区二区三区丝袜|
久久久久亚洲av无码专区喷水|
免费a级毛片在线观看|
最新69国产成人精品免费视频动漫
|
亚洲日本va一区二区三区|
亚洲成a人片77777群色|
色噜噜综合亚洲av中文无码|
日本二区免费一片黄2019|
亚洲精品无码久久千人斩|
免费视频精品一区二区|
亚洲av无码国产精品色在线看不卡
|
亚洲国产精品专区在线观看|
91香蕉视频免费|
妞干网免费视频观看|
四色在线精品免费观看|
国产jizzjizz免费视频|
亚洲色偷拍另类无码专区|
亚洲国产三级在线观看|
亚洲欧洲高清有无|
亚洲影院天堂中文av色|
成人在线免费视频|
18禁在线无遮挡免费观看网站|
亚洲免费福利视频|
亚洲国模精品一区|
亚洲AV午夜福利精品一区二区|
永久免费的网站在线观看|
国产精品玖玖美女张开腿让男人桶爽免费看
|
亚洲不卡AV影片在线播放|
伊人婷婷综合缴情亚洲五月|
亚洲国产激情在线一区|
免费一级毛片在线播放视频|
日韩激情淫片免费看|
亚洲精品白色在线发布|
国产免费内射又粗又爽密桃视频|
好吊妞在线成人免费|