Re: [問題] 常態pdf積分已回收
十分感激 問題已解決
還有另一類似問題想請教:
首先I1的部分同樣是對常態分配的PDF作積分
積分下界為負無窮 而積分上界為(a+sqrt(ro)*z)/sqrt(1-ro) 為一變數z的function
而val1為I1的k-1次方;
接著I2的部分則是對常態分配的PDF作積分 變數仍為z 積分下界為負無窮 而積分上界
為(g+sqrt(ro)*z)/sqrt(1-ro);
接著I3的部分則是val1*I2*normpdf(z,0,1) 變數仍為z 積分上下界分別正負無窮
這積分出來的結果應為一D的Function 但答案卻是以下的情況:
numeric::int((642775217703596110216784836936465041008881197513117134
1205504*2^(1/2)*pi^(3/2)*(erf((2^(1/2)*5^(1/2)*14^(1/2)*((3*14^(1/2)*z)/14 -
(5*D)/28..............以下省略
答案中同時包含了z和D 很明顯是不對的 請問問題出在哪呢?
先謝過了
a = 2;
b = 2;
u = 3;
syms D;
g = -D*b+a;
syms x z;
I1 = int(normpdf(x,0,1), x, -Inf, (a+sqrt(ro)*z)/sqrt(1-ro));
val1 = I1^(k-1);
I2 = int(normpdf(x,0,1), x, -Inf, (g+sqrt(ro)*z)/sqrt(1-ro));
I3 = int(val1*I2*normpdf(z,0,1), z, -Inf, Inf);
vpa(I3)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.96.50
→
04/09 17:29, , 1F
04/09 17:29, 1F
→
04/09 18:10, , 2F
04/09 18:10, 2F
討論串 (同標題文章)
完整討論串 (本文為第 3 之 3 篇)：
問題
0
4