李彩萍,李樂(lè)生
(太原理工大學(xué)電氣與動(dòng)力工程學(xué)院,山西 太原 030024)
節(jié)點(diǎn)分析法和網(wǎng)孔分析法是《電路分析基礎(chǔ)》課程中兩種主要的分析方法,可用于直流電阻電路及正弦穩(wěn)態(tài)電路的分析求解。其中,網(wǎng)孔分析法適用于網(wǎng)孔數(shù)少于節(jié)點(diǎn)數(shù)的平面電路的分析,其實(shí)質(zhì)是列電路中各網(wǎng)孔的KVL方程;節(jié)點(diǎn)分析法則適用于節(jié)點(diǎn)數(shù)少于網(wǎng)孔數(shù)的平面及非平面電路的分析求解,其實(shí)質(zhì)是列電路中各獨(dú)立節(jié)點(diǎn)的KCL方程。對(duì)平面電路來(lái)說(shuō),網(wǎng)孔分析法和節(jié)點(diǎn)分析法是兩種對(duì)偶方法,因此,本文針對(duì)節(jié)點(diǎn)分析法在授課過(guò)程中學(xué)生較難掌握以及較易出錯(cuò)的知識(shí)點(diǎn)進(jìn)行歸納總結(jié),最后給出了用MATLAB軟件對(duì)節(jié)點(diǎn)方程進(jìn)行輔助求解的過(guò)程,并舉例進(jìn)行了詳細(xì)說(shuō)明。
節(jié)點(diǎn)分析法是以節(jié)點(diǎn)電壓為未知量,聯(lián)立方程,求解各節(jié)點(diǎn)電壓值,然后進(jìn)一步求出待求量的分析方法[1]。這要求同一電路各節(jié)點(diǎn)電壓必須是一組獨(dú)立和完備的變量。以圖1直流電阻電路為例,列出其節(jié)點(diǎn)電壓方程如方程組(1)所示:
圖1 節(jié)點(diǎn)分析法示例1
式中,G11、G22分別為節(jié)點(diǎn)1、節(jié)點(diǎn)2的自電導(dǎo),節(jié)點(diǎn)的自電導(dǎo)等于與該節(jié)點(diǎn)相連的各電阻支路電導(dǎo)之和。G12、G21分別為節(jié)點(diǎn)1與節(jié)點(diǎn)2、節(jié)點(diǎn)2與節(jié)點(diǎn)1的互電導(dǎo),互電導(dǎo)是兩個(gè)相鄰節(jié)點(diǎn)公共電阻支路的電導(dǎo)。方程右端為流向各節(jié)點(diǎn)的理想電流源電流的代數(shù)和,用isii表示;其中流入節(jié)點(diǎn)電流取“+”,流出節(jié)點(diǎn)電流取“-”。因此將圖1中各元件參數(shù)代入方程組(1)可得方程組(2)。
對(duì)有N個(gè)節(jié)點(diǎn)的電路,同理可列出N-1個(gè)獨(dú)立節(jié)點(diǎn)方程。
(1)電路中如果含有理想(或受控)電壓源串聯(lián)電阻支路,則根據(jù)等效變換,可將其等效變換為電流源并聯(lián)電阻支路。電路中如果含有理想(或受控)電流源串聯(lián)電阻支路,則根據(jù)替換定理,可將理想(或受控)電流源串聯(lián)電阻支路等效變換為一個(gè)理想(或受控)電流源支路。對(duì)電流源串聯(lián)電阻支路,學(xué)生在分析時(shí)容易忽視這種等效變換,在一個(gè)方程的左端既考慮了串聯(lián)電阻,在同一個(gè)方程的右端又考慮了理想(或受控)電流源的電流。避免出現(xiàn)這種錯(cuò)誤的有效手段,除了將理想(或受控)電流源串聯(lián)電阻支路等效變換為一個(gè)理想(或受控)電流源支路方法外,還可以在理想(或受控)電流源和其串聯(lián)電阻之間增加一個(gè)節(jié)點(diǎn),將理想(或受控)電流源和其串聯(lián)電阻看成兩條支路來(lái)進(jìn)行處理。
(2)參考節(jié)點(diǎn)的選取。由于節(jié)點(diǎn)分析法實(shí)質(zhì)是列各個(gè)節(jié)點(diǎn)的KCL方程,這樣對(duì)有n個(gè)節(jié)點(diǎn)的電路,獨(dú)立節(jié)點(diǎn)數(shù)為n-1,根據(jù)n-1個(gè)獨(dú)立節(jié)點(diǎn)列出的KCL方程恰好是一組獨(dú)立的方程。因此在用結(jié)點(diǎn)分析法求解時(shí),需要首先選定其中一個(gè)節(jié)點(diǎn)作為參考節(jié)點(diǎn)。參考節(jié)點(diǎn)的選取應(yīng)該遵循以下幾個(gè)原則:①由于人們一般習(xí)慣認(rèn)為下面是地,因此在選取參考節(jié)點(diǎn)時(shí),可根據(jù)人們的習(xí)慣選電路下面的節(jié)點(diǎn)作為參考節(jié)點(diǎn)。如圖1中節(jié)點(diǎn)3恰好位于電路的下部。如果電路下方有兩個(gè)或兩個(gè)以上節(jié)點(diǎn),則選擇連接支路比較多的節(jié)點(diǎn)作為參考節(jié)點(diǎn)。②如果電路中含有理想(或受控)電壓源,通常選取電壓源的負(fù)極作為參考節(jié)點(diǎn)。這樣,電壓源正極對(duì)應(yīng)節(jié)點(diǎn)的節(jié)點(diǎn)電壓就是該電壓源電壓。因此不必再列寫(xiě)該節(jié)點(diǎn)電壓方程。部分同學(xué)易犯的錯(cuò)誤是不進(jìn)行參考節(jié)點(diǎn)的選取,給所有節(jié)點(diǎn)都列寫(xiě)節(jié)點(diǎn)方程。
(3)如果兩個(gè)獨(dú)立節(jié)點(diǎn)之間包含一個(gè)理想(或受控)電壓源支路,則在列寫(xiě)這兩個(gè)獨(dú)立節(jié)點(diǎn)的節(jié)點(diǎn)方程時(shí),需要將該理想(或受控)電壓源所在支路電流也考慮進(jìn)去。如果將該理想(或受控)電壓源所在支路電流放在方程左端,則流出節(jié)點(diǎn)電流前取“+”號(hào),否則取“-”號(hào)。如果將該理想(或受控)電壓源支路電流放在方程右端,則恰好相反。學(xué)生們易犯的錯(cuò)誤是忽視了該理想(或受控)電壓源支路流過(guò)的電流。因此授課時(shí),應(yīng)特別強(qiáng)調(diào)節(jié)點(diǎn)分析法實(shí)質(zhì)是列各個(gè)節(jié)點(diǎn)的KCL方程,與每個(gè)和獨(dú)立節(jié)點(diǎn)相連的所有支路電流都應(yīng)該考慮。
下面以圖(2)電路圖為例,進(jìn)一步說(shuō)明以上幾個(gè)難點(diǎn)。
首先,在該電路圖中含有一個(gè)1V理想電壓源串聯(lián)1個(gè)3S電導(dǎo),為簡(jiǎn)化分析,應(yīng)將其變換為一個(gè)電流源并聯(lián)電導(dǎo)的組合。
圖2 節(jié)點(diǎn)分析法示例2
圖3 節(jié)點(diǎn)分析法示例2分析圖1
圖4 節(jié)點(diǎn)分析法示例2分析圖2
其次,在該電路圖中含有一個(gè)8A理想電流源串聯(lián)1個(gè)2S電導(dǎo),此時(shí)可采用兩種方法處理,一種方法是在8A理想電流源和2S電導(dǎo)之間增加一個(gè)節(jié)點(diǎn)4,見(jiàn)圖3,則電路的獨(dú)立節(jié)點(diǎn)數(shù)為4個(gè)。另一種方法是將8A理想電流源串聯(lián)2S電導(dǎo)的組合等效變換為一個(gè)8A的電流源,見(jiàn)圖4,此時(shí)電路的獨(dú)立節(jié)點(diǎn)數(shù)為3個(gè)。比較簡(jiǎn)單的方法是采用圖4的簡(jiǎn)化電路,這樣可以減少一個(gè)節(jié)點(diǎn)方程。
另外,該電路如果選用節(jié)點(diǎn)0作為參考節(jié)點(diǎn),比較符合人們的思維習(xí)慣。但此時(shí)節(jié)點(diǎn)2和節(jié)點(diǎn)3之間有一個(gè)22 V獨(dú)立電壓源,因此需要在電路圖上標(biāo)注流經(jīng)22 V獨(dú)立電壓源電流及其參考方向。根據(jù)圖 4 列寫(xiě)的節(jié)點(diǎn)方程見(jiàn)(3)、(4)、(5)和(6)式:
如果根據(jù)圖(3)列寫(xiě)節(jié)點(diǎn)方程,則需要再增加方程(7):
MATLAB語(yǔ)言是一種廣泛應(yīng)用于工程計(jì)算及數(shù)值分析領(lǐng)域的新型高級(jí)語(yǔ)言,其強(qiáng)大的數(shù)值計(jì)算功能、圖形顯示功能、簡(jiǎn)潔的編程語(yǔ)言及可交互的集成環(huán)境,深受科研工作者及工程技術(shù)人員的歡迎。
用MATLAB進(jìn)行電路理論的輔助分析,主要包括3種方法[2]:第一種是先通過(guò)手工建立電路方程,再借助MATLAB強(qiáng)大的數(shù)值計(jì)算功能求解已建立的電路方程。第二種是利用節(jié)點(diǎn)支路關(guān)聯(lián)矩陣和支路VAR[3],通過(guò)矩陣運(yùn)算推導(dǎo)出節(jié)點(diǎn)方程,但這種方法需要利用電路的拓?fù)鋱D,因此建立節(jié)點(diǎn)方程的過(guò)程比較繁瑣。第三種方法是采用符號(hào)電路分析,可用于電路設(shè)計(jì)時(shí)確定元件的參數(shù)。下面以第一種方法為例,說(shuō)明采用MATLAB求解采用前述方法已建立的電路方程的詳細(xì)過(guò)程。
設(shè)已建立的電路方程如式(8)
其中:x是由電路中的一些電壓或和電流變量構(gòu)成的列向量;A為系數(shù)矩陣,取決于電路元件的值;B為右端列向量,其元素與電壓源的電壓和電流源的電流有關(guān)。則用MATLAB求解該方程的指令為
下面給出求解節(jié)點(diǎn)電壓方程(3)、(4)、(5)、(6)中未知變量U1、U2、U3及I的程序:%node.m;求解節(jié)點(diǎn)電壓方程的node.m程序U=zeros(4,1);%對(duì)代求變量進(jìn)行初始化,其中代求變量為U1、U2、U3、I。
A=[7-3-4 0;-3 4 0 1;-4 0 9-1;0-1 1 0];%建立系數(shù)矩陣A。
B=[-11;3;25;22];%根據(jù)方程右端參數(shù)建立矩陣 B。
U=AB
輸入程序后點(diǎn)擊運(yùn)行按鈕,則在MATLAB COMMAND WINDOW界面顯示如下結(jié)果:
則 U1=-4.5 V、U2=-15.5 V、U3=6.5 V、I=51.5 V。
用MATLAB輔助求解線性代數(shù)方程的優(yōu)勢(shì)不是很明顯,但在求解以復(fù)數(shù)為變量且方程前面的系數(shù)也是復(fù)數(shù)的節(jié)點(diǎn)電壓相量方程時(shí)的優(yōu)勢(shì)就比較顯著。以圖5為例進(jìn)行說(shuō)明。
圖5 正弦穩(wěn)態(tài)電路的節(jié)點(diǎn)分析法示例
圖6 圖5的等效相量模型
圖6選下面節(jié)點(diǎn)0作為參考節(jié)點(diǎn),則在節(jié)點(diǎn)2、3之間有一個(gè)受控電壓源,因此需要在圖中標(biāo)注出流經(jīng)該受控電壓源的電流及其參考方向。根據(jù)圖6列寫(xiě)相量形式的節(jié)點(diǎn)電壓方程:
觀察發(fā)現(xiàn),5個(gè)方程中有6個(gè)未知變量,這正是由于電路中含有受控電壓源,對(duì)含有受控電壓(電流)源的電路,一般需要根據(jù)受控電壓(電流)源控制量所在支路再增加一個(gè)方程,該電路中控制量為為流經(jīng) 2Ω 電阻電流,因此根據(jù)電阻 VAR可得式(15):
聯(lián)立(10)-(15)式可求解出各節(jié)點(diǎn)電壓的相量形式。如果直接用手工求解該復(fù)數(shù)方程組,其計(jì)算量是比較大的,且易出錯(cuò)。下面給出用MATLAB輔助求解該方程組的程序:%node1.m;求解節(jié)點(diǎn)電壓方程的node1.m程序U=zeros(6,1);%對(duì)代求變量進(jìn)行初始化,其中代求變量為
a11=1;a21=-1/20;a22=1/20+1/2+1/(-i×2);a23=1/(i×2);a24=1;a32=1/(i×2);
a33=1/2;a34=-1/2;a35=-1;a44=1;a52=1;a53=-1;a56=-5;a62=1/2;a66=-1;
A=[a11 0 0 0 0 0;a21 a22 a23 a24 0 0;0 a32 a33 a34 a35 0;0 0 0 a44 0 0;0 a52 a53 0 0 a56;0 a62 0 0 0 a66];%建立復(fù)數(shù)系數(shù)矩陣A
B=[2;0;1;4×i;0;0];%建立復(fù)數(shù)矩陣 B
U=AB
輸入程序后點(diǎn)擊運(yùn)行按鈕,則在MATLAB COMMAND WINDOW界面顯示如下結(jié)果:
節(jié)點(diǎn)分析法是電路理論中的一種重要的分析方法,本文首先介紹了節(jié)點(diǎn)分析法的原理,然后分析了用節(jié)點(diǎn)分析法分析電路的重點(diǎn)和難點(diǎn),最后通過(guò)示例給出了用MATLAB輔助求解節(jié)點(diǎn)電壓方程的詳細(xì)過(guò)程。
[1]史健芳,陳惠英.電路基礎(chǔ).北京:人民郵電出版社,2006.
[2]趙錄懷,楊育霞,張震.電路與系統(tǒng)分析-使用MATLAB.北京:高等教育出版社,2004.
[3]于舒娟,史學(xué)軍.Matlab應(yīng)用于大規(guī)模電網(wǎng)絡(luò)的節(jié)點(diǎn)分析法.現(xiàn)代電子技術(shù),2009:17(304),P121-123.