Solve Equation numerically for variable

The following error occurred converting from sym to double: Error using symengine (line 59) DOUBLE cannot convert the input expression into a double array. If the input expression contains a symbolic variable, use VPA.

Error in Post_Final (line 40) HourAngleSunrise_Degrees(i) = rad2deg* solve('A + B*cos(w)+C*sin(w)')

A, B, C are values calculated from matrixes, so A = sin(matrix(i)) etc

Philip Hoskinson 2016-7-21

在 MATLAB Online 中打开

Original question had an error that makes this more complicated: Actual equation:

 0 = A + Bcos(w) + *Csin(w)*
or with identity:
 0 = A = Bcos(w) + Csqrt(1-cos^2(w))

请先登录，再进行评论。

Answer 2

Star Strider 2016-7-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/296665-solve-equation-numerically-for-variable#answer_229387

在 MATLAB Online 中打开

If ‘A’, ‘B’ and ‘C’ are matrices, then:

w = acos(-(B+C)\A);

assuming that the matrices have the appropriate dimensions, and the resulting argument matrix ‘(-(B+C)\A)’ elements are all between -1 and +1.

2 个评论
显示无隐藏无

Philip Hoskinson 2016-7-21

在 MATLAB Online 中打开

Original question had an error that makes this more complicated: Actual equation:

 0 = A + Bcos(w) + *Csin(w)*
or with identity:
 0 = A = Bcos(w) + Csqrt(1-cos^2(w))

Star Strider 2016-7-21

在 MATLAB Online 中打开

The Symbolic Math Toolbox comes up with these solutions for ‘w’ that would also work for matrices:

syms A B C w
w_sol = solve(A == B*cos(w) + C*sin(w), w);
w_sol = simplify(w_sol, 'steps',10)
   -log((A + (A^2 - B^2 - C^2)^(1/2))/(B - C*1i))*1i
   -log((A - (A^2 - B^2 - C^2)^(1/2))/(B - C*1i))*1i

Use the sqrtm and logm functions. Also see the funm function. If you know that your matrices will always be real or positive, you can tell the Symbolic Math Toolbox to assume that. It may simplify the result.

请先登录，再进行评论。