请老师们帮忙看一下我这个二重积分计算不出来 NaN

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

这是一个计算二元t-copula 分布函数的程序：

syms x y
ro=1;v=13.47;
a=tinv(0.204,v);
b=tinv(0.211,v);
fun=@(x,y) 1/(2*pi*sqrt(1-ro.^2))*(1+(x.^2-2*ro*x.*y+y.^2)./(v*(1-ro^2))).^(-(v+2)/2);
y =integral2(fun,-inf,a,-inf,b)

结果： NaN   警告: 遇到无限值或非数字值。
> In integralCalc/iterateScalarValued (line 349)
  In integralCalc/vadapt (line 132)
  In integralCalc (line 91)
  In integral2Calc>@(xi,y1i,y2i)integralCalc(@(y)fun(xi*ones(size(y)),y),y1i,y2i,opstruct.integralOptions)
  In integral2Calc>@(x)arrayfun(@(xi,y1i,y2i)integralCalc(@(y)fun(xi*ones(size(y)),y),y1i,y2i,opstruct.integralOptions),x,ymin(x),ymax(x)) (line 17)
  In integralCalc/iterateScalarValued (line 314)
  In integralCalc/vadapt (line 132)
  In integralCalc (line 91)
  In integral2Calc>integral2i (line 20)
  In integral2Calc (line 7)
  In integral2 (line 106)
  In Jt (line 6)
警告: 积分未成功。
> In integral2 (line 109)
  In Jt (line 6)

学胜刚接触matlab不久，请各位老师指点一下！ 感谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：二重积分 Integral arrayfun options struct

请问你解决了吗我也出现同样的问题了，可以一起交流下吗