MATLAB牛顿拉夫逊法求潮流计算

程序如下：

clc

clear all

n=5;%节点数

%输入原始数据，求节点导纳矩阵

y(1,2)=1/(0.06+0.18i);

y(1,3)=1/(0.06+0.18i);

y(1,4)=1/(0.04+0.12i);

y(1,5)=1/(0.02+0.06i);

y(2,3)=1/(0.01+0.03i);

y(2,5)=1/(0.08+0.24i);

y(3,4)=1/(0.08+0.24i);

y(4,5)=0;

for i=1:5

for j=i:5

y(j,i)=y(i,j);

end

end

Y=0;

%求互导纳

for i=1:5

for j=i:5

if i~=j

Y(i,j)=-y(i,j);

end

end

end