1、<p>　　信息工程學(xué)院課程設(shè)計(jì)任務(wù)書(shū)</p><p>　　年 月 日 </p><p><b>　　摘 要</b></p><p>　　潮流計(jì)算是電力網(wǎng)絡(luò)設(shè)計(jì)及運(yùn)行中最基本的計(jì)算，對(duì)電力網(wǎng)絡(luò)的各種設(shè)計(jì)方案及各種運(yùn)行方式進(jìn)行潮流計(jì)算，可以得到各種電網(wǎng)各節(jié)點(diǎn)的電壓，并求得網(wǎng)絡(luò)的潮流及網(wǎng)絡(luò)中各元

2、件的電力損耗，進(jìn)而求得電能損耗。本文以一個(gè)具體例子分析潮流計(jì)算的具體方法，并給出MATLAB仿真。</p><p>　　關(guān)鍵詞： 電力系統(tǒng) 潮流計(jì)算 Mmatlab仿真</p><p><b>　　目 錄</b></p><p><b>　　1 需求分析6</b></p><p>　　1

3、.1 所給電路圖及設(shè)計(jì)要求6</p><p>　　1.2 潮流計(jì)算的意義6</p><p><b>　　2 總體設(shè)計(jì)8</b></p><p>　　2.1 等值電路圖8</p><p>　　2.2 設(shè)計(jì)思路及方法8</p><p><b>　　3 詳細(xì)設(shè)計(jì)9</b>

4、;</p><p><b>　　3.1參數(shù)計(jì)算9</b></p><p>　　3.2潮流計(jì)算12</p><p>　　3.2 MATLAB仿真 14 </p><p><b>　　4 總結(jié)23</

5、b></p><p><b>　　參考文獻(xiàn)24</b></p><p><b>　　1 需求分析</b></p><p><b>　　1.1 二級(jí)標(biāo)題</b></p><p>　　1.2 潮流計(jì)算的意義</p><p>　　潮流計(jì)算是研究電力系統(tǒng)

6、穩(wěn)態(tài)運(yùn)行情況的一種基本電氣計(jì)算，常規(guī)潮流計(jì)算的任務(wù)是根據(jù)給定的運(yùn)行條件和網(wǎng)路結(jié)構(gòu)確定整個(gè)系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。潮流計(jì)算的結(jié)果是電力系統(tǒng)穩(wěn)定計(jì)算和故障分析的基礎(chǔ)。通過(guò)潮流計(jì)算可以判斷電網(wǎng)母線電壓、支路電流和功率是否越限，如果有越限，就應(yīng)采取措施，調(diào)整運(yùn)行方式。對(duì)于正在規(guī)劃的電力系統(tǒng)，通過(guò)潮流計(jì)算，可以為選擇電網(wǎng)供電方案和電氣設(shè)備提供依據(jù)。潮流計(jì)算還可以為繼電保護(hù)和自動(dòng)裝置定整計(jì)算、

7、電力系統(tǒng)故障計(jì)算和穩(wěn)定計(jì)算等提供原始數(shù)據(jù)。</p><p>　　具體表現(xiàn)在以下方面：</p><p>　　(1)在電網(wǎng)規(guī)劃階段,通過(guò)潮流計(jì)算,合理規(guī)劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無(wú)功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。</p><p>　　(2)在編制年運(yùn)行方式時(shí),在預(yù)計(jì)負(fù)荷增長(zhǎng)及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計(jì)

8、算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對(duì)規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。</p><p>　　(3)正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開(kāi)機(jī)方式,有功、無(wú)功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。</p><p>　　(4)預(yù)想事故、設(shè)備退出運(yùn)行對(duì)靜態(tài)安全的影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案。&l

