My non-linear equations are—
cos(a)-cos(b)+cos(c)=0.725cos(5a)-cos(5b)+cos(5c)=0.5cos(7a)-cos(7b)+cos(7c)=0.5MATLAB: How to solve three non-linear equations having 3 unknowns in MATLAB
non-linear equation
Related Solutions
fsolve is a numerical solver - use vpasolve instead:
syms A B C t0eq1 = A*sin(2*B*pi*0.2 + 2*B*pi*(0 - t0)^2)+C - sin(2*pi*0.2); eq2 = A*sin(2*B*pi*0.4 + 2*B*pi*(0 - t0)^2)+C - sin(2*pi*0.4);eq3 = A*sin(2*B*pi*0.5 + 2*B*pi*(0 - t0)^2)+C - sin(2*pi*0.5);eq4 = A*sin(2*B*pi*0.7 + 2*B*pi*(-t0)^2)+C - sin(2*pi*0.7);result = vpasolve([eq1,eq2,eq3,eq4]);A = result.AB = result.BC = result.Ct0 = result.t0I allowed myself to replace 3.1416 by pi
It looks like the expected solution can be obtained with some work. Maybe someone else can get there faster.
>> syms x real>> A = sym([1 -2*cos(x) 0; -sin(x)/cos(x)*sin(x) 1 cos(x); 1 -2*sin(x)/cos(x)*sin(x) 2]) A = [ 1, -2*cos(x), 0][ -sin(x)^2/cos(x), 1, cos(x)][ 1, -(2*sin(x)^2)/cos(x), 2]>> [V,D]=eig(A) V = [ ((cos(x)^4 + 3*cos(x)^2*sin(x)^2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3) + 4/3))/(cos(x)^4 - 2*sin(x)^6) - (2*(cos(x)^4 + cos(x)^2*sin(x)^4 + cos(x)^2*sin(x)^2))/(cos(x)^4 - 2*sin(x)^6) - (cos(x)^2*sin(x)^2*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3) + 4/3)^2)/(cos(x)^4 - 2*sin(x)^6), - (2*(cos(x)^4 + cos(x)^2*sin(x)^4 + cos(x)^2*sin(x)^2))/(cos(x)^4 - 2*sin(x)^6) - ((cos(x)^4 + 3*cos(x)^2*sin(x)^2)*(1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 + (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 4/3))/(cos(x)^4 - 2*sin(x)^6) - (cos(x)^2*sin(x)^2*(1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 + (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 4/3)^2)/(cos(x)^4 - 2*sin(x)^6), - (2*(cos(x)^4 + cos(x)^2*sin(x)^4 + cos(x)^2*sin(x)^2))/(cos(x)^4 - 2*sin(x)^6) - ((cos(x)^4 + 3*cos(x)^2*sin(x)^2)*(1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 - (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 4/3))/(cos(x)^4 - 2*sin(x)^6) - (cos(x)^2*sin(x)^2*(1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 - (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 4/3)^2)/(cos(x)^4 - 2*sin(x)^6)][ ((3*cos(x)^3 + 2*cos(x)*sin(x)^4)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3) + 4/3))/(2*(cos(x)^4 - 2*sin(x)^6)) - (cos(x)^3 + 2*cos(x)*sin(x)^4 + cos(x)^3*sin(x)^2)/(cos(x)^4 - 2*sin(x)^6) - (cos(x)^3*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3) + 4/3)^2)/(2*(cos(x)^4 - 2*sin(x)^6)), - (cos(x)^3 + 2*cos(x)*sin(x)^4 + cos(x)^3*sin(x)^2)/(cos(x)^4 - 2*sin(x)^6) - ((3*cos(x)^3 + 2*cos(x)*sin(x)^4)*(1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 + (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 4/3))/(2*(cos(x)^4 - 2*sin(x)^6)) - (cos(x)^3*(1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 + (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 4/3)^2)/(2*(cos(x)^4 - 2*sin(x)^6)), - (cos(x)^3 + 2*cos(x)*sin(x)^4 + cos(x)^3*sin(x)^2)/(cos(x)^4 - 2*sin(x)^6) - ((3*cos(x)^3 + 2*cos(x)*sin(x)^4)*(1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 - (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 4/3))/(2*(cos(x)^4 - 2*sin(x)^6)) - (cos(x)^3*(1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 - (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 4/3)^2)/(2*(cos(x)^4 - 2*sin(x)^6))][ 1, 1, 1] D = [ 1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) + (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3) + 4/3, 0, 0][ 0, 4/3 - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 - (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)), 0][ 0, 0, 4/3 - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)/2 + (3^(1/2)*(1/(9*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3)) - (((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))*1i)/2 - 1/(18*(((cos(x)^2 + sin(x)^2 - 1/27)^2 - 1/729)^(1/2) - cos(x)^2 - sin(x)^2 + 1/27)^(1/3))] >> D=simplify(diag(D),300) % unclear why D(1) doesn't simplify to 2 - 1i
D = 2 - (2*(-1)^(1/3)*(2*3^(1/2) - 3)*(- 15*3^(1/2) - 26)^(1/3))/3 - 3^(1/2)/3 0 2 + 1i >> V2=V(:,2); % pick the eigenvector associated with lamda = 0
>> simplify(imag(V2),100) % verify the imaginary part is 0
ans = 0 0 0 >> V2=real(V2); % grab the real part
>> V2=simplify(V2,100) V2 = - 1/cos(2*x) - 1 -cos(x)/(2*cos(x)^2 - 1) 1 >> simplify(V2*(1-2*sin(x)^2),100) ans = -2*cos(x)^2 -cos(x) cos(2*x)Related Question
- How initial value of solution can be provided for system solving by fsolve
- How to solve this integration ???
- Solve command
- I am writing a code to solve 5 simultaneous equations with 5 unknowns. I am using the function vpasolve, however the code takes 50 minutes to run. Is there a quicker way of solving the equations
- Matrix differentiation
Best Answer