matlab應(yīng)用課程設(shè)計(jì)--matlab運(yùn)算與戴維南定理應(yīng)用分析

1、<p>　　Matlab應(yīng)用課程設(shè)計(jì)任務(wù)書</p><p>　　學(xué)生姓名:__ _ 專業(yè)班級(jí): </p><p>　　指導(dǎo)教師: 工作單位: </p><p>　　題 目:Matlab運(yùn)算與戴維南定

2、理應(yīng)用分析</p><p><b>　　初始條件:</b></p><p>　?。?）提供實(shí)驗(yàn)室機(jī)房及其Matlab6.5以上版本軟件； </p><p>　?。?）《MATLAB教程》學(xué)習(xí)</p><p>　　要求完成的主要任務(wù):(包括課程設(shè)計(jì)工作量及其技術(shù)要求,以及說明書撰寫等具體要求):</p>&l

3、t;p>　　選擇一本《MATLAB教程》認(rèn)真學(xué)習(xí)該教程的全部內(nèi)容，包括基本使用方法、數(shù)組運(yùn)</p><p>　　算、矩陣運(yùn)算、數(shù)學(xué)運(yùn)算、程序設(shè)計(jì)、符號(hào)計(jì)算、圖形繪制、GUI設(shè)計(jì)等內(nèi)容；</p><p>　　對(duì)該套綜合題的10道題，進(jìn)行理論分析，針對(duì)具體設(shè)計(jì)部分的原理分析、建模、必 </p><p>　　要的推導(dǎo)和可行性分析，畫出程序設(shè)計(jì)框圖，編寫程序代碼（含注

4、釋），上機(jī)調(diào)試運(yùn) </p><p>　　行程序，記錄實(shí)驗(yàn)結(jié)果（含計(jì)算結(jié)果和圖表）。</p><p>　?。?）對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行分析和總結(jié)；</p><p>　?。?）要求閱讀相關(guān)參考文獻(xiàn)不少于5篇；</p><p>　?。?）根據(jù)課程設(shè)計(jì)有關(guān)規(guī)范，按時(shí)、獨(dú)立完成課程設(shè)計(jì)說明書。</p><p><b>　　時(shí)間

5、安排:</b></p><p>　?。?）布置課程設(shè)計(jì)任務(wù)，查閱資料，學(xué)習(xí)《MATLAB教程》 十周；</p><p>　　（2）進(jìn)行編程設(shè)計(jì) 一周；</p><p>　?。?）完成課程設(shè)計(jì)報(bào)告書

6、一周；</p><p>　　指導(dǎo)教師簽名: 年 月 日</p><p>　　系主任（或責(zé)任教師）簽名: 年 月 日</p><p><b>　　目錄</b></p><p>　　1、MA

7、TLAB概述1</p><p>　　1.1MATLAB簡介1</p><p>　　1.2MATLAB的功能1</p><p>　　1.3MATLAB的典型應(yīng)用3</p><p><b>　　2、設(shè)計(jì)題目4</b></p><p><b>　　3、設(shè)計(jì)內(nèi)容6</b>

8、</p><p><b>　　第一題6</b></p><p><b>　　第二題7</b></p><p><b>　　第三題8</b></p><p><b>　　第四題10</b></p><p><b>　

9、　第五題14</b></p><p><b>　　第六題15</b></p><p><b>　　第七題16</b></p><p><b>　　第八題17</b></p><p><b>　　第九題18</b></p>

10、<p><b>　　第十題19</b></p><p><b>　　第十一題20</b></p><p>　　4、 心得體會(huì)22</p><p><b>　　摘要</b></p><p>　　MATLAB是一種高性能的工程計(jì)算語言。它集計(jì)算，可視化，可編程于一

11、體，并在一個(gè)易于使用的環(huán)境下用熟悉的數(shù)學(xué)符號(hào)解決問題。典型用途包括數(shù)學(xué)計(jì)算，算法，開發(fā)，日期，采集，建模，仿真，可視化，科學(xué)和工程圖形，包括圖形界面的應(yīng)用開發(fā)。</p><p>　　MATLAB有可以用于曲線擬合的內(nèi)建函數(shù)。MathWorks公式也提供了很多工具箱可以用于曲線擬合。這些方法可以用來做線性或者非線性曲線擬合。MATLAB也有一個(gè)開放的工具箱――曲線擬合工具箱（Curve Fitting Toolbo

