利用Matlab剖析交流电路.docVIP

利用Matlab剖析交流电路 　　摘 要：探讨利用Matlab软件辅助分析交流电路的方法。以Matlab软件为辅助工具，具体分析了动态电路、正弦稳态电路及非正弦交流电路。结果表明，利用Matlab软件辅助分析交流电路问题，逻辑清晰、程序简短、结果精确。因此，Matlab软件是分析和研究电路的有效工具。?? 　　关键词：动态电路；正弦稳态电路；非正弦交流电路；Matlab?? 　　中图分类号：TP31 文献标识码：A 文章编号：1672-7800（2011）10-0016-02?お? 　　?? 　　基金项目：广西教育厅科研项目(200911LX467)?? 　　作者简介：陆朝华(1964-)，男，广西来宾人，柳州师范高等专科学校物理与信息科学系讲师，研究方向为电路理论、计算机仿真；蓝海江(1963-)，男，广西柳城人，柳州师范高等专科学校物理与信息科学系教授，研究方向为量子光学、计算机仿真。 　　 　　0 引言?? 　　 Matlab软件是实现科学计算、符号运算和图形处理的高效工具，是教学和科研的得力助手。比如，其在电路分析及计算机仿真等方面的应用，目前仍然是人们研究的热点。交流电路是生产设备的重要组成部分，其结构复杂，分析和计算过程烦琐而费时，引入Matlab软件进行辅助分析和研究是化解这一难题的有效手段。?? 　　 本文通过实例，探讨利用Matlab软件辅助分析动态电路、正弦稳态电路及非正弦交流电路的方法。?? 　　1 利用Matlab分析动态电路?? 　　 例1 RLC与??U-S串联的电路如图1所示。已知，R=100Ω，L=0.5H，C=500μF，U-S=10V。求开关S闭合后电容C的端电压u-c??。?? 　　 　　 建立数学模型。依据基尔霍夫电压定律，电路的微分方程为: 　　?И?[SX(]d+2u-c[]dt+2[SX)]+[SX(]R[]L[SX)][SX(]du-c[]dt[SX)]+[SX(]1[]LC[SX)]u-c=[SX(]U-s[]LC[SX)] [JY](1)?И? 　　 电路的初始条件为??u-c(0)=0，du-c(0)/dt=0??。把已知条件及初始条件代入式(1)，得 　　?И?[SX(]d+2u-c[]dt+2[SX)]+200[SX(]du-c[]dt[SX)]+4000u-c=40000(2)?И? 　　 式(2)是二阶微分方程，利用数学方法求解，烦琐而复杂。而利用Matlab软件求解，只需1条dsolve语句即可完成。程序如下：?? 　　 %分析动态电路程序program1.m?? 　　 us=dsolve(D2y+200*Dy+4000*y=40000,Dy(0)=0,y(0)=0) %求解微分方程?? 　　 运行程序program1.m，即可得到电压 　　?И?u-c=10-11.4550e+{-22.5403t}+1.4550e+{-177.4597t}?И? 　　 由此可知，利用Matlab软件求解动态电路微分方程，程序简短、简单快捷。?? 　　2 利用Matlab分析正弦稳态电路?? 　　 例2 正弦稳态电路如图2所示。已知，??u-S=10sin(500t+2π/3)V，i-S=5sin(500t+π/6)A，R-1=2.5Ω，R-2=2Ω，R-3=1Ω，R-4=3Ω，L=0.50mH及 C=2000μF??。求各支路的电流。?? 　　 　　 建立数学模型。依据基尔霍夫电流定律，得 　　?И?[JB(]-[AKI#8226;]-1-[AKI#8226;]-4-[AKI#8226;]-6=[AKI#8226;]-S?? 　　-[AKI#8226;]-2-[AKI#8226;]-3+[AKI#8226;]-4=0?? 　　[AKI#8226;]-3-[AKI#8226;]-5+[AKI#8226;]-6=-[AKI#8226;]-S[JB)}][JY](3)?И? 　　 依据基尔霍夫电压定律，得 　　?И?[JB(]-R-1[AKI#8226;]-1+(R-2+jX-L)[AKI#8226;]-2+R-4[AKI#8226;]-4=0?? 　　-(R-2+jX-L)[AKI#8226;]-2+R-3[AKI#8226;]-3=[AKU#8226;]-S?? 　　-R-3[AKI#8226;]-3-R-4[AKI#8226;]-4-jX-C[AKI#8226;]-6=0[JB)}][JY](4)?И? 　　 　　 式(3)、(4)的矩阵形式为 　　?И?[JB([][HL(6]-1[]0[]0[]-1[]0[]-1?? 　　0[]-1[]-1[]1[]0[]0?? 　　0[]0[]1[]0[]-1