电力系统稳态实验报告(9)

120. deta_s=zeros(5,5); 121. p=zeros(5,1); 122. q=zeros(5,1); 123. for i=1:5 124. for j=1:5

125. p(i)=p(i)+e(i)*(G(i,j)*e(j)-B(i,j)*f(j))+f(i)*(G(i,j)*f(j)+B(i,j)*e(j)); 126. end 127. end 128. for i=1:5 129. for j=1:5

130. q(i)=q(i)+f(i)*(G(i,j)*e(j)-B(i,j)*f(j))-e(i)*(G(i,j)*f(j)+B(i,j)*e(j)); 131. end 132. end 133. for i=1:5 134. for j=1:5

135. deta_s(i,j)=s(i,j)+s(j,i); 136. end 137. end 138. disp(s); 139. disp(deta_s); 140. disp(p);

disp(q);

2.PQ解耦法

clear; clc;

yb=zeros(5,5);

yb(1,1)=(1.37874-6.26166i)/2;yb(1,2)=-0.62402+3.90015i;yb(1,3)=-0.75471+2.64150i; yb(2,2)=(1.45390-66.98082i)/2;yb(2,3)=-0.82987+3.11203i;yb(2,4)=63.49206i; yb(3,3)=(1.58459-35.73786i)/2;yb(3,5)=31.74603i; yb(4,4)=(-66.66667i)/2; yb(5,5)=(-33.33333i)/2; yb=yb+conj(yb'); k=1; kp=0; kq=0; eps1=10^-5; G=real(yb); B=imag(yb); B1=B(1:4,1:4); B2=B(1:3,1:3);

U=[1; 1 ;1 ;1.05 ;1.05]; theta=zeros(5,1); deta_theta=zeros(4,1);