單機系統(tǒng)潮流的經(jīng)濟分布控制 (1)

1、I遼遼 寧寧 工工 業(yè)業(yè) 大大 學學電力系統(tǒng)計算電力系統(tǒng)計算課程設(shè)計（論文）課程設(shè)計（論文）題目：題目： 單機系統(tǒng)潮流的經(jīng)濟分布控制（單機系統(tǒng)潮流的經(jīng)濟分布控制（4） 院（系）：院（系）： 電電 氣氣 工工 程程 學學 院院 專業(yè)班級：專業(yè)班級： 電氣電氣 082 學學 號：號： 學生姓名：學生姓名： 指導(dǎo)教師：指導(dǎo)教師： 教師職稱：教師職稱： 起止時間：起止時間：11-06-2711-06-27 至至 11-07-0811-07-08本科生課程設(shè)計（論文）II課程設(shè)計（論文）任務(wù)及評語課程設(shè)計（論文）任務(wù)及評語院（系）：電氣工程學院 教研室：電氣工程及其自動化課程設(shè)計（論文）任務(wù)原始資料：系

2、統(tǒng)如圖G 2 L24 4 S4 L12 L14 1 S1 各元件參數(shù)標幺值如下：L14：電阻 0.1，電抗 0.35，對地容納 0.04；L12: 電阻 0.1，電抗 0.28，對地容納 0.03；L24: 電阻 0.2，電抗 0.4，對地容納 0.035；給定功率及電壓如下：S1=0.3+j0.1， S4=0.35+j0.1， V2=1.05.任務(wù)要求：1 對系統(tǒng)進行潮流計算，計算支路功率分布；2 計算在 L14 支路中應(yīng)加入的附加電勢，實現(xiàn)支路功率經(jīng)濟分布；3 設(shè)計串聯(lián)補償器硬件電路；4 采用適當?shù)乃惴?、編寫程序，實現(xiàn)支路功率的經(jīng)濟分布。指導(dǎo)教師評語及成績平時考核： 設(shè)計質(zhì)量： 論文格式：

3、 總成績： 指導(dǎo)教師簽字： 年 月 日本科生課程設(shè)計（論文）III摘 要本課程設(shè)計主要研究某一電力系統(tǒng)網(wǎng)絡(luò)，對此網(wǎng)絡(luò)進行潮流計算。電力系統(tǒng)的潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。本文研究此電力網(wǎng)的支路功率分布，并完成對一支路加入附加電勢，使系統(tǒng)實現(xiàn)支路功率經(jīng)濟分布的效果。在電力網(wǎng)中會出現(xiàn)功率損失，為是電能質(zhì)量可靠，保證電能的質(zhì)量要求，則需對系統(tǒng)進行功率補償，本此設(shè)計對串聯(lián)補償器進行設(shè)計，并設(shè)計其硬件電路圖，以滿足電力用戶需要的電能質(zhì)量。并通過牛頓拉夫遜潮流計算使系統(tǒng)實現(xiàn)支路功率的經(jīng)濟分布，并用MATLAB 編寫程序。關(guān)鍵詞：潮流計算；經(jīng)

4、濟分布；補償器； 本科生課程設(shè)計（論文）IV目 錄第 1 章 緒論 .11.1 潮流計算概況 .11.2 本文研究內(nèi)容 .1第 2 章 潮流計算 .22.1 潮流計算 .22.1.1 潮流計算及其作用.22.1.2 系統(tǒng)的潮流計算.32.2 支路功率分布計算 .4第 3 章 功率經(jīng)濟分布 .73.1 功率經(jīng)濟分布 .73.2 系統(tǒng)的附加電勢計算 .7第 4 章 串聯(lián)補償器設(shè)計 .94.1 串聯(lián)補償作用 .94.2 電壓和電流的采集 .10第 5 章 牛頓拉夫遜潮流計算 .125.1 潮流計算的基本原理 .125.2 牛頓拉夫遜潮流計算 .125.2.1 功率方程.125.2.2 雅可比矩陣各元

5、素的表達式.125.3 編寫程序 .13第 6 章 課程設(shè)計總結(jié) .17參考文獻 .18本科生課程設(shè)計（論文）1第 1 章 緒論1.1 潮流計算概況電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，它根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)：各母線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃的設(shè)計和現(xiàn)有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性.可靠性和經(jīng)濟性。此外，電力系統(tǒng)潮流計算也是計算系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。所以潮流計算是研究電力系統(tǒng)的一種很重要和很基礎(chǔ)的計算。電力系統(tǒng)潮流計算也分為離線計算和

