| 查看: 663 | 回复: 5 | |||
| 当前主题已经存档。 | |||
H_LG银虫 (正式写手)
|
[交流]
【求助】帮忙看个程序怎么实现，谢谢
|
||
|
draw j0 at random with probability proportional to exp{}c^2/2sigma^2*|dj'r^(k-1)|^2/||dj||. 如上，是什么意思呢？用matlab 如何实现选择这个j0？ 请高手解答，多谢！！ |
回复此楼» 猜你喜欢
国自然评审意见
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
2026年叶企孙基金
已经有7人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
H_LG
银虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 157.8
- 散金: 228
- 帖子: 334
- 在线: 82.3小时
- 虫号: 579678
- 注册: 2008-07-11
- 性别: GG
- 专业: 计算数学与科学工程计算
高手来吧
2楼2009-11-07 11:51:02
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
★
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
★ ★ ★ ★
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
