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

证明：

因为 + = +

=

=

( )2= ( + ) = +

所以 + =Re + = = 附录： 1. 牛拉法

2. clear; 3. clc;

4. yb=zeros(5,5);

5. yb(1,1)=(1.37874-6.26166i)/2;yb(1,2)=-0.62402+3.90015i;yb(1,3)=-0.75471+2.64150i; 6. yb(2,2)=(1.45390-66.98082i)/2;yb(2,3)=-0.82987+3.11203i;yb(2,4)=63.49206i; 7. yb(3,3)=(1.58459-35.73786i)/2;yb(3,5)=31.74603i; 8. yb(4,4)=(-66.66667i)/2; 9. yb(5,5)=(-33.33333i)/2; 10. yb=yb+conj(yb'); 11. k=0;

12. eps1=10^-4; 13. jeps=1; 14. G=real(yb); 15. B=imag(yb);

16. e=[1; 1 ;1 ;1.05 ;1.05]; 17. f=zeros(5,1);

18. pis=[-1.6; -2 ;-3.7; 5]; 19. qis=[-0.8; -1; -1.3]; 20. deta_p=zeros(4,1); 21. deta_q=zeros(3,1); 22. deta=zeros(8,1); 23. deta_ef=zeros(8,1); 24. deta_e=zeros(4,1); 25. deta_f=zeros(4,1); 26. U=zeros(5,1); 27. while (jeps>eps1) 28. p=zeros(4,1); 29. q=zeros(3,1); 30. for i=1:4 31. for j=1:5

32. 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));

2

2