6、在線計算兩種，前者主要用于系統(tǒng)規(guī)劃設(shè)計和安排系統(tǒng)的運行方式，后者則用于正在運行系統(tǒng)的經(jīng)常監(jiān)視及實時控制。電力系統(tǒng)潮流計算問題在數(shù)學上是一組多元非線性方程式求解問題，其解法都離不開迭代。因此，對潮流計算方法，首先要求它能可靠地收斂，并給出正確答案。由于電力系統(tǒng)結(jié)構(gòu)及參數(shù)的一些特點，并且隨著電力系統(tǒng)不斷擴大，潮流問題的方程式階數(shù)越來越高，對這樣的方程式并不是任何數(shù)學方法都能保證給出正確答案的。這種情況成為促使電力系統(tǒng)計算人員不斷尋求新的更可靠方法的重要因素。1.2 本文研究內(nèi)容本文通過對電力系統(tǒng)網(wǎng)絡(luò)的潮流計算，究此電力網(wǎng)的支路功率分布，并完成對一支路加入附加電勢，使系統(tǒng)實現(xiàn)支路功率經(jīng)濟分布的效果。

7、在電力網(wǎng)中會出現(xiàn)功率損失，則需對系統(tǒng)進行功率補償，本此設(shè)計對串聯(lián)補償器進行設(shè)計，并設(shè)計其硬件電路圖，并通過牛頓拉夫遜潮流計算使系統(tǒng)實現(xiàn)支路功率的經(jīng)濟分布，并用 MATLAB 編寫程序。本科生課程設(shè)計（論文）2本科生課程設(shè)計（論文）3第 2 章 潮流計算2.1 潮流計算2.1.1 潮流計算及其作用潮流計算是電力系統(tǒng)分析中的一種最基本的計算，它的任務(wù)是對給定的運行條件確定系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角） 、網(wǎng)絡(luò)中的功率分布及功率損耗等。電力系統(tǒng)潮流計算的結(jié)果是電力系統(tǒng)穩(wěn)定計算和故障分析的基礎(chǔ)。潮流計算根據(jù)給定的電網(wǎng)結(jié)構(gòu)、參數(shù)和發(fā)電機、負荷等元件的運行條件，確定電力系統(tǒng)各部分穩(wěn)態(tài)運行狀

8、態(tài)參數(shù)的計算。通常給定的運行條件有系統(tǒng)中各電源和負荷點的功率、樞紐點電壓、平衡點的電壓和相位角。待求的運行狀態(tài)參量包括電網(wǎng)各母線節(jié)點的電壓幅值和相角，以及各支路的功率分布、網(wǎng)絡(luò)的功率損耗等。潮流計算是電力系統(tǒng)非常重要的分析計算，用以研究系統(tǒng)規(guī)劃和運行中提出的各種問題。對規(guī)劃中的電力系統(tǒng)，通過潮流計算可以檢驗所提出的電力系統(tǒng)規(guī)劃方案能否滿足各種運行方式的要求；對運行的電力系統(tǒng)，通過潮流計算可以預(yù)知各種負荷變化和網(wǎng)絡(luò)結(jié)構(gòu)的改變會不會危及系統(tǒng)的安全，系統(tǒng)中所有母線的電壓是否在允許的范圍以內(nèi)，系統(tǒng)中各種元件（線路、變壓器等）是否會出現(xiàn)過負荷，以及可能出現(xiàn)過負荷時應(yīng)事先采取哪些預(yù)防措施等。系統(tǒng)電路圖如圖

