基于牛顿拉夫逊法潮流计算的matlab实验报告(6)

-4.000000000099142 -0.000018005879814 -0.000000000557833

ans =

第6次迭代结果：

ans =

1.000000000000000 0.000000000028884 0 5.000000000000000 0.000000000090221 0 -4.000000000000000 -0.000000000099142 0

>> x =

1 5 -4

function y=F(x)

y(1)=x(1)-5*x(2)^2+7*x(3)^2+12; y(2)=3*x(1)*x(2)+x(1)*x(3)-11*x(1); y(3)=2*x(2)*x(3)+40*x(1); y=[y(1);y(2);y(3)]; end

function J=Fd(x)

J=[1,-10*x(2),14*x(3); 3*x(2)+x(3)-11, 3*x(1),x(1); 40,2*x(3),2*x(2)]; end