Re: [問題] 常態pdf積分已回收
※ 引述《xisthebest (hitori)》之銘言：
: 不好意思 我有積分的問題想請教
: 首先I1的部分是對常態分配的PDF作積分
: 積分下界為負無窮 而積分上界為(a+sqrt(ro)*z)/sqrt(1-ro) 為一變數z的function
: 而val1為I1的k次方
: 接著I2的部分則是val1*normpdf(z,0,1) 變數仍為z 積分上下界分別正負無窮
: 照理說這東西積出來的值應為0.93656
: 但答案卻不對
: 能否請板友看看問題出在哪裡呢?? 感激~
: k =3;
: ro = 0.25;
: a = 2;
: syms z;
: I1 = int(normpdf(z,0,1),z,-1e+10,(a+sqrt(ro)*z)/sqrt(1-ro));
: val1 = I1^k;
: I2 = int(val1*normpdf(z,0,1),z,-1e+10,1e+10);
: vpa(I2)
k =3;
ro = 0.25;
a = 2;
syms x z;
I1 = int(normpdf(x,0,1), x, -Inf, (a+sqrt(ro)*z)/sqrt(1-ro));
val1 = I1^k;
I2 = int(val1*normpdf(z,0,1), z, -Inf, Inf);
vpa(I2)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.94.230
討論串 (同標題文章)