電氣的地工程及其自動(dòng)化專(zhuān)業(yè)綜合實(shí)訓(xùn)一{新}

1、實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 電氣工程及其自動(dòng)化專(zhuān)業(yè)綜合實(shí)訓(xùn)一 原始資料 網(wǎng)絡(luò)接線(xiàn)如圖，各支路阻抗和各節(jié)點(diǎn)功率均以標(biāo)幺值標(biāo)于圖中， 其中節(jié)點(diǎn) 2 2 連接的實(shí)際是發(fā)額定功率的發(fā)電廠(chǎng)。設(shè)節(jié)點(diǎn) 1 1 的電壓 保持為 1.061.06，試計(jì)算圖示系統(tǒng)中的潮流分布。（計(jì)算精度 1010-5） 二、等值電路圖 電力系統(tǒng)潮流計(jì)算 二、報(bào)告要求： 1 1、 寫(xiě)出潮流計(jì)算的原理 2 2、 寫(xiě)出功率方程（極坐標(biāo)或直角坐標(biāo)） 3 3、 寫(xiě)出雅克比矩陣元素的表達(dá)式或 PQPQ 分解法中的系數(shù)矩陣 4 4、 利用 MATLABLEMATLABLE 語(yǔ)言或其他語(yǔ)言編程調(diào)試 5 5、 心得或體會(huì) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 6

2、6、 附程序代碼，打印 潮流計(jì)算步驟： 1、寫(xiě)出潮流計(jì)算的原理 潮流是指在發(fā)電機(jī)母線(xiàn)上功率被注入網(wǎng)絡(luò)，而在變（配）電站的 母線(xiàn)上接入負(fù)荷，其間，功率在網(wǎng)絡(luò)中流動(dòng)。對(duì)于這種流動(dòng)的功率， 電力生產(chǎn)部門(mén)稱(chēng)之為潮流。以電力網(wǎng)絡(luò)潮流、電壓計(jì)算為主要內(nèi)容的 電力網(wǎng)絡(luò)穩(wěn)態(tài)行為特性計(jì)算的目的在于估計(jì)對(duì)用戶(hù)電力供應(yīng)的質(zhì)量以 及為電力網(wǎng)運(yùn)行的安全性與經(jīng)濟(jì)性評(píng)估提供基礎(chǔ)數(shù)據(jù)。配電網(wǎng)潮流計(jì) 算是配電網(wǎng)絡(luò)分析的基礎(chǔ)，配電網(wǎng)的網(wǎng)絡(luò)重構(gòu)、無(wú)功功率優(yōu)化、狀態(tài) 估計(jì)和故障處理都需要用到配電網(wǎng)潮流數(shù)據(jù)。 電力系統(tǒng)穩(wěn)態(tài)運(yùn)行應(yīng)滿(mǎn)足以下要求： 1 1） 滿(mǎn)足系統(tǒng)經(jīng)濟(jì)性運(yùn)行的要求，每一臺(tái)發(fā)電機(jī)的輸出必須接近于 預(yù)先設(shè)定值； 2 2） 必

3、須確保聯(lián)絡(luò)線(xiàn)潮流低于線(xiàn)路熱極限和電力系統(tǒng)穩(wěn)定極限； 3 3） 必須保持某些中樞點(diǎn)母線(xiàn)上的電壓水平在容許范圍內(nèi)， 必要時(shí) 用無(wú)功功率補(bǔ)償計(jì)劃來(lái)達(dá)到； 4 4） 區(qū)域電網(wǎng)是互聯(lián)系統(tǒng)的一部分，必須執(zhí)行合同規(guī)定的輸送至鄰 網(wǎng)的聯(lián)絡(luò)線(xiàn)功率計(jì)劃； 5 5 ）用故障前的潮流控制策略使事故擾動(dòng)效應(yīng)最小化。 通常情況下，輸電線(xiàn)路電壓在輕載時(shí)會(huì)較高，重載時(shí)會(huì)較低，電 壓調(diào)整是指在負(fù)載由輕載到滿(mǎn)載變化過(guò)程中實(shí)時(shí)調(diào)整線(xiàn)路電壓滿(mǎn)足運(yùn) 行要求；對(duì)于超高壓輸電線(xiàn)路，線(xiàn)路電壓維持在額定電壓的士 5%5%之內(nèi), , 實(shí)際運(yùn)行時(shí)，通常電壓調(diào)整約為 10%10%。對(duì)于低壓輸電線(xiàn)路，電壓調(diào) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 整數(shù)值為 10%1

4、0%，包含了變壓器本身的電壓降落。 .1 潮流計(jì)算的基本物理量 潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，它的任務(wù)是對(duì) 給定的運(yùn)行條件確定系統(tǒng)的運(yùn)行狀態(tài)，就是在三相平衡穩(wěn)態(tài)狀態(tài)下計(jì) 算電力系統(tǒng)中每條母線(xiàn)的電壓幅值和相角，其中每一設(shè)備如傳輸線(xiàn)和 變壓器中的有功和無(wú)功潮流，以及各設(shè)備的損耗都需要計(jì)算出來(lái)。 潮流計(jì)算采用電力系統(tǒng)的單線(xiàn)圖，對(duì)于任意一條母線(xiàn) i i,需要以下 四個(gè)變量描述：電壓幅值 UiUi、相角，電網(wǎng)供給母線(xiàn)的有功 PiPi、無(wú)功 QiQi。若某一電力系統(tǒng)有 N N 個(gè)節(jié)點(diǎn)，則共有 4N4N 個(gè)變量，對(duì)于每條母線(xiàn)， 這些變量中的兩個(gè)指定為輸入數(shù)據(jù)，其它的兩個(gè)是潮流程

