Matlab Vs Mathematica using numeric solver

Question

0 个投票

eq.m

Hi all I am using Mathematica tool to solve polynomial equations. But now I want to solve the equations with matlab. I have a system contains 5 equations which are defined in the attached file. With Mathematica I used Nsolve to get the following solutions:

    {{q1 -> -0.0458716, q2 -> -0.710509, q3 -> -0.0525018, 
  q4 -> -0.700226, lamda -> -2579.73}, {q1 -> -0.0525018, 
  q2 -> 0.700226, q3 -> 0.0458716, q4 -> -0.710509, 
  lamda -> -2579.73}, {q1 -> 0.0525018, q2 -> -0.700226, 
  q3 -> -0.0458716, q4 -> 0.710509, 
  lamda -> -2579.73}, {q1 -> -0.837359, q2 -> 0.0330304, 
  q3 -> -0.542385, q4 -> 0.0596414, 
  lamda -> -2252.12}, {q1 -> 0.837359, q2 -> -0.0330304, 
  q3 -> 0.542385, q4 -> -0.0596414, 
  lamda -> -2252.12}, {q1 -> 0.542385, q2 -> 0.0596414, 
  q3 -> -0.837359, q4 -> -0.0330304, 
  lamda -> -2252.12}, {q1 -> -0.542385, q2 -> -0.0596414, 
  q3 -> 0.837359, q4 -> 0.0330304, 
  lamda -> -2252.12}, {q1 -> 0.208883, q2 -> -0.0654466, 
  q3 -> 0.975643, q4 -> -0.0143001, 
  lamda -> -270.132}, {q1 -> 0.975643, q2 -> 0.0143001, 
  q3 -> -0.208883, q4 -> -0.0654466, 
  lamda -> -270.132}, {q1 -> -0.208883, q2 -> 0.0654466, 
  q3 -> -0.975643, q4 -> 0.0143001, 
  lamda -> -270.132}, {q1 -> -0.975643, q2 -> -0.0143001, 
  q3 -> 0.208883, q4 -> 0.0654466, 
  lamda -> -270.132}, {q1 -> -0.208883, q2 -> 0.0654466, 
  q3 -> -0.975643, q4 -> 0.0143001, 
  lamda -> -270.132}, {q1 -> -0.590465, q2 -> 0.508565, 
  q3 -> -0.431456, q4 -> -0.454487, 
  lamda -> -2.15417}, {q1 -> 0.590465, q2 -> -0.508565, 
  q3 -> 0.431456, q4 -> 0.454487, 
  lamda -> -2.15417}, {q1 -> -0.431456, q2 -> 0.454487, 
  q3 -> 0.590465, q4 -> 0.508565, lamda -> -2.15417}, {q1 -> 0.431456,
   q2 -> -0.454487, q3 -> -0.590465, q4 -> -0.508565, 
  lamda -> -2.15417}, {q1 -> 0.495178, q2 -> 0.518363, 
  q3 -> -0.518277, q4 -> 0.466356, 
  lamda -> -3.76404}, {q1 -> -0.495178, q2 -> -0.518363, 
  q3 -> 0.518277, q4 -> -0.466356, 
  lamda -> -3.76404}, {q1 -> 0.518277, q2 -> 0.466356, q3 -> 0.495178,
   q4 -> -0.518363, lamda -> -3.76404}, {q1 -> -0.518277, 
  q2 -> -0.466356, q3 -> -0.495178, q4 -> 0.518363, 
  lamda -> -3.76404}, {q1 -> 0.00101825, q2 -> -0.997746, 
  q3 -> -0.0666776, q4 -> 0.00750728, 
  lamda -> -0.263298}, {q1 -> -0.00101825, q2 -> 0.997746, 
  q3 -> 0.0666776, q4 -> -0.00750728, 
  lamda -> -0.263298}, {q1 -> -0.0666776, q2 -> -0.00750728, 
  q3 -> -0.00101825, q4 -> -0.997746, 
  lamda -> -0.263298}, {q1 -> 0.0666776, q2 -> 0.00750728, 
  q3 -> 0.00101825, q4 -> 0.997746, 
  lamda -> -0.263298}, {q1 -> 0.168571, q2 -> 0.985018, 
  q3 -> 0.0312651, q4 -> -0.018617, 
  lamda -> -0.00357358}, {q1 -> -0.168571, q2 -> -0.985018, 
  q3 -> -0.0312651, q4 -> 0.018617, 
  lamda -> -0.00357358}, {q1 -> -0.0312651, q2 -> -0.018617, 
  q3 -> 0.168571, q4 -> -0.985018, 
  lamda -> -0.00357358}, {q1 -> 0.0312651, q2 -> 0.018617, 
  q3 -> -0.168571, q4 -> 0.985018, 
  lamda -> -0.00357358}, {q1 -> -0.171986, q2 -> 0.982162, 
  q3 -> 0.0760095, q4 -> 0.000512586, 
  lamda -> -0.00391533}, {q1 -> 0.171986, q2 -> -0.982162, 
  q3 -> -0.0760095, q4 -> -0.000512586, 
  lamda -> -0.00391533}, {q1 -> -0.0760095, q2 -> 0.000512586, 
  q3 -> -0.171986, q4 -> -0.982162, 
  lamda -> -0.00391533}, {q1 -> 0.0760095, q2 -> -0.000512586, 
  q3 -> 0.171986, q4 -> 0.982162, lamda -> -0.00391533}}

Now, With matlab I am trying with vpasolve but I got just one solution and it takes a lot of time.

q1 = sym('q1','real');
q2 = sym('q2','real');
q3 = sym('q3','real');
q4 = sym('q4','real');
lamda = sym('lamda','real');
S = vpasolve([ eq1==0 , eq2==0 , eq3==0 , eq4==0 , eq5==0],[q1, q2, q3, q4, lamda])