1、題目:Matlab運算與受控源電阻電路分析初始條件：1.Matlab7.0以上版本軟件；2.課程設計輔導資料：“Matlab語言基礎及使用入門”、“Matlab及在電子信息課程中的應用”、線性代數及相關書籍等；3.先修課程：高等數學、線性代數、電路、Matlab應用實踐及信號處理類相關課程等。要求完成的主要任務:（包括課程設計工作量及其技術要求，以及說明書撰寫等具體要求）1.課程設計時間：1.5周，分散；2.課程設計內容：根據指導老師給定的7套題目，按規定選擇其中1套完成；3.本課程設計統一技術要求：研讀輔導資料對應章節，對選定的設計題目進行理論分析，針對具體設計部分的原理分析、建模、必要的推

2、導和可行性分析，畫出程序設計框圖，編寫程序代碼（含注釋），上機調試運行程序，記錄實驗結果（含計算結果和圖表等）并對實驗結果進行分析和總結。具體設計要求包括：1初步了解Matlab、熟悉Matlab界面、進行簡單操作；2MATLAB的數值計算：創建矩陣、矩陣運算、多項式運算、線性方程組、數值統計；3基本繪圖函數：plot,plot3,mesh,surf等，要求掌握以上繪圖函數的用法、簡單圖形標注、簡單顏色設定等；4使用文本編輯器編輯m文件，函數調用；5能完成簡單的電路Matlab編程分析；6按要求參加課程設計實驗演示和答辯等。7按要求參加課程設計實驗演示和答辯等。4.課程設計說明書按學校“課程設

3、計工作規范”中的“統一書寫格式”撰寫，具體包括：1目錄；2與設計題目相關的理論分析、歸納和總結；3與設計內容相關的原理分析、建模、推導、可行性分析；4程序設計框圖、程序代碼（含注釋）、程序運行結果和圖表、實驗結果分析和總結；5課程設計的心得體會（至少500字）；6參考文獻（不少于5篇）；7其它必要內容等。時間安排：1.5周（分散進行）指導教師簽名：年月日系主任（或責任教師）簽名：年月日MatlabMatlab應用課程設計目錄1.1簡介11.2主要功能11.3典型應用21.4基本命令32 2.設計題目：MatlabMatlab運算與應用設計套題一63 3 . .任務完成過程83.1題一83.2題

4、二83.3題三93.4題四103.5題五153.6題六163.7題七173.8題八183.9題九193.10題十193.11題十一203.12題十二224 4 . .心得體會245 5 . .參考文獻251、Matlab概述1.1簡介Matlab是美國MathWorks公司出品的商業數學軟件，用于算法開發、數據可視化、數據分析以及數值計算的高級技術計算語言和交互式環境，主要包括MATLAB和Simulink兩大部分Matlab主要用于數值運算，同時還能用于控制系統設計與分析、圖像處理、信號處理與通訊、金融建模和分析等。Simulink部分提供了一個可視化開發環境，常用于系統模擬、動態/嵌入式系

5、統開發等方面。如今，Matlab不僅成為各大公司和科研機構的專用軟件，在大學校園也得到了普及，許多本科和專科的學生借助它來學習大學數學和計算方法等課程，而碩士生和博士生在做科學研究時， 也經常要用到MATLAB進行數值計算和圖形處理。 可以說，MATLAB已經深入到了各個專業的很多學科。1.2 Matlab的功能1.2.1強大的科學計算機數據處理能力Matlab是一個包含大量計算算法的集合。其擁有600多個工程中要用到的數學運算函數，可以方便的實現用戶所需的各種計算功能。函數中所使用的算法都是科研和工程計算中的最新研究成果，而前經過了各種優化和容錯處理。在通常情況下，可以用它來代替底層編程語言

