【正文】 m=p+1。J(p,q)=X6。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 p=2*i1。 m=p+1。J(p,q)=X6。 % N=N0+1 即 N=2*n+1 擴展列 △ P、△ Q for k2=k1:N1 % 擴展列 △ P、 △ Q J(k,k2)=J(k,k2)./J(k,k)。%消去運算 end J(k3,k)=0。 end 17 else for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 %修改節(jié)點電壓實部 k1=k+1。 %不滿足要求的節(jié)點數(shù)加 1 end end ICT2(a)=IT2。)。)。 E(k)=e(k)+f(k)*j。 disp(E)。39。)。 disp(39。各節(jié)點的電壓相角 sida 為 (節(jié)點號從小到大排列 )： 39。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。)。)。39。)。 Siz(i)=Si(p,q)。 SSi(p,q)=Si(p,q)*SB。,39。 disp(ZF)。)。 for i=1:nl p=B1(i,1)。 else Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。 ZF=[39。,num2str(p),39。 %disp(SSj(q,p))。 end disp(39。q=B1(i,2)。 ZF=[39。,num2str(q),39。 %disp(DDS(i))。 end 運行結果 21 22 23 6 心得體會 這次課程設計，雖然這次課程設計沒有想象中的那么簡單，花費了很多的時間來不斷地實驗，雖說花費了大量的時間但是最后終于完成了預定的目 的，收獲頗豐，也很好的鍛煉了自己運用知識的能力。同時由于求解過程中用到求節(jié)點導鈉矩陣，求矩陣的逆等等，又使我對以前所學的知識有了一次很好的溫習。同時在設計過程中要求我們有足夠的耐心。利用牛頓 拉夫遜法、采用 Matlab 語言進行潮流計算，提高了準確性，使用方便，已得到成功應用，而它在我們以后的專業(yè)課程的學習中應用也是很廣泛的，掌握好了 MATLAB，這為以后的專業(yè)課程的學習將有很大的幫助。從中收獲知識 、 收獲閱歷 、 收獲成熟，在此過程中，通過查找大量資料，以及不懈的努力，不僅培養(yǎng)了獨立思考、動手操作的能力，在各種其它能力上也都有了提高。水利電力出版社 ,1990. [5]陶海英 .基于 MATLAB 的電力系統(tǒng)穩(wěn)態(tài)仿真分析 [J].江西電力職業(yè)技術學院學報 ,2021(3):8082. [6]徐敏 ,彭瑜 .MATLAB 在《電力系統(tǒng)分析》教學中的應用 [J].電力系統(tǒng)及其自動化學報 ,2021,(7):152155. 。只有不斷的測試自己，挑戰(zhàn)自己，才能擁有更多的成功和快樂！享受過程，而不是結果！認真對待每一件事，珍惜每一分一秒，學到最多的知識和方法，鍛煉自己的能力，這個是我課程設計中學到的 最重要的東西，也是以后都將受益。 這次課程設計，雖然沒有想象中的那么簡單，花費了很多的時間來不斷地實驗，但是最后終于完成了預定的目的，收獲頗豐，也很好的鍛煉了自己運用知識的能力。傳統(tǒng)的手工計算 工作量大，易出錯，潮流計算費時費力。尤其是在剛開始的時候，不知道該從哪里入手，更是意識到這一點的重要性。 此次課程設計首先讓我明白了要使電力系統(tǒng)運行的穩(wěn)定，必須經(jīng)過精密的設 24 計和計算。39。,num2str(DDS(i))]。,num2str(p),39。 disp(DS(i))。)。39。,num2str(SSj(q,p))]。,num2str(q),39。 end disp(Sj(q,p))。 if B1(i,6)==0 Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))conj(E(p)))*conj(1./( 20 B1(i,3)*B1(i,5))))。各條支路的末端功率 Sj 為 (順序同您輸入 B1 時一致 )： 39。 disp(39。)=39。S(39。 Siz(i)=Si(p,q)。q=B1(i,2)。 disp(39。 disp(SS)。 disp(39。 end disp(39。 disp(sida)。)。 VV=V*UB。 disp(39。 disp(EE)。各節(jié)點的實際電壓標幺值 E 為 (節(jié)點號從小到大排列 )： 39。 for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 18 disp(39。 end %用高斯消去法解 w=J*V disp(39。 %修改節(jié)點電壓虛部 end %修改節(jié)點電壓 for k=3:N0 DET=abs(J(k,N))。 end end end %====上面是用 線性變換方式將 Jacobi 矩陣化成單位矩陣 ===== for k=3:2:N01 L=(k+1)./2。 end %========================================== for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 % 對角元規(guī)格化 if k~=3 % 不是第三行 %============================================================ k4=k1。 end end end end end %========= 以上為求雅可比矩陣的各個元素 ===================== 16 for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點） k1=k+1。J(m,N)=DP。J(p,q)=X5。% dP/df X5=2*e(i)。 elseif j1==iamp。J(m,N)=DP。J(p,q)=X5。 % dP/df X5=0。 % PV 節(jié)點電壓誤差 for j1=1:n if j1~=isbamp。J(m,N)=DP。%擴展列 △ Q m=p+1。% dQ/df p=2*i1。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 J(p,q)=X4。m=p+1。 % X1=dP/de X4=dQ/df p=2*i1。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 %電壓模平方 %========= 以下針對非 PV 節(jié)點來求取功率差及 Jacobi 矩陣元素 ========= if B2(i,6)~=3 %非 PV 節(jié)點 DP=P(i)P1。% 節(jié) 點 功 率 P 計算eiΣ(Gij*ejBij*fj)+fiΣ(Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。D(i)=0。a=0。 ICT1=0。 %PV 節(jié)點電壓給定模值 end for i=1:n %給定各節(jié)點注入功率 S(i)=B2(i,1)B2(i,2)。B=imag(Y)。 %對角元 1 側 end %求導納矩陣 disp(39。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。 disp(B2) end