应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 请教向量Q满足高斯分布，用matlab如何表示？
121/2返回列表12下一页
查看: 2566  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

li_clifff

银虫 (正式写手)

[求助] 请教向量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发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

绿遍山原

铜虫 (小有名气)
★ ★
余泽成: 金币+2, 谢谢参与应助！ 2012-03-30 14:09:40
内容已删除
赞一下(1人) 回复此楼
要夢遊，不要催眠。
5楼2012-03-25 09:41:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绿遍山原

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
li_clifff: 金币+10, ★★★很有帮助, 空间向量是不是会特殊点，一般的其他意义下的向量，它的分量或许是有关联的， 2012-03-28 12:36:46
余泽成: 金币+1, 谢谢参与应助！ 2012-03-30 14:10:47
http://www.mathworks.cn/help/toolbox/stats/mvnrnd.html
这个是产生多变量gaussian分布的函数
赞一下(1人) 回复此楼
要夢遊，不要催眠。
6楼2012-03-25 09:49:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

锐利的碎片

木虫 (正式写手)

star watcher

【答案】应助回帖

★
感谢参与，应助指数 +1
余泽成: 金币+1, 谢谢参与应助！ 2012-03-30 14:09:55
这样只是产生三个独立的gauss随机数，表示的矢量不满足gauss分布。我觉得你要把问题说得更清楚点。
赞一下 回复此楼
2楼2012-03-25 08:45:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sogooda

金虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
余泽成: 金币+3, 谢谢参与应助！ 2012-03-30 14:10:16
randn就是生成0均值、1方差的高斯分布随机数。
你确信要用三维的高斯随机分布矩阵（我感觉这个好像不常用啊，不知是不是我们理解有偏差）？
randn([m,n,p]) 将生成m*n*p的三维随机数组。

我看你帖子里的意思似乎是想生成一个N*3的二位随机数组。如果只是做这件事的话，
randn(m,n)
就可以。它将生成m*n的二维数组(矩阵)。也就是你帖子里用的
NQ = randn(N,3);
就可以。
赞一下 回复此楼
3楼2012-03-25 08:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sogooda

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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随机数，表示的矢量不满足gauss分布。我觉得你要把问题说得更清楚点。

不是的。
采用NQ = randn(N,3);这种方式生成的伪随机数组，在统计意义上就是满足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发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
引用回帖:
6楼: Originally posted by 绿遍山原 at 2012-03-25 09:49:08:
http://www.mathworks.cn/help/toolbox/stats/mvnrnd.html
这个是产生多变量gaussian分布的函数

我在想，二维是不是特例，刚好可以用分量来表示，
比如  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发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
引用回帖:
4楼: Originally posted by sogooda at 2012-03-25 08:59:33:
不是的。
采用NQ = randn(N,3);这种方式生成的伪随机数组，在统计意义上就是满足Gauss分布的。当N=100时，我们可以做如下简单的验证：

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

谢谢您的验证，对于满足标准高斯分布的向量，我想应该是可以写成它的分量形式，每一个都是独立同分布的，好像二维三维都可以，你觉得呢
赞一下 回复此楼
8楼2012-03-25 11:34:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绿遍山原

铜虫 (小有名气)
★ ★
余泽成: 金币+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发消息 送TA红花 TA的回帖

sogooda

金虫 (正式写手)
引用回帖:
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发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 li_clifff 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026 WR青拔 +3 liuchb715 2026-04-09 4/200 2026-04-13 08:13 by kudofaye
[考研] 297工科，求调剂? +10 河南农业大学-能 2026-04-12 10/500 2026-04-12 22:44 by lbsjt
[考研] 296求调剂 +14 汪！？！ 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 0854调剂 +8 长弓傲 2026-04-11 9/450 2026-04-11 21:44 by 蓝云思雨
[考研] 材料与化工调剂 10+11 下一站上岸@ 2026-04-10 36/1800 2026-04-11 10:26 by 89436494
[考研] 282，求调剂 +12 jggshjkkm 2026-04-09 14/700 2026-04-11 09:39 by 猪会飞
[考研] 346，工科0854求调剂，专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 273求调剂 +51 麦小叮当 2026-04-06 58/2900 2026-04-10 15:54 by jiajinhpu
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 招收有机化学、化工，药学，食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 085404，334分，求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 270求调剂 +3 031127 2026-04-06 4/200 2026-04-08 21:00 by 逆水乘风
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
信息提示
关闭
请填处理意见
关闭