Matlab隨機(jī)數(shù)產(chǎn)生的問題

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

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

3、然而不知道前面這一點小偏差在后面會造成什么結(jié)果。j,m,n變成什么都一樣。但如果隨后再生成一個randn . /quote請你再試試，種子不同時。得到的結(jié)果肯定是不同的>>   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隨機(jī)數(shù)的產(chǎn)生就是根據(jù)種子（用來決定初

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

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

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