matlab编程问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问matlab软件中，怎样产生均值为5，协方差阵是1  0.5  0.5  0.5  0.5      的一组数
                                                                                    0.5  1  0.5  0.5  0.5
                                                                                    0.5 0.5   1  0.5  0.5
                                                                                    0.5 0.5  0.5   1  0.5
                                                                                    0.5 0.5  0.5  0.5   1

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：MATLAB编程 MATLAB atlab matla Mat MATLAB

这是蒙特卡洛基础知识中的“生成随机数”问题。
但是，你又没有指定这些数的分布，所以还没发生成。

5+A^0.5*randn(n,5)，其中，A为矩阵，n为数列长度，正态分布。

真牛这里的人，经济行，编程又行

很简单；
A=one(5,5)\2+eye(5,5)\2;

看看除法错误没，基本就这样，就可以产生对角元素为 1，其他元素为 0.5  的一个 5*5 的矩阵。

A=ones(5,5)/2+eye(5,5)/2;
这个可以生成你要的那个协方差矩阵。

你说的要均值为5，然后其他的，这个我做不了。我在MATLAB里试过了，不知道其他人怎么做。

惨不忍睹的实验结果。呵呵

>> A=ones(5,5)/2+eye(5,5)/2;
>> s=5+randn(5+A^0.5*randn(100,5)
??? s=5+randn(5+A^0.5*randn(100,5)
                                  |
Error: Expression or statement is incorrect--possibly unbalanced (, {, or [.

>> s=5+A^0.5*randn(100,5)
??? Error using ==> mtimes
Inner matrix dimensions must agree.

>> s=5*ones(5,5)+A^0.5*randn(5,5)

s =

    5.3755    4.1730    4.6841    6.1727    4.7111
    4.4260    3.9711    5.1235    6.0819    4.0557
    4.6161    5.1142    3.1196    5.2683    4.3939
    4.4241    4.4237    5.2623    5.5658    4.5194
    5.3062    4.9788    6.2436    5.6105    4.3031

>> m=mean(s)

m =

    4.8296    4.5322    4.8866    5.7398    4.3966

大家看看这个做法如何：

1， 定义A，然后生产随机数向量。
A=ones(5,5)/2+eye(5)/2;
result=5+randn(1,5)*chol(A);

2, 生成很多（1e6）个随机数向量，然后检验是否符合相关性质
clear data;for h=1e5:-1:1;data(h,1:5)=5+randn(1,5)*chol(A);end
>> mean(data)%均值性质
ans =
       5.0053       4.9988       5.0077       5.0065       5.0027
>> cov(data)  %协方差矩阵性质
ans =
       1.0042      0.50162      0.49954      0.50216      0.50143
      0.50162      0.99988       0.5016       0.4986      0.50117
      0.49954       0.5016      0.99769      0.49795      0.50243
      0.50216       0.4986      0.49795      0.99382      0.50154
      0.50143      0.50117      0.50243      0.50154       1.0039

R = MVNRND(MU,SIGMA,N) returns a N-by-D matrix R of random vectors
    chosen from the multivariate normal distribution with 1-by-D mean
    vector MU, and D-by-D covariance matrix SIGMA.

汗  MATLAB最强大之处是最带函数    。。。。。