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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】Matlab随机数生成问题【已完成】
111/2返回列表12下一页
查看: 850  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 中青基了要发朋友圈吗？ +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