【正文】 end J(k3,k)=0。 end end end %====上面是用 線性變換方式將 Jacobi 矩陣化成單位矩陣 ===== for k=3:2:N01 L=(k+1)./2。 e(L)=e(L)J(k,N)。 %修改節(jié)點電壓實部 k1=k+1。 f(L)=f(L)J(k1,N)。 %修改節(jié)點電壓虛部 end %修改節(jié)點電壓 for k=3:N0 DET=abs(J(k,N))。 if DET=pr %電壓偏差量是否滿足要求 IT2=IT2+1。 %不滿足要求的節(jié)點數(shù)加 1 end end ICT2(a)=IT2。 ICT1=ICT1+1。 end %用高斯消去法解 w=J*V disp(39。迭代次數(shù)： 39。)。 disp(ICT1)。 18 disp(39。沒有達(dá)到精度要求的個數(shù)： 39。)。 disp(ICT2)。 for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 sida(k)=atan(f(k)./e(k))*180./pi。 E(k)=e(k)+f(k)*j。 end %=============== 計算各輸出量 =========================== disp(39。各節(jié)點的實際電壓標(biāo)幺值 E 為 (節(jié)點號從小到大排列 )： 39。)。 disp(E)。 EE=E*UB。 disp(EE)。 disp(39。39。)。 disp(39。各節(jié)點的電 壓大小 V 為 (節(jié)點號從小到大排列 )： 39。)。 disp(V)。 VV=V*UB。 disp(VV)。 disp(39。39。)。 disp(39。各節(jié)點的電壓相角 sida 為 (節(jié)點號從小到大排列 )： 39。)。 disp(sida)。 for p=1:n C(p)=0。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。 end S(p)=E(p)*C(p)。 end disp(39。各節(jié)點的功率 S 為 (節(jié)點號從小到大排列 )： 39。)。 disp(S)。 disp(39。~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~39。)。 19 SS=S*SB。 disp(SS)。 disp(39。39。)。 disp(39。各條支路的首端功率 Si 為 (順序同您輸入 B1 時一致 )： 39。)。 for i=1:nl p=B1(i,1)。q=B1(i,2)。 if B1(i,6)==0 Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。 Siz(i)=Si(p,q)。 else Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。 Siz(i)=Si(p,q)。 end disp(Si(p,q))。 SSi(p,q)=Si(p,q)*SB。 ZF=[39。S(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(SSi(p,q))]。 disp(ZF)。 %disp(SSi(p,q))。 disp(39。39。)。 end disp(39。各條支路的末端功率 Sj 為 (順序同您輸入 B1 時一致 )： 39。)。 for i=1:nl p=B1(i,1)。q=B1(i,2)。 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))))。 Sjy(i)=Sj(q,p)。 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))))。 Sjy(i)=Sj(q,p)。 end disp(Sj(q,p))。 SSj(q,p)=Sj(q,p)*SB。 ZF=[39。S(39。,num2str(q),39。,39。,num2str(p),39。)=39。,num2str(SSj(q,p))]。 disp(ZF)。 %disp(SSj(q,p))。 disp(39。39。)。 end disp(39。各條支路的功率損耗 DS 為 (順序同您輸入 B1 時一致 )： 39。)。 for i=1:nl p=B1(i,1)。q=B1(i,2)。 DS(i)=Si(p,q)+Sj(q,p)。 disp(DS(i))。 DDS(i)=DS(i)*SB。 ZF=[39。DS(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(DDS(i))]。 disp(ZF)。 %disp(DDS(i))。 disp(39。39。)。 end 運行結(jié)果 21 22 23 6 心得體會 這次課程設(shè)計，雖然這次課程設(shè)計沒有想象中的那么簡單，花費了很多的時間來不斷地實驗，雖說花費了大量的時間但是最后終于完成了預(yù)定的目 的，收獲頗豐，也很好的鍛煉了自己運用知識的能力。重新鞏固了有關(guān)穩(wěn)態(tài)部分的潮流計算的相關(guān)知識，培養(yǎng)了自己設(shè)計思維能力。 此次課程設(shè)計首先讓我明白了要使電力系統(tǒng)運行的穩(wěn)定，必須經(jīng)過精密的設(shè) 24 計和計算。在進(jìn)行課題設(shè)計的過程中，加深了我對潮流計算的認(rèn)識，尤其是對牛頓拉夫遜潮流計算的求解思路有了比較透徹的理解。同時由于求解過程中用到求節(jié)點導(dǎo)鈉矩陣，求矩陣的逆等等，又使我對以前所學(xué)的知識有了一次很好的溫習(xí)。同時也看到了研究性學(xué)習(xí)的效果，從研究中去學(xué)習(xí)，理論結(jié)合實際，將理論運用到實際，同時在實踐中發(fā)現(xiàn)問題，然后解決 問題，而且在此次課程設(shè)計中，讓我對課本上的知識有了更加深入的了解，特別是對系統(tǒng)設(shè)計部分的理論知識，在認(rèn)真閱讀課本的同時還要去圖書館查閱一些系統(tǒng)設(shè)計方面的資料，并且要考慮很多東西，更要學(xué)習(xí)一些其他以前未曾接觸的知識。尤其是在剛開始的時候，不知道該從哪里入手，更是意識到這一點的重要性。在這個設(shè)計中最重要的部分是程序圖的設(shè)計，只有在程序圖的基礎(chǔ)上，其他一切問題都好辦，最終才能得到正確的仿真結(jié)果。同時在設(shè)計過程中要求我們有足夠的耐心。 這段時間的不斷實驗，更是讓進(jìn)一步學(xué)習(xí)了 MATLAB 軟件的相關(guān)應(yīng)用。傳統(tǒng)的手工計算 工作量大，易出錯，潮流計算費時費力。目前，計算機已成為潮流計算的重要工具。利用牛頓 拉夫遜法、采用 Matlab 語言進(jìn)行潮流計算，提高了準(zhǔn)確性，使用方便，已得到成功應(yīng)用，而它在我們以后的專業(yè)課程的學(xué)習(xí)中應(yīng)用也是很廣泛的，掌握好了 MATLAB，這為以后的專業(yè)課程的學(xué)習(xí)將有很大的幫助。在以后的學(xué)習(xí)中，自己還要慢慢學(xué)習(xí)，慢慢探索，更好的掌握這個軟件的運用。 這次課程設(shè)計，雖然沒有想象中的那么簡單，花費了很多的時間來不斷地實驗，但是最后終于完成了預(yù)定的目的，收獲頗豐，也很好的鍛煉了自己運用知識的能力。重新鞏固 了有關(guān)穩(wěn)態(tài)部分的潮流計算的相關(guān)知識，培養(yǎng)了自己設(shè)計思維能力。從中收獲知識 、 收獲閱歷 、 收獲成熟，在此過程中，通過查找大量資料，以及不懈的努力，不僅培養(yǎng)了獨立思考、動手操作的能力，在各種其它能力上也都有了提高。努力的去彌補自己的缺點，發(fā)展自己的優(yōu)點，去充實自己，只有在了解了自己的長短之后，則會更加珍惜我所擁有的，更加努力的去完善它，增進(jìn)它。只有不斷的測試自己，挑戰(zhàn)自己，才能擁有更多的成功和快樂！享受過程，而不是結(jié)果！認(rèn)真對待每一件事，珍惜每一分一秒，學(xué)到最多的知識和方法，鍛煉自己的能力，這個是我課程設(shè)計中學(xué)到的 最重要的東西，也是以后都將受益。 參考文獻(xiàn) 25 [1] 李光琦 .電力系統(tǒng)穩(wěn)態(tài)分析 (第 2 版 ) [M].北京 :水利電力出版社 ,1995. [2] 吳國炎 .電力系統(tǒng)分析 [M].杭州 :浙江大學(xué)出版社 ,1993. [3]牛煥娜 ,井天軍 ,李漢成 ,等 .基于回路分析的含分布式電源配電網(wǎng)簡化潮流計算[J].電網(wǎng)技術(shù) ,2021(6):10331038. [4] 陸敏政 .電力系統(tǒng)分析 [M].北京 。水利電力出版社 ,1990. [5]陶海英 .基于 MATLAB 的電力系統(tǒng)穩(wěn)態(tài)仿真分析 [J].江西電力職業(yè)技術(shù)學(xué)院學(xué)報 ,2021(3):8082. [6]徐敏 ,彭瑜 .MATLAB 在《電力系統(tǒng)分析》教學(xué)中的應(yīng)用 [J].電力系統(tǒng)及其自動化學(xué)報 ,2021,(7):152155.