9、t;/p><p>　　總結(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式時(shí)，采用離線潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線潮流計(jì)算。 </p><p

10、><b>　　2 總體設(shè)計(jì)</b></p><p><b>　　2.1 等值電路圖</b></p><p><b>　　設(shè)計(jì)思路及方法</b></p><p>　　潮流計(jì)算（load flow calculation）根據(jù)電力系統(tǒng)接線方式、參數(shù)和運(yùn)行條件計(jì)算電力系統(tǒng)穩(wěn)態(tài)運(yùn)行狀態(tài)下的電氣量。通常給

11、定的運(yùn)行條件有電源和負(fù)荷節(jié)點(diǎn)的功率、樞紐點(diǎn)電壓、平衡節(jié)點(diǎn)的電壓和相位角。待求的運(yùn)行狀態(tài)量包括各節(jié)點(diǎn)電壓及其相位角和各支路（元件）通過(guò)的電流（功率）、網(wǎng)絡(luò)的功率損耗等。</p><p>　　此題先按電網(wǎng)平均額定電壓VN計(jì)算功率損耗，再利用上面求得的功率分布，從電源點(diǎn)開(kāi)始，順著功率的傳送方向，依次計(jì)算各段線路的電壓降落，求出各節(jié)點(diǎn)電壓。對(duì)于首端電壓和功率已知的線路，只需順著電路的方向算各個(gè)變壓器及線路的電壓降落和功率

12、損耗，以此求得各節(jié)點(diǎn)的電壓和功率。</p><p><b>　　3 詳細(xì)設(shè)計(jì)</b></p><p><b>　　3.1參數(shù)計(jì)算</b></p><p><b>　　變壓器1：</b></p><p><b>　　=</b></p><

13、p><b>　　勵(lì)磁功率 ：</b></p><p><b>　　線路一</b></p><p><b>　　設(shè)，，</b></p><p><b>　　產(chǎn)生的充電功率：</b></p><p><b>　　同理可得</b>&l

14、t;/p><p>　　變壓器2（2臺(tái)變壓器并聯(lián)運(yùn)行）：</p><p><b>　　;</b></p><p><b>　　線路2：</b></p><p><b>　　產(chǎn)生的充電功率：</b></p><p><b>　　變壓器3：</b&

15、gt;</p><p><b>　　線路3（雙回路）：</b></p><p><b>　　產(chǎn)生的充電功率：</b></p><p>　　變壓器4（兩臺(tái)變壓器并聯(lián)）：</p><p><b>　　;</b></p><p><b>　　線路4：

16、</b></p><p><b>　　產(chǎn)生的充電功率：</b></p><p>　　變壓器5（兩臺(tái)三繞組變壓器并聯(lián)）：</p><p><b>　　線路5（雙回路）：</b></p><p><b>　　產(chǎn)生的充電功率：</b></p><p&g

17、t;　　變壓器6（兩臺(tái)變壓器并聯(lián)）：</p><p><b>　　;</b></p><p><b>　　發(fā)電機(jī)A：</b></p><p><b>　　發(fā)電機(jī)B：</b></p><p><b>　　發(fā)電機(jī)C：</b></p><p&

18、gt;<b>　　發(fā)電機(jī)D：</b></p><p><b>　　發(fā)電機(jī)E：</b></p><p><b>　　3.2潮流計(jì)算</b></p><p>　　變壓器1的功率損耗：</p><p><b>　　a點(diǎn)功率</b></p><

19、p><b>　　線路1的功率損耗：</b></p><p>　　變壓器T2的功率損耗：</p><p><b>　　所以b點(diǎn)功率為</b></p><p><b>　　c點(diǎn)功率：</b></p><p>　　計(jì)算線路I的電壓分布：</p><p>

20、<b>　　VA=115V;</b></p><p>　　變壓器2中的電壓降落：</p><p>　　線路I中的電壓降落：</p><p><b>　　線路II中：</b></p><p>　　線路III的功率損耗和電壓降落：</p><p><b>　　線路IV的

21、損耗：</b></p><p><b>　　點(diǎn)f的功率和電壓：</b></p><p>　　變壓器6的功率和損耗：</p><p><b>　　點(diǎn)g的功率及電壓：</b></p><p>　　線路5 的功率損耗及電壓降落：</p><p><b>　　點(diǎn)

22、h的功率：</b></p><p>　　變壓器5-2的功率損耗和電壓降落：</p><p><b>　　點(diǎn)m的功率和電壓：</b></p><p>　　變壓器5-3的功率損耗和電壓降落：</p><p><b>　　點(diǎn)n的功率和電壓：</b></p><p>　　

23、3.3MATLAB仿真（由陳志彪完成）</p><p>　　相關(guān)的原始數(shù)據(jù)輸入格式如下：</p><p>　　1、B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點(diǎn)編號(hào)。節(jié)點(diǎn)編號(hào)由小到大編寫(xiě)。</p><p>　　2、對(duì)于含有變壓器的支路，第一列為低壓側(cè)節(jié)點(diǎn)編號(hào)，第二列為高壓側(cè)節(jié)點(diǎn)編號(hào)，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理，第三列為支路的串列阻抗參數(shù)，第四列為支路的對(duì)地導(dǎo)納參數(shù)，

24、第五烈為含變壓器支路的變壓器的變比，第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，“0”為不含有變壓器。</p><p>　　3、B2為節(jié)點(diǎn)參數(shù)矩陣，其中第一列為節(jié)點(diǎn)注入發(fā)電功率參數(shù)；第二列為節(jié)點(diǎn)負(fù)荷功率參數(shù)；第三列為節(jié)點(diǎn)電壓參數(shù)；第六列為節(jié)點(diǎn)類型參數(shù)，其中“1”為平衡節(jié)點(diǎn)，“2”為PQ節(jié)點(diǎn)，“3”為PV節(jié)點(diǎn)參數(shù)。</p><p>　　4、X為節(jié)點(diǎn)號(hào)和對(duì)地參數(shù)矩陣。其中第

25、一列為節(jié)點(diǎn)編號(hào)，第二列為節(jié)點(diǎn)對(duì)地參數(shù)。</p><p><b>　　程序如下：</b></p><p><b>　　close all</b></p><p><b>　　clear all</b></p><p>　　n=input('請(qǐng)輸入節(jié)點(diǎn)數(shù):n=');

26、</p><p>　　n1=input('請(qǐng)輸入支路數(shù):n1=');</p><p>　　isb=input('請(qǐng)輸入平衡節(jié)點(diǎn)號(hào):isb=');</p><p>　　pr=input('請(qǐng)輸入誤差精度:pr=');</p><p>　　B1=input('請(qǐng)輸入支路參數(shù):B1='

27、);</p><p>　　B2=input('請(qǐng)輸入節(jié)點(diǎn)參數(shù):B2=');</p><p>　　X=input('節(jié)點(diǎn)號(hào)和對(duì)地參數(shù):X=');</p><p>　　Y=zeros(n);</p><p>　　Times=1; %置迭代次數(shù)為初始值</p><p&

28、gt;<b>　　%創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣</b></p><p>　　for i=1:n1</p><p>　　if B1(i,6)==0 %不含變壓器的支路</p><p>　　p=B1(i,1);</p><p>　　q=B1(i,2);</p><p>　　Y(p,q)=Y(p,q)

29、-1/B1(i,3);</p><p>　　Y(q,p)=Y(p,q);</p><p>　　Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);</p><p>　　Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);</p><p>　　else %含有變壓器的支路&l

