请问用matlab怎么解这个定积分方程组
发布网友
共2个回答
热心网友
这个积分比较复杂，得不到解析解的表达式
所以用符号运算给不出解
可以考虑用数值办法求解
B=15;A=0.2;cv=5;r=0.07;
f=@(t,x) 1/B*log(B./(cv+x*exp(r*t)));
f1=@(T,x) quad(@(t) f(t,x),0,T)-A;
f2=@(T,x) T-1/r*log((B-cv)/x);
fun=@(X) [f1(X(1),X(2));f2(X(1),X(2))];
options=optimset('TolFun',1e-10,'TolX',1e-10);
[s v]=fsolve(fun,[1;1],options)
结果是
s =
11.9392
4.3355
v =
1.0e-016 *
-0.5551
0
也即是
T=11.9391571887040
x=4.335530988603708
将方程等式移项表示为
取这两个值的时候上边两个表达式的值都小于1e-16
热心网友
第二行不是int\r\n是syms追问我是刚开始用matlab的， 因为论文需要解这个方程组，不是特别清楚你说的是什么意思。 能不能麻烦更加具体点？可以的话能不能写出代码?
