MATLAB的答案

多元正态随机数与正态分布的随机数

16个视图(30天)
所罗门罗斯
所罗门罗斯 2021年5月28日
回答: 史蒂文的主 2021年5月28日
我使用两个样本抽样随机数字 多元正态分布 和样本 正态分布与指定的均值和方差 。
我有一个σμ平均向量和协方差矩阵:
μ= (-0.25,0.03,0.01)'
σ=眼睛(长度(μ))
使用多元正态随机数,我得到以下结果:
rng默认的%的再现性
R = mvnrnd(μ、σ3)
R =
0.2877 0.8922 -0.4236
1.5839 0.3488 0.3526
-2.5088 -1.2777 3.5884
利用正态分布的随机数与特定的意思(b)和方差(a),我得到以下结果:
= 1;%性病
b =μ(1);%的意思
rng默认的
y =。* randn (3,1) + b
y =
0.2877 - 0.5677
1.5839 - 1.8639
-2.5088 - -2.2288
创建第二列,我使用b =μ(2)。
我想知道为什么是不同的在这种情况下的结果。第一列的结果是相同的,但从第二列变化。难道这是相同的吗?如果不是,什么是正确的方式画随机数的平均向量μ。

在回答这个问题。

接受的答案

史蒂文的主
史蒂文的主 2021年5月28日
跑:
所以你这样做?
μ= (-0.25,0.03,0.01)';
= 1;%性病
b =μ(1);%的意思
rng默认的
y (: 1) = a。* randn (3,1) + b;
rng默认的
y (:, 2) =。* randnμ(3,1)+ (2)
y = 3×2
0.2877 0.5677 1.5839 1.8639 -2.5088 -2.2288
第二列不是完全随机的。考虑到第一列和μ向量我可以告诉你 完全 第二列将是什么。没有任何额外的我甚至可以计算它 randn 调用。
y (:, 3) = y(: 1) +(μμ(2)(1)
y = 3×3
0.2877 0.5677 0.5677 1.5839 1.8639 1.8639 -2.5088 -2.2288 -2.2288
规范(y (:, 3) - y (:, 2))%有效平等
ans = 4.4409 e-16
如果你现在离开,第二个调用 rng违约 或书面所以第二列生成使用发电机的状态,第一个“离开”:
rng默认的
z (: 1) =。* randn(3,1) +μ(1);
z (:, 2) = a。* randnμ(3,1)+ (2)
z = 3×2
0.2877 0.8922 1.5839 0.3488 -2.5088 -1.2777
rng默认的
z =。* randn(3 2) +μ(1:2)。%使用隐式扩张
z = 3×2
0.2877 0.8922 1.5839 0.3488 -2.5088 -1.2777

答案(1)

保罗
保罗 2021年5月28日
跑:
我不太确定你要创建y,由于代码生成y作为一个单独的列中。
它看起来像y是由:
μ= (-0.25,0.03,0.01)';
σ=眼睛(长度(μ));
y = 0 (3,2);
rng默认的
y (: 1) = randn (3,1) + mu (1);
rng默认的
y (:, 2) = randnμ(3,1)+ (2)
y = 3×2
0.2877 0.5677 1.5839 1.8639 -2.5088 -2.2288
然而,这不是正确的,因为y(:, 2)需要从样本生成独立的旅游房车的用于生成y (: 1)。
mvnrnd的输出()在这个例子中可以使用9重建如下样品randn(三维随机向量的三个样本):
rng默认的%的再现性
R = mvnrnd(μ、σ3)
R = 3×3
0.2877 0.8922 -0.4236 1.5839 0.3488 0.3526 -2.5088 -1.2777 3.5884
rng默认的
R1 = randn(3,3) +μ。
R1 = 3×3
0.2877 0.8922 -0.4236 1.5839 0.3488 0.3526 -2.5088 -1.2777 3.5884
R - R1
ans = 3×3
0 0 0 0 0 0 0 0 0

类别

找到更多的在随机数生成在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!