30、t;/p><p>　　p=B1(i,1);</p><p>　　q=B1(i,2);</p><p>　　Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));</p><p>　　Y(q,p)=Y(p,q);</p><p>　　Y(p,p)=Y(p,p)+1/B1(i,3);</p><

31、;p>　　Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　Y</b></p><p>　　OrgS=zeros(2*n-2,1);

32、</p><p>　　DetaS=zeros(2*n-2,1); %將OrgS、DetaS初始化</p><p>　　%創(chuàng)建OrgS，用于存儲(chǔ)初始功率參數(shù)</p><p><b>　　h=0;</b></p><p><b>　　j=0;</b></p><p>　　fo

33、r i=1:n %對(duì)PQ節(jié)點(diǎn)的處理</p><p>　　if i~=isb&B2(i,6)==2 </p><p><b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p><p>　　OrgS(2*h-1,1)=OrgS(2

34、*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p>　　OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j

35、))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>　　end</b></p><p><b>　　end</b></p><

36、;p><b>　　end</b></p><p>　　for i=1:n %對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0</p><p>　　if i~=isb&B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b&g

37、t;　　for j=1:n</b></p><p>　　OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p

38、><p>　　OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>　　end<

39、/b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　OrgS</b></p><p>　　%創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓</p><p>　　PVU=zeros(n-

40、h-1,1);</p><p><b>　　t=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if B2(i,6)==3</p><p><b>　　t=t+1;</b></p><p>　　PVU(t,1)=

41、B2(i,3);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　PVU</b></p><p>　　%創(chuàng)建DetaS，用于存儲(chǔ)有功功率、無(wú)功功率和電壓幅值的不平衡量</p><p>&

42、lt;b>　　h=0;</b></p><p>　　for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理</p><p>　　if i~=isb&B2(i,6)==2</p><p><b>　　h=h+1;</b></p><p>　　DetaS(2*h-1,1)=real(B2(i,2

43、))-OrgS(2*h-1,1);</p><p>　　DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　t=0;</b></

