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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】Matlab随机数生成问题【已完成】
111/2返回列表12下一页
查看: 740  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今年审到国自然15份，谈谈感受 +17 国自然国社科中 2026-05-17 18/900 2026-05-20 00:03 by kksky
[基金申请] 评审有感 +11 popular289 2026-05-18 19/950 2026-05-19 23:09 by popular289
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 11/550 2026-05-19 22:26 by 而立得乐
[教师之家] 上海大学实验技术岗位非升即走 +11 嘻嘻哈哈乐呵呵 2026-05-15 12/600 2026-05-19 22:17 by otani
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 7/350 2026-05-19 19:23 by bencaidai
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +3 光伏爱好者 2026-05-17 5/250 2026-05-19 19:13 by 光伏爱好者
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[有机交流] 反应很差，大量原料没有反应 5+3 Mr.Zot 2026-05-19 3/150 2026-05-19 13:56 by xtlilibin
[论文投稿] Sci. Bull. 悲剧经验 +4 jyang1999 2026-05-16 4/200 2026-05-19 10:56 by maolC
[基金申请] 同样的基金本子，换个专家直接从C变A！ +3 国自然国社科中 2026-05-19 3/150 2026-05-19 08:50 by Equinoxhua
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[基金申请] 青C资助名额大幅增加！ +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[考博] 光量子物理方向 博士招生 1人（2026.09） +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[有机交流] 求有机合成大神指点三硫酸乙烯酯（CAS：2793408-99-6）的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
关闭
请填处理意见
关闭