9、 2.1 所示G 2 L24 4 S4 L12 L14 1 S1 本科生課程設(shè)計（論文）4圖 2.1 系統(tǒng)電網(wǎng)2.1.2 系統(tǒng)的潮流計算潮流計算所涉及的公式：充電功率為： =1/2* （2-BQCB2NV1）線路的功率損耗為： =/* （2-LS2S2NUZ2）線路中電壓降降落的縱、橫分量分別為：=()/（2-3）LVLLPRQXAV =()/ （2-LVLLPXQRAV4）B 點電壓為： = （2-BV22()ALLVVV5）根據(jù)系統(tǒng)圖知： =0.1+j0.2812LZ=0.1+j0.3514LZ=0.2+j0.424LZ=0.3+j0.11S=0.35+j0.14S=1.052V系統(tǒng)的等效

10、電路圖如圖 2.2 所示：圖 2.2 系統(tǒng)的等效電路圖線路的充電功率為：= 0.5*0.03*1.05*1.05= 0.016512BQ本科生課程設(shè)計（論文）5= 0.5*0.04*1.05*1.05= 0.022114BQ= 0.5*0.035*1.05*1.05= 0.019324BQb 節(jié)點運算負荷為：=+=0.3+j0.1-j0.0165-j0.0221=0.3+j0.0614bS1S12BQ14BQc 節(jié)點運算負荷為：=+=0.35+j0.1-j0.0221-j0.0193=0.35+j0.0586cS4S14BQ24BQ線路損耗為：14L=0.0114+j0.0414LS2220.

11、350.0586(0.10.35)1.05j3 點的運算負荷為：=+=0.3614+j0.09843ScS14LS2 點的運算負荷為：=+=0.3614+j0.05982S3S14BQ12BQ線路損耗為：12L=0.0122+j0.034112LS2220.36140.0598(0.10.28)1.05j1 點的運算負荷為：=+=0.3736+j0.09391S2S12LSA 點運算負荷為：=+=0.3736+j0.0581AS1S12BQ24BQ線路中電壓降落的縱、橫分量分別為：12L=0.0511LV0.3736*0.1 0.0581*0.281.05=0.0941LV0.3736*0.2

12、80.0581*0.11.05b 點電壓為：=1.0033bV22(1.050.0511)0.0941線路中電壓降落的縱、橫分量分別為：14L=0.0420LV0.3614*0.1 0.0598*0.351.0033=0.1201LV0.3614*0.350.0598*0.11.0033c 點電壓為：=0.9688cV22(1.00330.0420)0.1201本科生課程設(shè)計（論文）62.2 支路功率分布計算根據(jù)系統(tǒng)提供的數(shù)據(jù)：=0.1+j0.2812LZ=0.1+j0.3514LZ=0.2+j0.424LZ=0.3+j0.11S=0.35+j0.14S=1.052V系統(tǒng)的等效電路圖如圖 2.

13、3 所示圖 2.3 系統(tǒng)的等效電路圖線路的充電功率為：= 0.5*0.03*1.05*1.05= 0.016512BQ= 0.5*0.04*1.05*1.05= 0.022114BQ= 0.5*0.035*1.05*1.05= 0.019324BQb 節(jié)點運算負荷為：=+=0.3+j0.1-j0.0165-j0.0221=0.3+j0.0614bS1S12BQ14BQc 節(jié)點運算負荷為：=+=0.35+j0.1-j0.0221-j0.0193=0.35+j0.0586cS4S14BQ24BQ網(wǎng)絡(luò)中的功率分布為：=12LS142424121424()bLLcLLLLS ZZS ZZZZ(0.30

14、.0614)*(0.10.350.20.4)(0.350.0586)*(0.20.4)0.10.280.10.350.20.4jjjjjjjj=0.3630+j0.0531=24LS121412121424()cLLbLLLLS ZZS ZZZZ本科生課程設(shè)計（論文）7(0.350.0586)*(0.10.280.20.35)(0.30.0614)*(0.10.28)0.10.280.10.350.20.4jjjjjjjj=0.2870+j0.0669驗算：+=0.65+j0.1212LS24LS +=065+j0.12bScS可見，計算結(jié)果誤差非常小，無需重算取=0.2870+j0.0669

15、24LS=-=0.35+j0.0586-(0.2870+j0.0669)=0.036-j0.008314LScS24LS其功率分布圖如圖 2.4 所示 圖 2.4 功率分布圖本科生課程設(shè)計（論文）8本科生課程設(shè)計（論文）9第 3 章 功率經(jīng)濟分布3.1 功率經(jīng)濟分布在環(huán)網(wǎng)中引入環(huán)路電勢使產(chǎn)生循環(huán)功率，是對環(huán)網(wǎng)進行潮流控制和改善功率分布的有效手段。功率在環(huán)形網(wǎng)絡(luò)中是與阻抗成反比分布的，這種分布稱為功率的自然分布。經(jīng)濟分布：功率在環(huán)式網(wǎng)中是按與電阻成反比分布時，其有功功率損耗最小。只有每段電力線路阻抗的 R/X 比值都相等的均一網(wǎng)絡(luò)中，功率的自然分布才與經(jīng)濟分布相符。在一般情況下，這兩者是有差別的

16、。各段線路的不均一程度越大，功率損耗的差別就越大。如果在環(huán)網(wǎng)中引入附加電勢，假定其產(chǎn)生與同方向的循環(huán)功率，且滿E1S足條件 += （3-1ScirS1ecS1）就可以使功率分布符合經(jīng)濟分布的要求。由此可得所要求的循環(huán)功率為 =（）+j()=+j （3-cirS1ecS1S1ecP1P1ecQ1QcirPcirQ2）為產(chǎn)生此循環(huán)功率所需的附加電勢則為 =*/=+j=+j （3-EZcirSNVcircirNP RQ XVcircirNP XQ RVXEYE3）式中，為環(huán)網(wǎng)的總阻抗，為網(wǎng)絡(luò)的額定電壓。ZNV3.2 系統(tǒng)的附加電勢計算系統(tǒng)電網(wǎng)如圖3.1 所示G 2 L24 4 S4 L12 L14