44、p><p>　　for i=1:n %對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0</p><p>　　if i~=isb&B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b>　　t=t+1;</b></p><

45、;p>　　DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);</p><p>　　DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;</p><p><b>　　end</b></p><p>&

46、lt;b>　　end</b></p><p><b>　　DetaS</b></p><p>　　%創(chuàng)建I，用于存儲(chǔ)節(jié)點(diǎn)電流參數(shù)</p><p>　　i=zeros(n-1,1);</p><p><b>　　h=0;</b></p><p><b&g

47、t;　　for i=1:n</b></p><p><b>　　if i~=isb</b></p><p><b>　　h=h+1;</b></p><p>　　I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));</p><p

48、><b>　　end</b></p><p><b>　　end</b></p><p><b>　　I</b></p><p>　　%創(chuàng)建Jacbi(雅可比矩陣)</p><p>　　Jacbi=zeros(2*n-2);</p><p><

49、;b>　　h=0;</b></p><p><b>　　k=0;</b></p><p>　　for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理</p><p>　　if B2(i,6)==2</p><p><b>　　h=h+1;</b></p><p>

50、;<b>　　for j=1:n</b></p><p><b>　　if j~=isb</b></p><p><b>　　k=k+1;</b></p><p>　　if i==j %對(duì)角元素的處理</p><p>　　Jacbi(2*h-1,2*k-1)=-imag

51、(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>　　Jacbi(2*h,2*k-1)=-Jacbi(2*h

52、-1,2*k)+2*real(I(h,1));</p><p>　　Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));</p><p>　　else %非對(duì)角元素的處理</p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j)

53、)*imag(B2(i,3));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);</p><p>　　Jacbi(2*h,2*k)=Jacbi(2*h-1,2*

54、k-1);</p><p><b>　　end</b></p><p>　　if k==(n-1) %將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行</p><p><b>　　k=0;</b></p><p><b>　　end</b></p><p&

55、gt;<b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　k=0;</b></p>

56、<p>　　for i=1:n %對(duì)PV節(jié)點(diǎn)的處理</p><p>　　if B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p><p><b>　　if j~=isb</b&g

57、t;</p><p><b>　　k=k+1;</b></p><p>　　if i==j %對(duì)角元素的處理</p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p

58、>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>　　Jacbi(2*h,2*k-1)=2*imag(B2(i,3));</p><p>　　Jacbi(2*h,2*k)=2*real(B2(i,3));</p><p>

59、;　　else %非對(duì)角元素的處理</p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</

60、p><p>　　Jacbi(2*h,2*k-1)=0;</p><p>　　Jacbi(2*h,2*k)=0;</p><p><b>　　end</b></p><p>　　if k==(n-1) %將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行</p><p><b>　　k

61、=0;</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b&

62、gt;　　end</b></p><p><b>　　Jacbi</b></p><p>　　%求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量</p><p>　　DetaU=zeros(2*n-2,1);</p><p>　　DetaU=inv(Jacbi)*DetaS;</p><p>&l

63、t;b>　　DetaU</b></p><p><b>　　%修正節(jié)點(diǎn)電壓</b></p><p><b>　　j=0;</b></p><p>　　for i=1:n %對(duì)PQ節(jié)點(diǎn)處理</p><p>　　if B2(i,6)==2</p><p&

64、gt;<b>　　j=j+1;</b></p><p>　　B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　f

65、or i=1:n %對(duì)PV節(jié)點(diǎn)的處理</p><p>　　if B2(i,6)==3</p><p><b>　　j=j+1;</b></p><p>　　B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>　　end&

66、lt;/b></p><p><b>　　end</b></p><p><b>　　B2</b></p><p>　　%開(kāi)始循環(huán)**********************************************************************</p><p>　　w

67、hile abs(max(DetaU))>pr</p><p>　　OrgS=zeros(2*n-2,1); %!!!初始功率參數(shù)在迭代過(guò)程中是不累加的，所以在這里必須將其初始化為零矩陣</p><p><b>　　h=0;</b></p><p><b>　　j=0;</b></p>&

68、lt;p><b>　　for i=1:n</b></p><p>　　if i~=isb&B2(i,6)==2</p><p><b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p><p>　　OrgS(2*h-1,1)