5、序所要計(jì)算 的未知量。為方便起見(jiàn)，在圖 3.13.1 中傳送給母線(xiàn) i i的功率可分為發(fā)電機(jī) 發(fā)出和負(fù)載吸收兩部分。也就是 Pi Pi = = PGi PGi - PLiPLi Qi Qi = = QGi QGi - QLiQLi I 5 圖 3 3- -1 1 節(jié)點(diǎn)的變量 每條母線(xiàn)被歸分為以下三種母線(xiàn)類(lèi)型中的某一類(lèi): 1 1） 平衡節(jié)點(diǎn)，一般一個(gè)系統(tǒng)只有一個(gè)平衡節(jié)點(diǎn)。在潮流分布 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 算出以前，網(wǎng)絡(luò)中的功率損耗是未知的，因此，至少有一個(gè)節(jié)點(diǎn)的 有功功率 P P 和無(wú)功功率 Q Q 不能給定。另外必須選定一個(gè)節(jié)點(diǎn)，制 定其電壓相角為零，作為其它節(jié)點(diǎn)電壓相位的參考，這個(gè)節(jié)點(diǎn)叫基

6、 準(zhǔn)節(jié)點(diǎn)。為了計(jì)算方便，常將平衡節(jié)點(diǎn)和基準(zhǔn)節(jié)點(diǎn)設(shè)在同一個(gè)節(jié)點(diǎn) 上。為方便起見(jiàn)在本書(shū)中把它標(biāo)號(hào)為母線(xiàn) 1 1。平衡節(jié)點(diǎn)是電壓參考 節(jié)點(diǎn)，該母線(xiàn)的-是給定值，作為輸入數(shù)據(jù)，典型取標(biāo)幺值 。 潮流程序計(jì)算 P1P1 和 Q1Q1。因?yàn)槠胶夤?jié)點(diǎn)的 P P、Q Q 事先無(wú)法確定， 為使潮流計(jì)算結(jié)果符合實(shí)際，常把平衡節(jié)點(diǎn)選在有較大調(diào)節(jié)裕量的 發(fā)電機(jī)節(jié)點(diǎn)，潮流計(jì)算結(jié)束時(shí)若平衡節(jié)點(diǎn)的有功功率、 無(wú)功功率和 實(shí)際情況不符，就要調(diào)整其他節(jié)點(diǎn)的邊界條件以使平衡節(jié)點(diǎn)的功率 滿(mǎn)足實(shí)際允許范圍。 2 2） PQPQ 節(jié)點(diǎn)，PiPi 和 QiQi 是輸入數(shù)據(jù)。這類(lèi)節(jié)點(diǎn)的有功功率 Pi Pi 和無(wú)功功率 QiQi是給定的，潮

7、流計(jì)算程序計(jì)算節(jié)點(diǎn)電壓幅值 UiUi 和相 角送。負(fù)荷節(jié)點(diǎn)和無(wú)功功率注入的聯(lián)絡(luò)節(jié)點(diǎn)都屬于這類(lèi)節(jié)點(diǎn)。有些 情況下，系統(tǒng)中某些發(fā)電廠(chǎng)送出的功率在一定時(shí)間內(nèi)為固定時(shí)，該 發(fā)電廠(chǎng)母線(xiàn)也可以作為 PQPQ 節(jié)點(diǎn)。在一個(gè)典型的潮流程序中絕大多 數(shù)母線(xiàn)作為 PQPQ 節(jié)點(diǎn)。 3 3） PUPU 節(jié)點(diǎn)（電壓控制母線(xiàn)），PiPi 和 UiUi 是輸入數(shù)據(jù)。這類(lèi)節(jié) 點(diǎn)的有功功率 PiPi和節(jié)點(diǎn)電壓幅值 UiUi 是給定的，潮流程序計(jì)算節(jié)點(diǎn) 的無(wú)功功率 QiQi 和電壓相角-。這類(lèi)節(jié)點(diǎn)必須具有足夠的無(wú)功可調(diào) 容量，用以保持給定的節(jié)點(diǎn)電壓幅值。在電力系統(tǒng)中這類(lèi)節(jié)點(diǎn)的數(shù) 目較少。例如與發(fā)電機(jī)、并聯(lián)補(bǔ)償電容器或者靜止無(wú)功

8、系統(tǒng)相連的 母線(xiàn)。設(shè)備無(wú)功功率最大值 QGimaxQGimax 和最小值 QGiminQGimin 都是輸入 數(shù)據(jù)。另一個(gè)例子是與抽頭可調(diào)節(jié)實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 變壓器相連的母線(xiàn)， 用潮流程序 計(jì)算抽頭的位置。 注意，當(dāng)母線(xiàn) i i 是無(wú)發(fā)電機(jī)相連接的負(fù)載母線(xiàn)時(shí)，Pi=Pi=- -PLiPLi 為 負(fù)值；也就是說(shuō)，在圖 3.13.1 中給母線(xiàn) i i 提供的有功為負(fù)值。如果負(fù) 荷是感性的，QiQi 二-QLiQLi 為負(fù)值。 綜上所述，若系統(tǒng)中有 n n 個(gè)節(jié)點(diǎn)，n n 為網(wǎng)絡(luò)中除去參考節(jié)點(diǎn)外 的節(jié)點(diǎn)數(shù)，本書(shū)中以大地為參考節(jié)點(diǎn)，選第 1 1 個(gè)節(jié)點(diǎn)為平衡節(jié)點(diǎn)， 剩下的 n n- -1 1 個(gè)

9、節(jié)點(diǎn)中有 r r 個(gè) PUPU 節(jié)點(diǎn)，則有 n n- -r r- -1 1 個(gè) PQPQ 節(jié)點(diǎn)。因 此，除了平衡節(jié)點(diǎn)外，有 n n- -1 1 個(gè)節(jié)點(diǎn)的注入有功功率、n n- -r r- -1 1 個(gè) PQ PQ 節(jié)點(diǎn)的注入無(wú)功功率和 r r 個(gè) PUPU 節(jié)點(diǎn)電壓幅值為已知量。 .2 潮流計(jì)算的數(shù)學(xué)模型 在穩(wěn)態(tài)潮流計(jì)算中，電力系統(tǒng)各元件（參數(shù)）等效成一個(gè)有源 網(wǎng)絡(luò)。將發(fā)電機(jī)和負(fù)荷用無(wú)阻抗線(xiàn)從網(wǎng)絡(luò)中抽出，剩下的是由接地 和不接地支路組成的無(wú)源線(xiàn)性網(wǎng)絡(luò)，可以用導(dǎo)納矩陣（ Y Y）或阻抗 矩陣（Z Z）來(lái)描述。 采用導(dǎo)納矩陣時(shí)，節(jié)點(diǎn)電流和節(jié)點(diǎn)電壓構(gòu)成以下方程： T T （ 3 3-