17、1 S1 本科生課程設(shè)計（論文）10圖 3.1 系統(tǒng)電網(wǎng)功率分布為：=0.3630+j0.053112LS142424121424()bLLcLLLLS ZZS ZZZZ=0.2870+j0.066924LS121412121424()cLLbLLLLS ZZS ZZZZ因為：=0.3+j0.0614bS所以：=0.3bP = j0.0614bQ功率損耗為：= =0.31251ecP232123()bcP RRPRRRR0.3*(0.1 0.2)0.35*0.10.1 0.1 0.2= =j0.06071ecQ232123()bcQ RRQ RRRR0.0614*(0.1 0.2)0.0586

18、*0.10.1 0.1 0.2循環(huán)功率為：=（）+j()cirS1ecS1S1ecP1P1ecQ1Q=(0.3125-0.2870)+j(0.0607-0.0669)=0.0255-j0.0012cirS附加電勢為：=*/=EZcirSNV2(0.41.03)*(0.02550.0012)1.05jj=0.0104+j0.0234E本科生課程設(shè)計（論文）11第 4 章 串聯(lián)補償器設(shè)計4.1 串聯(lián)補償作用在電力系統(tǒng)中，無功功率既有其有益的一面，也有其有害的一面。有益的是，電力網(wǎng)在運行時，電源供給的無功功率是電能轉(zhuǎn)換為其它形式能量的前提。它為電能的輸送、轉(zhuǎn)換創(chuàng)造了條件。而有害的是，在電力系統(tǒng)輸送電

19、能的過程中，無功功率不足，將使系統(tǒng)中輸送的總電流增加、使變壓器的出力減少、供電線路及系統(tǒng)設(shè)備有功功率損耗增大、線路末端電壓下降；同時，對于電力用戶來說，過多的從電力網(wǎng)中吸取無功，不僅使電網(wǎng)電能質(zhì)量下降，也影響自身的用電和生產(chǎn)。因此，為了減少無功的損失和避免其在電網(wǎng)中的不當流動，必須進行無功補償。串聯(lián)補償可改變傳輸線路的等效阻容或在線路中串入補償電壓，因此，通過串聯(lián)補償可方便地調(diào)節(jié)系統(tǒng)的有功和無功潮流，從而有效控制電力系統(tǒng)的電壓水平和功率平衡。在配電網(wǎng)中，大部分的電壓降是由線路電感引起的，串聯(lián)補償可使負載波動引起的電壓下降減少。由于串聯(lián)電容其實是一個無源的電路元件，其響應(yīng)自發(fā)而迅速，因此，串聯(lián)補