12、x），它可以用于參數(shù)擬合，也可以用于非參數(shù)擬合。</p><p>　　在MATLAB中，多項(xiàng)式擬合函數(shù)polyfit可以用來計(jì)算擬合的得到多項(xiàng)式系數(shù)。默認(rèn)的擬合目標(biāo)是最小方差最小，即最小二乘法擬合數(shù)據(jù)。此時(shí)，判斷的依據(jù)過擬合曲線的道德原始數(shù)據(jù)之間的平均誤差是否達(dá)到最小。</p><p>　　關(guān)鍵詞：工程語言；數(shù)學(xué)計(jì)算 ；曲線擬合</p><p><b>　

13、　Abstract</b></p><p>　　Matlab is a high-performance language for technical computing.gratinates computation,visualization,and programming in an easy-to-environment where problems and solutions are expr

14、essed in familiar mathematical notation.Typical uses include Math and computation,Algorithm development,Data acquisition,Modeling,simulation,and prototyping,Data analysis,exploration,and visualization Scientific and engi

15、neering graphics,Application development, including graphical user interface building.</p><p>　　Matlab can be used for curve fitting of the built-in function.Math Works formula also provides a lot of toolbox

16、 can used to curve fitting.MATLAB also has an open Toolbox-matlab curve fitting toolbox matlab(Curve Fitting Toolbox),it can be used for parameter fitting,can also be used for no parameter fitting. </p><p>　

17、　In Matlab, polynomial function poly fit can be used to calculate the fitting polynomial coefficients are obtained. The default of fitting target is the minimum variance minimum, is the least squares fitting data. At thi

18、s point, the judge basis of fitting curve of morality between raw data whether the minimum average error. </p><p>　　Keywords：Engineering language computing Curve Fitting </p><p>　　1、MATLAB概

19、述</p><p>　　1.1MATLAB簡介</p><p>　　MATLAB軟件由美國Math Works公司于1984年推出，經(jīng)過不斷的發(fā)展和完善，如今己成為覆蓋多個(gè)學(xué)科的國際公認(rèn)的最優(yōu)秀的數(shù)值計(jì)算仿真軟件。MATLAB具備強(qiáng)大的數(shù)值計(jì)算能力，許多復(fù)雜的計(jì)算問題只需短短幾行代碼就可在MATLAB中實(shí)現(xiàn)。作為一個(gè)跨平臺(tái)的軟件，MATLAB已推出Unix、Windows、Linux和Ma

20、c等十多種操作系統(tǒng)下的版本，大大方便了在不同操作系統(tǒng)平臺(tái)下的研究工作。</p><p>　　MATLAB軟件具有很強(qiáng)的開放性和適應(yīng)性。在保持內(nèi)核不變的情況下，MATLAB可以針對(duì)不同的應(yīng)用學(xué)科推出相應(yīng)的工具箱(toolbox)，目前己經(jīng)推出了圖象處理工具箱、信號(hào)處理工具箱、小波工具箱、神經(jīng)網(wǎng)絡(luò)工具箱以及通信工具箱等多個(gè)學(xué)科的專用工具箱，極大地方便了不同學(xué)科的研究工作。國內(nèi)已有越來越多的科研和技術(shù)人員認(rèn)識(shí)到MATL

21、AB的強(qiáng)大作用，并在不同的領(lǐng)域內(nèi)使用MATLAB來快速實(shí)現(xiàn)科研構(gòu)想和提高工作效率。</p><p>　　MATLAB提供了20類圖像處理函數(shù),涵蓋了圖像處理的包括近期研究成果在內(nèi)的幾乎所有的技術(shù)方法,是學(xué)習(xí)和研究圖像處理的人員難得的寶貴資料和加工工具箱。這些函數(shù)按其功能可分為:圖像顯示;圖像文件I/O;圖像算術(shù)運(yùn)算;幾何變換;圖像登記;像素值與統(tǒng)計(jì);圖像分析;圖像增強(qiáng);線性濾波;線性二元濾波設(shè)計(jì);圖像去模糊;圖像

22、變換;鄰域與塊處理;灰度與二值圖像的形態(tài)學(xué)運(yùn)算;結(jié)構(gòu)元素創(chuàng)建與處理;基于邊緣的處理;色彩映射表操作;色彩空間變換;圖像類型與類型轉(zhuǎn)換。</p><p>　　1.2MATLAB的功能</p><p>　　(1)強(qiáng)大的科學(xué)計(jì)算機(jī)數(shù)據(jù)處理能力</p><p>　　MATLAB是一個(gè)包含大量計(jì)算算法的集合。其擁有600多個(gè)工程中要用到的數(shù)學(xué)運(yùn)算函數(shù)，可以方便的實(shí)現(xiàn)用戶所需的

