%形成节点导纳矩阵
B1=[1 5 0.02 0.06;
1 2 0.08 0.24;
3 5 0.06 0.18;
4 5 0.02 0.18;
2 5 0.08 0.12;
3 4 0.01 0.03;
2 4 0.03 0.24;];
y=zeros(5);Y=zeros(5);
for i=1:7
p=B1(i,1);
q=B1(i,2);
y(p,q)=1.00/(B1(i,3)+B1(i,4)*sqrt(-1));
end
for b=1:5
for d=(b+1):5
y(d,b)=y(b,d);
end
end
m=5;n=5;
for i=1:n
for j=1:n
if(i==j)
for m=1:n
Y(i,i)=Y(i,i)+y(i,m);
end
else
Y(i,j)=-y(i,j);
end
end
end
%赋初始值
e(1)=1.00;e(3)=1.00;e(4)=1.00;e(2)=1.00;
f(1)=0;f(3)=0;f(4)=0;f(2)=0;
G=real(Y);B=imag(Y);
S(2)=(-67.5-41.83i)/100.00;S(1)=(97.5+60i)/100.00;S(3)=(-90.0-
55.77i)/100.00;S(4)=(-60.0-37.18i)/100.00;
P=real(S);Q=imag(S);k=0;p=1;N1=4;
%设置精度和最大循环次数
while p>0.00001&k<100
e(5)=1.05;
f(5)=0;
for m=1:N1
for n=1:N1+1 Pt(n)=(e(m)*(G(m,n)*e(n)-B(m,n)*f(n))
资源评论
