求教：simulation of a complex distribution - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Computation - 科学计算

求教：simulation of a complex distribution

求教：simulation of a complex distribution# Computation - 科学计算

t*t2010-01-09 08:01

1 楼

有人买过Norhtstar 5 listening and speaking 这书得第三版吗？想参加个ESL的培训
，5个礼拜level 5的培训，不知道效果如何，有上过类似这样的课的兄弟姐妹吗，对口
语提高有多少帮助。有用过的书可以卖给我，呵呵，跪求。。站内邮件等联系。

b*e2010-01-09 08:01

2 楼

按照信息产业部的说法
在中华人民共和国境内提供非经营性互联网信息服务，应当办理备案！未经备案，不得
在中华人民共和国境内从事非经营性互联网信息服务。而对于没有备案的网站将予以罚
款和关闭。
但是我糊涂了。是服务器在中国的都要备案，还是通过中国机构申请域名时候要备案，
还是原则上讲面向中国用户的都要备案？
我想建一个主要面向国内的站点。但是服务器想用在美国的。域名也想在这边直接申请
一个.com的。这样如果我不备案的话，即便内容不敏感，中国用户访问的时候也有可
能在某一天被墙挡住对吗？

q*z2010-01-09 08:01

3 楼

z*h2010-01-09 08:01

4 楼

f(x) is proptional to
exp{-sum(v_i*(x-b_i))-(x-c)^2/(2s)}/product{1+exp(b_i-x)}
i=1 to 20, v_i, b_i, c, s 已知。sum表示和，product表示积。
请问有什莫方法可以模拟这种复杂的分布？可以用ratio of uniforms吗？
多谢您的任何帮助！！！

W*k2010-01-09 08:01

5 楼

acceptance-rejection

【在 z*****h 的大作中提到】

: f(x) is proptional to
: exp{-sum(v_i*(x-b_i))-(x-c)^2/(2s)}/product{1+exp(b_i-x)}
: i=1 to 20, v_i, b_i, c, s 已知。sum表示和，product表示积。
: 请问有什莫方法可以模拟这种复杂的分布？可以用ratio of uniforms吗？
: 多谢您的任何帮助！！！

z*h2010-01-09 08:01

6 楼

多谢！是不是可以把它看作两部分：Normal+其他. 其他小于1。只要generate Normal
and Uniform? Thank tremendous!!!

z*h2010-01-09 08:01

7 楼

Sorry. I think the problem is that the function is not normalized and it is
impossible to calculate the normalization constant. Can I still use accept-
reject? Any help will be high appreciated,

t*s2010-01-09 08:01

8 楼

accept/reject 方法不需要做归一化
问题是你的PDF含有指数，如果最大值和最小值差别很大的话，用accept/reject方法效
率很低。可以考虑用monte carlo

is

【在 z*****h 的大作中提到】

: Sorry. I think the problem is that the function is not normalized and it is
: impossible to calculate the normalization constant. Can I still use accept-
: reject? Any help will be high appreciated,

z*h2010-01-09 08:01

9 楼

f f(x) need not to be normalized, can I use the following algorithm?
1. Generate X~N(c,s), U~U(0,1)
2. Accept Y=X if U<=1(or sqrt(2pi*s)?)
3. Return to 1, otherwise
I'm confused about how to determine g(x)/f(x). If f(x) is density function,
there will be no confusion. But here f(x) is proptional to ...
Any help will be highly appreciated,