產(chǎn)生高斯隨機(jī)數(shù)

1、一、產(chǎn)生0-均1勻分布數(shù)米用混合同余法,x=(ax,c)modM之間的隨機(jī)數(shù)。通常情況卞，選取a=y=x/Myiii2045,c=1,M=220即為一個(gè)在,我們?nèi)我饨o定一個(gè)初始值x，即可得到一個(gè)隨機(jī)數(shù)，在程序中采用指針變量使得每次產(chǎn)生的隨機(jī)數(shù)都不一樣，只要改變循環(huán)的次數(shù)，就可以得到任意多個(gè)0隨1機(jī)數(shù)。程序?qū)崿F(xiàn)代碼如下doublex1(doublea,doubleb,int*s)doublet;*s=2045*(*s)+1;*s=*s-(*s/1048576)*1048576;t=(*s)/1048576.0;return(t);二、產(chǎn)生標(biāo)準(zhǔn)高斯分布數(shù)標(biāo)準(zhǔn)的高斯分布均值卩，方差1根據(jù)中心極限定理

2、有z=y0.6，根據(jù)此式可以由個(gè)分布數(shù)產(chǎn)生一個(gè)高斯數(shù)。改變ii循環(huán)次數(shù)可以得到任意多個(gè)標(biāo)準(zhǔn)高斯分布數(shù)。程序?qū)崿F(xiàn)代碼如下：intmain()doublemean,sigma;doublex100;inti,j;ints;mean=0.0;sigma=1.0;s=19444;cout產(chǎn)生個(gè)高斯數(shù)：endl;for(j=0;j200;j+)xj=0;for(i=0;i12;i+)xj=xj+x1(mean,sigma,&s);xj=xj-6.0;printf(%13.7f,xj);if(j+1)%5=0)coutendl;coutendl;產(chǎn)生一百個(gè)高斯數(shù)：嚴(yán)生i酣個(gè)島斯數(shù)t-0.66422461.

3、53079800.9023342-0.12383460.06557270.2713375-1.5182590-0.1274357-9.1929111-0.16390420.9786640.3181820-0.1896763-1.5249271-0.19928930.4630184-0.6747227-1.56173131.0492737-9.44292641.6019497-0.0243168-0.3334446-9.1496525-1.19965920.3373165-0.07044410.5928492-9.5595493-0.22683141.07927510.0345516-0.497

4、7297-1.4667696-9.63428691.4254S1S-1.6741734-0.13247112.O1563070.161092S1.0266372-1.6512547-0.94509160.25897760.25S7643-0.45466041.7935848-9.1332188-1.1842899-1.12134740.4S1390Q1.23729369.9468746-0.4013157-0.63158610.61934470.9922579-0.94456521.79465660.26320460.93166000.1969343-0.2598820-0.58121680.

5、28568Q8-0.42090800724798213368886一078627970.34599880.908697122650967一00340214074962430.5418148-0.5441685079245570.0399685187415120.15828510.9431515-0.53296850155706416224575-0.33193401.28Q8132-03635197一090165141.2171993-0.26868630.5664730一01640415一015944861.8685329-0.804315600852871-02118778一0457529

6、1一022588541.0963345-0.6900883140312770.5517635一18888994-1.3440800Q.90Q5Q32-0.221368814035854007614711.28459740.704717621997890一01794071-1.19458961.580Q228-1.242288601392574一02870579一034545330.3273525-0.71785931.2571926一08217106115871240.49924281.688161804012508一0498209018485636-0.74415020.57343101.6

7、0658840.15610310.2102566-0.05517011.7231045-0.90413861.80138210.76544760.10133931.1098385-2.7207737-0.71471601.9912930-0.55196570.8937893-0.2456608-0.3570347-1.63955121.89507101.4226131-0.69364360.09708210.5330715-1.4598942-1.7685337-0.09206580.05779080.35681721.6682949-0.95699500.71922870.1227474-1

