【正文】 MATLAB電路仿真實驗報告武漢大學電氣工程學院MATLAB電路仿真實驗報告班級：0810 學號：2008302540299姓名：李德澳2010年7月目錄實驗一 直流電路（1） 2實驗二 直流電路（2） 8實驗三 正弦穩(wěn)態(tài) 17實驗四 交流分析和網(wǎng)絡函數(shù) 26實驗五 動態(tài)電路 31實驗六 頻率響應 43實驗一 直流電路（1）一 實驗目的1 加深對直流電路的節(jié)點電壓法和網(wǎng)孔電流法的理解2 學習使用MATLAB的矩陣運算的方法二 實驗示例1節(jié)點分析電路如圖所示（見書本12頁），求節(jié)點電壓V1,V2,V3.根據(jù)電路圖得到矩陣方程，根據(jù)矩陣方程使用matlab命令為Y = 節(jié)點v1，v2和v3： v = 2 回路分析電路如圖所示（見書本13頁），使用解析分析得到同過電阻RB的電流，另外求10V電壓源的輸出功率。分析電路得到節(jié)點方程，根據(jù)節(jié)點方程得到矩陣方程，根據(jù)矩陣方程，使用matlab的命令為z=[40,10,30。 10,30,5。 30,5,65]。v=[10,0,0]39。I=inv(z)*v。IRB=I(3)I(2)。fprintf(39。the current through R is % Amps \n39。,IRB)ps=I(1)*10。fprintf(39。the power supplied by 10v source is % watts\n39。,ps)結果為：the current through R is Ampsthe power supplied by 10V source is watts三 實驗內(nèi)容1 根據(jù)書本15頁電路圖，求解電阻電路，已知：R1=2Ω,R2=6Ω,R3=12Ω,R4=8Ω,R5=12Ω,R6=4Ω,R7=2Ω(1) 如果Us=10V,求i3,u4,u7(2) 如果U4=4V,求Us,i3,i7使用matlab命令為clear% 初始化阻抗矩陣Z=[20 12 0。 12 32 12。 0 12 18]。% 初始化電壓矩陣V=[10 0 0]39。% 解答回路電流I=inv(Z)*V。% I3的計算I3=I(1)I(2)。fprintf(39。the current I3 is % Amps\n39。,I3)% U4的計算U4=8*I(2)。fprintf(39。the voltage U4 is % Vmps\n39。,U4)% U7的計算U7=2*I(3)。fprintf(39。the voltage U7 is % Vmps\n39。,U7)結果the current I3 is Ampsthe voltage U4 is Vmpsthe voltage U7 is Vmpsclear% 初始化矩陣XX=[20 1 0。 12 0 12。 0 0 18]。% 初始化矩陣YY=[6 16 6]39。% 進行解答A=inv(X)*Y。% 計算各要求量Us=A(2)I3=A(1)I7=A(3)結果Us = I3 = I7 =2 求解電路里的電壓如圖14（書本16頁），求解V1,V2,V3,V4,V5使用matlab命令為clear% 初始化節(jié)點電壓方程矩陣Z=[ 5 。 0 0。 0 。 1 0 1 1 0。 0 0 1]。I=[0 6 5 0 0]39。% 解答節(jié)點電壓U1，U3，U4與Vb，IaA=inv(Z)*I。% 最終各電壓計算V1=A(1)V2=A(1)10*A(5)V3=A(2)V4=A(3)V5=24結果V1 =V2 = V3 =V4 =V5 = 243 如圖15（書本16頁），已知R1=R2=R3=4Ω,R4=2Ω,控制常數(shù)k1=，k2=4，is=2A，求i1和i2.使用matlab命令為clear% 初始化節(jié)點電壓方程矩陣Z=[ 0 。 1 1 。 0 0 1。 1 1 4 0]。I=[2 0 0 0]39。% 解答節(jié)點電壓V1，V2及電流I1，I2A=inv(Z)*I。% 計算未知數(shù)V1=A(1)V2=A(2)I1=A(3)I2=A(4)結果如下：V1 =6V2 =2I1 = 1I2 =1實驗二 直流電路（2）一 實驗目的1 加深多戴維南定律，等效變換等的了解2 進一步了解matlab在直流電路中的作用二實驗示例如圖所示（圖見書本17頁21），分析并使用matlab命令求解為clear,format pactR1=4。R2=2。R3=4。R4=8。is1=2。is2=。a11=1/R1+1/R4。a12=1/R1。a13=1/R4。a21=1/R1。a22=1/R1+1/R2+1/R3。a23=1/R3。a31=1/R4。a32=1/R3。a33=1/R3+1/R4。A=[a11,a12,a13。a21,a22,a23。a31,a32,a33]。B=[1,1,0。0,0,0。0,1,1]。X1=A\B*[is1。is2。0]。uoc=X1(3)。X2=A\B*[0。0。1]。Req=X2(3)。RL=Req。P=uoc^2*RL/(Req+RL)^2。RL=0:10,p=(RL*uoc./(Req+RL)).*uoc./(Req+RL),figure(1),plot(RL,p),gridfor k=1:21ia(k)=(k1)*。X=A\B*[is1。is2。ia(k)]。u(k)=X(3)。endfigure(2),plot(ia,u,39。x39。),gridc=polyfit(ia,u,1)。%ua=c(2)*ia=c(1) , 用擬合函數(shù)術,c(1),c(2)uoc=c(1),Req=c(2)RL = 0 1 2 3 4 5 6 7 8 9