23、各種計(jì)算功能。函數(shù)中所使用的算法都是科研和工程計(jì)算中的最新研究成果，而前經(jīng)過了各種優(yōu)化和容錯(cuò)處理。在通常情況下，可以用它來代替底層編程語言，如C和C++ 。在計(jì)算要求相同的情況下，使用MATLAB的編程工作量會(huì)大大減少。MATLAB的這些函數(shù)集包括從最簡單最基本的函數(shù)到諸如矩陣，特征向量、快速傅立葉變換的復(fù)雜函數(shù)。函數(shù)所能解決的問題其大致包括矩陣運(yùn)算和線性方程組的求解、微分方程及偏微分方程的組的求解、符號(hào)運(yùn)算、傅立葉變換和數(shù)據(jù)的統(tǒng)計(jì)分析

24、、工程中的優(yōu)化問題、稀疏矩陣運(yùn)算、復(fù)數(shù)的各種運(yùn)算、三角函數(shù)和其他初等數(shù)學(xué)運(yùn)算、多維數(shù)組操作以及建模動(dòng)態(tài)仿真等。</p><p>　　(2)出色的圖形處理功能</p><p>　　圖形處理功能MATLAB自產(chǎn)生之日起就具有方便的數(shù)據(jù)可視化功能，以將向量和矩陣用圖形表現(xiàn)出來，并且可以對(duì)圖形進(jìn)行標(biāo)注和打印。高層次的作圖包括二維和三維的可視化、圖象處理、動(dòng)畫和表達(dá)式作圖?？捎糜诳茖W(xué)計(jì)算和工程繪圖。

25、新版本的MATLAB對(duì)整個(gè)圖形處理功能作了很大的改進(jìn)和完善，使它不僅在一般數(shù)據(jù)可視化軟件都具有的功能（例如二維曲線和三維曲面的繪制和處理等）方面更加完善，而且對(duì)于一些其他軟件所沒有的功能（例如圖形的光照處理、色度處理以及四維數(shù)據(jù)的表現(xiàn)等），MATLAB同樣表現(xiàn)了出色的處理能力。同時(shí)對(duì)一些特殊的可視化要求，例如圖形對(duì)話等，MATLAB也有相應(yīng)的功能函數(shù)，保證了用戶不同層次的要求。另外新版本的MATLAB還著重在圖形用戶界面（GUI）的制作

26、上作了很大的改善，對(duì)這方面有特殊要求的用戶也可以得到滿足。</p><p>　　(3)應(yīng)用廣泛的模塊集合工具箱</p><p>　　MATLAB對(duì)許多專門的領(lǐng)域都開發(fā)了功能強(qiáng)大的模塊集和工具箱。一般來說，它們都是由特定領(lǐng)域的專家開發(fā)的，用戶可以直接使用工具箱學(xué)習(xí)、應(yīng)用和評(píng)估不同的方法而不需要自己編寫代碼。目前，MATLAB已經(jīng)把工具箱延伸到了科學(xué)研究和工程應(yīng)用的諸多領(lǐng)域，諸如數(shù)據(jù)采集、數(shù)據(jù)

27、庫接口、概率統(tǒng)計(jì)、樣條擬合、優(yōu)化算法、偏微分方程求解、神經(jīng)網(wǎng)絡(luò)、小波分析、信號(hào)處理、圖像處理、系統(tǒng)辨識(shí)、控制系統(tǒng)設(shè)計(jì)、LMI控制、魯棒控制、模型預(yù)測、模糊邏輯、金融分析、地圖工具、非線性控制設(shè)計(jì)、實(shí)時(shí)快速原型及半物理仿真、嵌入式系統(tǒng)開發(fā)、定點(diǎn)仿真、DSP與通訊、電力系統(tǒng)仿真等，都在工具箱（Toolbox）家族中有了自己的一席之地。</p><p>　　(4)實(shí)用的程序接口和發(fā)布平臺(tái)</p><