6、，如C和C+。在計算要求相同的情況下，使用MATLAB的編程工作量會大大減少。MATLAB的這些函數集包括從最簡單最基本的函數到諸如矩陣，特征向量、快速傅立葉變換的復雜函數。函數所能解決的問題其大致包括矩陣運算和線性方程組的求解、微分方程及偏微分方程的組的求解、符號運算、傅立葉變換和數據的統計分析、工程中的優化問題、稀疏矩陣運算、復數的各種運算、三角函數和其他初等數學運算、多維數組操作以及建模動態仿真等。1.2.2出色的圖形處理功能圖形處理功能Matlab自產生之日起就具有方便的數據可視化功能， 以將向量和矩陣用圖形表現出來，并且可以對圖形進行標注和打印。高層次的作圖包括二維和三維的可視化、圖

7、象處理、動畫和表達式作圖。可用于科學計算和工程繪圖。新版本的Matlab對整個圖形處理功能作了很大的改進和完善，使它不僅在一般數據可視化軟件都具有的功能方面更加完善，而且對于一些其他軟件所沒有的功能，Matlab同樣表現了出色的處理能力。同時對一些特殊的可視化要求，例如圖形對話等，Matlab也有相應的功能函數，保證了用戶不同層次的要求。另外新版本的Matlab還著重在圖形用戶界面（GUI）的制作上作了很大的改善，對這方面有特殊要求的用戶也可以得到滿足。1.2.3應用廣泛的模塊集合工具箱Matlab對許多專門的領域都開發了功能強大的模塊集和工具箱。一般來說，它們都是由特定領域的專家開發的，用戶

8、可以直接使用工具箱學習、應用和評估不同的方法而不需要自己編寫代碼。目前，Matlab已經把工具箱延伸到了科學研究和工程應用的諸多領域，如數據采集、數據庫接口、概率統計、樣條擬合、優化算法、偏微分方程求解、神經網絡、小波分析、信號處理、圖像處理、系統辨識、控制系統設計、LMI控制、魯棒控制、模型預測、模糊邏輯、金融分析、地圖工具、非線性控制設計、實時快速原型及半物理仿真、嵌入式系統開發、定點仿真、DSP與通訊、電力系統仿真等，在工具箱（Toolbox）都有。1.2.4實用的程序接口和發布平臺新版本的Matlab可以利用Matlab編譯器和C/C+數學庫以及圖形庫，將自己的Matlab程序自動轉換

9、為獨立于Matlab運行的C和C+的代碼。 允許用戶編寫可以和Matlab進行交互的C或C+語言程序。 另外，Matlab網頁服務程序還容許在Web應用中使用自己的Matlab數學和圖形程序。Matlab的一個重要特色就是具有一套程序擴展系統和一組稱之為工具箱的特殊應用子程序。工具箱是Matlab函數的子程序庫，每一個工具箱都是為某一類學科專業和應用而定制的，主要包括信號處理、控制系統、神經網絡、模糊邏輯、小波分析和系統仿真等方面的應用。1.2.5應用軟件開發（包括用戶界面）在開發環境中，使用戶更方便地控制多個文件和圖形窗口；在編程方面支持了函數嵌套，有條件中斷等；在圖形化方面，有了更強大的圖

10、形標注和處理功能，包括對性對起連接注釋等；在輸入輸出方面，可以直接向Excel和HDF5進行連接1.3 Matlab的典型應用Matlab的應用領域十分廣闊，典型的應用舉例如下:(1)數據分析(2)數值與符號計算；(3)工程與科學繪圖；(4)控制系統設計；(5)航天工業；(6)汽車工業；(7)生物醫學工程；(8)語音處理；(9)圖像與數字信號處理；(10)財務、金融分析；(11)建模、仿真及樣機開發；(12)新算法研究開發；(13)圖形用戶界面設計。1.4基本命令1.4.1變量賦值Matlab是矩陣運算，因此Matlab運算的基本變量通常都是矩陣(1)直接賦值對2X2的矩陣A進行賦值A=13;

