【正文】 分之一取負(fù)值，若有變壓器，Yij為線路阻抗乘以KT后分之一再取負(fù)值。②求導(dǎo)納矩陣Y中對角元元素Yii，無變壓器時Yii為Yij加上線路對地電導(dǎo)的一半乘j，有變壓器時，對角元元素就與所輸入的折算到哪一側(cè)有關(guān)，如果支路起始端處于高壓側(cè)，支路起始節(jié)點(diǎn)的自導(dǎo)納中要加上變壓器等值導(dǎo)納模型的對地支路的(1KT)/KT^2倍，支路終止節(jié)點(diǎn)的自導(dǎo)納要加上變壓器等值導(dǎo)納模型的對地支路的(KT1)/KT倍，如果支路起始端處于低壓側(cè)，情況正好相反支路起始節(jié)點(diǎn)的自導(dǎo)納中要加上變壓器等值導(dǎo)納模型的對地支路的(KT1)/KT倍，支路終止節(jié)點(diǎn)的自導(dǎo)納要加上變壓器等值導(dǎo)納模型的對地支路的(1KT)/KT^2倍。 導(dǎo)納矩陣Y為 雅克比矩陣J（k=0） 算例運(yùn)行結(jié)果輸出程序段本程序在結(jié)束時，會通過fopen這個函數(shù)，在MATLAB運(yùn)行程序所儲存的文件夾內(nèi)生成一個TXT文本文檔輸出滿足精度迭代后的節(jié)點(diǎn)電壓、相角、功率分布等運(yùn)行結(jié)果，清晰直觀。以下是程序運(yùn)行結(jié)果輸出程序段：fid=fopen(39。39。,39。wt39。)。fprintf(fid,39。******************************電力系統(tǒng)08 李進(jìn)才 121408131******************************\n39。)。fprintf(fid,39。*****************潮流計(jì)算輸出結(jié)果*************\n39。)。fprintf(fid,39。迭代次數(shù)k為: %d \n39。,k)。fprintf(fid,39。================================================\n39。)。fprintf(fid,39。各節(jié)點(diǎn)的實(shí)際電壓標(biāo)么值U為(節(jié)點(diǎn)號從小到大排列): \n39。)。for m=1:nfprintf(fid,39。第%d個實(shí)際電壓標(biāo)么值U:%f+i*(%f)\n39。,m,real(U(m)), imag(U(m))39。 )。endfprintf(fid,39。================================================\n39。)。fprintf(fid,39。各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): \n39。)。for m=1:n fprintf(fid,39。第%d個節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): %f\n39。,m,V(m))。endfprintf(fid,39。================================================\n39。)。fprintf(fid,39。各節(jié)點(diǎn)的電壓相角shita為(節(jié)點(diǎn)號從小到大排列): \n39。)。for m=1:n fprintf(fid,39。第%d個節(jié)點(diǎn)的電壓相角時shita為(節(jié)點(diǎn)號從小到大排列): %f\n39。,m,shita(m))。endfprintf(fid,39。================================================\n39。)。fprintf(fid,39。各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列): \n39。)。for m=1:nfprintf(fid,39。第%d個的復(fù)功率: %f+i*(%f)\n39。,m,real(S(m)),imag(S(m)) 39。)。endfprintf(fid,39。================================================\n39。)。fprintf(fid,39。節(jié)點(diǎn)有功功率P為: \n39。)。for m=1:n fprintf(fid,39。第%d個節(jié)點(diǎn)有功功率: %f\n39。,m,real(S(m)))。endfprintf(fid,39。================================================\n39。)。fprintf(fid,39。節(jié)點(diǎn)無功功率Q為: \n39。)。for m=1:n fprintf(fid,39。第%d個節(jié)點(diǎn)無功功率Q: %f\n39。,m,imag(S(m)))。endfprintf(fid,39。================================================\n39。)。fprintf(fid,39。各條支路的功率損耗DS為: \n39。)。n1=6。for m=1:n1fprintf(fid,39。各條支路的功率損耗為:%f+i*(%f)\n’,real(DS(m)), imag(DS (m)))。endfclose(fid)。 %關(guān)閉文件,程序結(jié)束 算例運(yùn)行結(jié)果輸出程序運(yùn)行結(jié)束后生成的TXT文檔如下所示：******************************電力系統(tǒng)08 李進(jìn)才 121408131*************************************************潮流計(jì)算輸出結(jié)果*************迭代次數(shù)k為: 101 ================================================各節(jié)點(diǎn)的實(shí)際電壓標(biāo)么值U為(節(jié)點(diǎn)號從小到大排列): 第1個實(shí)際電壓標(biāo)么值U: +i*()第2個實(shí)際電壓標(biāo)么值U: +i*()第3個實(shí)際電壓標(biāo)么值U: +i*()第4個實(shí)際電壓標(biāo)么值U: +i*()第5個實(shí)際電壓標(biāo)么值U: +i*()第6個實(shí)際電壓標(biāo)么值U: +i*()================================================各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): 第1個節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): 第2個節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): 第3個節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): 第4個節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): 第5個節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): 第6個節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列): ================================================各節(jié)點(diǎn)的電壓相角shita為(節(jié)點(diǎn)號從小到大排列): 第1個節(jié)點(diǎn)的電壓相角shita為(節(jié)點(diǎn)號從小到大排列): 第2個節(jié)點(diǎn)的電壓相角shita為(節(jié)點(diǎn)號從小到大排列): 第3個節(jié)點(diǎn)的電壓相角shita為(節(jié)點(diǎn)號從小到大排列): 第4個節(jié)點(diǎn)的電壓相角shita為(節(jié)點(diǎn)號從小到大排列): 第5個節(jié)點(diǎn)的電壓相角shita為(節(jié)點(diǎn)號從小到大排列): 第6個節(jié)點(diǎn)的電壓相角shita為(節(jié)點(diǎn)號從小到大排列): ================================================各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列): 第1個節(jié)點(diǎn)復(fù)功率: +i*()第2個節(jié)點(diǎn)復(fù)功率: +i*()第3個節(jié)點(diǎn)復(fù)功率: +i*()第4個節(jié)點(diǎn)復(fù)功率: +i*()第5個節(jié)點(diǎn)復(fù)功率: +i*()第6個節(jié)點(diǎn)復(fù)功率: +i*()================================================節(jié)點(diǎn)有功功率P為: 第1個節(jié)點(diǎn)有功功率: 第2個節(jié)點(diǎn)有功功率: 第3個節(jié)點(diǎn)有功功率: 第4個節(jié)點(diǎn)有功功率: 第5個節(jié)點(diǎn)有功功率: 第6個節(jié)點(diǎn)有功功率: ================================================節(jié)點(diǎn)無功功率Q為: 第1個節(jié)點(diǎn)無功功率Q: 第2個節(jié)點(diǎn)無功功率Q: 第3個節(jié)點(diǎn)無功功率Q: 第4個節(jié)點(diǎn)無功功率Q: 第5個節(jié)點(diǎn)無功功率Q: 第6個節(jié)點(diǎn)無功功率Q: ================================================各條支路的功率損耗DS為: 各條支路的功率損耗為: +i*() 各條支路的功率損耗為: +i*() 各條支路的功率損耗為: +i*() 各條支路的功率損耗為: +i*() 各條支路的功率損耗為: +i*() 各條支路的功率損耗為: +i*() 潮流計(jì)算GUI界面我們在matlab中繪制出的圖形是由一些圖形對象的基本圖形對象組成，各個圖形對象的具體顯示效果，是依舊所具有的能控制其特性的屬性來決定的。 用戶和windows應(yīng)用程序簡的信息交換是通過控件實(shí)現(xiàn)的，控件不僅是一般窗口的主要組成部分，也是構(gòu)成對話框的主要部件?？丶泻芏鄬傩?，以便用戶根據(jù)需要設(shè)計(jì)不同種類的控件。matlab提供10種控件分別為：命令按鈕（push button）；單選按鈕(radio button)；檢查框(cheek bos)；列表框(list box)；下拉菜單(popup menu)；滑塊(slider)；編輯框(edit box)；靜態(tài)文字(static text)；框架(frame)。 首先介紹MATLAB的工作界面點(diǎn)擊MATLAB的桌面快捷方式圖標(biāo)，直接啟動MATLAB軟件，或者運(yùn)行MATLAB的安裝目錄的快捷起動圖標(biāo)，啟動MATLAB軟件。啟動后的MATLAB操作界面的默認(rèn)情況（Default Desktop Layout）下有3個上層窗口：指令（命令）窗口（Command Window），工作空間瀏覽器窗口（Workspace Browser）和歷史指令（命令）窗口（Command Hisdory），如下圖所示：圖54 MATLAB的工作界面在指令（命令）窗口（Command Window）內(nèi)輸入GUIDE命令，按下回車鍵，打開如下對話框：圖55 新建GUI界面點(diǎn)擊OK，即可打開GUIDE的空白操作界面。圖56 空白GUI界面在GUIDE的設(shè)計(jì)編輯區(qū)layout area，拖入一個static text鍵到設(shè)計(jì)區(qū)，選中static text鍵，拖拽其邊沿覆蓋整個設(shè)計(jì)區(qū)layout area，點(diǎn)擊右鍵，選中彈出菜單中的send to back設(shè)為低層，雙擊打開其property inspector屬性檢查器，將string欄修改為空白，BackgroundColor背景顏色欄修改為淡綠色（可以根據(jù)個人愛好，修改為任意顏色），其它欄的屬性不變，這樣就完成背景顏色的設(shè)置。再拖入一個static text鍵到設(shè)計(jì)區(qū)，將其BackgroundColor背景顏色欄修改為灰色，在ForegroundColor欄把字體的顏色改為黃色，在FontSize字體大小欄把字號改為20，將string欄修改為“電力系統(tǒng)潮流計(jì)算”， 其它欄的屬性不變。其它各文本框的設(shè)置都仿照此方法編輯。在拖入一個push button鍵，將其BackgroundColor背景顏色欄修改為藍(lán)色，在ForegroundColor欄把字體的顏色改為黑色，在FontSize字體大小欄把字號改為20，將string欄修改為“開始計(jì)算”，其它欄的屬性不變。下面，著重介紹返回函數(shù)call back欄的設(shè)置。點(diǎn)擊GUIDE的設(shè)計(jì)編輯區(qū)layout area內(nèi)的push button鍵，即“開始計(jì)算”，右鍵選中View Callback,左鍵點(diǎn)擊Callback，這就進(jìn)入了push button1回調(diào)函數(shù)的編輯界面。在下面插入潮流計(jì)算的迭代程序，并添加以下程序代碼（結(jié)果顯示）：set(,39。String39。,e(1))。 set(,39。String39。,e(4))。set(,39。String39。,f(1))。 set(,39。String39。,f(4))。set(,39。String39。,e(2))。 set(,39。String39。,e(5))。set(,39。String39。,f(2))。 set(,39。String39。,f(5))。set(,39。String39。,e(3))。 set(,39。String39。,e(6))。set(,39。String39。,f(3))。 set(,39。String39。,f(6))。set(,39。String39。,V(1))。 set(,39。String39。,V(4))。set(,39。String39。,V(2))。 set(,39。String39。,V(5))。set(,39。String39。,V(3))。 set(,39。String39。,V(6))。set(,39。String39。,shita(1))。 set(,39。String39。,shita(4))。set(,39。String39。,shita(2))。 set(,39。String39。,shita(5))。set(,39。String39。,shita(3))。 set(,39。String39。,shita(6))。set(,39。String39。,real(S(1)))。 set(,39。String39。,real(S(4)))。set(,39。String39。,imag(S(1)))。 set(,39。String39。,imag(S(4)))。set(,39。String39。,real(S(2)))。 set(,39。String39。,real(S(5)))。set(,39。String39。,imag(S(2)))。 set(,39。String39。,imag(S(5)))。set(,39。String39。,real(S(3)))。 set(,39。String39。,real(S(6)))。set(,39。String39。,imag(S(3)))。 set(,39。String39。,imag(S(6)))。set(,39。String39。,real(DS(1)))。 set(,39。String39。,real(D