How to use mvnrnd function ?

Question

Betha Shirisha 2015-4-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/196744-how-to-use-mvnrnd-function

评论： Betha Shirisha 2015-4-10

I need to generate a matrix whose elements follow normal distribution with zero mean and Q covariance matrix . I know Q (some positive definite matrix)

I'm using the following code

Q size is (4,4) MU=zeros(1,4); noise=mvnrnd(MU,Q,10)+1j*mvnrnd(MU,Q,10);

MY noise size should be (4,10) but i m getting size as (10,4)

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Alka Nair 2015-4-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/196744-how-to-use-mvnrnd-function#answer_174671

The documentation of MVNRND suggests that the syntax is

>> r = mvnrnd(MU,SIGMA,cases) and r = mvnrnd(MU,SIGMA,cases) returns a cases-by-d matrix R of random vectors chosen from the multivariate normal distribution with a common 1-by-d mean vector MU, and a common d-by-d covariance matrix SIGMA. As cases is 10 in your case and d is 4, you are observing the dimension 10*4 for 'noise'.

Please refer to the documentation at:

http://www.mathworks.com/help/stats/mvnrnd.html

3 个评论
显示 1更早的评论隐藏 1更早的评论

John D'Errico 2015-4-10

编辑：John D'Errico 2015-4-10

No, no no.

cases, the third argument is the NUMBER of samples to be generated. Those samples are independent of each other.

The covariance matrix is the covariance of those random variables, but not from sample to sample. Again, the samples are INDEPENDENT.

So if the random variable lives in a d-dimensional space, then you will need to pass in a vector of length d for mu, and a matrix of size dxd for sigma. Then expect that the result will be of size (cases,d).

Betha Shirisha 2015-4-10