Unable to find explicit solution

조회 수: 214 (최근 30일)
Berkay Can Tuncel
Berkay Can Tuncel 2022년 1월 19일
댓글: Yongjian Feng 2022년 1월 20일
Hello,
I am trying to solve an equation for one variable but I get this error message:
Warning: Unable to find explicit solution. For options, see help.
> In sym/solve (line 317)
In untitled (line 15)
sol =
Empty sym: 0-by-1
>>
I searched but could'nt find any solution. I couldn't even understand what the problem is. You can see my code below.
Thanks in advance.
syms t
b=2*t+14.3;
w=3*t+21.45;
A=(w*b)-(w-t)*(b-2*t);
rc=(w*b)*450-(w-t)*(b-2*t)*(450+t/2)/A;
rn=A/(b*log((450-w/2+t)/(450-w/2))+2*w*log((450+w/2)/(450-w/2+t)));
E=rc-rn;
ro=450+w/2;
ri=450-w/2;
Ay=-40260.604300;
M1=56124.523326;
sol=solve((Ay/A)+(M1*(ro-rn))/(A*E*ro)==116, t)
  댓글 수: 1
Torsten
Torsten 2022년 1월 19일
편집: Torsten 2022년 1월 19일
Not all equations can explicitly be solved.
Use "vpasolve" instead of "solve".
If this doesn't yield a result, plot
(Ay/A)+(M1*(ro-rn))/(A*E*ro)-116
for a certain t-interval to see whether a zero really exists.

댓글을 달려면 로그인하십시오.

답변 (3개)

Yongjian Feng
Yongjian Feng 2022년 1월 19일
This equation might not even have a solution. Try to plot it:
syms t
b=2*t+14.3;
w=3*t+21.45;
A=(w*b)-(w-t)*(b-2*t);
rc=(w*b)*450-(w-t)*(b-2*t)*(450+t/2)/A;
rn=A/(b*log((450-w/2+t)/(450-w/2))+2*w*log((450+w/2)/(450-w/2+t)));
E=rc-rn;
ro=450+w/2;
ri=450-w/2;
Ay=-40260.604300;
M1=56124.523326;
x = 1:100;
y = feval(matlabFunction((Ay/A)+(M1*(ro-rn))/(A*E*ro)), x);
plot(x, y)
It looks like when t goes up, (Ay/A)+(M1*(ro-rn))/(A*E*ro)) only goes up from negative to 0. Not sure it will ever reach 116.
  댓글 수: 2
Yongjian Feng
Yongjian Feng 2022년 1월 19일
Try vpasolve as recommended by Torsten above, your equation results in empty solution. Most likely there is not solution for it.
If you plot, you can see (Ay/A)+(M1*(ro-rn))/(A*E*ro)) is always negative. It might never reach 116 as you want.
sol=vpasolve((Ay/A)+(M1*(ro-rn))/(A*E*ro)==0, t) % this gives empty solution. This equation never reaches positive
sol=vpasolve((Ay/A)+(M1*(ro-rn))/(A*E*ro)==-1, t) % this gives 77.287
Yongjian Feng
Yongjian Feng 2022년 1월 20일
Walter is right. I forgot about t<0 range.
If you plot -100:100, you can see most likely there shall be only one solution.

댓글을 달려면 로그인하십시오.


Walter Roberson
Walter Roberson 2022년 1월 19일
편집: Walter Roberson 2022년 1월 19일
If you look carefully at the graph, you can see a couple of discontinuities at negative t values.
At first I thought it was a simple case of the system having non-zero imaginary components in some ranges, but it turns out that there are some steep +/- infinities being generated -- and each of those is an opportunity for a solution.
I saw at least one other candidate stretch for a solution between t = 0 and t = -10
syms t
b=2*t+14.3;
w=3*t+21.45;
A=(w*b)-(w-t)*(b-2*t);
rc=(w*b)*450-(w-t)*(b-2*t)*(450+t/2)/A;
rn=A/(b*log((450-w/2+t)/(450-w/2))+2*w*log((450+w/2)/(450-w/2+t)));
E=rc-rn;
ro=450+w/2;
ri=450-w/2;
Ay=-40260.604300;
M1=56124.523326;
eqn = (Ay/A)+(M1*(ro-rn))/(A*E*ro) == 116
eqn = 
sol=solve(eqn, t)
Warning: Unable to find explicit solution. For options, see help.
sol = Empty sym: 0-by-1
string(lhs(eqn) - rhs(eqn))
ans = "5533375321142303/(137438953472*((143*t)/5 - (2*t + 143/10)*(3*t + 429/20) + 61347/200)) - ((23141087250120879*t)/274877906944 + (7713695750040293*((143*t)/5 - (2*t + 143/10)*(3*t + 429/20) + 61347/200))/(137438953472*(log((t/2 - 17571/40)/((3*t)/2 - 17571/40))*(2*t + 143/10) + log(-((3*t)/2 + 18429/40)/(t/2 - 17571/40))*(6*t + 429/10))) + 142155698977492559697/5497558138880)/(((3*t)/2 + 18429/40)*(((143*t)/5 - (2*t + 143/10)*(3*t + 429/20) + 61347/200)/(log((t/2 - 17571/40)/((3*t)/2 - 17571/40))*(2*t + 143/10) + log(-((3*t)/2 + 18429/40)/(t/2 - 17571/40))*(6*t + 429/10)) + 450*(2*t + 143/10)*(3*t + 429/20) + ((t/2 + 450)*((143*t)/5 + 61347/200))/((143*t)/5 - (2*t + 143/10)*(3*t + 429/20) + 61347/200))*((143*t)/5 - (2*t + 143/10)*(3*t + 429/20) + 61347/200)) - 116"
vpasolve(eqn, [-7.985, -7.98])
ans = 

Alex Sha
Alex Sha 2022년 1월 20일
One solution seems to be:
t: -8.33315732250827

카테고리

Help CenterFile Exchange에서 Calculus에 대해 자세히 알아보기

제품


릴리스

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by