28、p>　　新版本的MATLAB可以利用MATLAB編譯器和C/C++數(shù)學(xué)庫以及圖形庫，將自己的MATLAB程序自動(dòng)轉(zhuǎn)換為獨(dú)立于MATLAB運(yùn)行的C和C++的代碼。允許用戶編寫可以和MATLAB進(jìn)行交互的C或C++語言程序。另外，MATLAB網(wǎng)頁服務(wù)程序還容許在Web應(yīng)用中使用自己的MATLAB數(shù)學(xué)和圖形程序。MATLAB的一個(gè)重要特色就是具有一套程序擴(kuò)展系統(tǒng)和一組稱之為工具箱的特殊應(yīng)用子程序。工具箱是MATLAB函數(shù)的子程序庫，每一

29、個(gè)工具箱都是為某一類學(xué)科專業(yè)和應(yīng)用而定制的，主要包括信號(hào)處理、控制系統(tǒng)、神經(jīng)網(wǎng)絡(luò)、模糊邏輯、小波分析和系統(tǒng)仿真等方面的應(yīng)用。</p><p>　　(5)應(yīng)用軟件開發(fā)（包括用戶界面）</p><p>　　在開發(fā)環(huán)境中，使用戶更方便地控制多個(gè)文件和圖形窗口；在編程方面支持了函數(shù)嵌套，有條件中斷等；在圖形化方面，有了更強(qiáng)大的圖形標(biāo)注和處理功能，包括對(duì)性對(duì)起連接注釋等；在輸入輸出方面，可以直接向E

30、xcel和HDF5進(jìn)行連接。</p><p>　　1.3MATLAB的典型應(yīng)用</p><p>　　MATLAB的應(yīng)用領(lǐng)域十分廣闊，典型的應(yīng)用舉例如下：</p><p><b>　　(1) 數(shù)據(jù)分析</b></p><p>　　(2) 數(shù)值與符號(hào)計(jì)算；</p><p>　　(3) 工程與科學(xué)繪圖；

31、</p><p>　　(4) 控制系統(tǒng)設(shè)計(jì)；</p><p><b>　　(5) 航天工業(yè)；</b></p><p><b>　　(6) 汽車工業(yè)；</b></p><p>　　(7) 生物醫(yī)學(xué)工程；</p><p><b>　　(8) 語音處理；</b>

32、;</p><p>　　(9) 圖像與數(shù)字信號(hào)處理；</p><p>　　(10) 財(cái)務(wù)、金融分析；</p><p>　　(11) 建模、仿真及樣機(jī)開發(fā)；</p><p>　　(12) 新算法研究開發(fā)；</p><p>　　(13) 圖形用戶界面設(shè)計(jì)。</p><p><b>　　2、

33、設(shè)計(jì)題目</b></p><p>　　計(jì)算y1=和y2=；</p><p>　　畫出衰減震蕩曲線及其他的包絡(luò)線,t的取值范圍是[0,6].</p><p>　　畫出所表示的三維曲面。x,y的取值范圍是[-9，9]。</p><p>　　分析下面每條指令的功能并運(yùn)行，觀察執(zhí)行結(jié)果。</p><p>　　X=0

34、:0.2:2;</p><p>　　Y=X.*exp(-X);</p><p>　　plot(X,Y),xlabel(‘x’), ylabel(‘y’),title(‘y=x*exp(-x)’);</p><p>　　A=zeros(3,4)</p><p><b>　　A(:)=-4:7</b></p>

35、<p>　　L=abs(A)>3</p><p>　　islogical(L)</p><p><b>　　X=A(L)</b></p><p>　　A=[1:5;6:10]</p><p><b>　　pow2(A)</b></p><p>　　A=zero

36、s(3,5)</p><p><b>　　A(:)=1:15</b></p><p><b>　　A=A*(1+i)</b></p><p><b>　　A1=A.’;</b></p><p><b>　　B1=A’;</b></p><

37、;p>　　A=ones(2,3)</p><p><b>　　B=ones(2)</b></p><p><b>　　C=eye(4)</b></p><p><b>　　D=diag(C)</b></p><p>　　E=repmat(C,1,3)</p>

38、<p>　　計(jì)算在x=0.1與10處的值。</p><p>　　求函數(shù), n=100的值</p><p>　　求1000個(gè)元素的隨機(jī)數(shù)向量A中大于0.5的元素個(gè)數(shù)。</p><p>　　8.求線性方程組AX=b的根。</p><p>　　9.畫出y=|1000sin(4x)|+1的y軸為對(duì)數(shù)圖。</p><p&