8、.1331577-1.2477055-0.23261450.0878963-0.92289161.28580280.45226100.50226400.54909321.89352990.02385520.6158504-1.4672031-1.1745243-0.2678318-1.32134440.27821920.3316402-0.17280010.94067960.0353603-1.33797650.5589504-0.34807781.0542202-0.93337442.17741970.0623837-0.4152012-0.2045536-1.06739230.422063

9、8-1.12290380.09848591.07451441.9809628三、由得到的高斯分布數(shù)繪制概率密度函數(shù)采用API圖形處理中的Rectangle矩陣?yán)L制函數(shù)繪制圖形，將得到的高斯數(shù)分成一個(gè)個(gè)小區(qū)間，計(jì)算每個(gè)區(qū)間內(nèi)的高斯數(shù)的個(gè)數(shù)，用這些個(gè)數(shù)作為矩陣的高。程序代碼如下：#defineNN1000#definePI3.14159826#include#include/#include#include#include#include#include#include#includelongWINAPIWndProc(HWNDhWnd,UINTiMessage,UINTwParam,LONGlP

10、aram);BOOLInitWindowsClass(HINSTANCEhInstance);BOOLInitWindows(HINSTANCEhInstance,intnCmdShow);intWINAPIWinMain(HINSTANCEhInstance,HINSTANCEhPrevInstance,LPSTRlpCmdLine,intnCmdShow)MSGMessage;if(!InitWindowsClass(hInstance)returnFALSE;if(!InitWindows(hInstance,nCmdShow)returnFALSE;while(GetMessage(&

11、Message,0,0,0)TranslateMessage(&Message);DispatchMessage(&Message);returnMessage.wParam;longWINAPIWndProc(HWNDhWnd,UINTiMessage,UINTwParam,LONGlParam)HDChDC;HBRUSHhBrush;HPENhPen;PAINTSTRUCTPtStr;/doublex1(doublea,doubleb,int*s);doublemean=0,sigma=1.0;srand(unsigned)time(NULL);doublex10000,b10000=0,

12、temp;inti,j;ints1;intk;s1=19444;for(j=0;j10000;j+)xj=0;for(i=0;i12;i+)xj=xj+x1(mean,sigma,&s1);xj=xj-6.0;for(i=0;i10000;i+)temp=xi;/printf(%lft,temp);for(j=0;j10000;j+)if(-3+0.2*jtemp)bj=bj+1;/switch(iMessage)caseWM_PAINT:hDC=BeginPaint(hWnd,&PtStr);SetMapMode(hDC,MM_ANISOTROPIC);hPen=(HPEN)GetStock

13、Object(BLACK_PEN);hBrush=(HBRUSH)GetStockObject(DKGRAY_BRUSH);SelectObject(hDC,hBrush);SelectObject(hDC,hPen);/for(k=0;k100;k+)Rectangle(hDC,20+10*k,400-0.3*bk,30+10*k,400);EndPaint(hWnd,&PtStr);return0;caseWM_DESTROY:PostQuitMessage(0);return0;default:return(DefWindowProc(hWnd,iMessage,wParam,lPara

14、m);BOOLInitWindows(HINSTANCEhInstance,intnCmdShow)HWNDhWnd;hWnd=CreateWindow(WinFill,高斯分布,WS_OVERLAPPEDWINDOW,CW_USEDEFAULT,0,CW_USEDEFAULT,0,NULL,NULL,hInstance,NULL);if(!hWnd)returnFALSE;ShowWindow(hWnd,nCmdShow);UpdateWindow(hWnd);returnTRUE;BOOLInitWindowsClass(HINSTANCEhInstance)WNDCLASSWndClass;WndClass.cbClsExtra=0;WndClass.cbWndExtra=0;WndClass.hbrBackground=(HBRUSH)(GetStockObject(WHITE_BRUSH);WndClass.hCursor=LoadCursor(NULL,IDC_ARROW);WndClass.hIcon=LoadIcon(NULL,END);WndClass.hInstance=hInstance;WndClass.lpfnWndProc=WndProc;WndClass