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

76. jacobi(8,7)=-2*e(4); 77. jacobi(8,8)=-2*f(4); 78. for i=1:2:7

79. deta(i)=deta_p((i+1)/2); 80. end

81. for i=2:2:6

82. deta(i)=deta_q(i/2); 83. end

84. deta(8)=deta_UU;

85. deta_ef=-inv(jacobi)*deta; 86. for i=1:4

87. deta_e(i)=deta_ef(2*i-1); 88. deta_f(i)=deta_ef(2*i); 89. end 90. for i=1:4

91. e(i)=e(i)+deta_e(i); 92. f(i)=f(i)+deta_f(i); 93. end

94. jeps=max(max(abs(deta_ef)),max(abs(deta_p))); 95. jeps=max(jeps,max(abs(deta_q))); 96. k=k+1; 97. end 98. for i=1:5

99. U(i)=e(i)+sqrt(-1)*f(i); 100. end

101. disp('µçѹµÄ·ùÖµ£º'); 102. disp(abs(U)); 103. disp('µçѹ£º'); 104. disp(U);

105. theta=zeros(5,1); 106. for i=1:5

107. theta(i)=angle(U(i)); 108. end

109. disp(theta);

110. y0=[-4i;-8i;-4i;0;0]; 111. s=zeros(5,5); 112. for i=1:5 113. for j=1:5

114. s(i,j)=U(i)*(conj(U(i))*conj(y0(i))+(conj(U(i))-conj(U(j)))*conj(-yb(i,j))); 115. end 116. end 117. for i=1:5 118. s(i,i)=0;