matlab在電路分析中的應(yīng)用研究10036new

《matlab在電路分析中的應(yīng)用研究10036new》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。

1、第27卷第1期沈陽(yáng)師范大學(xué)學(xué)報(bào)(自然科學(xué)版)Vol127,No.12009年1月JournalofShenyangNormalUniversity(NaturalScience)Jan.2009文章編號(hào):1673-5862(2009)01-0059-04MATLAB在電路分析中的應(yīng)用研究1222張亞琴,趙宏音,吳建軍,陳嵐峰(1.沈陽(yáng)師范大學(xué)實(shí)驗(yàn)中心,遼寧沈陽(yáng)110034;2.沈陽(yáng)師范大學(xué)物理科學(xué)與技術(shù)學(xué)院,遼寧沈陽(yáng)110034)摘要:MATLAB具有強(qiáng)大的矩陣運(yùn)算能力、簡(jiǎn)便的繪圖功能、豐富的算法工具箱,以及高效、可視化、推理能力強(qiáng)等特點(diǎn),

2、是目前理論研究和工程應(yīng)用中流行最廣的科學(xué)計(jì)算語(yǔ)言.結(jié)合電路分析的具體實(shí)例,介紹了應(yīng)用MATLAB語(yǔ)言編程對(duì)電路進(jìn)行輔助計(jì)算和輔助分析的方法.將MATLAB軟件引入到電路分析中,大大地提高了計(jì)算精度和工作效率,為電路分析提供了一個(gè)有效的輔助工具.關(guān)鍵詞:MATLAB;電路分析;輔助計(jì)算;輔助分析中圖分類號(hào):TP312文獻(xiàn)標(biāo)識(shí)碼:A0引言電路理論系統(tǒng)地介紹了電路的基本理論、概念、定理和定律,并討論各種電路的分析與計(jì)算方法,內(nèi)容多,涉及范圍廣.而電路的分析計(jì)算問(wèn)題往往最終歸結(jié)為有關(guān)數(shù)學(xué)模型的求解,如代數(shù)方程、微分方程、復(fù)數(shù)代數(shù)方程、矩陣運(yùn)算、拉普

3、拉斯變換等的數(shù)值計(jì)算和符號(hào)計(jì)算問(wèn)題,還涉及大量的時(shí)域和頻域的繪圖工作.手工求解相當(dāng)繁瑣,不僅費(fèi)時(shí)而且容易出錯(cuò).MATLAB具有強(qiáng)大的矩陣運(yùn)算能力、簡(jiǎn)便的繪圖功能、豐富的算法工具箱、高效、可視化及推理能力強(qiáng)等特點(diǎn),是目前理論研究和工程應(yīng)用中流行最廣的科學(xué)計(jì)算語(yǔ)言,已成為科研和工程技術(shù)人員的有力開(kāi)發(fā)工具.在這極具挑戰(zhàn)的時(shí)代,將計(jì)算機(jī)充分運(yùn)用到實(shí)際的工程計(jì)算,MATLAB軟件引入到電路分析中進(jìn)行輔助分析與計(jì)算,具有極其重要的意義.下面以一些典型電路為例,研究MATLAB軟件在電路分析中的具體運(yùn)用.1MATLAB輔助計(jì)算1.1方程組的求解穩(wěn)態(tài)電路的

4、核心是基爾霍夫方程,n個(gè)節(jié)點(diǎn)(或回路)就有n個(gè)代數(shù)方程,交流電路則是復(fù)數(shù)代數(shù)方程.這類方程組常用矩陣求解,對(duì)于四階以上的矩陣方程,在計(jì)算矩陣的行列式時(shí),就要展開(kāi)為代數(shù)余子式之和,展開(kāi)式很長(zhǎng),階數(shù)越高式子越龐大,手工操作應(yīng)該說(shuō)理論上可行,而實(shí)際上幾乎不可能,復(fù)數(shù)代數(shù)方程更是如此.而借助于MATLAB軟件,只需幾條語(yǔ)句問(wèn)題就得以解決.下面以簡(jiǎn)單例子來(lái)說(shuō)明問(wèn)題,復(fù)雜問(wèn)題可依此類推.例1求解代數(shù)方程組1.2u1-u2-0.1u4=1-u1+2.5u2-0.5u3=-0.5-0.5u2+1.25u3-0.25u4=0.5-0.1u1-0.25u3+0

5、.6u4=0解MATLAB程序:收稿日期:2008209226基金項(xiàng)目:遼寧省教育科學(xué)“十一五”規(guī)劃項(xiàng)目(職35223).作者簡(jiǎn)介:張亞琴(1962-),女,遼寧沈陽(yáng)人,沈陽(yáng)師范大學(xué)高級(jí)實(shí)驗(yàn)師.60沈陽(yáng)師范大學(xué)學(xué)報(bào)(自然科學(xué)版)第27卷a=[1.2-10-0.1;-12.5-0.50;0-0.51.25-0.25;-0.10-0.250.6];b=[1-0.50.50]′;mformatbank;mu=ab輸出結(jié)果為:u=1.230.420.670.48例2求解復(fù)數(shù)代數(shù)方程組(3+j4)?I1-j4?I2=10(2-j4)?I1+j2?I

6、2=0解a=[3+4j-4j;2-4j2j];b=[100]′;formatbank;I=ab;%在MATLAB程序中,變量都看做復(fù)數(shù),并且不能在變量上方加點(diǎn)L=abs(I),Q=degree(I)%degree為自編的求相角度數(shù)的函數(shù)輸出結(jié)果為:L=1.242.77Q=29.7456.311.2拉普拉斯變換的求解例3求u(t)=40sin(3t)V的拉普拉斯變換解symst;%指定t為符號(hào)變量ut=403sin(33t);us=laplace(ut)輸出結(jié)果為:us=120/(s^2+9)例4求us=120/(s^2+9)的逆拉普拉斯變

7、換解symss;us=120/(s^2+9);ut=ilaplace(us)輸出結(jié)果為:ut=403sin(33t)1.3常微分方程的符號(hào)解9y例5求解微分方程=t-2yy(0)=0dt解dsolve(’Dy=t-23y’,’y(0)=0’)輸出結(jié)果為:ans=1/23t-1/4+1/43exp(-23t)2MATLAB輔助分析2.1建模分析法2.1.1公式法例6正弦穩(wěn)態(tài)電路電路如圖1所示,求各電流相量.解1)建模:Z2+Z4Z=Z1+Z3;?I=U?S/Z;Z3+Z2+Z4?IZ3Z2+Z4?IC=;?IL=?IZ3+Z2+Z4Z2+Z2

8、+Z4圖1例6的電路圖2)MATLAB程序:Z1=1.5;Z2=1;Z3=j;Z4=-2j;Us=43exp(0j3pi/180);Z=Z1+Z33(Z2+Z4)/(Z3+Z2+