11、24A=1 324各行元素由分號分隔，同行元素由空格、逗號分隔(2)冒號賦值A=S1:S2:S3其中S1為起始值，S2為步長，S3為終止值。A=1:2:10A=13579步長S2的默認值為1。A=1:6A=123456(3)函數賦值n返回變量列表=函數名(輸入變量列表)x=1y=sin(x)y=0.84151.4.2矩陣運算Matlab的運算都是以矩陣的方式進行矩陣與矩陣的運算：+、-：相同維數的矩陣才能進行加減運算.*、./:相同維數的矩陣對應元素的乘除運算*：滿足矩陣乘法的規則A：方陣才能進行A的幕次運算n矩陣與標量的運算+、-：矩陣各元素都與標量進行加、減*、/:矩陣各元素都與標量進行乘

12、、除矩陣元素標識：如果要得到矩陣中某一元素或者一組元素，可采用以下方法：A(m,n)：第m行、第n列的元素A(：,n)：第n列的所有元素組成的向量A(m,：):第m行的所有元素元素組成的向量A(m1:m2,n1:n2):從m1行到m2行，以及從n1列到n2列的子矩陣1.4.3繪圖命令Plot線性X-Y坐標圖Stem柱狀圖subplot在一個figure中子圖繪制Loglog雙對數坐標圖figure打開圖形窗口hold是否允許新圖覆蓋舊圖Title畫題頭Xlabelx軸標注Ylabely軸標注Text任意定位的標注grid圖形網格1.4.4命令文件和函數文件M文件有兩種形式： 命令文件(Scri

13、ptFile)和函數文件(FunctionFile)。這兩種文件的擴展名，均為“.mo對于一些比較簡單的問題，在指令窗中直接輸入指令計算；對于復雜計算，采用命令文件(Scriptfile)最為合適。Matlab只是按文件所寫的指令執行。命令文件的構成比較簡單，只是一串按用戶意圖排列而成的MATLAB指令集合。命令文件運行后，所產生的所有變量都駐留在Matlab基本工作空間(Baseworkspace中。只要用戶不使用清除指令(clear),Matlab指令窗不關閉，這些變量將一直保存在基本工作空間中。與命令文件不同，函數文件可以傳遞參數，命令文件不具備參數傳遞的功能。Matlab提供的函數指令

14、大部分都是由函數文件定義的。從形式上看，與命令文件不同，函數文件的笫一行總是以“function引導的函數中明行”；從運行上看，與命令文件運行不同，每當函數文件運行，MATLAB就會專門為它開辟一個臨時工作空間，稱為函數工作空間。當執行文件最后一條指令時，就結束該函數文件的運行，同時該臨時函數空間及其所有的中間變量就立即被清除。(1)Matlab的程序結構1)順序結構2)循環結構：forend循環while-end循環。3)選擇結才ifelseendswitch-case結構。4)異常處理機制：try-catch結構。(2)程序流控制常用指令：return,echo,input,pause,k

15、eyboard,break警示令：error,warning(3)Matlab語言的函數的基本結構Function返回變量列表二函數名(輸入變量列表)2.設計題目Matlab運算與受控源電阻電路分析“2sin(0.3二)2cos(0.3二)1 .計算y1=)和y2=)；151.52.畫出衰減震蕩曲線y=e(t/3)sin3t及其他的包絡線y0=e3islogical(L)X=A(L)(3)A=1:4;5:8pow2(A)(4)A=zeros(2,3)A(:)=1:6A=A*(1+i)A1=A.?B1=A?(5)A=ones(1,2)B=ones(2)C=eye(3)D=diag(C)E=rep

16、mat(C,1,3)545.計算y=x5+0.1x4-100在x=0.1與10處的值。6.求函數1+1/2+1/3+-+1/n,n=100的值。7.求500個元素的隨機數向量A中大于0.5的元素個數。8.求線性方程組AX=b的根。0.41560.23190.20240.29300.369010.32060.46500.56620.18120.4631A=0.51460.37340.83150.35750.62780.71920.05000.42260.93410.94450.84200.07390.83950.58940.4058-0.05050.4449b=0.5959卜.74830.864

