电力系统三节点潮流算例（An example of power system three - point tidal current）.doc

电力系统三节点潮流算例（An example of power system three - point tidal current） 格式长 0057 =[1 4 0 6 0; 2 7 - 0 - 0062 - 5; 三9 - 0 - 6 - 0 0058。 0085—0088年4 5 0 01; 4 0017 0092 0079; - 5 0032 0161 0153; 六9 0039 .17 0039; 7 8 0008 0072 0072 5; 8 - 0011 0100 0011); y = 0(9); k = 1:9 m = 1(k); n =(k,2); (m m)= y(m m)+ 1 /((k,3)+ 1 *(4 k)+ 1 *(k,5); (n,n)= y(n,n)+ 1 /((k,3)+ 1 *(4 k)+ 1 *(k,5); (m,n)= y(m(1)(3)(k + 1 *(k),4); (n))= y(n m)1 /((k,3)+ 1 *(k),4); 世界尽头 U =[1,1025,1025,1,1,1,1,1,1); =[0,0,0,0,0,0,0,0,0); P =[1,1或清漆,再0 0 1 = 9 0 0 1 0); Q =[0,0,0,0,0,0 . 3,0,0,0); dP =[0,0,0,0,0,0,0,0,0); dQ =[0,0,0,0,0,0,0,0,0); 真实(y);B G = = y(上面); 分= 0(1 - 9); Qt = 0(9); H = 0 8(8); 8 N = 0(6); M = 0(8); L = 0(6,6); 哦,= 0(9); id = 0(9); 名= 0(14); I = 0(9); k = 0;精密= 1; 虽然精度 0000 - 0001 二是m =也 n = 1:9 第(n)= U(m)*(n *(G(m,n * cos(w)-(n))+ B(m,n)* sin(w)她(n))); 世界尽头 P(m)= P(m)-sum(pt); 世界尽头 m = 4:9 n = 1:9 qt(n)= U(m)*(n *(G(m,n * sin(w)- m,n(n)b()* cos(w)她(n))); 世界尽头 dQ(m)= Q(m)-sum qt(); 世界尽头 1米= 8 是1 n = 8 如果m = =( 其他的 H(m,n)= U(m + 1)*(n + 1)*(G(m + 1,n + 1)*(((m + 1)(n + 1)b(m + 1,n + 1)* cos(m + 1)(n + 1))); 世界尽头 世界尽头 世界尽头 1米= 8 n = 1:6 如果m = =( 其他的 N(N)= U(m + 1)*(N + 3)*(G(m + 1,N + 3)* cos(m + 1)(N + 3)+ B(m + 1,N + 3)*(((m + 1)(N + 3))); 世界尽头 世界尽头 世界尽头 m = 1:6 是1 n = 8 如果m = =( 其他的 米(M,n)= U(M + 3)*(n + 1)*(G(M + 3 n + 1)* cos(M + 3)她(n + 1)n + B(M + 3,+ 1)*(((M + 3)她(n + 1))); 世界尽头 世界尽头 结束 m = 1:6 为n = 1:6 如果m = = n 其他的 L(m,n)= - U(m + 3)* U(n + 3)*(G(m + 3,n + 3)*罪((m + 3)——(n + 3))- b(m + 3,n + 3)* cos((m + 3)——(n + 3))); 结束 结束 结束 %雅克比我不等于j m = 1:8 为n = 1:9 人工智能(n)= U(n)*(G(m + 1,n)* sin((m + 1)——(n))- b(m + 1,n)* cos((m + 1)——(n))); 结束 H(m m)= U(m + 1)*(m + 1)* B(m + 1,m + 1)+ U(m + 1)*(Ai)总和; 结束 m = 1:6 为n = 1:9 Bi(n)= U(n)*(G(m + 1,n)* cos((m + 1)——(n))+ B(m + 1,n)* sin((m + 1)——(n))); 结束 N(m m)= - U(m + 1)* U(m + 1)* G(m + 1,m + 1)- U(m + 1)*(Bi)总和; 结束 m = 1:6 为n = 1:9 Bi(n)= U(n)*(G(m + 3,n)* cos((m + 3)——(n))+ B(m + 3,n)* sin((m + 3)——(n))); 结束 米(M M)= U(M + 3)*(M + 3)* G