69、=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p>　　OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(rea

70、l(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>　　end</b></p><p><b>　　end</b></p&g

71、t;<p><b>　　end</b></p><p><b>　　for i=1:n</b></p><p>　　if i~=isb&B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b>　　for j

72、=1:n</b></p><p>　　OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><

73、;p>　　OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>　　end</b><

74、;/p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　OrgS</b></p><p><b>　　%創(chuàng)建DetaS</b></p><p><b>　　h=0;&

75、lt;/b></p><p><b>　　for i=1:n</b></p><p>　　if i~=isb&B2(i,6)==2</p><p><b>　　h=h+1;</b></p><p>　　DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1)

76、;</p><p>　　DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　t=0;</b></p><p>

77、<b>　　for i=1:n</b></p><p>　　if i~=isb&B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b>　　t=t+1;</b></p><p>　　DetaS(2*h-1,1)=real(B2(i

78、,2))-OrgS(2*h-1,1);</p><p>　　DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;</p><p><b>　　end</b></p><p><b>　　end</b></p>

79、;<p><b>　　DetaS</b></p><p><b>　　%創(chuàng)建I</b></p><p>　　i=zeros(n-1,1);</p><p><b>　　h=0;</b></p><p><b>　　for i=1:n</b>&

80、lt;/p><p><b>　　if i~=isb</b></p><p><b>　　h=h+1;</b></p><p>　　I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));</p><p><b>　　end</b

81、></p><p><b>　　end</b></p><p><b>　　I</b></p><p><b>　　%創(chuàng)建Jacbi</b></p><p>　　Jacbi=zeros(2*n-2);</p><p><b>　　h=0;

82、</b></p><p><b>　　k=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if B2(i,6)==2</p><p><b>　　h=h+1;</b></p><p><b>

83、;　　for j=1:n</b></p><p><b>　　if j~=isb</b></p><p><b>　　k=k+1;</b></p><p><b>　　if i==j</b></p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(

84、i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>　　Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,

85、2*k)+2*real(I(h,1));</p><p>　　Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));</p><p><b>　　else</b></p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))

86、*imag(B2(i,3));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);</p><p>　　Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k

87、-1);</p><p><b>　　end</b></p><p>　　if k==(n-1)</p><p><b>　　k=0;</b></p><p><b>　　end</b></p><p><b>　　end</b>&

88、lt;/p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　k=0;</b></p><p><b>　　for i=1:

89、n</b></p><p>　　if B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p><p><b>　　if j~=isb</b></p><p>&l