20、償有助于電壓調(diào)節(jié)，可較完美地解決電壓閃變問題。4.2 串聯(lián)補償電路的硬件結(jié)構(gòu)該裝置核心硬件采用 AT89C51 單片機實現(xiàn)，AT89C51 單片機是現(xiàn)在人們用的最多的單片機種類之一，其結(jié)構(gòu)簡單、成本低廉、功能較全且便于擴展，雖然處理速度上趕不上現(xiàn)在大多數(shù)的微處理器，如 DSP、嵌入式微處理器 ARM9 等，但是在此設(shè)計中完全可以滿足和實現(xiàn)全自動補償裝置的要求。實時參數(shù)分別采集電網(wǎng)側(cè)電壓、電流，經(jīng)過電壓形成電路分別對電壓電流進行變換。同時由 CPU 中的定時器定時控制對經(jīng)過低通濾波出來的信號進行采樣并保持，再通過多路轉(zhuǎn)換開關(guān)，分別經(jīng)高速 A/D 轉(zhuǎn)換單元進行模數(shù)轉(zhuǎn)換，并最后將各種數(shù)據(jù)實時顯示在數(shù)

21、碼管上，同時根據(jù)存儲器在 RAM 中的功率因素表判定該投切哪幾組電容器，并最終通過 CPU 驅(qū)動 SSR 固態(tài)繼電器實現(xiàn)電容器的自動投切。裝置硬件框圖如圖 4.1 所示本科生課程設(shè)計（論文）12圖4.1 裝置硬件框圖4.2 電壓和電流的采集電網(wǎng)側(cè)的電壓和電流分別經(jīng)電壓互感器和電流互感器再到電壓形成電路產(chǎn)生微機可以利用的信號。其中交流電壓信號采用電壓變換器即可滿足要求，而交流電流信號可以采用電抗變換器和電流變換器，在此選用電流變換器對電流信號進行變換。此外，為保證后面的電路不受干擾，兩個經(jīng)變換器出來的信號還需經(jīng)低通濾波環(huán)節(jié)進行濾波才允許輸入到采樣保持電路，供 CPU 發(fā)出信號進行采樣。低通濾波環(huán)

22、節(jié)可由電阻及電容組合構(gòu)成一模擬低通濾波器。而采樣保持環(huán)節(jié)則選擇較常用的 LF-398 采樣芯片。其邏輯輸入端 S/H 由 CPU 中的定時器按一定時序控制，從而控制何時采樣、何時保持。當邏輯輸入端即 S/H=1，即電平為高時，為采樣功能；當邏輯輸入端即 S/H=0，即電平為低時，為保持功能。低通濾波及采樣保持電路如圖 4.2 所示本科生課程設(shè)計（論文）13圖 4.2 低通濾波及采樣保持電路經(jīng)采樣保持電路出來的信號需通過多路選擇開關(guān)（MPX）來選擇哪一個信號需要 A/D 轉(zhuǎn)換。當電壓信號需 A/D 轉(zhuǎn)換時，MPX 選擇到電壓信號檔進而通過A/D 芯片進行模數(shù)轉(zhuǎn)換，當電流信號需 A/D 轉(zhuǎn)換時，M

23、PX 選擇到電流信號檔進而通過 A/D 芯片進行模數(shù)轉(zhuǎn)換。在此設(shè)計中，MPX 選擇 CD4051，而 A/D 轉(zhuǎn)換芯片則選擇 MAX521。二者與單片機連接的電路原理圖如圖 4.3 所示。圖 4.3 MPX 和 A/D 轉(zhuǎn)換芯片與單片機連接的電路原理圖本科生課程設(shè)計（論文）14第 5 章 牛頓拉夫遜潮流計算5.1 潮流計算的基本原理電力系統(tǒng)潮流計算是電力系統(tǒng)運行和規(guī)劃中最基本和最經(jīng)常的計算，其任務(wù)是在已知某些運行參數(shù)的情況下，計算出系統(tǒng)中全部的運行參數(shù)，一般來說，各個母線所提供負荷的功率是已知的，各個節(jié)點電壓是未知的（平衡節(jié)點除外） ，可以根據(jù)網(wǎng)絡(luò)結(jié)構(gòu)形成節(jié)點導(dǎo)納矩陣，然后由節(jié)點導(dǎo)納矩陣和網(wǎng)絡(luò)

24、拓撲結(jié)構(gòu)寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉(zhuǎn)化為求解非線性方程組的問題了。為了便于用迭代法解方程組，需要將上述功率方程改寫為成功率平衡方程，并對功率平衡方程求偏導(dǎo)，得出對應(yīng)的雅可比矩陣，給未知節(jié)點賦電壓初值，一般為額定電壓，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點電壓不平衡量（未知的）構(gòu)成了誤差方程，解誤差方程，得到節(jié)點電壓不平衡量，節(jié)點電壓不平衡量構(gòu)成新的節(jié)點電壓初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。5.

25、2 牛頓拉夫遜潮流計算5.2.1 功率方程=0 (5-1)1P11111111()()gljjjjjjjjj ij iPPeG eB ffG fB e=0 (5-2)1Q1111111()()glijjjjjjjjj ij iQQfG eB feG fB e=0 (5-3)222222222()()gljjjjjjjjj ij iPPPeG eBffGfB e=0 (5-4)222222222()()gljjjjjjjjj ij iQQQfG eBfeGfB e =0 (5-5)333333333()()gljjjjjjjjj ij iPPPeG eB ffGfB e =0 (5-6)3333

26、33333()()gljjjjjjjjj ij iQQQfG eB feGfB e5.2.2 雅可比矩陣各元素的表達式當 ji 時：本科生課程設(shè)計（論文）15 (5-7)iijijjijijPHB eG ff (5-8)iijijjijijPNG eB fe (5-9)iijijiijjjQMB fG ef (5-10)iijijiijjjQLG fB ee 當 j=i 時：- (5-11)iijijjijijPHB eG ffib- (5-12)iijijjijijPNG eB fe ia- (5-13)iijijiijjjQMB fG efia+ (5-14)iijijiijjjQLG f

27、B ee ib其中， = (5-15)ia1()nijjijiiG eB f = (5-16)ib1()nijiijjiG fB e5.3 編寫程序e=1.05 1.05 1.05f=0 0 0=0.3;2p=0.0614;2q=0.35;3p=0.0586;3qG=0,0,0;B=-0.0165,-0.0221,-0.0193;K=0;Maxp=0;Maxq=0I=0,0,0;本科生課程設(shè)計（論文）16for v=1:15for n=1:3;(1,1)(1,1)(1, )* ( )(1, )*( )IIGne nBnf n;(1,2)(1,2)(1, )*( )(1, )* ( )IIGnf

28、 nBne nendfor n=1:3;(2,1)(2,1)(2, )* ( )(2, )*( )IIGne nBnf n;(2,2)(2,2)(2, )*( )(2, )* ( )IIGnf nBne nendfor n=1:3;(3,1)(3,1)(3, )* ( )(3, )*( )IIGne nBnf n;(3,2)(3,2)(3, )*( )(3, )* ( )IIGnf nBne nendH=;N=;M=;L=;J=;20.3(2)* (2,1)(2)* (2,2)PeIfI;20.0614(2)* (2,1)(2)* (2,2)QfIeI;30.35(3)* (3,1)(3)*

29、(3,2)PeIfI;30.0531(3)* (3,1)(3)* (3,2)QfIeIfor m=2:3for n=2:3if(m=n) ;H(m,m)=B(m,m)* ( )( ,)*( )( ,2)e mG m mf mI m ;N(m,m)=-G(m,m)* ( )( ,)*( )( ,1)e mB m mf mI m ;M(m,m)=G(m,m)* ( )( ,)*( )( ,1)e mB m mf mI m ;L(m,m)=B(m,m)* ( )( ,)*( )( ,2)e mG m mf mI melse ;( , )( , )* ( )( , )*( )H m nB m ne m

30、G m nf m ;( , )( , )* ( )( , )*( )N m nG m ne mB m nf m ;( , )( , )M m nN m n ;( , )( , )L m nH m n本科生課程設(shè)計（論文）17endendendJ=H(2,2),N(2,2),H(2,3),N(2,3),;M(2,2),L(2,2),M(2,3),L(2,3);H(3,2),N(3,2),H(3,3)N(3,3);M(3,2),L(3,2),M(3,3),L(3,3);A=;C=P2;Q2;P3;Q3A=-inv(J)*Cmaxp=C(1);maxq=C(2);for d=1:2:8if(C(d)maxp) maxp=C(d);endendfor t=2:2:8if(C(t)maxq) maxq=C(t);endendif(maxp0.00001&maxq0.00001)e(2)=e(2)+A(2)e(3)=e(3)+A(3)f(2)=f(2)+A(1)f(3)=f(3)+A(3)endendS1=(1.05+j*0)*(G(1,1)+j*B(1,1)*(1.05+j*0)+ (G(1,2)