| 查看: 489 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
H_LG银虫 (正式写手)
|
[交流]
【求助】帮忙看个程序怎么实现，谢谢
|
||
|
draw j0 at random with probability proportional to exp{}c^2/2sigma^2*|dj'r^(k-1)|^2/||dj||. 如上，是什么意思呢？用matlab 如何实现选择这个j0？ 请高手解答，多谢！！ |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
★
H_LG(金币+1,VIP+0):这个很接近了，能否提供一个matlab实现呢？ 11-11 09:55
H_LG(金币+1,VIP+0):这个很接近了，能否提供一个matlab实现呢？ 11-11 09:55
|
这是用已知概率密度函数p产生一组样本使其服从这个概率分布吧? 如果是可以这么求: 1.确定样本取值范围及数值精度.如0~255,精度为1.也就是样本都是整数. 2.根据已知的密度函数p确定每个精度范围内[n, n+1]的分布的比例f_n,即求区间分布函数,也就是将密度函数p在[n,n+1]上积分. 做计算的时候用梯形公式逼近. 由于概率密度函数在其定义域内积分为1, 所以算的这些分布比例f_n加和为1.当然计算上近似为1. 3. 取定要取的样本个数m, 确定取值为n的数据的个数:m*f_n; 4. 将3中产生的数合在一起, 对其空间位置进行随机扰动(matlab用函数randperm可以实现). 这样就产生了m个样本满足概率密度函数p. 如果是数学上一些典型的分布,在matlab中可由函数random加一些参数直接产生. e的右大括号在哪里? 如果在最后, 括号里应该还有个负号. |
5楼2009-11-10 12:35:53
3楼2009-11-07 14:06:48
H_LG
银虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 157.8
- 散金: 228
- 帖子: 334
- 在线: 82.3小时
- 虫号: 579678
- 注册: 2008-07-11
- 性别: GG
- 专业: 计算数学与科学工程计算
4楼2009-11-09 16:17:14
★ ★ ★ ★
bslt(金币+3,VIP+0):感谢热心回复，欢迎常来！ 11-14 10:09
H_LG(金币+1,VIP+0):这个是我真正想要的，非常感谢，我回去测试下 11-14 11:49
bslt(金币+3,VIP+0):感谢热心回复，欢迎常来！ 11-14 10:09
H_LG(金币+1,VIP+0):这个是我真正想要的，非常感谢，我回去测试下 11-14 11:49
|
%%This is a Matlab programme which can produce a sample %%with any given probability density function (PDF). clear; clc; %some parameters. sam_min=0; sam_max=255; pre=1.0; sam_num=50000; sam_x=sam_min:pre:sam_max; %two parameters for the PDF. sig=[15 20]; mea=[80 150]; %PDF. e.g. Gaussian mixture. pro=0.3*1/sqrt(2*pi*sig(1)^2).*exp(-(sam_x-mea(1)).^2/(2*sig(1)^2))+... 0.7*1/sqrt(2*pi*sig(2)^2).*exp(-(sam_x-mea(2)).^2/(2*sig(2)^2)); %distribution function F=pro*pre; F_num=round(F*sam_num); %producing the sample. sam=[]; for ind=1:length(sam_x) sam=[sam repmat(sam_x(ind),[1 F_num(ind)])]; end sam(randperm(length(sam)))=sam(1:end); %test. figure, his=hist(sam,sam_x); plot(sam_x,his/(pre*sum(his)),'*'); hold on; plot(sam_x,pro,'r'); legend('PDF of the sample','Theoretic PDF'); |
6楼2009-11-11 12:55:03