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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】Matlab随机数生成问题【已完成】
111/2返回列表12下一页
查看: 624  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

woshilsh

荣誉版主 (职业作家)

优秀版主

[交流] 【求助】Matlab随机数生成问题【已完成】

我想生成一个三维的数据，

数据服从双峰高斯分布，要400个点，中心在（1，0，0）和（-1，0，0）

试了几次，搞不定，那位大侠知道，交流下！

[ Last edited by nono2009 on 2009-9-23 at 22:10 ]
回复此楼

» 猜你喜欢
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了！[/img][/url][/center]
1楼 2009-07-01 19:21:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mininurse

金虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
woshilsh(金币+3,VIP+0):感谢，但是你生成的只是随机数，不一定复合双峰高斯分布的，期待继续关注 7-2 12:36
不知道这样做对不对

%%%% generate the random x, y, z dataset for the first Gaussian
A=randn(400,1,3);
X1=A(:,:,1)+1; %% x axis dataset
Y1=A(:,:,2);  %% y axis dataset
Z1=A(:,:,3);  %% z axis dataset

%%%% generate the random x, y, z dataset for the second Gaussian
B=randn(400,1,3);
X2=B(:,:,1)-1; %% x axis dataset
Y2=B(:,:,2);  %% y axis dataset
Z2=B(:,:,3);  %% z axis dataset

%% the outcome

newX=(X1+X2)/2; %% x dataset for the new mutual gaussian distribution
newY=(Y1+Y2)/2; %% y dataset for the new mutual gaussian distribution
newZ=(Z1+Z2)/2; %% z dataset for the new mutual gaussian distribution
引用回帖:
Originally posted by woshilsh at 2009-7-1 19:21:
我想生成一个三维的数据，

数据服从双峰高斯分布，要400个点，中心在（1，0，0）和（-1，0，0）

试了几次，搞不定，那位大侠知道，交流下！

[ Last edited by mininurse on 2009-7-2 at 09:06 ]
赞一下(21人) 回复此楼
2楼2009-07-02 08:55:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主
引用回帖:
Originally posted by mininurse at 2009-7-2 08:55:
不知道这样做对不对

%%%% generate the random x, y, z dataset for the first Gaussian
A=randn(400,1,3);
X1=A(:,:,1)+1; %% x axis dataset
Y1=A(:,:,2);  %% y axis dataset
Z1=A(:,:,3);  %% z axis ...

偶需要的是服从双峰高斯分布的数据，你这个思路对了，可是不一定服从正太分布
赞一下 回复此楼
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了！[/img][/url][/center]
3楼2009-07-02 12:36:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
woshilsh(金币+5,VIP+0):期待张兄有空再看看，三维数据咋弄，如何得到 7-2 19:54
引用回帖:
Originally posted by woshilsh at 2009-7-1 19:21:
我想生成一个三维的数据，

数据服从双峰高斯分布，要400个点，中心在（1，0，0）和（-1，0，0）

试了几次，搞不定，那位大侠知道，交流下！

[ Last edited by woshilsh on 2009-7-2 at 16:34 ]

>> MU = [-1 0;1 0];
>> SIGMA = cat(3,[.5 0;0 .5],[.5 0;0 .5]);
>> p = ones(1,2)/2;
>> obj = gmdistribution(MU,SIGMA,p);
>> ezsurf(@(x,y)pdf(obj,[x y]),[-2:.1:2],[-2:.1:2])
>>
赞一下(21人) 回复此楼
4楼2009-07-02 18:26:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mininurse

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
传说中的高手终于出现了，不是一般的高啊
赞一下(1人) 回复此楼
5楼2009-07-02 19:18:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主
引用回帖:
Originally posted by hitzhang at 2009-7-2 18:26:


>> MU = [-1 0;1 0];
>> SIGMA = cat(3,[.5 0;0 .5],[.5 0;0 .5]);
>> p = ones(1,2)/2;
>> obj = gmdistribution(MU,SIGMA,p);
>> ezsurf(@(x,y)pdf(obj,[x y]),[-2:.1:2 ...

呵呵，强大，

我的语句就复杂多了，都是C语言惹的祸啊

clear all;
mu = [1 0]; Sigma = [.9 .4; .4 .3];mu2=[0,1];
[X1,X2,X3] = meshgrid(linspace(-1,2,20)', linspace(-1,2,20)');
%%
X = [X1( : ) X2( : )];
p1= mvnpdf(X, mu, Sigma);
p2=mvnpdf(X,mu2,Sigma);
p=0.2*p2+0.8*p1;
surf(X1,X2,reshape(p,20,20));
figure(2);
subplot(1,2,1);plot(X(:,1),p,'-');
subplot(1,2,2);plot(X(:,2),p,'-');
figure(3);
plot3(X(:,1),X(:,2),p,'-');


但现在新问题出来了，我想要的三维数据呢？？？？

按理说横纵坐标都出来，然后把密度值作为Z轴之后，可以写为[x y z]就是得到的三维数据，可是仿真结果出来不对啊！

麻烦高手们再弄弄！
赞一下 回复此楼
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了！[/img][/url][/center]
6楼2009-07-02 19:53:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主
★ ★
sunxiao(金币+2,VIP+0):友情鼓励一下，呵呵 7-4 01:24
引用回帖:
Originally posted by hitzhang at 2009-7-2 18:26:


>> MU = [-1 0;1 0];
>> SIGMA = cat(3,[.5 0;0 .5],[.5 0;0 .5]);
>> p = ones(1,2)/2;
>> obj = gmdistribution(MU,SIGMA,p);
>> ezsurf(@(x,y)pdf(obj,[x y]),[-2:.1:2 ...

按照张兄的程序是不是应该这就是三维数据？

y=pdf(obj,[x y]);
[[-2:.1:2]' [-2:.1:2]' y']?

要不还是二维数据啊

[ Last edited by woshilsh on 2009-7-2 at 20:11 ]
赞一下(10人) 回复此楼
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了！[/img][/url][/center]
7楼2009-07-02 20:02:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
woshilsh(金币+2,VIP+0):感谢张兄百忙中抽空来看看，呵呵 7-2 21:57
楼上正解，这个就是熟能生巧
赞一下(21人) 回复此楼
8楼2009-07-02 21:47:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mininurse

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
大侠，为什么你用的一些函数我这边都没有呢？你的是MATLAB什么版本？
赞一下(1人) 回复此楼
9楼2009-07-03 09:16:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
★
kuhailangyu(金币+1,VIP+0):欢迎参与交流 7-3 16:32
2008a
回复此楼
10楼2009-07-03 10:31:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 woshilsh 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