请问用matlab怎么解这个定积分方程组只有两个未知数T和π0, β=15 A=0.2 c+v=5 r=0.07.下面是我写的代码,提示没有解.matlab新手,还请懂的人麻烦写一下具体的代码,谢谢! syms T x>> [T,x]=solve('0.2=1

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/24 02:56:57
请问用matlab怎么解这个定积分方程组只有两个未知数T和π0, β=15 A=0.2 c+v=5 r=0.07.下面是我写的代码,提示没有解.matlab新手,还请懂的人麻烦写一下具体的代码,谢谢! syms T x>> [T,x]=solve('0.2=1

请问用matlab怎么解这个定积分方程组只有两个未知数T和π0, β=15 A=0.2 c+v=5 r=0.07.下面是我写的代码,提示没有解.matlab新手,还请懂的人麻烦写一下具体的代码,谢谢! syms T x>> [T,x]=solve('0.2=1
请问用matlab怎么解这个定积分方程组

只有两个未知数T和π0, β=15 A=0.2 c+v=5 r=0.07.

下面是我写的代码,提示没有解.matlab新手,还请懂的人麻烦写一下具体的代码,谢谢!

 syms T x

>> [T,x]=solve('0.2=1/3*exp(-7/100*T)*log(5)+1/3*exp(-7/100*T)*log(1/(5+exp(7/100*T)*T))+1/15*exp(-7/100*T)*exp(7/100*T)*T*log(3)+1/15*exp(-7/100*T)*exp(7/100*T)*T*log(5)+1/15*exp(-7/100*T)*exp(7/100*T)*T*log(1/(5+exp(7/100*T)*T))+1/15*exp(-7/100*T)*exp(7/100*T)*T','T=(log((15-5)/x))/0.07')

 

 

Warning: Explicit solution could not be found.

> In solve at 140

 

T =

 

[ empty sym ]

 

 


x =


     []


悬赏10分是因为只有那么多了。。如果有可以解答的本人愿意支付宝转账50RMB

请问用matlab怎么解这个定积分方程组只有两个未知数T和π0, β=15 A=0.2 c+v=5 r=0.07.下面是我写的代码,提示没有解.matlab新手,还请懂的人麻烦写一下具体的代码,谢谢! syms T x>> [T,x]=solve('0.2=1
这个积分比较复杂,得不到解析解的表达式
所以用符号运算给不出解
可以考虑用数值办法求解


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.939157164887040 
x=4.335530988603708
将方程等式移项表示为

取这两个值的时候上边两个表达式的值都小于1e-16

第二行不是int\r\n是syms