Orbit Determination with Matlab Revisited

[stupid6625]

I has posted on here about 6 months ago about trying to create a MATLAB code for determining the orbital elements of a body about Earth from three separate observing times . I am having trouble with some of the code. The very last command is the "solve" command and it is supposed to solve the two equations in parenthesis, but it returns an error saying "unable to find closed form solution"... any ideas why? or any suggestions for another approach to solve these two equations??--Thanks

R2Magnitude=(sqrt(sum(R2.^2)))^2;
variable=transpose(u2)*R2;
[rho2,rr]=solve('(A+(B/rr^3))-rho2=0','((rho2)^2)+(R2Magnitude)-((2*(rho2))*(variable))-(rr^2)=0');

 

[Aaron Barnard]

The error you are getting is likely because the equations are not linear, and MATLAB's "solve" command only works for linear equations. If you want to solve these equations, you could use one of the numerical solvers in MATLAB such as fsolve or fzero. These numerical solvers require an initial guess for the solution, so you will need to determine a reasonable starting point for your equations.

 

[tacman]

It is difficult to give a specific answer without seeing the code and understanding the context of the problem. However, here are a few suggestions that may help:

1. Check your equations and make sure they are correct. Double-check the syntax and variables used in the "solve" command.

2. Make sure you have defined all the necessary variables and parameters before running the code.

3. Consider using a different solver method, such as "fsolve", which may be more suitable for your equations.

4. Break down the problem into smaller steps and check the results at each step. This can help identify where the error may be occurring.

5. Consider using a different approach altogether, such as using a different software or algorithm for orbit determination.

In general, it is always a good idea to carefully review your code and make sure all the inputs and outputs are correct before running it. Good luck with your project!