我现在只能求解有限个数的方程组,在此我令N=9,其中 ,x(10)对应V,x(j)对应fai(j)
function y=fc(x,f1,f2,f3,f4,f5,f6,f7,f8,f9)
N=3*3;
y(1)=x(1)-atan(imag(x(10)/abs(x(10))*exp(1i*f1))/real(x(10)/abs(x(10))*exp(1i*f1)));
y(2)=x(2)-atan(imag(x(10)/abs(x(10))*exp(1i*f2))/real(x(10)/abs(x(10))*exp(1i*f2)));
y(3)=x(3)-atan(imag(x(10)/abs(x(10))*exp(1i*f3))/real(x(10)/abs(x(10))*exp(1i*f3)));
y(4)=x(4)-atan(imag(x(10)/abs(x(10))*exp(1i*f4))/real(x(10)/abs(x(10))*exp(1i*f4)));
y(5)=x(5)-atan(imag(x(10)/abs(x(10))*exp(1i*f5))/real(x(10)/abs(x(10))*exp(1i*f5)));
y(6)=x(6)-atan(imag(x(10)/abs(x(10))*exp(1i*f6))/real(x(10)/abs(x(10))*exp(1i*f6)));
y(7)=x(7)-atan(imag(x(10)/abs(x(10))*exp(1i*f7))/real(x(10)/abs(x(10))*exp(1i*f7)));
y(8)=x(8)-atan(imag(x(10)/abs(x(10))*exp(1i*f8))/real(x(10)/abs(x(10))*exp(1i*