39、gt;　　10．用FFT直接計(jì)算x=sin(2*pi*40*t)+2cos(2*pi*120*t)+w(t)的功率譜。</p><p>　　11. 在如圖 所示的電路中，已知，，，，，負(fù)載電阻可變。</p><p>　　（1）負(fù)載電阻為何值時(shí)可吸收最大功率？求此最大功率。</p><p>　?。?）研究在0~10Ω范圍內(nèi)變化時(shí)，其吸收功率的情況。</p>

40、<p><b>　　3、設(shè)計(jì)內(nèi)容</b></p><p><b>　　第一題</b></p><p>　　計(jì)算y1=和y2=；</p><p>　　分析：直接運(yùn)行即可。</p><p>　　>> y1=2*sin(0.5*pi)/(1+sqrt(6))</p>

41、<p>　　y2=3*cos(0.5*pi)/(1+sqrt(7))</p><p><b>　　y1 =</b></p><p><b>　　0.5798</b></p><p><b>　　y2 =</b></p><p>　　5.0387e-017

42、 </p><p><b>　　圖1操作界面</b></p><p><b>　　第二題</b></p><p>　　畫出衰減震蕩曲線及其他的包絡(luò)線,t的取值范圍是[0,6].</p><p>　　分析：先聲明t的范圍，寫出y、y0關(guān)于t的方程，運(yùn)行程序。</p>

43、<p>　　>> t=0:6*pi; %設(shè)置t的范圍</p><p>　　y=exp(-3*t).*cos(t/5); %生成y的數(shù)據(jù)</p><p><b>　　plot(t,y)</b></p><p><b>　　圖2-1函

44、數(shù)曲線</b></p><p>　　>> t=0:6*pi;</p><p>　　y=exp(-3*t); </p><p><b>　　plot(t,y)</b></p><p>　　圖2-1函數(shù)曲線的包絡(luò)線</p><p><b>　　第三題</b&g

45、t;</p><p>　　畫出所表示的三維曲面。x,y的取值范圍是[-9，9]。</p><p>　　分析：先設(shè)置x的范圍，用mesh函數(shù)繪制給定的三維曲面。</p><p>　　>> [x,y]=meshgrid(-9:0.5:9); %設(shè)置x，y的范圍并生成平面網(wǎng)格數(shù)據(jù)</p><p>　　z=cos(

46、sqrt(2*x.^2+2*y.^2))./sqrt(x.^2+y.^2+eps));</p><p>　　subplot(2,2,1);</p><p>　　mesh(x,y,z);</p><p>　　title(mesh(x,y,z))</p><p>　　subplot(2,2,2);</p><p>　　me

47、shc(x,y,z);</p><p>　　title(meshc(x,y,z))</p><p>　　subplot(2,2,3);</p><p>　　meshz(x,y,z)</p><p>　　title(meshz(x,y,z))</p><p>　　subplot(2,2,4);</p>&l

48、t;p>　　surf(x,y,z); %繪制該曲面 </p><p>　　Title(surf(x,y,z))</p><p>　　圖3-1 三維曲面 </p><p>　　圖3-2 操作界面</p><p><b&

49、gt;　　第四題</b></p><p>　　4.分析下面每條指令的功能并運(yùn)行，觀察執(zhí)行結(jié)果。</p><p>　　X=0:0.2:2; %X取值從0 到2 每隔0.2 取一個(gè)數(shù)</p><p>　　Y=X.*exp(-X); %Y對(duì)于X 的表達(dá)式,求出X 對(duì)應(yīng)下的Y 值</p><p>　　plot(X,Y),x

50、label(‘x’), ylabel(‘y’),title(‘y=x*exp(-x)’);</p><p>　　圖4-1題4.1執(zhí)行結(jié)果</p><p>　　分析：畫x與e的-x次方的乘積的圖像</p><p>　?。?）A=zeros(3,4) %產(chǎn)生3 x 4 的全零矩陣A</p><p>　　A(:)

51、=-4:7 %將-4 到7 依次賦值給A</p><p>　　L=abs(A)>3 %判斷A 內(nèi)各元素絕對(duì)值是否大于3,大于返回1,小于返回0,并將返回 值構(gòu)成與A相同行列的矩陣L</p><p>　　islogical(L) %判斷矩陣L 內(nèi)是否為邏輯值,是返回1,不是則返回0</p

