24小时热门版块排行榜    

Znn3bq.jpeg
查看: 6683  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +19 硕星赴 2026-04-09 20/1000 2026-04-10 05:27 by chenzhimin
[考研] 求调剂 +6 翩翩一书生 2026-04-09 6/300 2026-04-10 05:14 by chenzhimin
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 一志愿中国科学院上海有机所,有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 086000调剂 +4 十七sa 2026-04-07 4/200 2026-04-09 09:27 by 壹往無前
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 081700,311,求调剂 +17 冬十三 2026-04-04 18/900 2026-04-07 12:50 by Sammy2
[考研] 材料工程310专硕调剂 +14 捞捞我…. 2026-04-04 15/750 2026-04-06 14:18 by lqwchd
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 334求调剂 +8 曾仰之 2026-04-03 8/400 2026-04-04 11:16 by w_xuqing
[考研] 26调剂 086003 +6 失活的细胞 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[考研] 357求调剂 +13 1050389037 2026-04-03 13/650 2026-04-03 22:27 by 无际的草原
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
信息提示
请填处理意见