MATLAB: How to solve the symbolic nonlinear equation

equationsymbolic

the system of equation has two nonlinear equations in two variables. The equations are:

F(1) = (x(1)-b2*cos(theta2)+a)^2+(x(2)-b2*sin(theta2))^2-c2^2;
F(2) = (x(1)-b1*cos(theta1)+a)^2+(x(2)-b1*sin(theta1))^2-c1^2;

Thanks!

Best Answer

syms  b2 a c2 c1 b1 theta2 theta1 x1 x2   % x1 is x(1) so is x2
F=zeros(2,1);
e1=F(1) == (x1-b2*cos(theta2)+a)^2+(x2-b2*sin(theta2))^2-c2^2;
e2=F(2) == (x1-b1*cos(theta1)+a)^2+(x2-b1*sin(theta1))^2-c1^2;
[x1,x2]=solve(e1,e2,x1,x2)

Related Solutions

MATLAB: Solve many independent equation in a run

syms x1 x2 x3 x4
sol = vpasolve( [ a1 + log(x1 - d1)*(c1 - x1) + x1*log(x1 - b1), a2 + log(x2 - d2)*(c2 - x2) + x2*log(x2 - b2), a3 + log(x3 - d3)*(c3 - x3) + x3*log(x3 - b3), a4 + log(x4 - d4)*(c4 - x4) + x4*log(x4 - b4)], [x1, x2, x3, x4])

MATLAB: How to find the jacobian of the system of nonlinear equations

Hi,.

try:

syms x1 x2 x3
h=4;
gma=18.4;
ka1=0.2;
kp1=8.76
ka2=0.2;
kp2=8.76;
sindel=0.437;
cosdel=0.9;
pa1=ka1*gma*(h*x1+0.5*x1^2);
pp1=kp1*gma*0.5*x1^2;
pa2=ka2*gma*(x1*x2*+0.5*x2^2);
pp2=kp2*(x2*(h+x1)+0.5*x2^2);
za1=(0.5*h*x1^2+(x1^3)/6)/(h*x1+0.5*x1^2);
zp2=(0.5*(h+x1)*x2^2+((x2^3)/3))/((h+x1)*x2+0.5*x2^2);
za2=(0.5*x1*x2^2+((x2^3)/3))+(x1*x2+0.5*x2^2);
e1=pp1*sindel-pa1*sindel-pp2*sindel-pa2*sindel;
e2=pp1*cosdel+pa2*cosdel-pa1*cosdel-pp2*cosdel-x3;
e3=pp1*cosdel*(x1/3)+pp2*cosdel*zp2-pa1*cosdel*za1-pa2*cosdel*za2-x3*(x1+(h/3));
g=[e1; e2; e3];
J=jacobian([e1, e2, e3], [x1, x2, x3])

results in:

J =
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                     (1075457*x1)/15625 - (95703*x2)/25000 - (10051*x2^3)/12500 - 20102/3125,                                                                                                                                                                                                                                                                                                                                              - (95703*x1)/25000 - (95703*x2)/25000 - (30153*x1*x2^2)/12500 - 95703/6250,          0]
[                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  (207*x2^3)/125 - (1971*x2)/250 + (88596*x1)/625 - 1656/125,                                                                                                                                                                                                                                                                                                                                                            (621*x1*x2^2)/125 - (1971*x2)/250 - (1971*x1)/250 - 3942/125,         -1]
[ (44298*x1^2)/625 - x3 - (207*x2^3*(x1*x2 + (x1*x2^2)/2 + x2^2/2 + x2^3/3))/125 - (1656*x1)/125 + (x2^2*((1971*x2*(x1 + 4))/250 + (1971*x2^2)/500))/(2*(x2^2/2 + (x1 + 4)*x2)) - (((414*x1)/125 + 1656/125)*(x1^3/6 + 2*x1^2))/(x1^2/2 + 4*x1) + (1971*x2*(x2^3/3 + (x1/2 + 2)*x2^2))/(250*(x2^2/2 + (x1 + 4)*x2)) - (207*x1*x2^3*(x2^2/2 + x2))/125 + (((207*x1^2)/125 + (1656*x1)/125)*(x1^3/6 + 2*x1^2)*(x1 + 4))/(x1^2/2 + 4*x1)^2 - (x2*(x2^3/3 + (x1/2 + 2)*x2^2)*((1971*x2*(x1 + 4))/250 + (1971*x2^2)/500))/(x2^2/2 + (x1 + 4)*x2)^2, ((x2^3/3 + (x1/2 + 2)*x2^2)*((1971*x1)/250 + (1971*x2)/250 + 3942/125))/(x2^2/2 + (x1 + 4)*x2) - (207*x1*x2^3*(x1 + x2 + x1*x2 + x2^2))/125 - (621*x1*x2^2*(x1*x2 + (x1*x2^2)/2 + x2^2/2 + x2^3/3))/125 + ((2*x2*(x1/2 + 2) + x2^2)*((1971*x2*(x1 + 4))/250 + (1971*x2^2)/500))/(x2^2/2 + (x1 + 4)*x2) - ((x2^3/3 + (x1/2 + 2)*x2^2)*((1971*x2*(x1 + 4))/250 + (1971*x2^2)/500)*(x1 + x2 + 4))/(x2^2/2 + (x1 + 4)*x2)^2, - x1 - 4/3]

Best regards

Stephan

Best Answer

Related Solutions

MATLAB: Solve many independent equation in a run

MATLAB: How to find the jacobian of the system of nonlinear equations

Related Question