Solve issue in Matlab R2018b

9 次查看（过去 30 天）
John Jarvis
John Jarvis 2019-3-28
评论： Star Strider 2019-3-28
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 个评论
显示 1更早的评论
madhan ravi
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
John Jarvis 2019-3-28
编辑：John Jarvis 2019-3-28
@Image Analyst:
The exact errors are as follows (all the red text)
Error in solve>getEqns (line 429)
[eqns, vars] = sym.getEqnsVars(argv{:});
Error in solve (line 226)
[eqns,vars,options] = getEqns(varargin{:});
Error in SolutionsPlot (line 11)
sol = solve(eqn,a,[0 pi]);

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Star Strider
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 个评论
显示 2更早的评论
John Jarvis
John Jarvis 2019-3-28
@Star Strider, Many thanks for your elaboration.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by