24小时热门版块排行榜    

查看: 6652  |  回复: 9

serpenttom

新虫 (小有名气)

[交流] 如何生成多为柯西分布的随机数

概率密度函数如下:
其中m是维数
1/(1+x_1^2+...+x_m^2)^{(m+1)/2}

当m=1时,
CODE:
x=atan(rand(1)*pi-pi/2)

(Matlab代码)可以生成一维的柯西分布随机数,再后来,就傻眼了......望高人指点
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
https://www.mathworks.com/matlab ... r/view_thread/20751

r = rand(1);
x=atan(rand(1)*pi-pi/2)
2楼2016-08-15 13:06:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

serpenttom

新虫 (小有名气)

引用回帖:
2楼: Originally posted by FMStation at 2016-08-15 13:06:55
https://www.mathworks.com/matlabcentral/newsreader/view_thread/20751

r = rand(1);
x=atan(rand(1)*pi-pi/2)

首先谢谢您啦,那个帖子我是看过的,但它只是生成的多个一维柯西的随机数,而不是一个多维柯西的随机数

发自小木虫Android客户端
3楼2016-08-15 14:19:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
>> original_x=rand(1,100000);
>> cauchy_x=tan((original_x-1/2)*pi);
4楼2016-08-15 16:46:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主的分布密度函数并不是Cauchy分布;如果问的的确是Cauchy分布,那么一维情况下参数为时的概率密度函数



二维情况下参数时,如果新的的概率密度函数可以表示成两者的乘积:



那么,从直接思路的角度是求二维分布函数的反函数即可。这在一般的数学类工具软件中有标准的实现函数。关键是看到底要什么样的二维。如果两个维数并不独立,实现稍微复杂些。
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
5楼2016-08-15 16:58:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

serpenttom

新虫 (小有名气)

6楼2016-08-15 19:07:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

serpenttom

新虫 (小有名气)

引用回帖:
5楼: Originally posted by Mr__Right at 2016-08-15 16:58:29
楼主的分布密度函数并不是Cauchy分布;如果问的的确是Cauchy分布,那么一维情况下参数为a,b时的概率密度函数

\dfrac{1}{\pi\,b\left(\dfrac{(x-a)^2}{b^2}+1\right)}

二维情况下参数a_1,b_1;\quad a_2,b_2时, ...

首先十分感谢您的回答。你看一下柯西高维的概率分布密度函数,我发在了你的楼下,这是维基网的图片,向量各分量之间是不存在独立与否,两个随机变量才会涉及是否独立的问题,同时多变量函数是没有反函数可言不具有一一对应关系,不过还是感谢你的热情

发自小木虫Android客户端
7楼2016-08-15 19:12:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by serpenttom at 2016-08-15 19:07:44

Note that in this example, even though there is no analogue to a covariance matrix, and are not statistically independent

即使严格意义上不是,也算一种类比吧。“反函数”应该是不对的,不过实现的方式类似。
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
8楼2016-08-15 19:36:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

serpenttom

新虫 (小有名气)

9楼2016-10-07 18:38:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

丶米国总理

禁虫 (初入文坛)

本帖内容被屏蔽

10楼2018-05-22 18:33:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 serpenttom 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285化工学硕求调剂(081700) +11 柴郡猫_ 2026-03-12 11/550 2026-03-19 09:37 by laoshidan
[考研] 274求调剂 +6 S.H1 2026-03-18 6/300 2026-03-19 09:34 by 花店相见
[基金申请] 国自科面上基金字体 +7 iwuli 2026-03-12 8/400 2026-03-19 08:55 by sjtu2012
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
请填处理意见