【正文】 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 擴(kuò)展列 △ P、△ Q for k2=k1:N1 % 擴(kuò)展列 △ P、 △ Q J(k,k2)=J(k,k2)./J(k,k)。%消去運(yùn)算 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é)點(diǎn)電壓實(shí)部 k1=k+1。 %不滿足要求的節(jié)點(diǎn)數(shù)加 1 end end ICT2(a)=IT2。)。)。 E(k)=e(k)+f(k)*j。 disp(E)。39。)。 disp(39。各節(jié)點(diǎn)的電壓相角 sida 為 (節(jié)點(diǎn)號從小到大排列 )： 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 運(yùn)行結(jié)果 21 22 23 6 心得體會(huì) 這次課程設(shè)計(jì)，雖然這次課程設(shè)計(jì)沒有想象中的那么簡單，花費(fèi)了很多的時(shí)間來不斷地實(shí)驗(yàn)，雖說花費(fèi)了大量的時(shí)間但是最后終于完成了預(yù)定的目 的，收獲頗豐，也很好的鍛煉了自己運(yùn)用知識的能力。同時(shí)由于求解過程中用到求節(jié)點(diǎn)導(dǎo)鈉矩陣，求矩陣的逆等等，又使我對以前所學(xué)的知識有了一次很好的溫習(xí)。同時(shí)在設(shè)計(jì)過程中要求我們有足夠的耐心。利用牛頓 拉夫遜法、采用 Matlab 語言進(jìn)行潮流計(jì)算，提高了準(zhǔn)確性，使用方便，已得到成功應(yīng)用，而它在我們以后的專業(yè)課程的學(xué)習(xí)中應(yīng)用也是很廣泛的，掌握好了 MATLAB，這為以后的專業(yè)課程的學(xué)習(xí)將有很大的幫助。從中收獲知識 、 收獲閱歷 、 收獲成熟，在此過程中，通過查找大量資料，以及不懈的努力，不僅培養(yǎng)了獨(dú)立思考、動(dòng)手操作的能力，在各種其它能力上也都有了提高。水利電力出版社 ,1990. [5]陶海英 .基于 MATLAB 的電力系統(tǒng)穩(wěn)態(tài)仿真分析 [J].江西電力職業(yè)技術(shù)學(xué)院學(xué)報(bào) ,2021(3):8082. [6]徐敏 ,彭瑜 .MATLAB 在《電力系統(tǒng)分析》教學(xué)中的應(yīng)用 [J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào) ,2021,(7):152155. 。只有不斷的測試自己，挑戰(zhàn)自己，才能擁有更多的成功和快樂！享受過程，而不是結(jié)果！認(rèn)真對待每一件事，珍惜每一分一秒，學(xué)到最多的知識和方法，鍛煉自己的能力，這個(gè)是我課程設(shè)計(jì)中學(xué)到的 最重要的東西，也是以后都將受益。 這次課程設(shè)計(jì)，雖然沒有想象中的那么簡單，花費(fèi)了很多的時(shí)間來不斷地實(shí)驗(yàn)，但是最后終于完成了預(yù)定的目的，收獲頗豐，也很好的鍛煉了自己運(yùn)用知識的能力。傳統(tǒng)的手工計(jì)算 工作量大，易出錯(cuò)，潮流計(jì)算費(fèi)時(shí)費(fèi)力。尤其是在剛開始的時(shí)候，不知道該從哪里入手，更是意識到這一點(diǎn)的重要性。 此次課程設(shè)計(jì)首先讓我明白了要使電力系統(tǒng)運(yùn)行的穩(wěn)定，必須經(jīng)過精密的設(shè) 24 計(jì)和計(jì)算。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 時(shí)一致 )： 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é)點(diǎn)的實(shí)際電壓標(biāo)幺值 E 為 (節(jié)點(diǎn)號從小到大排列 )： 39。 for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 18 disp(39。 end %用高斯消去法解 w=J*V disp(39。 %修改節(jié)點(diǎn)電壓虛部 end %修改節(jié)點(diǎn)電壓 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 %========= 以上為求雅可比矩陣的各個(gè)元素 ===================== 16 for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點(diǎn)） 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é)點(diǎn)電壓誤差 for j1=1:n if j1~=isbamp。J(m,N)=DP。%擴(kuò)展列 △ 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é)點(diǎn)來求取功率差及 Jacobi 矩陣元素 ========= if B2(i,6)~=3 %非 PV 節(jié)點(diǎn) DP=P(i)P1。% 節(jié) 點(diǎn) 功 率 P 計(jì)算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é)點(diǎn)電壓給定模值 end for i=1:n %給定各節(jié)點(diǎn)注入功率 S(i)=B2(i,1)B2(i,2)。B=imag(Y)。 %對角元 1 側(cè) end %求導(dǎo)納矩陣 disp(39。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。 disp(B2) end