52、><p>　　X=A(L) %將A 對(duì)應(yīng)L 內(nèi)為真的元素取出來,構(gòu)成列向量X,即絕對(duì)值大于3的值取出</p><p>　　圖4-2 題4.2執(zhí)行結(jié)果</p><p>　?。?）A=[1:5;6:10]</p><p>　　pow2(A) %將A 內(nèi)各元素冪以2 為底的結(jié)果</p><p>

53、;　　圖4-3題4.3執(zhí)行結(jié)果</p><p>　?。?）A=zeros(3,5) %產(chǎn)生3 x 5 的全零矩陣A</p><p>　　A(:)=1:15 %將1 到15 依次賦值給各元素</p><p>　　A=A*(1+i) %將A 各元素乘以(1+i),成為復(fù)數(shù)

54、矩陣</p><p>　　A1=A.’; %令A(yù)1 成為A 的轉(zhuǎn)置非共軛矩陣</p><p>　　B1=A’; %令B1 成為A 的轉(zhuǎn)置共軛矩陣</p><p>　　圖4-4-1 題4.4執(zhí)行結(jié)果1</p><p>　　圖4-4-2 題4.4執(zhí)行結(jié)果2<

55、;/p><p>　　A=ones(2,3) %產(chǎn)生一個(gè)2 x 3 全1 矩陣A</p><p>　　B=ones(2) %產(chǎn)生一個(gè)行列為2 的全1 的方陣B</p><p>　　C=eye(4) %產(chǎn)生一個(gè)行列為4 單位方陣C</p><p>　　D=diag(C)

56、 %將C 的主對(duì)角元素取出形成列向量D</p><p>　　E=repmat(C,1,3) %將C 看成一個(gè)元素,重復(fù)C 構(gòu)成一個(gè)1 x 3 的矩陣E</p><p>　　圖4-5-1題4.5執(zhí)行結(jié)果1</p><p>　　圖4-5-2 題4.5執(zhí)行結(jié)果2</p><p><b>　　第五題</b&

57、gt;</p><p>　　計(jì)算在x=0.1與10處的值。</p><p>　　分析：此題可以直接創(chuàng)建一個(gè)匿名函數(shù)進(jìn)行運(yùn)算，也可以利用M文件計(jì)算。</p><p><b>　　>> x=0.1;</b></p><p>　　y=x^5+0.1*x^4-100</p><p><b

58、>　　x=10;</b></p><p>　　y=x^5+0.1*x^4-100</p><p><b>　　y =</b></p><p><b>　　-100.0000</b></p><p><b>　　y =</b></p><p&

59、gt;<b>　　100900</b></p><p><b>　　圖5執(zhí)行結(jié)果</b></p><p><b>　　第六題</b></p><p>　　求函數(shù), n=100的值。</p><p>　　>> syms n;s=symsum(1/(n^2),1,100

60、)</p><p><b>　　eval(s)</b></p><p><b>　　s = </b></p><p>　　1589508694133037873112297928517553859702383498543709859889432834803818131090369901/972186144434381030

61、589657976672623144161975583995746241782720354705517986165248000</p><p><b>　　ans =</b></p><p><b>　　1.6350</b></p><p><b>　　圖6執(zhí)行結(jié)果</b></p>&l

62、t;p><b>　　第七題</b></p><p>　　7.求1000個(gè)元素的隨機(jī)數(shù)向量A中大于0.5的元素個(gè)數(shù)。</p><p>　　>> A=rand(1,1000); %隨機(jī)產(chǎn)生1 x 1000 矩陣A</p><p>　　sum=0;

63、 %令次數(shù)sum為零</p><p>　　for ii=1:length(A)</p><p>　　if A(ii)>0.5</p><p>　　sum=sum+1; %循環(huán)比較,如果元素大于0.5,次數(shù)加1</p><p><b>　　end</b&g

64、t;</p><p><b>　　end</b></p><p>　　Sum %輸出次數(shù)times</p><p><b>　　sum =</b></p><p><b>　　509</b></p>&

65、lt;p><b>　　圖7執(zhí)行結(jié)果</b></p><p><b>　　第八題</b></p><p>　　8.求線性方程組AX=b的根。</p><p>　　分析：此題可用左除法求解線性方程的根。</p><p>　　>> a=[0.4156,0.2319,0.2024,0.27

66、30,0.3690;0.3206,0.3650,0.5662,0.1812,0.4631;</p><p>　　0.5146,0.3734,0.7315,0.3575,0.6278;0.7192,0.0500,0.4226,0.6341,0.9445;</p><p>　　0.8420,0.0739,0.8395,0.5894,0.4058];</p><p>　　

