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