24小时热门版块排行榜    

查看: 2344  |  回复: 20

wsry0112

银虫 (小有名气)

[求助] 二维正态概率密度函数

二维正态概率密度函数,而且变量之间相互独立,求概率,为什么mathmatic无法算出计算结果。这么经典的问题,为什么计算不出来呢
二维正态概率密度函数
计算公式


二维正态概率密度函数-1
计算结果
回复此楼
不是因为有希望而坚持,而是因为坚持才会有希望。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

feixiaolin

荣誉版主 (文坛精英)

优秀版主

mathmatics没有分离变量的智慧。
2楼2013-11-25 13:17:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wsry0112: 金币+30, ★★★很有帮助, 大哥,我真的没做出来,咱先预付这些,等出结果一定悉数奉上 2013-11-25 21:54:11
1) 直接计算
Integrate[exp[-(x-1)^2/2]*Integrate[ exp[-(y-1)^2/8], {y, 0, x} ]), {x, 0, Infinity} ]
结果出不来!

2)变两分离,分步计算
Integrate[ exp[-(y-1)^2/8], {y, 0, x}
结果为sqrt(2*pi)*(erf((x-1)/sqrt(8))+ erf(1/sqrt(8)))
进一步有:
sqrt(2*pi)*(erf((x-1)/ 2.828427125)+ 0.38292492255)

3)继续分步计算
a=Integrate[exp[-(x-1)^2/2]*sqrt(2*pi)*erf((x-1)/ 2.828427125)], {x, 0, Infinity} ]= 0.541176
b=Integrate[exp[-(x-1)^2/2]*sqrt(2*pi)* 0.38292492255], {x, 0, Infinity} ]= 2.0242655708

4)积分结果
Sum(a, b)/(4*pi)=0.204151353603134
6楼2013-11-25 17:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

hank612

至尊木虫 (著名写手)

引用回帖:
2楼: Originally posted by feixiaolin at 2013-11-25 13:17:43
mathmatics没有分离变量的智慧。

feixiaolin大神,这次你看走眼了,这个双重积分没法分离变量,积分区域是个三角形 0<y<x<Infinity
We_must_know. We_will_know.
3楼2013-11-25 13:57:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wsry0112

银虫 (小有名气)

引用回帖:
2楼: Originally posted by feixiaolin at 2013-11-25 13:17:43
mathmatics没有分离变量的智慧。

其他软件可以计算吗
不是因为有希望而坚持,而是因为坚持才会有希望。
4楼2013-11-25 14:28:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
3楼: Originally posted by hank612 at 2013-11-25 13:57:46
feixiaolin大神,这次你看走眼了,这个双重积分没法分离变量,积分区域是个三角形 0<y<x<Infinity...

不是大神,是大胆。
本身非数学专业,胆子大的过了一些。
5楼2013-11-25 15:55:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wsry0112

银虫 (小有名气)

引用回帖:
6楼: Originally posted by feixiaolin at 2013-11-25 17:12:56
1) 直接计算
Integrate[exp*Integrate), {x, 0, Infinity} ]
结果出不来!

2)变两分离,分步计算
Integrate[ exp, {y, 0, x}
结果为sqrt(2*pi)*(erf((x-1)/sqrt(8))+ erf(1/sqrt(8)))
进一步有:
sqrt(2* ...

a=Integrate[exp[-(x-1)^2/2]*sqrt(2*pi)*erf((x-1)/ 2.828427125)], {x, 0, Infinity} ]= 0.541176 这一步计算不出来啊,显示Warning: Explicit integral could not be found.  大侠,你是怎么计算出来的啊?
不是因为有希望而坚持,而是因为坚持才会有希望。
7楼2013-11-25 20:45:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

看图,就那样做出来滴。
记住付我报酬50金币。
二维正态概率密度函数-2
4U.JPG

8楼2013-11-25 21:30:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wsry0112

银虫 (小有名气)

引用回帖:
8楼: Originally posted by feixiaolin at 2013-11-25 21:30:53
看图,就那样做出来滴。
记住付我报酬50金币。

4U.JPG

大哥,你这是在哪个软件里算的?
不是因为有希望而坚持,而是因为坚持才会有希望。
9楼2013-11-25 21:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wsry0112

银虫 (小有名气)

引用回帖:
8楼: Originally posted by feixiaolin at 2013-11-25 21:30:53
看图,就那样做出来滴。
记住付我报酬50金币。

4U.JPG

大哥啊,我在mathmatic和matlab都试过了,还是不行啊
不是因为有希望而坚持,而是因为坚持才会有希望。
10楼2013-11-25 22:07:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wsry0112 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 05:51 by LH5NkK5Ud4bh
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:35 by LH5NkK5Ud4bh
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 04:47 by LH5NkK5Ud4bh
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:11 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-16 03:23 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +6 两三岁ss 2026-09-14 6/300 2026-09-15 20:17 by 只吃肉葡萄
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-15 18:11 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
信息提示
请填处理意见