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

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

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

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

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

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

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è)計題目4</b></p><p><b>　　3、設(shè)計內(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、 心得體會22</p><p><b>　　摘要</b></p><p>　　MATLAB是一種高性能的工程計算語言。它集計算，可視化，可編程于一

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

12、x），它可以用于參數(shù)擬合，也可以用于非參數(shù)擬合。</p><p>　　在MATLAB中，多項式擬合函數(shù)polyfit可以用來計算擬合的得到多項式系數(shù)。默認(rèn)的擬合目標(biāo)是最小方差最小，即最小二乘法擬合數(shù)據(jù)。此時，判斷的依據(jù)過擬合曲線的道德原始數(shù)據(jù)之間的平均誤差是否達(dá)到最小。</p><p>　　關(guān)鍵詞：工程語言；數(shù)學(xué)計算 ；曲線擬合</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ā)展和完善，如今己成為覆蓋多個學(xué)科的國際公認(rèn)的最優(yōu)秀的數(shù)值計算仿真軟件。MATLAB具備強大的數(shù)值計算能力，許多復(fù)雜的計算問題只需短短幾行代碼就可在MATLAB中實現(xiàn)。作為一個跨平臺的軟件，MATLAB已推出Unix、Windows、Linux和Ma

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

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

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

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

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

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

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

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

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

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

31、</p><p>　　(4) 控制系統(tǒng)設(shè)計；</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ù)字信號處理；</p><p>　　(10) 財務(wù)、金融分析；</p><p>　　(11) 建模、仿真及樣機(jī)開發(fā)；</p><p>　　(12) 新算法研究開發(fā)；</p><p>　　(13) 圖形用戶界面設(shè)計。</p><p><b>　　2、

33、設(shè)計題目</b></p><p>　　計算y1=和y2=；</p><p>　　畫出衰減震蕩曲線及其他的包絡(luò)線,t的取值范圍是[0,6].</p><p>　　畫出所表示的三維曲面。x,y的取值范圍是[-9，9]。</p><p>　　分析下面每條指令的功能并運行，觀察執(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>　　計算在x=0.1與10處的值。</p><p>　　求函數(shù), n=100的值</p><p>　　求1000個元素的隨機(jī)數(shù)向量A中大于0.5的元素個數(shù)。</p><p>　　8.求線性方程組AX=b的根。</p><p>　　9.畫出y=|1000sin(4x)|+1的y軸為對數(shù)圖。</p><p&

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

40、<p><b>　　3、設(shè)計內(nèi)容</b></p><p><b>　　第一題</b></p><p>　　計算y1=和y2=；</p><p>　　分析：直接運行即可。</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的方程，運行程序。</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.分析下面每條指令的功能并運行，觀察執(zhí)行結(jié)果。</p><p>　　X=0:0.2:2; %X取值從0 到2 每隔0.2 取一個數(shù)</p><p>　　Y=X.*exp(-X); %Y對于X 的表達(dá)式,求出X 對應(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)各元素絕對值是否大于3,大于返回1,小于返回0,并將返回 值構(gòu)成與A相同行列的矩陣L</p><p>　　islogical(L) %判斷矩陣L 內(nèi)是否為邏輯值,是返回1,不是則返回0</p

52、><p>　　X=A(L) %將A 對應(yīng)L 內(nèi)為真的元素取出來,構(gòu)成列向量X,即絕對值大于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)生一個2 x 3 全1 矩陣A</p><p>　　B=ones(2) %產(chǎn)生一個行列為2 的全1 的方陣B</p><p>　　C=eye(4) %產(chǎn)生一個行列為4 單位方陣C</p><p>　　D=diag(C)

56、 %將C 的主對角元素取出形成列向量D</p><p>　　E=repmat(C,1,3) %將C 看成一個元素,重復(fù)C 構(gòu)成一個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>　　計算在x=0.1與10處的值。</p><p>　　分析：此題可以直接創(chuàng)建一個匿名函數(shù)進(jìn)行運算，也可以利用M文件計算。</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個元素的隨機(jī)數(shù)向量A中大于0.5的元素個數(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軸為對數(shù)圖。</p><p>　　分析：此題關(guān)鍵是用semillogy(x,y)函數(shù)將y坐標(biāo)軸表示成對數(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)軸是對數(shù)坐標(biāo)系</p><p>　　grid on;

71、 %給圖形加上網(wǎng)格</p><p>　　title('y=|1000sin(4x)|+1的x軸為對數(shù)的圖'); %添加標(biāo)題</p><p><b>　　圖9對數(shù)圖像</b></p><p><b>　　第十題</b></p><p>　　用FFT直接計算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點數(shù)</p><p>　　X=fft(x,tfft); %對 x 進(jìn)行快速傅里葉變換</p><p&g

74、t;　　Pxx=abs(X).^2; %功率譜的量度為變換后幅度的平方 </p><p>　　m=0:round(tfft/2-1); %由于對稱性,取一半</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>　　（1）負(fù)載電阻為何值時可吸收最大功率？求此最大功率。</p><p>　?。?）研究在0~10Ω范圍內(nèi)變化時，其吸收功率的情況。</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) %計算最大功率Pmax </p><p>　　RL=0:0.1:10;P=Uoc^2*RL./((Req+RL).*(Req+R

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

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

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

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

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

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