function main()

clear;clc

z_ini=[1j*184.1 1j*67.6 1j*161.4];%初值

options=optimoptions(@lsqnonlin,'Algorithm','levenberg-marquardt');

[z,resnorm]=lsqnonlin(@eq17, z_ini,[],[],options);

fprintf('zxx=%gj   zxy=%fj   zyy=%fj   Resnorm=%g\n\n',imag(z(1)),imag(z(2)),imag(z(3)),resnorm)

F=eq17(z_ini)

function F=eq17(z0)

kk=[0.477 0.720 0.987 0.936 0.647];%kk就是原方程的x，tk是原方程的θ

tk=[0 pi/6 pi/3 pi/2 2*pi/3];

zxx=z0(1);zxy=z0(2);zyy=z0(3);

F(1)=(-1i*(377^2-zxy^2+zxx*zyy)+(-(377^2-zxy^2+zxx*zyy)^2+4*377^2*(zyy*(cos(tk(1)))^2-zxy*sin(2*tk(1))+zxx*(sin(tk(1)))^2)*(zxx*(cos(tk(1)))^2+zxy*sin(2*tk(1))+zyy*(sin(tk(1)))^2))^(1/2))*(2*377*(zyy*(cos(tk(1)))^2-zxy*sin(2*tk(1))+zxx*(sin(tk(1)))^2))^(-1)-kk(1);

F(2)=(-1i*(377^2-zxy^2+zxx*zyy)+(-(377^2-zxy^2+zxx*zyy)^2+4*377^2*(zyy*(cos(tk(2)))^2-zxy*sin(2*tk(2))+zxx*(sin(tk(2)))^2)*(zxx*(cos(tk(2)))^2+zxy*sin(2*tk(2))+zyy*(sin(tk(2)))^2))^(1/2))*(2*377*(zyy*(cos(tk(2)))^2-zxy*sin(2*tk(2))+zxx*(sin(tk(2)))^2))^(-1)-kk(2);

F(3)=(-1i*(377^2-zxy^2+zxx*zyy)+(-(377^2-zxy^2+zxx*zyy)^2+4*377^2*(zyy*(cos(tk(3)))^2-zxy*sin(2*tk(3))+zxx*(sin(tk(3)))^2)*(zxx*(cos(tk(3)))^2+zxy*sin(2*tk(3))+zyy*(sin(tk(3)))^2))^(1/2))*(2*377*(zyy*(cos(tk(3)))^2-zxy*sin(2*tk(3))+zxx*(sin(tk(3)))^2))^(-1)-kk(3);

F(4)=(-1i*(377^2-zxy^2+zxx*zyy)+(-(377^2-zxy^2+zxx*zyy)^2+4*377^2*(zyy*(cos(tk(4)))^2-zxy*sin(2*tk(4))+zxx*(sin(tk(4)))^2)*(zxx*(cos(tk(4)))^2+zxy*sin(2*tk(4))+zyy*(sin(tk(4)))^2))^(1/2))*(2*377*(zyy*(cos(tk(4)))^2-zxy*sin(2*tk(4))+zxx*(sin(tk(4)))^2))^(-1)-kk(4);

F(5)=(-1i*(377^2-zxy^2+zxx*zyy)+(-(377^2-zxy^2+zxx*zyy)^2+4*377^2*(zyy*(cos(tk(5)))^2-zxy*sin(2*tk(5))+zxx*(sin(tk(5)))^2)*(zxx*(cos(tk(5)))^2+zxy*sin(2*tk(5))+zyy*(sin(tk(5)))^2))^(1/2))*(2*377*(zyy*(cos(tk(5)))^2-zxy*sin(2*tk(5))+zxx*(sin(tk(5)))^2))^(-1)-kk(5);