| 查看: 1317 | 回复: 1 | |||||
leking137joey金虫 (正式写手)
材料小生
|
[求助]
Matlab的熟练用户或者大神看过来，两个随机变量存在相互关系的简单问题。重金求助！
|
|
大神们好，本人Matlab小白，今天才开始学。有一个问题我不知道怎么着手，诚心请指教。问题描述如下： 有一种针状纳米粒子，通过实验测得其具有一定的长度分布y【单位nm】，如下 y=[21.55, 24.05, 27.82, 29.26, 32.6, 34.24, 35.06, 36, 37.11, 38.08, 39.55, 39.63, 41.24, 42.47, 42.87, 44.66, 44.94, 46.18, 46.62, 46.93, 48.52, 49.14, 50.7, 52.65, 53.8, 55.66, 56.35, 56.65, 57.44, 58.87, 61.12, 62.92, 62.96, 63.55, 64.93, 66.83, 68.48, 69.68, 72.1, 75.02, 77.1, 77.49, 78.82, 79.26, 79.85, 81.14, 81.18, 82.67, 85.69, 88.75, 89.83, 90.13, 90.86, 91.14, 95.96, 96.47, 96.53, 96.53, 97.47, 104.99, 105.94, 108.76, 108.76, 109.07, 109.63, 111.36, 112.15, 112.75, 114.42, 116.13, 119.81, 123.18, 124.35, 125.3, 128.12, 130.39, 131.72, 133.74, 135.43, 135.47, 136.26, 136.93, 139.23, 144.48, 145.32, 150.83, 155.62, 158.91, 173.96, 182.82, 187.27, 187.67, 206.03, 225.03, 227.02, 230.33, 245.05, 246.98, 247.6, 276.26, 291.42, 302.53]; [F,X] = ksdensity(y) [N,C] = hist(y,15) plotyy(C,N,X,F,'bar','plot') 在matlab中运算后，得到附件所示的图片，曲线表示概率分布函数，直方图表示一定尺寸范围内的粒子数。 现在另一组电镜实验中某一区域【厚度t为302 nm】观测到282个粒子，由于观测方向平行于针的方向，因此在这个电镜图片中无法知道每一个粒子的长度。因此我打算通过分布函数来赋282个随机长度的y值，但是这些针状颗粒的中心点在厚度方向的坐标z与l，t必须满足以下关系：|z|-0.5*y<0.4*t=120.8，才能服从“在电镜下可见”这一事实。也就是在|z|-0.5*y<120.8的限制条件下，随机产生z值，并依据分布函数随机产生y值，z和y的组合数目总共为282个。请问这个在Matlab中如何实现？ 求大神们指点！如能给出运算式子，感激不尽！最佳答案100金币奉上，其他答案亦会追加金币作为感谢！  在线等  untitled.jpg |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
这样的filecode谁见过
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
关于代码变化问题，想知道的进来
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
mygt_hit
专家顾问 (职业作家)
-
专家经验: +362 - 应助: 438 (硕士)
- 贵宾: 0.019
- 金币: 19579.9
- 散金: 5130
- 红花: 135
- 沙发: 3
- 帖子: 4939
- 在线: 990.5小时
- 虫号: 1489764
- 注册: 2011-11-13
- 性别: GG
- 专业: 结构工程
- 管辖: 土木建筑
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
leking137joey: 金币+20, ★★★很有帮助, 谢谢，请问“判断|z|-0.5*y<120.8，如果符合就接受，不符合就继续对Y采样”应该用什么语句或者函数来实现呢？ 2015-12-20 14:03:07
leking137joey: 金币+100, ★★★★★最佳答案, 按照你的思路，问题解决了，金币奉上，小小意思，不成敬意。 2016-01-08 18:36:22
感谢参与，应助指数 +1
leking137joey: 金币+20, ★★★很有帮助, 谢谢，请问“判断|z|-0.5*y<120.8，如果符合就接受，不符合就继续对Y采样”应该用什么语句或者函数来实现呢？ 2015-12-20 14:03:07
leking137joey: 金币+100, ★★★★★最佳答案, 按照你的思路，问题解决了，金币奉上，小小意思，不成敬意。 2016-01-08 18:36:22
|
先生成Z的随机数，再生成Y的随机数，再加判断|z|-0.5*y<120.8，如果符合就接受，不符合就继续对Y采样，直到符合要求。这样循环282次，就可以得到所需的样本对(Z Y)了。 大致的想法，最好数学上推导一下，感觉像拒绝采样的内容。 |
2楼2015-12-19 11:58:20
100