用乘同余法產生隨機數(共7頁)

1、精選優(yōu)質文檔-傾情為你奉上例2.1 用乘同余法產生隨機數 (見光盤 FLch2sjxleg1.m)編程如下：A=6; N=100; %初始化；x0=1; M=255;for k=1:N %乘同余法遞推100次； x2=A*x0; %x2和x0分別表示xi和xi-1； x1=mod (x2,M); %將x2存儲器的數除以M，取余數放x1（xi）中； v1=x1/256; %將x1存儲器的數除以256得到小于1的隨機數放v1中； v(:,k)=v1; % 將v1中的數（）存放在矩陣存儲器v的第k列中，v(:,k)%表示行不變、列隨遞推循環(huán)次數變化； x0=x1; %xi-1= xi； v0=v1;

2、end %遞推100次結束；v2=v %該語句末無；，實現矩陣存儲器v中隨機數放在v2中，%且可直接顯示在MATLAB的window中；k1=k;%grapher %以下是繪圖程序；k=1:k1;plot(k,v,k,v,r);xlabel(k), ylabel(v);tktle(0-1)均勻分布的隨機序列) 程序運行結果如圖2.5所示。圖2.5 采用MATLAB產生的(0,1)均勻分布的隨機序列圖 產生的(0-1)均勻分布的隨機序列在程序運行結束后，產生的(0，1)均勻分布的隨機序列，直接從MATLAB的window界面中copy出來如下（v2中每行存6個隨機數）： v2 = 0.0234

3、0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.25780.5508 0.3164

4、0.9023 0.4336 0.6094 0.66800.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.31640.9023 0.4336 0.6094 0.6680 0.0234 0.14060.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.08200.4922 0.9609 0.7852 0.7

5、266 0.3750 0.25780.5508 0.3164 0.9023 0.4336 0.6094 0.66800.0234 0.1406 0.8438 0.0820第二章的白噪聲產生程序例2.2 用乘同余法產生(見光盤 FLch2bzsheg2.m) 編程如下：A=6; x0=1; M=255; f=2; N=100； %初始化；x0=1; M=255;for k=1: N %乘同余法遞推100次； x2=A*x0; %分別用x2和x0表示xi+1和xi-1； x1=mod (x2,M); %取x2存儲器的數除以M的余數放x1（xi）中； v1=x1/256; %將x1存儲器中的數除以2

6、56得到小于1的隨機數放v1中； v(:,k)=(v1-0.5 )*f; %將v1中的數（）減去0.5再乘以存儲器f中的系數，存放在矩陣存儲器v的第k列中，v(:,k)表示行不變、列隨遞推循環(huán)次數變化； x0=x1; % xi-1= xi； v0=v1;end %遞推100次結束；v2=v %該語句后無；，實現矩陣存儲器v中隨機數放在v2中，且可直接顯示在MATLAB的window中；k1=k;%grapher %以下是繪圖程序；k=1:k1;plot(k,v,k,v,r);xlabel(k), ylabel(v);tktle( (-1,+1)均勻分布的白噪聲) 程序運行結果如圖2.6所示。圖

7、2.6 采用MATLAB產生的(-1,+1)均勻分布的白噪聲序列 產生的(-1，1)均勻分布的白噪聲序列在程序運行結束后，產生的(-1，1)均勻分布的白噪聲序列，直接從MATLAB的window界面中copy出來如下（v2中每行存6個隨機數）：v2 = -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703

8、0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188

9、0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359 -0.0156 0.9219 0.5703 0.4531 -0.2500 -0.4844 0.1016 -0.3672 0.8047 -0.1328 0.2188 0.3359 -0.9531 -0.7188 0.6875 -0.8359*另外，書中圖2.3白噪聲的產生如下：顯然，只

10、要在例2.2程序的初始化部分中給N=300，f=6，運行程序就可以得到如圖2.3所示的（-3，3）的白噪聲過程.編程如下： A=6; x0=1; M=255; f=6; N=300； %初始化；x0=1; M=255;for k=1: N %乘同余法遞推100次； x2=A*x0; %分別用x2和x0表示xi+1和xi-1； x1=mod (x2,M); %取x2存儲器的數除以M的余數放x1（xi）中； v1=x1/256; %將x1存儲器中的數除以256得到小于1的隨機數放v1中； v(:,k)=(v1-0.5 )*f; %將v1中的數（）減去0.5再乘以存儲器f中的系數，存放在矩陣存儲器v

11、的第k列中，v(:,k)表示行不變、列隨遞推循環(huán)次數變化； x0=x1; % xi-1= xi； v0=v1;end %遞推100次結束；v2=v %該語句后無；，實現矩陣存儲器v中隨機數放在v2中，且可直接顯示在MATLAB的window中；k1=k;%grapher %以下是繪圖程序；k=1:k1;plot(k,v,k,v,r);xlabel(k), ylabel(v);tktle( (-1,+1)均勻分布的白噪聲) 程序運行結果如圖2.3所示。圖2.3 白噪聲過程第二章的M序列產生程序例2.3 用移位寄存器產生M序列的MATLAB軟件實現(見光盤 FLch2bzsheg3.m) 編程如下

12、： X1=1;X2=0;X3=1;X4=0; %移位寄存器輸入Xi初T態(tài)（0101）， Yi為移位寄存器各級輸出m=60; %置M序列總長度for i=1:m %1# Y4=X4; Y3=X3; Y2=X2; Y1=X1; X4=Y3; X3=Y2; X2=Y1; X1=xor(Y3,Y4); %異或運算 if Y4=0 U(i)=-1; else U(i)=Y4;end endM=U%繪圖i1=ik=1:1:i1;plot(k,U,k,U,rx)xlabel(k)ylabel(M序列)title(移位寄存器產生的M序列) 程序運行結果如圖2.8 所示。圖2.8 軟件實現的移位寄存器產生的M序列圖. 四級移位寄存器產生的M序列M = Columns 1 through 10 -1 1 -1 1 1 1 1 -1 -1 -1 Columns 11 through 20 1 -1 -1 1 1 -1 1 -1 1 1 Column