作者reddragon00 (00)
看板MATLAB
标题[问题]有关积分问题想请教大家一下
时间Tue Feb 26 09:58:13 2013
请教一下各位大大
刚接触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
1F：推 sunev:Nonscalar arrays of function handles are not allowed; 02/26 10:18
2F：→ sunev: use cell arrays instead. 02/26 10:19
3F：→ reddragon00:可以请教一下如何改善这个问题吗? 02/26 10:38