[問題]有關積分問題想請教大家一下
請教一下各位大大
剛接觸matlab不久,忽然要寫出這樣的function對我而言有些吃力
對一個連乘的標準常態(theta*x)*Gammapdf做0.001到T的積分
我寫的 function 如下,但是一直跑不出結果,請各位大大能為我指點一下
萬分感謝
function [Pd Ps]=defaul_survival_prob(T,nu,sigma,theta,)
N=length(theta);
for j=1:N
if j==1
phi(j)=@(x)normcdf((theta(j).*x)./(sigma(j).*sqrt(x)));
else j>1;
phi(j)=@(x)normcdf((theta(j).*x)./(sigma(j).*sqrt(x)))...
*normcdf((Li(j-1)-m(j-1)*T-VG_theta(j-1).*x)./(sigma(j-1).*sqrt(x)));
end
end
Fzd=@(phi)phi(N).*gampdf(T/nu,1/nu);
Pd=quad(Fzd,0.001,T);
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 220.132.65.140
推
02/26 10:18, , 1F
02/26 10:18, 1F
→
02/26 10:19, , 2F
02/26 10:19, 2F
→
02/26 10:38, , 3F
02/26 10:38, 3F