用隐式函数作为下界来评估积分

时间:2014-07-10 16:16:55

标签: matlab

我有这个功能:

 psi = @(x,y)(-(11*x^(6/5)*y^(1/5) - 19*x + 9)/(10*x^(1/10)*y^(1/10)))

解决x只能隐式完成:

 R_sp = @(y)solve(psi(x,y), x);

我需要用R_sp作为下限评估以下双积分:

 int(s * int(s*t, t, R_sp(s), 1), s, .3, .4)

评估时,matlab说:MuPAD命令出错:需要列表或条目表。 [(DOM ::矩阵(DOM :: ExpressionField()))::创建]

R_sp的值在[.3,.4]中有明确定义,可以很容易地计算出来:

R_sp(.3) = 0.845993307852830503940553192774
R_sp(.4) = 0.89641082924669996315283636152907

所以我假设matlab的int函数不能处理该隐式函数中的连续值,对吧?当然在这种情况下,我可以使用trapz函数通过计算for循环中被积函数的离散值对来近似上面的积分。但除了使用陷阱近似(使用连续值)之外,还有其他任何方法可以在matlab中评估这个积分吗?

0 个答案:

没有答案