li_clifff

銀蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 212.8
散金: 1
紅花: 1
帖子: 354
在線: 171.4小時
蟲號: 846028
注冊: 2009-09-12
專業(yè): 流體力學

[求助] 請教向量Q滿足高斯分布，用matlab如何表示？

請教個問題：

      三維的向量Q是高斯隨機變量，均值為零，方差為1，用matlab如何賦值?

1)  Q = randn(1,3)?
2) Q1 = randn(1); Q2=randn(1);Q3=randn(1); Q =[Q1;Q2;Q3];  ?

如果有N個Q這樣的向量，該怎么一起賦值呢？

NQ = randn(N,3); ?

回復此樓

1樓 2012-03-25 03:59:32

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

綠遍山原

銅蟲 (小有名氣)

應助: 64 (初中生)
金幣: 707.7
紅花: 5
帖子: 159
在線: 201.2小時
蟲號: 1128371
注冊: 2010-10-21
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

★ ★
余澤成: 金幣+2, 謝謝參與應助！ 2012-03-30 14:09:40

內容已刪除

贊一下(1人)

回復此樓

要夢遊，不要催眠。

5樓2012-03-25 09:41:48

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

綠遍山原

銅蟲 (小有名氣)

應助: 64 (初中生)
金幣: 707.7
紅花: 5
帖子: 159
在線: 201.2小時
蟲號: 1128371
注冊: 2010-10-21
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
li_clifff: 金幣+10, ★★★很有幫助, 空間向量是不是會特殊點，一般的其他意義下的向量，它的分量或許是有關聯(lián)的， 2012-03-28 12:36:46
余澤成: 金幣+1, 謝謝參與應助！ 2012-03-30 14:10:47

http://www.mathworks.cn/help/toolbox/stats/mvnrnd.html
這個是產生多變量gaussian分布的函數(shù)

贊一下(1人)

回復此樓

要夢遊，不要催眠。

6樓2012-03-25 09:49:08

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

銳利的碎片

木蟲 (正式寫手)

star watcher

應助: 136 (高中生)
金幣: 3637.1
散金: 252
紅花: 22
帖子: 988
在線: 1224.9小時
蟲號: 961933
注冊: 2010-03-05
專業(yè): 凝聚態(tài)物性 II ：電子結構

【答案】應助回帖

★
感謝參與，應助指數(shù) +1
余澤成: 金幣+1, 謝謝參與應助！ 2012-03-30 14:09:55

這樣只是產生三個獨立的gauss隨機數(shù)，表示的矢量不滿足gauss分布。我覺得你要把問題說得更清楚點。

贊一下

回復此樓

2樓2012-03-25 08:45:54

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

sogooda

金蟲 (正式寫手)

應助: 19 (小學生)
金幣: 1497.4
紅花: 2
帖子: 737
在線: 94.6小時
蟲號: 438143
注冊: 2007-09-16
專業(yè): 防災工程

【答案】應助回帖

★ ★ ★
感謝參與，應助指數(shù) +1
余澤成: 金幣+3, 謝謝參與應助！ 2012-03-30 14:10:16

randn就是生成0均值、1方差的高斯分布隨機數(shù)。
你確信要用三維的高斯隨機分布矩陣（我感覺這個好像不常用啊，不知是不是我們理解有偏差）？
randn([m,n,p]) 將生成m*n*p的三維隨機數(shù)組。

我看你帖子里的意思似乎是想生成一個N*3的二位隨機數(shù)組。如果只是做這件事的話，
randn(m,n)
就可以。它將生成m*n的二維數(shù)組(矩陣)。也就是你帖子里用的
NQ = randn(N,3);
就可以。

贊一下

回復此樓

3樓2012-03-25 08:48:31

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

sogooda

金蟲 (正式寫手)

應助: 19 (小學生)
金幣: 1497.4
紅花: 2
帖子: 737
在線: 94.6小時
蟲號: 438143
注冊: 2007-09-16
專業(yè): 防災工程

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
li_clifff: 金幣+10, ★★★很有幫助 2012-03-28 12:35:20
余澤成: 金幣+2, 謝謝參與應助！歡迎常來程序語言版！ 2012-03-30 14:10:33

引用回帖:

2樓: Originally posted by 銳利的碎片 at 2012-03-25 08:45:54:
這樣只是產生三個獨立的gauss隨機數(shù)，表示的矢量不滿足gauss分布。我覺得你要把問題說得更清楚點。

不是的。
采用NQ = randn(N,3);這種方式生成的偽隨機數(shù)組，在統(tǒng)計意義上就是滿足Gauss分布的。當N=100時，我們可以做如下簡單的驗證：

>>NQ=randn(100,3);
>> mean(NQ)
ans =
-0.1108 -0.0647 0.0155
>> std(NQ)
ans =
1.0163 0.9420 0.8736
>> hist(NQ(:,1),20)
>> hist(NQ(:,2),20)
>> hist(NQ(:,3),20)