90、t;b>　　k=k+1;</b></p><p><b>　　if i==j</b></p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>　　Jacbi(2*h-1,2

91、*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>　　Jacbi(2*h,2*k-1)=2*imag(B2(i,3));</p><p>　　Jacbi(2*h,2*k)=2*real(B2(i,3));</p><p><b>　　else<

92、;/b></p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>　

93、　Jacbi(2*h,2*k-1)=0;</p><p>　　Jacbi(2*h,2*k)=0;</p><p><b>　　end</b></p><p>　　if k==(n-1)</p><p><b>　　k=0;</b></p><p><b>　　end

94、</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>

95、　　Jacbi</b></p><p>　　DetaU=zeros(2*n-2,1);</p><p>　　DetaU=inv(Jacbi)*DetaS;</p><p><b>　　DetaU</b></p><p><b>　　%修正節(jié)點(diǎn)電壓</b></p><p

96、><b>　　j=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if B2(i,6)==2</p><p><b>　　j=j+1;</b></p><p>　　B2(i,3)=B2(i,3)+DetaU(2*j,1)+Deta

97、U(2*j-1,1)*sqrt(-1);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　for i=1:n</b></p><p>　　if B2(i,6)==3</p><p>&l

98、t;b>　　j=j+1;</b></p><p>　　B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>

99、;　　B2</b></p><p>　　Times=Times+1; %迭代次數(shù)加1</p><p><b>　　end</b></p><p><b>　　Times</b></p><p><b>　　4 總結(jié)</b></p><p

100、>　　透過(guò)這次電力系統(tǒng)分析課程設(shè)計(jì)，對(duì)自己以往的學(xué)習(xí)有了很深的認(rèn)識(shí)。電力系統(tǒng)分析專業(yè)課的系統(tǒng)學(xué)習(xí)已經(jīng)過(guò)去了三個(gè)多月，由于電力系統(tǒng)分析課程本身較難，在學(xué)習(xí)過(guò)程中沒(méi)有做到足夠的理解學(xué)習(xí)以及在以后的時(shí)間沒(méi)有做到充分的溫習(xí)，因此在這次設(shè)計(jì)作業(yè)的解題過(guò)程中感遇到了很大的阻力。</p><p>　　首先拿到圖就感到無(wú)從下手，已經(jīng)超出我們平時(shí)作業(yè)的水平，也超出了我們的能力范圍，以至好長(zhǎng)時(shí)間都沒(méi)有動(dòng)筆，只是每次拿出來(lái)看幾眼

101、。慢慢看出了些眉頭后，又用過(guò)請(qǐng)教同學(xué)和老師把部分不懂的細(xì)節(jié)弄清了。完成這個(gè)艱難的第一步后，對(duì)整個(gè)線路圖有了宏觀的把握，對(duì)解題方法也有了大致思路。第二步就應(yīng)該對(duì)電路圖做詳細(xì)分析，計(jì)算線路參數(shù)作出等值電路，這一步相對(duì)進(jìn)行的比較順利。接下來(lái)就是理清思路認(rèn)真做潮流計(jì)算，在做潮流計(jì)算的時(shí)候由于計(jì)算量比較大，線路參數(shù)比較多，所花的時(shí)間很多。完成整個(gè)潮流計(jì)算后，面臨這一個(gè)更艱巨的任務(wù)，就是手寫(xiě)稿轉(zhuǎn)化成電子稿。由于公式比較復(fù)雜，上下標(biāo)多，對(duì)公式編輯及畫(huà)

102、圖軟件應(yīng)用甚少，操作起來(lái)很費(fèi)時(shí)。總得來(lái)說(shuō)，還是基本上完成了任務(wù)，現(xiàn)將在此次設(shè)計(jì)的過(guò)程中遇到的問(wèn)題及收獲總結(jié)如下：</p><p><b>　　存在的問(wèn)題：</b></p><p>　　不能熟記電力系統(tǒng)課本中有關(guān)公式，應(yīng)用時(shí)需要查看課本。</p><p>　　對(duì)三繞組的變壓器應(yīng)用不深。</p><p>　　對(duì)潮流計(jì)算的幾種

103、模型及基本運(yùn)算了解不透，不能熟練應(yīng)用。</p><p><b>　　計(jì)算能力不強(qiáng)。</b></p><p>　　對(duì)公式編輯及畫(huà)圖軟件viso操作不熟。</p><p>　　由于在計(jì)算過(guò)程中忽略了變壓器的勵(lì)磁損耗，而且只做了一次迭代，沒(méi)有保證精度，與MATLAB仿真存在很一定的差異，有的甚至出現(xiàn)錯(cuò)誤的數(shù)據(jù)。</p><p>

104、;<b>　　收獲：</b></p><p>　　為了較好的完成此次作業(yè)，對(duì)主要的相關(guān)章節(jié)又進(jìn)行較全面的學(xué)習(xí)，尤其是加深了對(duì)有關(guān)公式的應(yīng)用。</p><p>　　提高了對(duì)有關(guān)軟件的操作能力。</p><p>　　由于此模板類似于畢業(yè)論文的模板，讓我們提前熟悉了畢業(yè)論文的要求和格式。這次的課程設(shè)計(jì)為明年的畢業(yè)論文打下了很好的基礎(chǔ)。</p&g

105、t;<p>　　同時(shí)，通過(guò)此次的課程設(shè)計(jì)，讓我認(rèn)識(shí)到一個(gè)很嚴(yán)重的問(wèn)題：我們的學(xué)習(xí)不扎實(shí)。考完之后或是經(jīng)過(guò)一段時(shí)間我們就對(duì)前面學(xué)過(guò)的知識(shí)沒(méi)有映像了，這主要是因?yàn)闆](méi)有理解性學(xué)習(xí)。所以，在以后的學(xué)習(xí)中，一定要學(xué)的扎實(shí)，學(xué)的牢靠，并且能學(xué)以致用，活學(xué)活用。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　1.何仰贊 溫増銀 《電力系統(tǒng)分析》第三

106、版</p><p>　　3.王守相 劉玉田 電力系統(tǒng)潮流計(jì)算研究現(xiàn)狀--《山東電力技術(shù)》1996年05期 </p><p>　　4.飛思科技產(chǎn)品研發(fā)中心.MATLAB基礎(chǔ)與提高.北京：電子工業(yè)出版社，2005.</p><p>　　5.Stephen J. Chapman. MATLAB Programming for Engineers Second Edi