matlab軟件潮流計算論文-資料下載頁
2025-06-19 14:22本頁面
【正文】 J(p,q)=X5。 % PV節(jié)點電壓誤差J(p,N)=DV。 m=p+1。 J(m,q)=X1。q=q+1。J(p,q)=X6。 % PV節(jié)點有功誤差J(m,N)=DP。 J(m,q)=X2。 end end end %(if B2(i,5)~=3 else) end %(if i~=isb) end %(for i=1:n)n個節(jié)點2n行(每節(jié)點兩個方程P和Q或U) JZ0=[39。形成的第(39。,num2str(a),39。)次Jacobi矩陣：39。]。 disp(JZ0)。disp(J)。 %=============================== 以上為形成完整的Jacobi矩陣 ================================ %====下面用高斯消去法對由Jacobi矩陣形成的修正方程進行求解(按列消去、回代) ========== for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點） for k1=k+1:N1 % 從k+1列的Jacobi元素到擴展列的△P、△Q 或 △U J(k,k1)=J(k,k1)./J(k,k)。% 用K行K列對角元素去除K行K列后的非對角元素進行規(guī)格化 end J(k,k)=1。 % 對角元規(guī)格化K行K列對角元素賦1 %================== 按列消去運算 ================================== for k2=k+1:N0 % 從k+1行到2*n最后一行 for k3=k+1:N1 %從k2+1列到擴展列消去k+1行后各行下三角元素 J(k2,k3)=J(k2,k3)J(k2,k)*J(k,k3)。%消去運算 end %用當前行K3列元素減去當前行K列元素乘以第k行K3列元素 J(k2,k)=0。 %當前行第k列元素已消為0 end end %JZ=[39。Jacobi矩陣前代第(39。,num2str(a),39。)次39。]。JZ1=[39。Jacobi矩陣回代第(39。,num2str(a),39。)次39。]。 %disp(JZ)。disp(J)。 %==================== 按列回代運算 ======================================= for k=N0:1:3 for k1=k1:1:3 J(k1,N1)=J(k1,N1)J(k1,k)*J(k,N1)。 J(k1,k)=0。 end end %for m=1:N0 % JJN1(m)=J(m,N1)。 %end %disp(JZ1)。disp(J)。disp(JJN1)。 %修改節(jié)點電壓 for k=3:2:N01 L=(k+1)./2。 e(L)=e(L)J(k,N1)。 %修改節(jié)點電壓實部 k1=k+1。 f(L)=f(L)J(k1,N1)。 %修改節(jié)點電壓虛部 end %============================== 結束一次迭代 ==============================end迭代次數： 4沒有達到精度要求的個數： 17 18 18 7 0各節(jié)點的電壓復數值E為(節(jié)點號從小到大排列)： +002 * Columns 1 through 5 + Columns 6 through 10 + + 各節(jié)點的電壓模值大小V為(節(jié)點號從小到大排列)： Columns 1 through 9 Column 10 各節(jié)點的電壓相角sida為(節(jié)點號從小到大排列)： Columns 1 through 9 0 Column 10 各節(jié)點的功率S為(節(jié)點號從小到大排列)： +002 * Columns 1 through 5 + + Columns 6 through 10 + 各條支路的首端功率Si為(順序同您輸入B1時一致)：S(1,3)=+S(1,5)=+S(1,7)=+S(2,9)=S(3,4)=+S(5,6)=+S(7,8)=+S(7,9)=+S(9,10)=+各條支路的末端功率Sj為(順序同您輸入B1時一致)：S(3,1)=S(5,1)=S(7,1)=S(9,2)=+S(4,3)=S(6,5)=S(8,7)=S(9,7)=S(10,9)=各條支路的功率損耗DS為(順序同您輸入B1時一致)：DS(1,3)=DS(1,5)=DS(1,7)=DS(2,9)=DS(3,4)=+DS(5,6)=+DS(7,8)=+DS(7,9)=DS(9,10)=+總網損S=39
外語相關推薦
鄂ICP備17016276號-1