【正文】 0 0 220 0 0 2。 200 0 229 229 0 3]。e=zeros(1,n)。V=zeros(1,n)。S1=zeros(nl)。UB=220。請(qǐng)輸入功率基準(zhǔn)值 :SB=39。 %UB=input(39。)。 BB1=B1。 for i=1:nl 12 B1(i,3)=B1(i,3)*YB。 end disp(39。)。 B2(i,2)=B2(i,2)./SB。 B2(i,4)=B2(i,4)./UB。 end disp(39。)。q=B1(i,2)。q=B1(i,1)。 %非對(duì)角元 Y(q,p)=Y(p,q)。 %對(duì)角元 K 側(cè) Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。導(dǎo)納矩陣 Y=39。 disp(Y) 13 % G=real(Y)。 %分解出導(dǎo)納陣的實(shí)部和虛部 for i=1:n %給定各節(jié)點(diǎn)初始電壓的實(shí)部和虛部 e(i)=real(B2(i,3))。 V(i)=B2(i,4)。 %i 節(jié)點(diǎn)注入功率 SGSL B(i,i)=B(i,i)+B2(i,5)。Q=imag(S)。IT2=1。N=N0+1。 while IT2~=0 IT2=0。 for i=1:n if i~=isb %非平衡節(jié)點(diǎn) C(i)=0。 for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。%Σ(Gij*fj+Bij*ej) end P1=C(i)*e(i)+f(i)*D(i)。% 節(jié)點(diǎn)功率 Q 計(jì)算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej) %求 P39。 14 V2=e(i)^2+f(i)^2。 %節(jié)點(diǎn)有功功率差 DQ=Q(i)Q1。j1~=i %非平衡節(jié)點(diǎn) amp。 % dP/de=dQ/df X2=B(i,j1)*e(i)G(i,j1)*f(i)。 % X2=dp/df X3=dQ/de X4=X1。q=2*j11。J(p,N)=DQ。 J(m,q)=X1。q=q+1。J(m,q)=X2。j1~=isb %非平衡節(jié)點(diǎn) amp。% dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 % dQ/de X4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。q=2*j11。J(p,N)=DQ。 J(m,q)=X1。J(p,q)=X4。%擴(kuò)展列 △ P J(m,q)=X2。 % PV 節(jié)點(diǎn)有功誤差 DV=V(i)^2V2。j1~=i %非平衡節(jié)點(diǎn) amp。 % dP/de X2=B(i,j1)*e(i)G(i,j1)*f(i)。X6=0。q=2*j11。J(p,N)=DV。 J(m,q)=X1。q=q+1。 J(m,q)=X2。j1~=isb %非平衡節(jié)點(diǎn) amp。% dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 X6=2*f(i)。q=2*j11。J(p,N)=DV。 J(m,q)=X1。q=q+1。 J(m,q)=X2。N1=N。 % 非對(duì)角元規(guī)格化 end J(k,k)=1。 for k3=3:k4 % 用 k3 行從第三行開始到當(dāng)前行前的 k4行消去 for k2=k1:N1 % k3 行后各行下三角元素 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end if k==N0 break。%消去運(yùn)算 end J(k3,k)=0。%消去運(yùn)算 end J(k3,k)=0。 e(L)=e(L)J(k,N)。 f(L)=f(L)J(k1,N)。 if DET=pr %電壓偏差量是否滿足要求 IT2=IT2+1。 ICT1=ICT1+1。迭代次數(shù)： 39。 disp(ICT1)。沒(méi)有達(dá)到精度要求的個(gè)數(shù)： 39。 disp(ICT2)。 sida(k)=atan(f(k)./e(k))*180./pi。 end %=============== 計(jì)算各輸出量 =========================== disp(39。)。 EE=E*UB。 disp(39。)。各節(jié)點(diǎn)的電 壓大小 V 為 (節(jié)點(diǎn)號(hào)從小到大排列 )： 39。 disp(V)。 disp(VV)。39。 disp(39。)。 for p=1:n C(p)=0。 end S(p)=E(p)*C(p)。各節(jié)點(diǎn)的功率 S 為 (節(jié)點(diǎn)號(hào)從小到大排列 )： 39。 disp(S)。~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~39。 19 SS=S*SB。 disp(39。)。各條支路的首端功率 Si 為 (順序同您輸入 B1 時(shí)一致 )： 39。 for i=1:nl p=B1(i,1)。 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))))。 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))))。 end disp(Si(p,q))。 ZF=[39。,num2str(p),39。,num2str(q),39。,num2str(SSi(p,q))]。 %disp(SSi(p,q))。39。 end disp(39。)。q=B1(i,2)。 Sjy(i)=Sj(q,p)。 Sjy(i)=Sj(q,p)。 SSj(q,p)=Sj(q,p)*SB。S(39。,39。)=39。 disp(ZF)。 disp(39。)。各條支路的功率損耗 DS 為 (順序同您輸入 B1 時(shí)一致 )： 39。 for i=1:nl p=B1(i,1)。 DS(i)=Si(p,q)+Sj(q,p)。 DDS(i)=DS(i)*SB。DS(39。,39。)=39。 disp(ZF)。 disp(39。)。重新鞏固了有關(guān)穩(wěn)態(tài)部分的潮流計(jì)算的相關(guān)知識(shí)，培養(yǎng)了自己設(shè)計(jì)思維能力。在進(jìn)行課題設(shè)計(jì)的過(guò)程中，加深了我對(duì)潮流計(jì)算的認(rèn)識(shí)，尤其是對(duì)牛頓拉夫遜潮流計(jì)算的求解思路有了比較透徹的理解。同時(shí)也看到了研究性學(xué)習(xí)的效果，從研究中去學(xué)習(xí)，理論結(jié)合實(shí)際，將理論運(yùn)用到實(shí)際，同時(shí)在實(shí)踐中發(fā)現(xiàn)問(wèn)題，然后解決 問(wèn)題，而且在此次課程設(shè)計(jì)中，讓我對(duì)課本上的知識(shí)有了更加深入的了解，特別是對(duì)系統(tǒng)設(shè)計(jì)部分的理論知識(shí)，在認(rèn)真閱讀課本的同時(shí)還要去圖書館查閱一些系統(tǒng)設(shè)計(jì)方面的資料，并且要考慮很多東西，更要學(xué)習(xí)一些其他以前未曾接觸的知識(shí)。在這個(gè)設(shè)計(jì)中最重要的部分是程序圖的設(shè)計(jì)，只有在程序圖的基礎(chǔ)上，其他一切問(wèn)題都好辦，最終才能得到正確的仿真結(jié)果。 這段時(shí)間的不斷實(shí)驗(yàn)，更是讓進(jìn)一步學(xué)習(xí)了 MATLAB 軟件的相關(guān)應(yīng)用。目前，計(jì)算機(jī)已成為潮流計(jì)算的重要工具。在以后的學(xué)習(xí)中，自己還要慢慢學(xué)習(xí)，慢慢探索，更好的掌握這個(gè)軟件的運(yùn)用。重新鞏固 了有關(guān)穩(wěn)態(tài)部分的潮流計(jì)算的相關(guān)知識(shí)，培養(yǎng)了自己設(shè)計(jì)思維能力。努力的去彌補(bǔ)自己的缺點(diǎn)，發(fā)展自己的優(yōu)點(diǎn)，去充實(shí)自己，只有在了解了自己的長(zhǎng)短之后，則會(huì)更加珍惜我所擁有的，更加努力的去完善它，增進(jìn)它。 參考文獻(xiàn) 25 [1] 李光琦 .電力系統(tǒng)穩(wěn)態(tài)分析 (第 2 版 ) [M].北京 :水利電力出版社 ,1995. [2] 吳國(guó)炎 .電力系統(tǒng)分析 [M].杭州 :浙江大學(xué)出版社 ,1993. [3]牛煥娜 ,井天軍 ,李漢成 ,等 .基于回路分析的含分布式電源配電網(wǎng)簡(jiǎn)化潮流計(jì)算[J].電網(wǎng)技術(shù) ,2021(6):10331038. [4] 陸敏政 .電力系統(tǒng)分析 [M].北