計算說明NQ的每一列都是0均值，1方差的。直方圖也可以看出，分布是形狀與Gauss分布接近。

贊一下

回復此樓

4樓2012-03-25 08:59:33

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

li_clifff

銀蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 212.8
散金: 1
紅花: 1
帖子: 354
在線: 171.4小時
蟲號: 846028
注冊: 2009-09-12
專業(yè): 流體力學

引用回帖:

6樓: Originally posted by 綠遍山原 at 2012-03-25 09:49:08:
http://www.mathworks.cn/help/toolbox/stats/mvnrnd.html
這個是產生多變量gaussian分布的函數(shù)

我在想，二維是不是特例，剛好可以用分量來表示，
比如 f(r)=f(x,y)= 1/(2pi) exp{- (x*x+y*y)/2 } = f(x) f(y),
其中x,y表示向量r的兩個坐標分量，這樣用兩獨立的一維標準高斯隨機變量表示是可以的，但是如果r是三維的呢？

贊一下

回復此樓

7樓2012-03-25 11:26:19

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

li_clifff

銀蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 212.8
散金: 1
紅花: 1
帖子: 354
在線: 171.4小時
蟲號: 846028
注冊: 2009-09-12
專業(yè): 流體力學

引用回帖:

4樓: Originally posted by sogooda at 2012-03-25 08:59:33:
不是的。
采用NQ = randn(N,3);這種方式生成的偽隨機數(shù)組，在統(tǒng)計意義上就是滿足Gauss分布的。當N=100時，我們可以做如下簡單的驗證：

>>NQ=randn(100,3);
>> mean(NQ)
ans =
-0.1108 - ...

謝謝您的驗證，對于滿足標準高斯分布的向量，我想應該是可以寫成它的分量形式，每一個都是獨立同分布的，好像二維三維都可以，你覺得呢

贊一下

回復此樓

8樓2012-03-25 11:34:55

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

綠遍山原

銅蟲 (小有名氣)

應助: 64 (初中生)
金幣: 707.7
紅花: 5
帖子: 159
在線: 201.2小時
蟲號: 1128371
注冊: 2010-10-21
性別: GG
專業(yè): 凝聚態(tài)物性I:結構、力學和

★ ★
余澤成: 金幣+2, 鼓勵討論！ 2012-03-30 14:11:00

引用回帖:

7樓: Originally posted by li_clifff at 2012-03-25 11:26:19:
我在想，二維是不是特例，剛好可以用分量來表示，
比如 f(r)=f(x,y)= 1/(2pi) exp{- (x*x+y*y)/2 } = f(x) f(y),
其中x,y表示向量r的兩個坐標分量，這樣用兩獨立的一維標準高斯隨機變量表示是可以的，但是如 ...

二維的分布也不是這個形式
自己看，http://202.117.96.226:8090/gltj/4/gltj04080101.htm
所以直接用獨立變量是不行的，三維的用mvnrnd可以求吧。

贊一下

回復此樓

要夢遊，不要催眠。

9樓2012-03-25 11:50:38

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

sogooda

金蟲 (正式寫手)

應助: 19 (小學生)
金幣: 1497.4
紅花: 2
帖子: 737
在線: 94.6小時
蟲號: 438143
注冊: 2007-09-16
專業(yè): 防災工程

引用回帖:

5樓: Originally posted by 綠遍山原 at 2012-03-25 09:41:48:
這個不一樣的，矢量滿足的是多變量高斯分布，而randn的算法只產生一維的高斯分布。你做的只是檢查了每個分量的均值。
參考這里http://en.wikipedia.org/wiki/Multivariate_normal_distribution#Drawing_values_ ...

恩，謝謝提醒。確實是這樣的，看來我對樓主的意思有誤解。
建議樓主還是用mnvrnd來生成吧，這樣比較嚴謹。

贊一下

回復此樓

10樓2012-03-25 16:51:09

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 【考研調劑】化學專業(yè) 281分，一志愿四川大學，誠心求調劑 +6	吃吃吃才有意義 2026-03-19	6/300	2026-03-20 10:47 by 盡舜堯1
[考研] 求調劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學校友
[考研] 材料學碩318求調劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 0703化學調劑 +10	妮妮ninicgb 2026-03-15	14/700	2026-03-19 22:59 by 學員8dgXkO
[考研] 320求調劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學員8dgXkO
[考研] 266求調劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調劑 +11	yangfz 2026-03-17	11/550	2026-03-19 15:06 by houyaoxu
[考研] 本科鄭州大學物理學院，一志愿華科070200學碩，346求調劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 330求調劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 311求調劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 328求調劑，英語六級551，有科研經歷 +3	生物工程調劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 化學工程321分求調劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 268求調劑 +6	簡單點0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 085601求調劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 278求調劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 11408 一志愿西電，277分求調劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 070300化學學碩求調劑 +6	太想進步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 0703 物理化學調劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 326求調劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 求老師收留調劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e