【正文】 /B1(i,5))conj(E(p)))*conj(1./( 20 B1(i,3)*B1(i,5))))。 disp(39。S(39。q=B1(i,2)。 disp(SS)。 end disp(39。)。 disp(39。各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值 E 為 (節(jié)點(diǎn)號(hào)從小到大排列 )： 39。 18 disp(39。 %修改節(jié)點(diǎn)電壓虛部 end %修改節(jié)點(diǎn)電壓 for k=3:N0 DET=abs(J(k,N))。 end %========================================== for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end end end end end %========= 以上為求雅可比矩陣的各個(gè)元素 ===================== 16 for k=3:N0 % N0=2*n （從第三行開(kāi)始，第一、二行是平衡節(jié)點(diǎn)） k1=k+1。J(p,q)=X5。 elseif j1==iamp。J(p,q)=X5。 % PV 節(jié)點(diǎn)電壓誤差 for j1=1:n if j1~=isbamp。%擴(kuò)展列 △ Q m=p+1。對(duì)角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。m=p+1。非對(duì)角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。% 節(jié) 點(diǎn) 功 率 P 計(jì)算eiΣ(Gij*ejBij*fj)+fiΣ(Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。a=0。 %PV 節(jié)點(diǎn)電壓給定模值 end for i=1:n %給定各節(jié)點(diǎn)注入功率 S(i)=B2(i,1)B2(i,2)。 %對(duì)角元 1 側(cè) end %求導(dǎo)納矩陣 disp(39。 disp(B2) end % % % for i=1:nl %支路數(shù) if B1(i,6)==0 %左節(jié)點(diǎn)處于低壓側(cè) p=B1(i,1)。 disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB。 YB=SB./UB./UB。 SB=100。 0 60+ 220 0 0 2。 B2=[0 0 229 229 0 1。 1 6 + 1 0。經(jīng)驗(yàn)證其他三路均滿足，調(diào)壓成功。 ??????? （ 17） 位置 4 由 G1 提供的電壓為 121 239。39。0????????????? （ 7）② T3變壓器損耗 MV AjjXRUSS TTNT)()(110 )( 2 2233223 ???????? （ 8）③ 100Km 與 70Km 線路交點(diǎn) 4 末端功率損耗 M V AjSSSS LDT ]()[(39。附加的工具箱（單獨(dú)提供的專(zhuān)用 MATLAB 函數(shù)集）擴(kuò)展了 MATLAB 環(huán)境，以解決這些應(yīng)用領(lǐng)域內(nèi)特定類(lèi)型的問(wèn) 題。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗 5 環(huán)境 中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語(yǔ)言（如C、 Fortran）的編輯模式，代表了當(dāng)今國(guó)際科學(xué)計(jì)算軟件的先進(jìn)水平。這個(gè)方法把一個(gè)大系統(tǒng)分割為幾個(gè)小的地區(qū)系統(tǒng)，在計(jì)算機(jī)內(nèi)只需要存儲(chǔ)各個(gè)地區(qū)系統(tǒng)的阻抗矩陣及它們之間聯(lián)絡(luò)的阻抗，這樣不僅大幅度的節(jié)省了內(nèi)存容量，同時(shí)也提高了計(jì)算速度。實(shí)際電力系統(tǒng)的潮流技術(shù)那主要采用牛頓 — 拉夫遜法。潮流計(jì)算是根據(jù)給定的電網(wǎng)結(jié)構(gòu)、參數(shù)和發(fā)電機(jī)、負(fù)荷等元件的運(yùn)行條件，確定電力系統(tǒng)各部分穩(wěn)態(tài)運(yùn)行狀態(tài)參數(shù)的計(jì)算。 設(shè)計(jì)題目 電力系統(tǒng)潮流計(jì)算（牛頓 拉夫遜法） 設(shè)計(jì)內(nèi)容 （ 1） 根據(jù)電力系統(tǒng)網(wǎng)絡(luò)推導(dǎo)電力網(wǎng)絡(luò)數(shù)學(xué)模型，寫(xiě)出節(jié)點(diǎn)導(dǎo)納矩陣； （ 2） 賦予各節(jié)點(diǎn)電壓變量（直角坐標(biāo)系形式）初值 后，求解不平衡量； （ 3） 形成雅可比矩陣； （ 4） 求解的電壓變量達(dá)到所要求的精度時(shí)，再計(jì)算各支路功率分布、功率損耗和平衡節(jié)點(diǎn)功率； 3 （ 5） 上機(jī)編程調(diào)試； （ 6） 計(jì)算分析給定系統(tǒng)潮流分析并與手工計(jì)算結(jié)果做比較分析； （ 7） 書(shū)寫(xiě)課程設(shè)計(jì)報(bào)告。電力系統(tǒng)穩(wěn)態(tài)分析包括潮流計(jì)算（或潮流分析 )和靜態(tài)安全分析。 設(shè)計(jì)題目 ........................................................................................... 錯(cuò)誤 !未定義書(shū)簽。 電力系統(tǒng)簡(jiǎn)介 ................................................................................... 錯(cuò)誤 !未定義書(shū)簽。原則上講，靜態(tài)安全分析也可 U 用潮流計(jì)算來(lái)代替。電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)穩(wěn) 定計(jì)算和故障分析的基礎(chǔ)。 潮流計(jì)算的意義及其發(fā)展 電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，是對(duì)復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運(yùn)行狀態(tài)的計(jì)算。 在用數(shù)字見(jiàn)算計(jì)算機(jī)解電力系統(tǒng)潮流問(wèn)題的開(kāi)始階段，普遍采取以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的逐次代入法。在解決電力系統(tǒng)潮流計(jì)算問(wèn)題時(shí)，是以導(dǎo)納矩陣為基礎(chǔ)的 ，因此，只要我們能在迭代過(guò)程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓法潮流程序的效率。 MATLAB 可以進(jìn)行 矩陣 運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、連接其他編程語(yǔ)言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號(hào)處理與通訊、 圖像處理 、 信號(hào)檢測(cè) 、金融建模設(shè)計(jì)與分析等領(lǐng)域。5﹪） /， ΔPs=52kW， ΔP0=, Us﹪ =， I0﹪ =； 變壓器 T4： 2SFL116000/110，（ 110177。39。39。（經(jīng)驗(yàn)證其他 3 路均不滿足關(guān)系），所以需要降低 121KV 端的輸出電壓和提高 110KV 的輸額定值。 nl=8。 4 6 + 1 0。 0 0 220 0 0 2。e=zeros(1,n)。請(qǐng)輸入功率基準(zhǔn)值 :SB=39。 for i=1:nl 12 B1(i,3)=B1(i,3)*YB。 B2(i,4)=B2(i,4)./UB。q=B1(i,1)。 disp(Y) 13 % G=real(Y)。Q=imag(S)。 for i=1:n if i~=isb %非平衡節(jié)點(diǎn) C(i)=0。 14 V2=e(i)^2+f(i)^2。 % X2=dp/df X3=dQ/de X4=X1。q=q+1。 % dQ/de X4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。J(p,q)=X4。 % dP/de X2=B(i,j1)*e(i)G(i,j1)*f(i)。 J(m,q)=X1。% dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 J(m,q)=X1。 % 非對(duì)角元規(guī)格化 end J(k,k)=1。%消去運(yùn)算 end J(k3,k)=0。 ICT1=ICT1+1。 disp(ICT2)。 EE=E*UB。 disp(V)。)。 disp(S)。)。 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))))。,num2str(q),39。 end disp(39。 Sjy(i)=Sj(q,p)。)=39。各條支路的功率損耗 DS 為 (順序同您輸入 B1 時(shí)一致 )： 39。DS(39。 disp(39。同時(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í)還要去圖書(shū)館查閱一些系統(tǒng)設(shè)計(jì)方面的資料，并且要考慮很多東西，更要學(xué)習(xí)一些其他以前未曾接觸的知識(shí)。在以后的學(xué)習(xí)中，自己還要慢慢學(xué)習(xí)，慢慢探索，更好的掌握這個(gè)軟件的運(yùn)