67、b=[0.0305;0.3449;0.5759;0.7483;0.4642];</p><p>　　x=inv(a)*b</p><p><b>　　x =</b></p><p><b>　　-2.6121</b></p><p><b>　　-0.0920</b></

68、p><p><b>　　0.8847</b></p><p><b>　　3.0216</b></p><p><b>　　0.3618</b></p><p><b>　　圖8執(zhí)行結(jié)果</b></p><p><b>　　第

69、九題</b></p><p>　　畫出y=|1000sin(4x)|+1的y軸為對(duì)數(shù)圖。</p><p>　　分析：此題關(guān)鍵是用semillogy(x,y)函數(shù)將y坐標(biāo)軸表示成對(duì)數(shù)坐標(biāo)系。</p><p>　　>> x=1:0.0001:10; %設(shè)置x的范圍</p><

70、;p>　　y=abs(1000*sin(4*x))+1; %生成y的數(shù)據(jù)</p><p>　　semilogx(x,y);xlabel('x'),ylabel('log(y)'); %表示x坐標(biāo)軸是對(duì)數(shù)坐標(biāo)系</p><p>　　grid on;

71、 %給圖形加上網(wǎng)格</p><p>　　title('y=|1000sin(4x)|+1的x軸為對(duì)數(shù)的圖'); %添加標(biāo)題</p><p><b>　　圖9對(duì)數(shù)圖像</b></p><p><b>　　第十題</b></p><p>　　用FFT直接計(jì)算x=sin(2*pi*40*

72、t)+2cos(2*pi*120*t)+w(t)的功率譜。</p><p>　　分析：利用傅里葉變換處理。</p><p>　　>> Fs=1000; %采樣頻率</p><p>　　t=0:1/Fs:1;</p><p>　　x=sin(2*pi*40*t)+2

73、*cos(2*pi*120*t)+rand(size(t)); %產(chǎn)生含有噪聲的序列</p><p>　　tfft=1024; %fft點(diǎn)數(shù)</p><p>　　X=fft(x,tfft); %對(duì) x 進(jìn)行快速傅里葉變換</p><p&g

74、t;　　Pxx=abs(X).^2; %功率譜的量度為變換后幅度的平方 </p><p>　　m=0:round(tfft/2-1); %由于對(duì)稱性,取一半</p><p>　　f=m*Fs/tfft; %使頻率變

75、化從零開始</p><p>　　P=10*log10(Pxx(m+1)); </p><p>　　plot(f,P),xlabel('頻率'),ylabel('功率') %畫出功率譜</p><p><b>　　圖10功率譜</b></p><p><b>

76、;　　第十一題</b></p><p>　　11. 在如圖 所示的電路中，已知，，，，，負(fù)載電阻可變。</p><p>　?。?）負(fù)載電阻為何值時(shí)可吸收最大功率？求此最大功率。</p><p>　?。?）研究在0~10Ω范圍內(nèi)變化時(shí)，其吸收功率的情況。</p><p>　　>> R1=2;R2=2;R3=4;K1=2

77、;K2=4; %設(shè)置元件參數(shù)</p><p>　　%按照A*X=B*u列寫此電路的矩陣方程，其中X=[Ua Ub I1]’;u=[Us Is]’</p><p>　　A=[1/R3 -1/R3 -K1/R3;-1/R3 1/R1+1/R2+1/R3 K1/R3-K2;0 1/R1 1];</p><p>　　B=[0 1;1/R1 0;1/R

78、1 0];</p><p>　　Us=6;Is=0; %令I(lǐng)s=0，求Us=X1(1)</p><p>　　X1=A\B*[Us Is]';Uoc=X1(1)</p><p>　　Us2=0;Is2=1; %令I(lǐng)s=1A，并將Us置零，求Req=Ua

79、/Is=X2(1)</p><p>　　X2=A\B*[Us2 Is2]';Req=X2(1) %解出Req</p><p>　　RL=Req,Pmax=Uoc^2/(4*Req) %計(jì)算最大功率Pmax </p><p>　　RL=0:0.1:10;P=Uoc^2*RL./((Req+RL).*(Req+R

80、L)); %設(shè)置RL序列，并計(jì)算相應(yīng)功率</p><p>　　figure(1),plot(RL,P); %繪制功率隨RL變化的曲線</p><p>　　xlabel('RL(Ω)'),ylabel('P(W)'),grid</p><p>　　圖11 P隨RL變化的曲線&

