【正文】 水利電力出版社 ,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é)到的 最重要的東西，也是以后都將受益。從中收獲知識 、 收獲閱歷 、 收獲成熟，在此過程中，通過查找大量資料，以及不懈的努力，不僅培養(yǎng)了獨(dú)立思考、動(dòng)手操作的能力，在各種其它能力上也都有了提高。 這次課程設(shè)計(jì)，雖然沒有想象中的那么簡單，花費(fèi)了很多的時(shí)間來不斷地實(shí)驗(yàn)，但是最后終于完成了預(yù)定的目的，收獲頗豐，也很好的鍛煉了自己運(yùn)用知識的能力。利用牛頓 拉夫遜法、采用 Matlab 語言進(jìn)行潮流計(jì)算，提高了準(zhǔn)確性，使用方便，已得到成功應(yīng)用，而它在我們以后的專業(yè)課程的學(xué)習(xí)中應(yīng)用也是很廣泛的，掌握好了 MATLAB，這為以后的專業(yè)課程的學(xué)習(xí)將有很大的幫助。傳統(tǒng)的手工計(jì)算 工作量大，易出錯(cuò)，潮流計(jì)算費(fèi)時(shí)費(fèi)力。同時(shí)在設(shè)計(jì)過程中要求我們有足夠的耐心。尤其是在剛開始的時(shí)候，不知道該從哪里入手，更是意識到這一點(diǎn)的重要性。同時(shí)由于求解過程中用到求節(jié)點(diǎn)導(dǎo)鈉矩陣，求矩陣的逆等等，又使我對以前所學(xué)的知識有了一次很好的溫習(xí)。 此次課程設(shè)計(jì)首先讓我明白了要使電力系統(tǒng)運(yùn)行的穩(wěn)定，必須經(jīng)過精密的設(shè) 24 計(jì)和計(jì)算。 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)用知識的能力。39。 %disp(DDS(i))。,num2str(DDS(i))]。,num2str(q),39。,num2str(p),39。 ZF=[39。 disp(DS(i))。q=B1(i,2)。)。 end disp(39。39。 %disp(SSj(q,p))。,num2str(SSj(q,p))]。,num2str(p),39。,num2str(q),39。 ZF=[39。 end disp(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))))。 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))))。 for i=1:nl p=B1(i,1)。各條支路的末端功率 Sj 為 (順序同您輸入 B1 時(shí)一致 )： 39。)。 disp(39。 disp(ZF)。)=39。,39。S(39。 SSi(p,q)=Si(p,q)*SB。 Siz(i)=Si(p,q)。 Siz(i)=Si(p,q)。q=B1(i,2)。)。 disp(39。39。 disp(SS)。)。 disp(39。)。 end disp(39。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。 disp(sida)。各節(jié)點(diǎn)的電壓相角 sida 為 (節(jié)點(diǎn)號從小到大排列 )： 39。)。 disp(39。 VV=V*UB。)。 disp(39。39。 disp(EE)。 disp(E)。各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值 E 為 (節(jié)點(diǎn)號從小到大排列 )： 39。 E(k)=e(k)+f(k)*j。 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)數(shù)加 1 end end ICT2(a)=IT2。 %修改節(jié)點(diǎn)電壓虛部 end %修改節(jié)點(diǎn)電壓 for k=3:N0 DET=abs(J(k,N))。 %修改節(jié)點(diǎn)電壓實(shí)部 k1=k+1。 end end end %====上面是用 線性變換方式將 Jacobi 矩陣化成單位矩陣 ===== for k=3:2:N01 L=(k+1)./2。 end 17 else for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end %========================================== for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。%消去運(yùn)算 end J(k3,k)=0。 % 對角元規(guī)格化 if k~=3 % 不是第三行 %============================================================ k4=k1。 % 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)。 end end end end end %========= 以上為求雅可比矩陣的各個(gè)元素 ===================== 16 for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點(diǎn)） k1=k+1。J(p,q)=X6。J(m,N)=DP。 m=p+1。J(p,q)=X5。 p=2*i1。% dP/df X5=2*e(i)。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 elseif j1==iamp。J(p,q)=X6。J(m,N)=DP。 m=p+1。J(p,q)=X5。 p=2*i1。 % dP/df X5=0。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 % PV 節(jié)點(diǎn)電壓誤差 for j1=1:n if j1~=isbamp。 end end 15 else %=============== 下面是針對 PV 節(jié)點(diǎn)來求取 Jacobi 矩陣的元素 =========== DP=P(i)P1。J(m,N)=DP。q=q+1。%擴(kuò)展列 △ Q m=p+1。J(p,q)=X3。% dQ/df p=2*i1。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 elseif j1==iamp。 J(p,q)=X4。J(m,N)=DP。m=p+1。 J(p,q)=X3。 % X1=dP/de X4=dQ/df p=2*i1。 % dP/df=dQ/de X3=X2。非對角元