10、 -1 1 ） 其中：丫為 n n xn n 階導(dǎo)納矩陣，其階數(shù) n n 為網(wǎng)絡(luò)中除去參考節(jié) 點(diǎn)外的節(jié)點(diǎn)數(shù)，如果不考慮網(wǎng)絡(luò)元件的非線(xiàn)性及變壓器的相位偏 移，Y Y 為對(duì)稱(chēng)矩陣，】為 n n X1 1 維節(jié)點(diǎn)注入電流列向量，在電力系統(tǒng) 計(jì)算中，節(jié)點(diǎn)注入電流可理解為該節(jié)點(diǎn)電源電流與負(fù)荷電流之和， 并規(guī)定流入節(jié)點(diǎn)電流為正。因此僅有負(fù)荷的節(jié)電電流就為負(fù)值，某 些僅起聯(lián)絡(luò)作用的節(jié)點(diǎn)，圖 3 3- -2 2中節(jié)點(diǎn) n=3n=3，其注入電流為零。 為 n n X1 1 維節(jié)點(diǎn)電壓列向量。網(wǎng)絡(luò)中有接地支路實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 時(shí)，節(jié)點(diǎn)電壓通常 指該節(jié)點(diǎn)的對(duì)地電壓，以大地作為參考節(jié)點(diǎn)，并規(guī)定其編號(hào)為零。 圖