17、2-10.用FFT直接計算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w11.如圖所示的電路，已知：R=2C,R=4C,R3=12Q,R4=陽，R5=12QRJR-2fl.如us=10V,求i3,u4,u7；(2)如已知U4=6V，求Us3,12.如圖所示的電路，各元件參數已在圖中標出，求i及U29.畫出y=|1000sin(4x)|+1的x軸為對數的圖。的功率譜。R2+U4R4R63.任務完成過程2sin(03)2cos(03)3計算丫仁1+J5和y2=VA/5；3.1.1分析：直接在主界面輸入程序代碼即可計算運行。3.1.2結果：如圖3.1.1和圖3.1.2所示。

18、yl=2*sin(0.3+pi)/(1+sqrtyl=2*sin(0.3+pi)/(1+sqrt(5)(5) y2=2*cos(0.3*pi)/(1+sqrty2=2*cos(0.3*pi)/(1+sqrt(5)(5) )yl=y2y2= =0.50000.5000I I36333633圖3.1.1題1.1程序結果圖圖3.1.2題1.2程序結果圖3.2畫出衰減震蕩曲線y=eY/3)/3)sin3t及其他的包絡線y0=十,t的取值范圍是0,4二.3.2.1分析：此類畫圖先聲明t范圍，再列出y和t、y0和t的關系式，最后畫圖3.2.2程序：t=0:pi/50:4*pi;holdon;y=exp(-

19、t/3).*sin(3*t);y0=exp(-t/3);plot(t,y,r);plot(t,y0,b);%t的取值范圍%允許在同一坐標系下繪制不同的圖形%y關于t的函數%y0關于t的函數%畫出y的二維圖形，紅色。%畫出y的二維圖形，藍色。x,y=meshgrid(-8:0.25:8);%在0,4pi04pi區域生成網格坐xlabel(時間t);ylabel(y軸)；標記x軸與y軸。legend(y=exp(-t./3).*sin(3.*t),y0=exp(-t./3);%添加文字標grid%畫出網格3.2.3結果：如圖3.2.1所示。圖3.2.1衰減震蕩曲線及包絡線圖sin、x2y23.3畫

20、出z=/22所表示的三維曲面。x,y的取值范圍是.xy-8,8。3.3.1分析：用meshgric數產生x,y平面區域內的網格坐標矩陣，再列出z與x,y的關系，然后再畫圖。3.3.2程序：z=sin(sqrt(x.*x+y.*y)./sqrt(x.*x+y.*y);%z關于x與y的函數surf(x,y,z)%畫出三維曲面3.3.3結果：如圖3.3.1所示。圖3.3.1三維曲面圖3.4分析下面每條指令的功能弁運行，觀察執行結果。3.4.1X=0:0.1:1;Y=X.*exp(-X);plot(X,Y),xlabel(,x?),ylabel(,y?),title(,y=x*exp(-x)?);分析

21、：X=0:0.1:1是以0.1為步長從0到1取實數，來確定X的取值。Y=X.*exp(-X)是Y關于X的函數。然后用plot(X,Y)畫出二維曲線，并用title(,y=x*exp(-x)?)為他加上圖形標題；以及用xlabel(,x?)為他加上x軸說明；和用ylabel(,y?)為他加上y軸說明。運行結果如圖3.4.1圖3.4.1題4.1運行圖3.4.2A=zeros(2,5)A(:)=-4:5L=abs(A)3islogical(L)X=A(L)分析：A=zeros(2,5)是創建一個2X5的零矩陣，A(:)=-4:5是給矩陣A從-4到5(步長為1)依次賦值。L=abs(A)3是判斷其絕對

22、值是否大于3,若是值為1,否則為0并把值賦給矩陣L。islogical(L)為判斷L是否為邏輯數組，若是值為1,否則為0。X是列向量，其數據為數組A對應邏輯數組L為1的位置的數。運行結果如圖3.4.2。ansans-4-44 45 5圖3.4.2題4.2運行圖3.4.3A=1:4;5:8pow2(A)分析：A=1:4;5:8是創建一個2季的矩陣， 其值是第一行為1到4,第二行為5到8（步長為1）。pow2（A）將矩陣A中的每個值都變為2Ax（x為矩陣中的每個值）并把該值賦給個缺省變量ans運行結果如圖3.4.312345678ansans248163264128256圖3.4.3題4.3運行圖

23、3.4.4A=zeros(2,3)A(:)=1:6A=A*(1+i)A1=A.?;B1=A?;分析：A=zeros(2,3捉創建一個2M的零矩陣。A(:)=1:6是從1到6(步長為1)給矩陣依次賦值。A=A*(1+i)是將A*(1+i)重新賦值給A。A1=A.?將人行列轉置后賦給A1;B1=A?將A轉置后再取其共腕復數將其賦給B1。運行結果如圖3.443.3. 00000000+ +3.OOOOi3.OOOOi5.5.00000000+ +5.5.00001000014.4. 0000+0000+4.4.OOOOiOOOOi6.6.00000000+6.+6.OOOOiOOOOiBl=Bl=1

24、.1.0000-1.OOOOi0000-1.OOOOi3.0000-3.OOOOi3.0000-3.OOOOi5.5.00000000-5.-5.OOOOiOOOOi圖3.4.4題4.4運行結果圖3.4.5A=ones(1,2)B=ones(2)C=eye(3)D=diag(C)E=repmat(C,1,3)1.1.0000+0000+1.OOOOi1.OOOOi2.2.0000+0000+2.OOOOi2.OOOOi1.1.00000000+ +1.OOOOi1.OOOOi3.00003.0000+ +3,OOOOi3,OOOOi5.00005.0000+ +6.OOOOi6.OOOOi2.

25、2.00000000+ +2.OOOOi2.OOOOi4.4.0000+4.OOOOi0000+4.OOOOi6.0000+6.6.0000+6.0 0000100012.0000-2.0000-2,2,OOOOiOOOOi4.4.00000000- -4.OOOOi4.OOOOi6.6.00000000- -6.OOOOi6.OOOOi分析：A=ones(1,2)創建一個1X2的全1矩陣；B=ones(2)& 建一個二階的全1方陣；C=eye(3)創建一個三階白單位矩陣；D=diag(C)取C矩陣主對角線上的元素構成向量D。E=repmat(C,1,3)將矩陣復制1毛塊，構成新的矩陣

26、E。運行結果如圖C C= =100000 010 0000010 00 010010010010010100 010 00 01 1圖3.4.5題4.5運行結果圖3.5計算y=x5+0.1x4-100在x=0.1與10處的值。3.5.1分析：我們可以寫一個函數，即編輯一個M文件，通過函數通過函數來計算y的值。3.5.2程序：functiony=T5(x)%編寫T5函數y=x.A5+0.1*x.A4-100%y關于x的函數end3.5.3結果：如圖3.5.1。T6CO.T6CO.1199-99-999999997 799999999949999999994zmszms= =99-99-99999

27、9997 799999999949999999994 T5T5(103(1031 10090000900100900100900圖3.5.1題五運行結果圖3.6求函數1+1/2+1/3+1/n,n=10的值。3.6.1分析：此題我們可以通過編寫M文件，利用循環機結構while語句或者3.6.2程序：functionsum=T6(n)%T6函數sum=0;i=1;while(iSUJUSUJU=T6=T6(100)(100)SUM=SUM=5.18745.1874圖3.6.2題六運行結果圖3.7求500個元素的隨機數向量A中大于0.5的元素個數。3.7.1分析：首先我們因該利用rand函數產生5

28、00個元素的隨機向量，我們可以將其定義為1X500的矩陣。統計其中大于0.5的元素個數時，我們可以借助于for循環語句將矩陣中的每一個元素與0.5一比較而統計個數。若其大于0.5,則N的值加一。最后輸出N的值即為我們所需要的結果。3.7.2程序：A=rand(1,500);%生成隨機數組An=0;for(i=1:1:500)%循環判斷A數組元素是否大于0.5ifA(1,i)0.5;n=n+1;endendfprintf(所產生的500個隨機數中大于0.5的有dj,n)%輸出n3.7.3結果：如圖3.7.1。T7所產生舸。幗機姊奸0摘26沙圖3.7.1題七運行結果圖3.8求線性方程組AX=b的根

29、。0.05050.4449b=0.59590.74830.8642-3.8.1分析：此題利用線性代數中所學的求解線性方程組根的求法，直接先求A的逆矩陣再乘以Bo3.8.2程序：A=0.4156,0.2319,0.2024,0.2930,0.3690;.嗨俞入A矩陣0.3206,0.4650,0.5662,0.1812,0.4631;.0.5146,0.3734,0.8315,0.3575,0.6278;.0.7192,0.0500,0.4226,0.9341,0.9445;.0.8420,0.0739,0.8395,0.5894,0.4058;B=0.0505;0.4449;0.5959;0.

30、7483;0.8642;%輸入B矩陣X=inv(A)*B%取人逆矩陣乘B3.8.3結果：如圖3.8.1-2.-2.274227421.43791,389038904.07104.0710-2.-2.199819980.41560.3206A=0.51460.71920.84200.23190.20240.2930 0.36900.46500.56620.1812 0.46310.37340.83150.3575 0.62780.05000.42260.9341 0.94450.07390.83950.5894 0.40583.9畫出y=|1000sin(4x)|+1的x軸為對數的圖。3.9.1

31、分析：做此題時，先定義x,然后寫出y關于x的函數，而后將x軸設為對數模式，最后畫圖。%畫出網格3.9.3結果：如圖3.9.1IZUUi75jrT T- -tt丁：.廠IUUUAnnouuii(也軸III-3皿=H:-1：HI500丁力三rynnJ1Ann4UU:11力於：h：=：；：200二二三二二二后301aq二產二產力廣廣 n 工w.1三三三-、 0 二rn 卜二產wq=二，三;w 口J4Ulin/：I：I：?：,：iri；u-1二31(101010圖3.9.1y=|1000sin(4x)|+1的x軸為對數圖3.10用FFT直接計算x=cos(2*pi*40*t)+2sin(2*pi*12

32、0*t)+w(t)的功率3.10.1分析：直接法又稱周期圖法，我們把隨機序列x(n)的N個觀測數據視為一能量有限的序列，直接計算x(n)的離散傅立葉變換，得X(k),然后再取其幅值的平方，并除以N,作為序列x(n)真實功率譜的估計。3.9.2程序：x=0:0.1:100;%定義x的范圍及步進長度y=abs(1000.*(sin(x)+1;semilogx(y);%y關于x函數%x軸取對數模式grid圖3.10.1功率譜圖3.11如圖所示的電路，已知：R=2C,R2=4C,R3=12R=3,R=4R7=*如Us=10V，求i3,U4,U7;如已知U4=6V，求Us,i3,U73.10.2程序:F

33、s=1000;%采樣頻率t=0:1/Fs:1;%定義t的范圍與刻度x=cos(2.*pi.*40.*t)+2*sin(2*pi*120*t)+rand(size(t);%產生含有噪聲的window=boxcar(length(x);nfft=1024;P,f=periodogram(x,window,nfft,Fs);plot(f,10*log10(P);xlabel(頻率),ylabel(功率)序列%矩形窗%fft點數%調用PscB數%畫出功率譜%給乂軸、y軸添加標簽3.10.3結果：如圖3.10.13.11.1分析：解此題時先對電路進行分析，用回路法列出方程，然后用線性代數的求解方程組的方

34、法進行求解即可。3.11.2程序：問題(1)：R1=2;R2=4;R3=12;R4=4;R5=12;R6=4;R7=2;Us=10;R11=R1+R2+R3;R12=-R3;R13=0;R21=-R3;R22=R3+R4+R5;R23=-R5;R31=0;R32=-R5;R33=R5+R6+R7;U1=Us;U2=0;U3=0;R=R11,R12,R13;R21,R22,R23;R31,R32,R33;U=U1;U2;U3;I=RU;Ia=I(1);Ib=I(2);Ic=I(3);i3=Ia-Ib,u4=R4*Ib,u7=R7*Ic問題(2):R1=2;R2=4;R3=12;R4=4;R5=1

35、2;R6=4;R7=2;u4=6;Ib=u4./R4;R11=R1+R2+R3;R12=0;Us1=-1;R21=R3;R22=R5;Us2=0;R31=0;R32=R5+R6+R7;Us3=0;U1=Ib*R3;U2=Ib*(R3+R4+R5);U3=Ib.*R5;A=R11,R12,Us1;R21,R22,Us2;R31,R32,Us3;B=U1;U2;U3;C=AB;Ia=C(1);IC=C(2);Us=C(3);Us,i3=Ia-Ib,u7=R7*Ib3.11.3結果：問題一如圖3.11.1,問題二如圖3.11.2Rl=2：R2=4;R3=12：R4=4;R6=12;R6=4;R7=2

36、;Us=10;KI1=R1+R2+K3;RI2=R3：KI3=0;K21=-R3;K2,2=RS-i-R4+R5;K23=-R5：R31=0;K32=-R5:：R33=R5+R6+R7;U1=U3:U2=0kU3=a.RnCRll.IUR113;R21.RW2.R23：R：3：LR32,R33:U=U1;U2;U3二I=RU;Ia=I(1);工匕=工C2)二Ta=T(3):u4=R4*Ibu7=R7*Ici3=O.3704u4-2.2222u7=CL7407圖3.11.2題11問題一結果及程序圖Rl=2;R2=4;R3=12;R4=4;R5=12;R6=4：R7=2;u4=6;Ib=u4-/R

37、4;RIi=RH-R2-kR3;：R12=0；Usl=-l;R21=R3;K22=R5,Us2=0;R31=0:R32=R5+R6+RT:Us3=0;Ul-Ib.*R3;U2-Ib.*(R3+R4+R5);U3=Ib.*R5A=RI1,R12.Umi;RZ1,底也Us2：R3i,R32,Us3;B=UI；U2；U3:C=AXB;Ia-C(1);IC=C(2);Usr=C(3);U5g3三工寸工6單u7=K7*IbUs=27i3=1 1u7u7三3 3圖3.11.2題11問題二結果及程序圖3.12如圖所示的電路，各元件參數已在圖中標出，求i i及弘3.12.1分析：解此題時先對電路進行分析，用網

38、孔法列出方程，通過整理列出線性方程組，然后用線性代數的求解方程組的方法進行求解即可。3.12.2程序：R1=1;R2=0.2;R3=0.5;R4=1;R5=1;R6=0.5;R11=R1+R2;R12=-R2;R13=0;R14=0;R21=-R2;R22=R2+R3;R23=-R3;R24=0;R31=0;R32=-R3;R33=R3+R4+R5;R34=-R4;R41=0.3;R42=-0.3;R43=-1;R44=1.5;U1=2;U2=1;U3=-3;U4=-1;R=R11,R12,R13,R14;R21,R22,R23,R24;R31,R3,R33,R34;R41,R42,R43,R

39、44;U=U1;U2;U3;U4;I=RU;i=I(2)-I(4),u2=(I(1)-I(2)*R23.12.3結果:如圖3.12.1 Rl=l;R2=0.2Rl=l;R2=0.2；R3=0.5;R4=l;R5=1;R5=0.R3=0.5;R4=l;R5=1;R5=0.5;5; RI1=*L-B2RI1=*L-B2：ltl2=-S2ltl2=-S2：R13=0;R14=0R13=0;R14=0; ; R21=-R2R21=-R2; ;R22=R2+R3R22=R2+R3：R23=-R3R23=-R3：R24=0;R24=0; R31=0R31=0：R32-R3;R32-R3;良33=R3+R4+R533=R3+R4+R5: :R34=-R4;R34=-R4; R41=0,3;R42=-0.3;3=-R41=0,3;R42=-0.3;3=-1 1;R44=1;R44=1. .5;5; U1=2U1=2：U2=1U2=1：U3=-3U3=-3：U4=-1U4=-1： h hRll,Rl2,R13,RiqRll,Rl2,R13,Riq；R21,K22,R23,R24;R31,R32,E33,R34;Ml,R42