Matlab隨機數產生的問題

1、Matlab隨機數產生的問題2009-02-13 08:132006/04/212006/04/21/672276.aspx/url/color/size/font/b最后歡迎各位版友發(fā)表自己的看法和問題，不討論也就不容易發(fā)現問題和錯誤i 本帖最后由 eight 于 2007-9-5 15:02 編輯 /i 不化頑石 發(fā)表于 2007-9-5 16:20關于這個語句： randn('state',2);randn(1,6)在我的電腦里面，也是一樣的結果。只要是輸入“randn('state',j); randn(m,n)"出來的結果一定是相同的。j,m

2、,n變成什么都一樣。但如果隨后再生成一個randn(m,n)，而不重新設定seed，就是與第一個序列不同的了，不過同一個seed每次生成的第二個序列也都是一樣的。這大約說明randn的偽隨機過程還是按照某種固定模式迭代的吧。另外，關于花兄說的”去均值后可以嚴格保證均值為0，但是個人覺得意義不大?！拔蚁胙a充一點：我現在正在做的一個工作是與Monte Carlo Simulation相關的，它要求生成標準白噪聲樣本，是因為在隨后的過程中，要用到“樣本標準差為1”來進一步計算。所以，如果在這種情況下嚴格地對已有序列標準化還是必須的，而不能完全相信randn函數，雖然它在方差上造成的偏差要遠小于均值，

3、然而不知道前面這一點小偏差在后面會造成什么結果。j,m,n變成什么都一樣。但如果隨后再生成一個randn . /quote請你再試試，種子不同時。得到的結果肯定是不同的>>   randn('state',2);randn(1,6)ans =1.7491 0.1326 0.3252 -0.7938 0.3149 -0.5273>> randn('state',3);randn(1,6)ans =0.9280 0.1733 -0.6916 -0.7230 -0.5744 -0.3077隨機數的產生就是根據種子（用來決定初

4、始值），然后按照一定的規(guī)則外推。具體的方法有很多，可以看我附的鏈接。而randn用的是什么方法，我也弄不清楚，因為沒法看到其真正的源文件9191-1-8.html高斯信號和非高斯信號/url”帖子里提到的那篇文章另外，關于花兄你的那篇貼子，我有一點疑問，就是下面兩句：subplot(121),plot(y);           %此處可以看到高斯白噪聲的頻譜還是白噪聲subplot(122),plot(real(Y);     %高斯函數的傅立葉變換還是高斯函數這里是不是有點問題？按說頻譜不是

5、這么看的。那個是有些不太準確，我已改過。用abs好些，不過結果基本是一樣的。那個結論有嚴格的數學證明，可以推導出來的另外，我對頻譜其實了解也不多少，只是看到別人都是這么用的。你覺得應該怎么看頻譜呢？現在用得比較多的，其實是功率譜。你所說的abs(Y)，應該是能量譜。我是做大氣的，接觸小波分析要比純粹的FFT多一些，在小波中，單看實部圖或者虛部圖，從原理上講，都是不能反映出序列的振蕩性質的。當然這部分并不是我的方向，所以了解的也許不是很透徹，期待再交流。rand產生的是0到1(不包括1)的隨機數.matlab的rand函數生的是偽隨機數,即由種子遞推出來的,相同的種子,生成相同的隨機數.matl

6、ab剛運行起來時,種子都為初始值,因此每次第一次執(zhí)行rand得到的隨機數都是相同的.1.多次運行,生成相同的隨機數方法:用rand('state',S)設定種子 S為35階向量，最簡單的設為0就好例:rand('state',0);rand(10)2. 任何生成相同的隨機數方法:試著產生和時間相關的隨機數,種子與當前時間有關.rand('state',sum(100*clock)即:rand('state',sum(100*clock) ;rand(10)只要執(zhí)行rand('state',sum(100*clock) ;的當前計算機時間不現,生成的隨機值就不現.也就是如果時間相同,生成的隨機數還是會相同.在你計算機速度足夠快的情況下,試運行一下:rand('state',sum(100*clock);A=rand(5,5);rand('state