11、3 3- -2 2 運(yùn)用節(jié)點(diǎn)電壓法時(shí)的電力網(wǎng)路等值電路 對(duì)于第 i i 個(gè)節(jié)點(diǎn)，展開(kāi)為如下形式： 河 (3 3- -2 2 ) 若米用阻抗矩陣可表示為： 展開(kāi)為： 川 (3 3- -3 3 ) 在潮流計(jì)算時(shí)一般以節(jié)點(diǎn)電壓方程進(jìn)行。節(jié)點(diǎn)導(dǎo)納矩陣與阻抗 矩陣互為逆陣，在短路計(jì)算時(shí)可直接利用導(dǎo)納矩陣求逆得到阻抗矩 陣以求得短路點(diǎn)的短路電流。 由于實(shí)際系統(tǒng)中一般不給出節(jié)點(diǎn)電流而是節(jié)點(diǎn)功率， 因此將式 (3 3- -2 2 )中的節(jié)點(diǎn)注入電流用節(jié)點(diǎn)注入功率來(lái)表示為： 如果節(jié)點(diǎn)電壓用極坐標(biāo)表示，令 Ui = U嚴(yán) 5 = U嚴(yán) rh Si=jQ = U=Ut LH j 丿=1憶，嚴(yán) (3(3- -4 4

12、) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 n n 個(gè)節(jié)點(diǎn)電力系統(tǒng)的潮流方程的一般形式是 $三百七2三5碼丿is）*1憶 丿 問(wèn) （3 3- -5 5） 臥耳十為呂沖百-毎-培） Q二q土誡應(yīng)他-闌 $=耳何口鮒-勺）+站啲-切 X Q二返巧6 血 一如耳噸 級(jí) ?-1 若米用直角坐標(biāo)系，節(jié)點(diǎn)電壓可表示為 導(dǎo)納矩陣元素可以表示為 玲-磅十為 將上述表達(dá)式帶入式（2 2- -8 8 ）的右端，展開(kāi)并分出實(shí)部和虛部, 便得 耳=兮工q號(hào)-馬 Q *力刀橢力-亠烏弓） /-I 2二總嗎巒恥）弋壬（釗十罷） z -1 可見(jiàn)，原來(lái)電流電壓的線(xiàn)性方程組變換為功率和電壓的非線(xiàn)性 方程組，式（3 3- -6 6）（ 3 3-

13、 -7 7 ）就是潮流計(jì)算的基本方程。它是一組共 有 n n 個(gè)非線(xiàn)性方程組成的復(fù)數(shù)方程組，如果把實(shí)部和虛部分開(kāi)便得 到 2n2n 個(gè)實(shí)數(shù)方程，由該方程組可解出 2n2n 個(gè)運(yùn)行參數(shù)。但是每一 個(gè)節(jié)點(diǎn)都有 P P、Q Q、U U、四個(gè)運(yùn)行變量，共有 4n4n 個(gè)運(yùn)行參數(shù)，所 以要事先給定其余(3(3- -6 6) (3(3- -7 7 ) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 2n2n 個(gè)參數(shù)。這就要根據(jù)節(jié)點(diǎn)的分類(lèi)，將每個(gè)節(jié) 點(diǎn)的 4 4 個(gè)運(yùn)行參數(shù)中的兩個(gè)作為原始數(shù)據(jù)，另外兩個(gè)作為待求量。 .3 潮流計(jì)算的約束條件 為了保證電力系統(tǒng)的正常運(yùn)行，潮流問(wèn)題中某些變量應(yīng)滿(mǎn)足一 定的約束條件，常用

14、的約束條件有： （1 1）所有節(jié)點(diǎn)電壓必須滿(mǎn)足 I I 入卩河込5 從保證電能質(zhì)量和供電安全的要求來(lái)看， 電力系統(tǒng)的所有電氣 設(shè)備都必須運(yùn)行在額定電壓附近。PUPU 節(jié)點(diǎn)的電壓幅值必須按上述 條件給定。因此，這一定約束條件主要是對(duì) PQPQ 節(jié)點(diǎn)而言 （2 2）所有電源節(jié)點(diǎn)的有功功率和無(wú)功功率必須滿(mǎn)足 PQPQ 節(jié)點(diǎn)的有功功率和無(wú)功功率以及 PUPU 節(jié)點(diǎn)有功功率，在給 定時(shí)就必須滿(mǎn)足此條件。因此，對(duì)平衡節(jié)點(diǎn)的 P P 和 Q Q 以及 PUPU 節(jié) 點(diǎn)的 Q Q 應(yīng)按此條件進(jìn)行檢驗(yàn)。 (3)(3) 某些節(jié)點(diǎn)之間電壓的相位差應(yīng)滿(mǎn)足 燉-即姬-孚 L L 為了保證系統(tǒng)運(yùn)行的穩(wěn)定性，要求某些輸電線(xiàn)

15、兩端的電壓相位差不 超過(guò)一定的數(shù)值。因此，潮流計(jì)算可以歸結(jié)為求解一組非線(xiàn)性方程 組，并使其解答滿(mǎn)足一定的約束條件。如果不能滿(mǎn)足，則應(yīng)修改某 些變量。甚至修改系統(tǒng)的運(yùn)行方式，重新進(jìn)行計(jì)算。 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 2、 寫(xiě)出功率方程(極坐標(biāo)或直角坐標(biāo)) 3、 寫(xiě)出雅克比矩陣元素的表達(dá)式或 PQ 分解法中的系數(shù) 矩陣 (2 (2 )、形成導(dǎo)納矩陣: fi.46-j6.3l -5.19+J2.7 -1 27tj3 54 0 o -5.19+J277 10.53j25.09 -l.72-hj4.31 -1.72+J4.31 2+je YB- -L27+j3.54 -1.27+J131 M4j25 叩

16、-6.70tjl7.24 0 0 *1.27+j4.31 -6.70+j 17.24 9.8-j25 4 3S+j34 0 -2+j 6 0 -1.3E+J3.4 3.38-j?.45 (3)(3)、B B ,形成以及其逆矩陣。 由題可知除 1 1 為平衡節(jié)點(diǎn)外，其他節(jié)點(diǎn)均為 PQPQ 節(jié)點(diǎn)，系數(shù) B B B B 階數(shù)相同，又因?qū)υ摰戎稻W(wǎng)絡(luò)不存在去除與有功功率和電壓或無(wú)功功 率和電壓大小關(guān)系較小因數(shù)的可能性， 這兩個(gè)矩陣 B B B B完全相同, 它們就由導(dǎo)納矩陣的虛數(shù)部分中除第一行和第二實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 行的各個(gè)元素組成。 由此可見(jiàn)，網(wǎng)絡(luò)的節(jié)點(diǎn)導(dǎo)納矩陣為奇異矩陣，但他的虛數(shù)部分的子 矩

17、陣 B B和 3 3，貝煜非記異矩陣，可以求逆矩陣，其逆矩陣： B1 B1 0.0768 0.0552 0.0611 0.0711 0.0552 0.1192 0.1017 0.0722 0.0611 0.1017 0.1327 0.0873 0.0711 0.0722 0.0873 0.1828 （4 4）、計(jì)算各節(jié)點(diǎn)有功功率不平衡量Vpi。 取 Ui 1.06,1 0U2（0）U3（0）U4（0）U 1.0；2（0）3（0）4（0）5（0）0 , ,按 下式計(jì)算各節(jié)點(diǎn)有功率不平衡： j n vp(0) Pi Ui(0)Uj(0) (GjjCOSjj BijSinjj) j i 25.092

18、5.09 4.314.31 4.314.31 6 6 4.314.31 25.0925.09 17.2417.24 0 0 B=B = = 4.314.31 17.2417.24 2525 3.453.45 6 6 0 0 3.453.45 9.459.45 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 Vp2(0) 0.251.0*1.06( 1.59o0 2.77*sir0 1.0*1.010.63Co0 25.9ir0 1.0*1.01.72Co0 4.31Sirt) 1.0*10 1.72Co0 4.3*sir0 1.0*1.*( 2cos06sin0) 0.584 相似的可得： Vp/0) 0.9Vp4(

19、0) 0.4Vp5(0) 0.6 ； (5) 、計(jì)算各節(jié)點(diǎn)電壓的相位角i : 由下列矩陣方程式： B 1(VP(0)/U(0) (U(0)V (0)實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 0.0768 0.0552 0.0611 0.0711 0.5814 1.0* (0) 2 0.0552 0.1192 0.1017 0.0722 0.9 /1.0 1.0* (0) 3 0.0611 0.1017 0.1327 0.0873 0.4 /1.0 = =1.0* (0) 4 0.0711 0.0722 0.0873 0.1823 0.6 /1.0 1.0* (0) 5 從而可得: (0) 2 0.0721 (1

20、) 2 (0) 2 0.0721 (0) 3 0.1592 (1) 3 (0) 3 0.1592 (0) 4 0.1615 ; (1) 4 (0) 4 0.1615 (0) 5 0.1679 (1) 5 (0) 5 0.1679 (6) 、計(jì)算各節(jié)點(diǎn)無(wú)功功率不平衡VQi。 按下列計(jì)算各節(jié)點(diǎn)無(wú)功功率不平衡： j n VQj(0) Qj Uj(0)U j(0) (Gjj sin cos J j i實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 VQ,(0) 0.25 1.0 1.06 5.19 sin( 0.07210) 2.77cos(0.07210) 1.0 1.0 (10.63sin(0)25.09cos0) 1.

21、0 1.0 1.72 sin(0.07210 0.15920)4.31 cos(0.072100.15920) 1.0 1.0 1.72 sin( 0.072100.1615)4.31 co$ 0.072100.16150) 1.0 1.0 2 sin(0.0721 0.16790)6 cos(0.0721 0.16790) 0.264267 相似可得: VQ(0) 0.46527V4(0) 0.77744VQ(0) 0.757824; (7(7)、計(jì)算各節(jié)點(diǎn)電壓的大小 5。 由下列矩陣方程式： IIII B 1 (VQ (0) / U (0) VU (0) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 0.076

22、80.05520.0611 0.711 7.264267也(0) 0.05520.11920.10170.0722 0.465275VU(0) 0.06110.10170.13270.08730.777443VU4(0) 0.07110.07220.08730.18280.757824VU5(0) VU(0) 0.684962 VU(0) 0.590229 V4(0) 0.660490 V5(0) 0.756140 U2 U)VU(0) 0.31503U U)U(0)U4 U)VJ4(0) 0.33951U U)U(0) 0.243896 4、利用 MATLABLE 語(yǔ)言或其他語(yǔ)言編程調(diào)試 P

23、Q 分解法潮流程序 /文件輸入格式：節(jié)點(diǎn)總數(shù) n（包括聯(lián)絡(luò)節(jié)點(diǎn)），支路數(shù) zls / /節(jié)點(diǎn)數(shù)（發(fā)電機(jī)和負(fù)荷）nb,接地電抗數(shù) mdk,迭代精度 eps / /考慮負(fù)荷靜特性標(biāo)志 kk2（0 考慮），平衡節(jié)點(diǎn)號(hào)，優(yōu)化標(biāo)志（0 不優(yōu)化） / /最大迭代次數(shù) it1,支路左右節(jié)點(diǎn)號(hào) izl,jzl,支路電阻 zr,電抗 zx / 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 / 支路容納 zyk,節(jié)點(diǎn)號(hào) nob及標(biāo)志 nobt（0-PQ,-1-PV） /發(fā)電機(jī)和負(fù)荷有功、無(wú)功 pg,qg,pl,ql / /電壓 vO（pv 節(jié)點(diǎn)輸入實(shí)際值，PQ 節(jié)點(diǎn)任輸入一值） / /電抗節(jié)點(diǎn)號(hào) idk,電抗值 dkk / #in e

24、lude math.h #i nclude stdio.h FILE *fp1, *fp2; char inn ame12, out name12; / fp1 輸入數(shù)據(jù)文件指針 fp2 輸出文件指針 / inname輸入數(shù)據(jù)文件名 outname輸出數(shù)據(jù)文件名 int n, zls, nb, mdk, mpj, bn sopt on, it1, dsd, kk2, n zls; /節(jié)點(diǎn)總數(shù) n（包括聯(lián)絡(luò)節(jié)點(diǎn)） 支路數(shù)（回路數(shù)）zls 節(jié)點(diǎn)數(shù) nb（發(fā)電機(jī)和負(fù)荷） #defi ne NS 2000 /最大節(jié)點(diǎn)數(shù) #defi ne NS2 NS * 2 #defi ne NS4 1000 /NS

25、4、NS 必須大于 #defi ne ZS 3000 /最大支路數(shù) #defi ne ZS2 ZS * 2 #defi ne DKS 200 /最大電抗器數(shù) #defi ne N2 ZS * 4 #defi ne N3 ZS * 8 + NS * 4 2*zls / 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 II 接地電抗數(shù) mdk 精度 eps 平衡節(jié)點(diǎn)號(hào) mpj / 節(jié)點(diǎn)優(yōu)化（標(biāo)志）bnsopton（=0 節(jié)點(diǎn)不優(yōu)化，!=0 節(jié)點(diǎn)優(yōu)化） II 最大迭代次數(shù) it1 最低電壓或最大功率誤差節(jié)點(diǎn)號(hào) dsd II 負(fù)荷靜特性標(biāo)志（=0 考慮負(fù)荷靜特性） II 支路數(shù)（雙回線(xiàn)算一條支路） int izlZS, jz

26、lZS, idkDKS, ydsNS, ydzNS, iyZS2; II izl,jzl,idk:分別存放左、右節(jié)點(diǎn)號(hào)和電抗器節(jié)點(diǎn)號(hào)。 II yds存放各行非零非對(duì)角元素的個(gè)數(shù)。 II ydzi是第 i 行第一個(gè)非零非對(duì)角元素的首地址， II 即在所有非零非對(duì)角元素中的次序號(hào) II iy存放列足碼。 int nnewNS4, oldNS, nobNS, nobtNS; II nn ew,old存放的是新、舊節(jié)點(diǎn)號(hào)。 II nnewi中為 i 對(duì)應(yīng)的新號(hào) II nob存放的是節(jié)點(diǎn)號(hào)。nobt存放的是節(jié)點(diǎn)類(lèi)型,0: pq 節(jié)點(diǎn),-1: pv 節(jié)點(diǎn) double eps, dsm, vmin, dp

27、h, dqh, af3; II eps 迭代收斂精度，dsm 最大功率誤差 II vmin:系統(tǒng)最低電壓值。dph,dqh:系統(tǒng)有、無(wú)功損耗。 II af0和 af1分別是負(fù)荷有功功率、無(wú)功功率靜態(tài)特性系數(shù)。 double v00; II v00:系統(tǒng)平均電壓 ci,cj 分別作為節(jié)點(diǎn) i,j 的電壓相角的臨時(shí)存儲(chǔ)單元。 double zrZS, zxZS, zykZS, dkkDKS, giiNS, biiNS, ygZS2, ybZS2; 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 double pgNS, qgNS, plNS, qlNS, vONS, vNS, vaNS; /支路電阻 zr 支路電抗 zx

28、輸電線(xiàn)路充電容納 zyk(y0/2) /接地電抗 dkk對(duì)角元實(shí)部 gii 對(duì)角元虛部 /非對(duì)角元實(shí)部 yg非對(duì)角元虛部 yb / pg,qg,pl,ql：發(fā)電機(jī)，負(fù)荷功率實(shí)、虛部 / v是電壓幅值,va是電壓相角。 double wNS2, kg3, bNS2; int n ewsortNS4; / newsorti存放 i 對(duì)應(yīng)的老號(hào) void in itial(); void pqflow(); void out(); void dataio(); void bn sopt(); void zlsort(i nt* nn ew); void prin to(); void prin ty

29、(); void y2(); void ya0(); void yzb(); void jdgl(int kq0); 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 void bbhl(int kqO); void calc(); int iabs(i nt a); void bran ch_output(); void n ewval(double* aa); void prin tc(); void iswap(); void swap(); void prin tf2(double* aa, double* bb, int n); void calc(i nt* iu, double* u, double* d

30、i, i nt* nfd, double* b); void printi(int* aa, int n); void prin tf1(double* aa, int n); int find(int k, int a, int* z); void yzb(i nt t, i nt* iu, double* u, double* di, i nt* n fd); int isg n(i nt a, int b); void yy1(); void y3(); void n ewtoold(); int main(v oid)實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 in itial(); /初始化 pqflo

31、w(); /pq 潮流計(jì)算 out(); /輸出節(jié)點(diǎn)和支路數(shù)據(jù) return 1; int isg n(i nt a, int b) /* 本函數(shù)功能返回值為 a 的絕對(duì)值 b 的符號(hào) */ /參數(shù) 1 提供值，參數(shù) 2 提供符號(hào)/ if (b 0) a = -a; return a; int find(int k, int a, int* z) /* 本函數(shù)查找 a中是否有 fabs(k)有則返回 0,無(wú)則返回 1 */ /參數(shù) 1 為待查找量，參數(shù) 2 待搜索數(shù)組，參數(shù) 3 返回 k 在 a中的次序號(hào)/ 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 int i; for (i = 1; i = n; i +)

32、if(iabs(k) = ai) *z = i; return 1; return 0; void oldt on ew() /* 本函數(shù)將輸入數(shù)據(jù)中的節(jié)點(diǎn)號(hào)變成從 1 開(kāi)始的連續(xù)節(jié)點(diǎn)號(hào) */ int i, j, k, ii1, ii2, zls2, k1, k2, k3, k4, ip; zls2 = zls + zls; for (i = 1; i = zls2; i +) n ewsorti = 0; ii1 = 0; for (i = 1; i = zls; i +) k = izli; if (!find(k, newsort, &ii2) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 ii1 +

33、; n ewsortii1 = iabs(k); k = jzli; if (!find(k, newsort, &ii2) ii1 +; n ewsortii1 = iabs(k); for (i = 1; i = ii1-1; i +) for (j = i+1; j newsortj) k = n ewsorti; n ewsorti = n ewsortj; n ewsortj = k; 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 for (i = 1; i = zls; i +) k = izli; if (find(k, newsort, &ii2) izli = isgn(ii2,

34、k); else prin tf(error!); k = jzli; if (find(k, newsort, &ii2) jzli = isgn(ii2, k); else prin tf(error!); prin tf(izl%d = %d, jzl%d = %dn, i, izli, i, jzli); for (i = 1; i = nb; i +) for (j = 1; j = n; j +) if (n obi = n ewsortj) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 n obi = j; break; printf(n ob%d = %dn, i, no bi); for

35、(j = 1; j = n; j +) if (mpj = n ewsortj) mpj = j; break; /電抗器節(jié)點(diǎn)號(hào)轉(zhuǎn)變 for (j = 1; j = mdk; j +) for (i = 1; i = n; i +) if (idkj = n ewsorti) idkj = i; break; 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 void in itial() /* 本函數(shù)進(jìn)行初始化工作 */ int i, k1; dataio();/輸入原始數(shù)據(jù) oldto new(); 轉(zhuǎn)化為新號(hào) if (bnsopton = 0) /節(jié)點(diǎn)不優(yōu)化，新節(jié)點(diǎn)號(hào)即為老節(jié)點(diǎn)號(hào)。 for (i = 1; i =

36、 n; i +) oldi = i; nn ewi = i; else bnsopt(); /節(jié)點(diǎn)優(yōu)化 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 mpj = nn ewmpj;/mpj: 平衡節(jié)點(diǎn) zlsort (nn ew); / sort the r,x and b for (i = 1; i = mdk; i +) k1 = idki; idki = nn ewk1; for (i = 1; i = n; i +) vi = vOO; vai = 0.0; / 所有節(jié)點(diǎn)的電壓幅值初值都為 1.000(v00)，電壓相角初值都為 0 / excha nge the node before and after

37、 sort for (i = 1; i = n; i +) ydsi = 0; / the immediate for (i = 1; i = nb; i +) k1 = nnewn obi; ydsk1 = nobti; for (i = 1; i = n; i +) n obti = ydsi; 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 n ewval(pg); n ewval(qg); n ewval(pl); n ewval(ql); n ewval(vO); for (i = 1; i = n; i +) / no bt is type of node if (n obti = -1) vi = v0

38、i; / n ob is serials nu mbe nobt = -1: pv 節(jié)點(diǎn)，v0存放的是最后一個(gè)節(jié)點(diǎn)數(shù)據(jù)， /對(duì)于 pv 節(jié)點(diǎn)，即為該點(diǎn)應(yīng)維持的電壓值。 nobt = 0: pq 節(jié)點(diǎn)，v0存放的是最后一個(gè)節(jié)點(diǎn)數(shù)據(jù), /對(duì)于 pq 節(jié)點(diǎn)，即為系統(tǒng)平均電壓值。 prin to(); /輸出 af、vOO 和節(jié)點(diǎn)排序后的支路、節(jié)點(diǎn)和 /接地電抗數(shù)據(jù)(僅僅查看中間結(jié)果) yaO(); 獲得 yds、ydz、列足碼 iy。( P407) void prin to() /* 輸出 af、v00 和節(jié)點(diǎn)排序后的支路、節(jié)點(diǎn)和接地電抗數(shù)據(jù) */ int i; 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 fprin

39、tf(fp2, n *AF AND V0 *n) fprin tf(fp2, n %7.3f%7.3f%7.3fn, af0, af1, v00); prin tc(-, 78); fprin tf(fp2, nn *ZLB*n) for (i = 1; i = zls; i +) fprin tf(fp2, n); fprin tf(fp2. %8d%8d%8d%8d, izli, jzli, oldabs(izli), oldabs(jzli); fprintf(fp2, %9.4f%9.4f%9.4f, zri, zxi, zyki); printc(-, 78); fprin tf(f

40、p2, nn*BUS*n); for (i = 1; i = nb; i +) fprin tf(fp2, n); fprintf(fp2, %8d%8d%8d, nobi,oldnobi, nobti); fprin tf(fp2, %9.4f%9.4f%9.4f%9.4f%9.4f, pgi, qgi, pli, qli, v0i); prin tc(-, 78); fprin tf(fp2,nn*DKK*n); for (i = 1; i = mdk; i +) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 fprin tf(fp2, n); fprintf(fp2, %8d%8d%7.4f, idki, o

41、ldidki, dkki); void dataio() II* 系統(tǒng)數(shù)據(jù)初始化 */ int i;實(shí)用標(biāo)準(zhǔn)文案 for (i = 1; i = nb; i +) 精彩文檔 afO = 0.6; af1 = 2.0;/af0 和 af1分別是負(fù)荷有功功率、無(wú)功功率靜態(tài)特性系數(shù)。 v00 = 1.000;/系統(tǒng)平均電壓 prin tf(nplease in put the n ame of data file n); sea nf(%s, inn ame); fp1 = fope n(i nn ame, r); prin tf(nplease output the n ame of data f

42、ile n); sea nf(%s, out name); fp2 = fope n(outn ame, w); fseanf(fp1, %d %d %d %d, &n, &zls, &nb, &mdk); / the nu mber of node ,bra nches, node fseanf(fp1, %lf %d %d %d %d, &eps, &kk2, &mpj, &bnsopton, &it1); /preeisi on, swing no de,sort the no de,iterati on nu mb

43、ers for (i = 1; i = zls; i +) fscanf(fp1, %d %d, &izli, &jzli); fscan f(fp1, %lf %lf %lf , &zri, & zxi, &zyki); 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 fscanf(fp1, %d %d, &nobi, &nobti); fscan f(fp1, %lf %lf %lf %lf %lf, & pgi, & qgi, & pli, &qli, &v0i); for (i = 1; i = mdk; i +)

44、 fscan f(fp1, %d %lf, & idki, & dkki); fclose(fp1); void pqflow() int kqO, iu1N2, nfd1NS, iu2N2, nfd2NS; int i, t; double u1N2, u2N2, di1NS, di2NS; yy1(); yzb(0, iu1, u1, di1, nfd1); /form the B matrix of P-0 iteration y2(); yzb(1, iu2, u2, di2, nfd2); /form the B matrix of Q-V iteration實(shí)用標(biāo)準(zhǔn)

45、文案 for (i = 1; i eps) kgkq0=1; if (kq0 = 0) calc(iu1, u1, di1, n fd1, b); if (kq0 = 1) calc(iu2, u2, di2, nfd2, b); if(kq0 = 0 )實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 vai = vai - bi / vOO; else vi = vi - bi; else kgkqO = 0； if(kqO = 0) kqO = 1; else kqO = O; t +; if(t it1) break; while(fabs(dsm) eps) | (kgkqO != 0); fprin tf(

46、fp2, n%s%d, times = , t); void out() 實(shí)用標(biāo)準(zhǔn)文案 for (i = 1; i = nb; i +) 精彩文檔 /* 本函數(shù)輸出節(jié)點(diǎn)和支路數(shù)據(jù)*/ zlsort(old); / recover the data if sorted / n ewtoold(); no de_output(); / node data bran ch_output(); /branch data prin tc(-, 78); prin tc(*, 78); fprin tf(fp2, n); void n ewval(double* aa) /* 本函數(shù)將舊號(hào)換成新號(hào) */

47、int i, k1; for (i = 1; i = n; i +) bi = 0.0; k1 = nnewn obi; bk1 = aai; for (i = 1; i = n; i +) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 aai = bi; void yzb(i nt t, i nt* iu, double* u, double* di, i nt* nfd) /* 本函數(shù)求因子表 */ /參數(shù) 1 為標(biāo)志(t=0 求 B,t=1 求 B)/ /參數(shù) 2 因子表上三角矩陣非零非對(duì)角元素的列足碼 /參數(shù) 3 因子表上三角矩陣非零非對(duì)角元素的數(shù)值 /參數(shù) 4 因子表上三角矩陣對(duì)角元素 /參數(shù) 5 因子表

48、上三角各行非零元素個(gè)數(shù) int i, j, k, i1, i2; int jj, jj1, jj2, im, x, fdNS; double ai, bNS; nfd1 = 1;實(shí)用標(biāo)準(zhǔn)文案 for (j = 1; j = n; j +) 精彩文檔 for (i = 1; i = n; i +) nobt 存放的是節(jié)點(diǎn)類(lèi)型，0: pq 節(jié)點(diǎn),-1: pv 節(jié)點(diǎn)。 if (t != 1) |(nobti != -1) & i != mpj) / -| / 1 for (j = i + 1; j = n; j +) / 1 bj = 0.0; / 1 bi = biii; / 1 if (

49、kk2 = 0) & (t = 1) & (nobti != -1)/ 存在(t = 1) 的情況，不多余。 bi = bi + af1 * qli / v0i / v0i;/af1 11 = ydzi; 12 = ydzi + 1 - 1; for (j = i1; j = i2; j +) k = iyj; bk = ybj; bmpj = 0.0; if (t = 1) if (no btj = -1) 實(shí)用標(biāo)準(zhǔn)文案 dii = 1.0 / bi; 精彩文檔 bj = 0.0； i1 = i - 1; for (im = 1; im = i1; im +) jj1 = n

50、 fdim; jj2 = nfdim + 1 - 1; for (jj = jj1; jj = jj2; jj +) if(iujj = i) ai = ujj / diim; for(k = jj; k = jj2; k +) j = iuk; bj = bj - ai * uk; break; x = n fdi;實(shí)用標(biāo)準(zhǔn)文案 n fdi+1 = n fdi + fdi; 精彩文檔 ai = dii; k = 0; i1 = i + 1; for (j = i1; j 1.0e-15) ux = bj * ai; iux = j; k+; x+; fdi = k; else fdi = 0

51、; dii = 0.0; fprin tf(fp2, n *U*); 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 for (i = 1; i = x; i +) if(i % 3 = 1) fprin tf(fp2, n); fprintf(fp2, %10.5f%8i, ui, iui); fprin tf(fp2, n * DI* *); prin tf1(di, n); void prin tf1(double* aa, int n) /* 本函數(shù)輸出 aai,i=1-n */ int i; for (i = 1; i = n; i +) if(i % 5 = 1) fprin tf(fp2, n); fp

52、rin tf(fp2, %9.5f, aai);fprin tf(fp2, nn); 實(shí)用標(biāo)準(zhǔn)文案 bi = bi * dii; 精彩文檔 void calc(i nt* iu, double* u, double* di, i nt* nfd, double* b) /* 本函數(shù)利用因子表解線(xiàn)形方程組。（P417 圖 F1-9） */ double bi; int i, j, k, i1, i2; for (i = 1; i = n; i +) / 前代過(guò)程 bi = bi; 11 = n fdi; 12 = n fdi + 1; for (j = i1; j = 1; i -) / 回代過(guò)

53、程 bi = bi; 11 = n fdi; 12 = n fdi + 1 - 1; for (j = i2; j = i1; j -) k = iuj; bi = bi - bk * uj; bi = bi; void zlsort(i nt* nn ew) Illi 。 II* 本函數(shù)進(jìn)行支路數(shù)據(jù)排序 */ II 小節(jié)點(diǎn)號(hào)放左邊，大節(jié)點(diǎn)號(hào)放右邊 II 左右皆按從小到大順序排列/ int ip, k1, k2, k3, k4; int i, j; 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 for (i = 1; i k2) izli = k4; jzli = k3; for (i = 1; i = zls -

54、 1; i +) ip = i;實(shí)用標(biāo)準(zhǔn)文案 swap(&zyki, &zykip); 精彩文檔 k1 = iabs(izli); k3 = iabs(jzli); for (j = i + 1; j = zls; j +) k2 = iabs(izlj); k4 = iabs(jzlj); if(k2 k1 | (k2 = k1 & k4 k3) ip = j; k1 = k2; k3 = k4; if(i != ip) isw ap(&i zli, & izlip); iswap(&jzli, &jzlip); swap(&z

55、ri, & zrip); swap(&zxi, &zxip); void bn sopt() 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 /* 節(jié)點(diǎn)優(yōu)化 */ int ii1, ii2, zls2, no max; int i, j, l, k1, k; int temp; zls2 = zls + zls; for (i = 1; i = zls2; i +) oldi = nnewi = 0;/ 先清零。由此可知：NS4、NS 必須大于 2*zls。 for (i = 1; i = zls; i +) oldi = iabs(izli); oldi + zls = iabs(jzli)

56、; /變壓器節(jié)點(diǎn)號(hào)由正變負(fù)，old前 zls 個(gè)為左節(jié)點(diǎn)號(hào)，后 zls 個(gè)為右節(jié)點(diǎn)號(hào)。 for (i = 1; i = zls2; i +)/ 冒泡法排序。 k1 = i + 1; for (j = k1; j oldj) isw ap(&oldi, & oldj); /交換整數(shù) oldi、oldj o小節(jié)點(diǎn)號(hào)排在支路左側(cè) 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 即是最大節(jié)點(diǎn)號(hào)。lee30.dat - 30 l = 1; for (i = 1; i = n; i +) ii1 = oldl; oldi = ii1; for (j = l; j = zls2; j +) ii2 = oldj;

57、if (ii1 != ii2) l = j; break; nn ewi +; no max = oldzls2; no max 實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 for (i = 1; i = n - 1; i +) for (j = i + 1; j = n; j +) if (nn ewj = nn ewi) if (nn ewj != nn ewi) | (oldj oldi) isw ap(&oldi, & oldj); iswap(&nn ewi, &nn ewj); for (i = 1; i = no max; i +) nn ewi = 0; for (

58、i = 1; i = n; i +) j = oldi; nn ewj = i; void yaO() /( P407 ) /* 本函數(shù)獲得 yds、ydz、列足碼 iy */ int i, j, I, II; 實(shí)用標(biāo)準(zhǔn)文案 ydz1 = 1; 精彩文檔 for (i = 1; i = n; i +) ydsi=0; / yds存放各行非零非對(duì)角元素的個(gè)數(shù)。 II = 1; for (I = 1; I = zIs; I +) i = iabs(izll); j = iabs(jzll); if (i = j) con ti nue; iyll = j; if(i != iabs(izll +

59、1) | (j != iabs(jzll + 1) II +; /II 統(tǒng)計(jì)總支路數(shù)(雙回線(xiàn)算一條支路) ydsi +; nzls = II - 1; /總支路數(shù)(雙回線(xiàn)算一條支路)實(shí)用標(biāo)準(zhǔn)文案 精彩文檔 for (i = 1; i = n - 1; i +) ydzi + 1 = ydsi + ydzi; / 由 yds得 ydz。 II ydzi是第 i 行第一個(gè)非零非對(duì)角元素的首地址， /即在所有非零非對(duì)角元素中的順序號(hào)。 prin ti(ydz, n); fprin tf(fp2, n *YDS*); prin ti(yds, n); fprin tf(fp2, nn); void p

60、rin ti(i nt* aa, int n) /* 本函數(shù)輸出 aa1-aan的值 */ int i; for(i = 1; i = n; i +) if(i % 10 = 1) fprin tf(fp2, n); fprintf(fp2, %5d, aai); fprin tf(fp2, n * YDZ* * ); 實(shí)用標(biāo)準(zhǔn)文案 ygi = 0.0; 精彩文檔 void y2() /* 本函數(shù)形成節(jié)點(diǎn)導(dǎo)納陣，一次形成，不分 intj1; double r, x, yk, zf, gij, bij; int i, j, i1, l, ll; for (i = 1; i = n ; i +) giii = 0.0; biii = 0.0; /導(dǎo)納陣對(duì)角元(與節(jié)點(diǎn)一一對(duì)應(yīng))先清零 for (i = 1; i = mdk; i +) j = idki; biij = -1.0 / dkki; /計(jì)算接地支路導(dǎo)納，只影響導(dǎo)納陣對(duì)角元 for (i = 1; i = zls; i +) B和 B */ 伯導(dǎo)納） 實(shí)用標(biāo)準(zhǔn)文案 giii =