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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 请教向量Q满足高斯分布，用matlab如何表示？
汕头大学海洋科学接受调剂
121/2返回列表12下一页
查看: 2573  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学308分求调剂（一志愿华东师大）接受跨专业 +7 相信必会光芒万� 2026-04-13 7/350 2026-04-13 19:49 by nxybio2007
[考研] 085600材料与化工349分求调剂 +11 李木子啊哈哈 2026-04-12 12/600 2026-04-13 19:28 by Kilig0317
[教师之家] 转长聘了 +3 简单化xn 2026-04-13 3/150 2026-04-13 14:18 by wwwkkk83
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen～
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 一志愿华东师范生物学326分，求调剂 +8 刘墨墨 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 一志愿华南师范大学0702物理学305调剂 +4 念常安 2026-04-07 6/300 2026-04-08 22:53 by bljnqdcc
[考研] 求调剂，现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
信息提示
关闭
请填处理意见
关闭