Solve issue in Matlab R2018b
显示 更早的评论
I was using Matlab R2014a earlier where I was able to solve this equation easily, but in Matlab R2018b it is returning the
error in sol = solve(eqn,a,[0 pi]); I couldn't figure out why. The working code in Matlab 2014a is
syms a T
v2=-2.3750
g=1;
b=0;
e2=0.5;
k=2.5;
w=-2*cos(k);
eqn = sin(3*k+a)./sin(2*k+a)==v2-w+(g.*T.^2)+(e2.*T.^2.*sin(k)^2)./(sin(2*k+a)^2+b*T.^2*sin(k).^2);
sol = solve(eqn,a,[0 pi]);
digits(5)
solutions = vpa(subs(sol),3)
Please note that "a" is to be bound to take values between 0 and pi.
3 个评论
Image Analyst
2019-3-28
What is the exact error (ALL the red text)? Have you tried it with R2019a?
madhan ravi
2019-3-28
Only you can specify the range in vpasolve() where only one parameter is not known but with solve you can't specify the bounds.
John Jarvis
2019-3-28
编辑：John Jarvis
2019-3-28
采纳的回答
Star Strider
2019-3-28
Restrict ‘a’ using an assume call.
Try this:
syms a T
assume(a >= 0 & a <= pi)
v2=-2.3750
g=1;
b=0;
e2=0.5;
k=2.5;
w=-2*cos(k);
eqn = sin(3*k+a)./sin(2*k+a)==v2-w+(g.*T.^2)+(e2.*T.^2.*sin(k)^2)./(sin(2*k+a)^2+b*T.^2*sin(k).^2);
[sol,prms,conds] = solve(eqn,a, 'ReturnConditions',true)
digits(5)
solutions = vpa(subs(sol),3)
vpaconds = vpa(conds)
4 个评论
John Jarvis
2019-3-28
Star Strider
2019-3-28
My pleasure.
The ‘conds’ result is likely as good as it gets. Note that it is a function of ‘T’, and with the inequality conditions stated in the vector, neither coeffs or any related function (such as numden or sym2poly) will produce a purely numeric result.
Adding:
Tsol(:,1) = solve(conds(1), T)
Tsol(:,2) = solve(conds(2), T)
doesn’t completely resolve the problem, since the result is now a funciton of ‘x’:
Tsol =
-2*exp(-5i/4)*(-((exp(x*2i)*exp(10i) - 1)*(4478027322974943*exp(5i/2) - 1125899906842624*exp(x*2i)*exp(15i) - 4478027322974943*exp(x*2i)*exp(25i/2) + 1125899906842624))/(4503599627370496 + 4503599627370496*exp(x*4i)*exp(20i) - 12233297969360201*exp(x*2i)*exp(10i)))^(1/2)
2*exp(-5i/4)*(-((exp(x*2i)*exp(10i) - 1)*(4478027322974943*exp(5i/2) - 1125899906842624*exp(x*2i)*exp(15i) - 4478027322974943*exp(x*2i)*exp(25i/2) + 1125899906842624))/(4503599627370496 + 4503599627370496*exp(x*4i)*exp(20i) - 12233297969360201*exp(x*2i)*exp(10i)))^(1/2)
The only other option I can provide is:
Tfcn = matlabFunction(Tsol)
that will produce an anonymous function to evaluate those.
John Jarvis
2019-3-28
Star Strider
2019-3-28
As always, my pleasure.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