81、lt;/p><p><b>　　心得體會(huì)</b></p><p>　　在本次的課設(shè)中我學(xué)到了很多，從開始的一無所知道現(xiàn)在能解決一點(diǎn)簡單的問題，以下是我的一點(diǎn)心得：</p><p>　　1.最好找一個(gè)熟悉編程的人來輔助你的學(xué)習(xí)，這就包括很多編程的技巧問題，程序的結(jié)構(gòu)設(shè)計(jì)問題，對(duì)于程序的運(yùn)行效率非常有幫助。有的時(shí)候，你編出來的程序，能夠運(yùn)行，但是耗時(shí)太

82、長，也就是說你的程序沒有錯(cuò)，但是不適合實(shí)際?；蛘哒f，對(duì)于規(guī)模小的問題能夠解決，但是規(guī)模大一點(diǎn)的問題就需要很長很長的時(shí)間，這就需要對(duì)程序的結(jié)構(gòu)和算法問題進(jìn)行改進(jìn)</p><p>　　2.需要找一本matlab的函數(shù)工具詞典，就像漢語詞典一樣，你要盡量多的熟悉matlab自帶的函數(shù)，及其作用，因?yàn)閙atlab的自帶函數(shù)特別多，基本上能夠滿足一般的數(shù)據(jù)和矩陣的計(jì)算，所以基本上不用你自己編函數(shù)（如vb中，大部分的函數(shù)都需

83、要自己編）。這一點(diǎn)對(duì)你的程序非常有幫助，可以使你的程序簡單，運(yùn)行效率高，可以節(jié)省很多時(shí)間。</p><p>　　3.你把基本的知識(shí)看過之后，就需要找一個(gè)實(shí)際的程序來動(dòng)手編一下，不要等所有的知識(shí)都學(xué)好之后再去編程，你要在編程的過程中學(xué)習(xí)，程序需要什么知識(shí)再去補(bǔ)充（這一條是別人教我的，很管用），編程是一點(diǎn)一點(diǎn)積累的，所以你要需做一些隨手筆記什么的。</p><p>　　4.編程問題最頭疼的不是

84、編程序，而是調(diào)程序，所以在你的程序編完之后，一定要進(jìn)行驗(yàn)證其正確性，你要盡量多的設(shè)想你的問題的復(fù)雜性，當(dāng)然，要一步一步復(fù)雜，這樣才能保證你的程序的適用性很強(qiáng)。</p><p>　　Matlab 確實(shí)是一個(gè)很強(qiáng)大的工具，這在我完成任務(wù)的時(shí)候逐步的認(rèn)識(shí)到了，以后還有多加以利用。這次課設(shè)不僅是一個(gè)任務(wù)，它不想僅讓我們學(xué)會(huì)了使用 Matlab，還讓我們有了一個(gè)自主學(xué)習(xí)與和別人交流的機(jī)會(huì)。而學(xué)會(huì)的東西都是對(duì)以后的專業(yè)學(xué)習(xí)很

85、有幫助的，相信在這個(gè)強(qiáng)大的工具下，讓以后的專業(yè)學(xué)習(xí)會(huì)更輕松。</p><p><b>　　5、參考文獻(xiàn)</b></p><p>　　[1].唐向宏，岳恒立，鄭雪峰,MATLAB及在電子信息類課程中的應(yīng)用,電子工業(yè)出版 社,2006.</p><p>　　[2].（美）穆爾,高會(huì)生,劉童娜,李聰聰,MATLAB實(shí)用教程(第二版),電

86、子工業(yè)出版 社,2010.</p><p>　　[3].陳杰,MATLAB寶典(第3版),電子工業(yè)出版社,2011.</p><p>　　[4].劉保柱,蘇彥華,張宏林.MATLAB 7.0從入門到精通(修訂版),人民郵電出版</p><p><b>　　社,2010.</b></p><p>　　[5].陳

87、懷琛,吳大正,高西全,MATLAB及在電子信息課程中的應(yīng)用,電子工業(yè)出版</p><p><b>　　社,2006.</b></p><p>　　[6].張琨.MATLAB 2010從入門到精通,電子工業(yè)出版社,2011.</p><p>　　[7].MATLAB基礎(chǔ)教程/薛山編著,清華大學(xué)出版社.2011.3</p><p