Question

我正在尝试在Matlab中解决x的等式，但不断得到错误：

空sym：0-by-1

该等式具有积分，其中x是上限，也是被积函数1的一部分。我使用的代码如下：

a = 0.2; b= 10; c = -10; d = 15; mu = 3; sig = 1;
syms x t
eqn = 0 == a + b*normcdf(x,mu,sig)+c*int( normcdf(d + x - t,mu,sig)*normpdf(t,mu,sig),t,0,x);
A = vpasolve(eqn,x)

关于我错在哪里的任何提示？

Answer 1

我认为符号工具箱可能不足以解决那个问题...也许某些assume或其他一些技巧可以完成这项工作，我个人无法找到方法。

但是，为了测试这是否可以解决，我尝试了Wolfram Alpha。它给出了一个可以使用的结果。

eq1=a + b*normcdf(x,mu,sig);
resint=c*(t^3*(d - t + x)*erfc((mu - x)/(sqrt(2)*sig)))/(4*sig*exp((-mu + x)^2/(2*sig^2))*sqrt(2*pi));
A=vpasolve(eq1+subs(resint,t,x)-subs(resint,t,0) ==0)

在我的电脑中提供1.285643225712432599485355373093。

在MATLAB中求解x的积分，其中x是绑定的并且是被积函数的一部分

1 个答案: