劉雄峰，姚思遠(yuǎn)

（三峽大學(xué)電氣與新能源學(xué)院，湖北宜昌，443002）

0 引言

在電路課程中經(jīng)常涉及到線性方程組的求解，具體的方法有精確法和迭代法，精確法求解不需要采取近似舍入，而是采用初等變換方法求出方程組的解；迭代法則是通過(guò)有限次的迭代，在允許的精度范圍內(nèi)求解方程組的近似解，精度要求設(shè)定越高，求解值越趨近與真實(shí)值。

1 MATLAB軟件和迭代算法簡(jiǎn)介

MATLAB是美國(guó)mathworks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開(kāi)發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級(jí)技術(shù)計(jì)算語(yǔ)言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分；MATLAB具有高效的數(shù)值計(jì)算及符號(hào)計(jì)算功能，能使用戶從繁雜的數(shù)學(xué)運(yùn)算分析中解脫出來(lái)。

迭代法也稱為轉(zhuǎn)輾法，是一種不斷用變量舊值遞推新值的過(guò)程，迭代算法是用計(jì)算機(jī)解決數(shù)學(xué)問(wèn)題的一種基本方法，它利用計(jì)算機(jī)運(yùn)算速度快、適合做重復(fù)性操作的特點(diǎn)，讓計(jì)算機(jī)對(duì)一組指令（或一定步驟）進(jìn)行重復(fù)執(zhí)行，在每次執(zhí)行這組指令或步驟時(shí)，都將從變量的原值推出它的一個(gè)新值[1]。

設(shè)有線性方程組Ax=b，其中

選取x0作為方程的初始近似解，以之迭代出一個(gè)方程組近似解收斂數(shù)組(xk),k≥0。[2]即可將式（1）變?yōu)椋?/p>

根據(jù)初始近似解x0，通過(guò)如下的迭代過(guò)程可以產(chǎn)生迭代根數(shù)組。

Pxk+1=Qxk+b,k=0,1,...，當(dāng)‖P-Q‖＜0時(shí)數(shù)組收斂。

線性方程組的解滿足：

將矩陣A作標(biāo)準(zhǔn)分解可得到：A=D+L+U

其中，D=diag(a11,a22,...,amm)

此時(shí)可將式（1）變?yōu)椋?/p>

此時(shí)，應(yīng)用雅克比迭代法可以得到：

當(dāng)應(yīng)用高斯-賽德?tīng)柕傻茫?/p>

式中 ：i=1,...,m；k=1,2,...

計(jì)算xi與k-1次迭代值的加權(quán)平均ω作為第k次的迭代值，即：

上式也可整理成：

其中ω稱為松弛因子，要求0＜ω＜2。當(dāng)ω＞1時(shí)，上式稱為逐次超松弛迭代法；當(dāng)ω=1時(shí)，上式即為高斯-賽德?tīng)柕?；?dāng)0＜ω＜1時(shí)，上式稱為低松弛迭代法[2]。

2 應(yīng)用案例

如圖1所示的電路圖中，已知Ue1=8 V,Ue2=4 V,Ue3=2 V,R1=R2=R4=1 Ω,R3=3 Ω,R5=2 Ω,R6=7Ω，利用迭代法求解Il1Il2Il3。

圖1 電路圖

本文將針對(duì)該應(yīng)用案例使用MATLAB軟件實(shí)現(xiàn)雅克比迭代法、高斯-賽德?tīng)柕ê椭鸫纬沙诘ㄔ谇蠼饩€性電路方程中的應(yīng)用[3]，并對(duì)運(yùn)行結(jié)果進(jìn)行分析。（計(jì)算精度控制在10-5）

根據(jù)電路等效分析可得：

根據(jù)戴維南定理可得：

代入數(shù)值：

上述方程組對(duì)應(yīng)于（1）式中的A、b矩陣：

3 迭代法的MATLAB實(shí)現(xiàn)

3.1 雅克比迭代法的MATLAB實(shí)現(xiàn)

Step1：編寫雅克比迭代法的程序

Step2：編寫驗(yàn)證程序

Step3：運(yùn)行驗(yàn)證程序jacobi.m，得到在10-5精度范圍內(nèi)的Il1Il2Il3值和計(jì)算機(jī)迭代次數(shù)如表1所示。

表 1

3.2 高斯-賽德?tīng)柕ǖ腗ATLAB實(shí)現(xiàn)

Step1：編寫高斯-賽德?tīng)柕ǖ某绦騕4]

Step2：編寫驗(yàn)證程序

Step3：運(yùn)行驗(yàn)證程序gauss.m，得到在10-5精度范圍內(nèi)的Il1Il2Il3值和計(jì)算機(jī)迭代次數(shù)如表2所示。

表2

3.3 逐次超松弛迭代法的MATLAB實(shí)現(xiàn)

Step1：編寫逐次超松弛迭代法的程序

Step2：編寫驗(yàn)證程序

Step3：運(yùn)行驗(yàn)證程序sor.m，得到在10-5精度范圍內(nèi)的Il1Il2Il3值和計(jì)算機(jī)迭代次數(shù)如表3所示。

表3

說(shuō)明：這里的最佳松弛因子w為1.15，可以通過(guò)w取不同值進(jìn)行驗(yàn)證。

4 總結(jié)

由表1-表3可以看出三種迭代法解線性方程組的收斂速度，三種迭代法達(dá)到要求的精度（例題要求精度：10-5）所需的迭代次數(shù)如下表所示：

表4

分析表4可以得到逐次超松弛迭代法的收斂速度最快，高斯-賽德?tīng)柕ù沃?，雅克比迭代法最慢，通過(guò)上面的應(yīng)用案例可以知道，高斯-賽德?tīng)柕煽醋魇茄趴吮鹊ǖ囊环N改進(jìn)，逐次超松弛迭代法是高斯-賽德?tīng)柕ǖ囊环N修正。鑒于程序的可移植性，在求解更為復(fù)雜的電路時(shí)只需要修改程序里的數(shù)據(jù)參數(shù)即可，電路越復(fù)雜越能體現(xiàn)計(jì)算機(jī)的高效。

本文使用了MATLAB軟件使電路分析的計(jì)算過(guò)程得到了極大的簡(jiǎn)化，同時(shí)也為廣大學(xué)生提供了一種新的解題思路，通過(guò)計(jì)算機(jī)高效的特點(diǎn)，激發(fā)學(xué)生利用編程軟件學(xué)習(xí)專業(yè)知識(shí)的興趣，具有一定的參考價(jià)值。