【正文】 t(39。)。 % 傳遞函數(shù) H(s)subplot(1,3,1),pzmap(sys)。 % 分母多項(xiàng)式系數(shù)a=[L*C R*C 1]。 % 階躍響應(yīng)s(t)。設(shè)采樣頻率f=1000Hz，信號(hào)頻率分布在50 ～ 100Hz，然后疊加隨機(jī)噪聲。 Pn=abs(Yn)。 不含噪聲信號(hào)39。 含噪聲信號(hào)39。 不含噪聲信號(hào)的幅值譜39。 含噪聲信號(hào)的幅值譜39。 檢測(cè)示例圖 MATLAB 在自動(dòng)控制原理中的應(yīng)用MATLAB 語(yǔ)言開(kāi)發(fā)了控制系統(tǒng)工具箱，即CONTROL 工具箱，并進(jìn)一步設(shè)計(jì)了控制系統(tǒng)結(jié)構(gòu)圖程序設(shè)計(jì)工具箱，即SIMULINK 工具箱，使得MATLAB 語(yǔ)言成為控制系統(tǒng)仿真分析的必要工具。 仿真結(jié)構(gòu)圖 仿真曲線圖 仿真曲線圖從上面例子可以看出，采用MATLAB 對(duì)自動(dòng)控制原理中的內(nèi)容進(jìn)行仿真是很方便的，不用搭建任何電子物理器件就可以進(jìn)行仿真，同時(shí)仿真結(jié)果以圖形的形式顯示。我也熟悉了MATLAB的工作環(huán)境，可以很熟練的對(duì)MATLAB進(jìn)行常規(guī)的操作，快速進(jìn)行程序編輯和仿真。畢業(yè)論文恰恰給我提供了一個(gè)應(yīng)用自己所學(xué)知識(shí)去驗(yàn)證總結(jié)一些東西。感謝這篇論文所涉及到的各位學(xué)者。寫(xiě)作畢業(yè)論文是一次再系統(tǒng)學(xué)習(xí)的過(guò)程，畢業(yè)論文的完成，同樣也意味著新的學(xué)習(xí)生活的開(kāi)始。最后再一次感謝所有在設(shè)計(jì)中曾經(jīng)幫助過(guò)我的良師益友和同學(xué)。相信畢業(yè)論文答辯過(guò)后，我的理論水平和實(shí)踐動(dòng)手能力都會(huì)有一個(gè)質(zhì)的飛躍，也為今后的工作和學(xué)習(xí)做了鋪墊。能夠順利完成該畢業(yè)論文，令我覺(jué)得受益匪淺，也是大學(xué)期間收獲最大的一次實(shí)踐。從上到下是解調(diào)波形、原始波形、調(diào)制波形。在SIMULINK 仿真界面上作仿真結(jié)構(gòu)圖（）時(shí)域仿真，仿真參數(shù)可以進(jìn)行設(shè)置修改，和現(xiàn)實(shí)用仿真儀器進(jìn)行參數(shù)調(diào)整是一樣的，這樣就可以得到在不同情況下的仿真曲線了。下面以邊緣檢測(cè)為例來(lái)說(shuō)明MATLAB 在數(shù)字圖像處理課程中的應(yīng)用。)。)。)。)。Yn=fft(sn,512)。通過(guò)上述語(yǔ)句，可見(jiàn)MATLAB 用于信號(hào)與系統(tǒng)的實(shí)驗(yàn)仿真中是很方便的[15]。 % 系統(tǒng)函數(shù)分母多項(xiàng)式系數(shù)a=[1 ]。 % 電感L=50uHC=470*(10^12)。 % 系統(tǒng)函數(shù)分子多項(xiàng)式系數(shù)a=[1 4 5]。 % 繪制幅頻特性曲線title(39。axis([20 20 0 ])。周期矩形脈沖的頻譜39。 % 給定頻譜的整數(shù)自變量,eps代表0fn=subs(fn,n,39。 信號(hào)調(diào)制波形如圖410所示周期矩形脈沖，試求其幅度譜。)。title(39。a=[1 1 ]。 卷積和（4）求z變換 程序如下： syms n z OMEGA % 定義符號(hào)變量fn=cos(n*OMEGA)。filled39。 % 確定h(n)序列值left=nf(1)+nh(1)。sinc(t)*cos(20*pi*t)39。 信號(hào)相加結(jié)果程序如下：syms t。若給定初始狀態(tài)，可由系統(tǒng)函數(shù)得齊次差分方程；進(jìn)一步求得。系統(tǒng)函數(shù)224。幾類(lèi)序列的ROC：有限長(zhǎng)序列，右邊序列，左邊序列，雙邊序列。求解零狀態(tài)響應(yīng)。3)FT的性質(zhì)線性、奇偶性、對(duì)稱(chēng)性、尺度變換、時(shí)移、頻移、卷積定理（時(shí)域、頻域）。2)沖激響應(yīng)定義，求解（經(jīng)典法），注意應(yīng)用LTI系統(tǒng)零狀態(tài)響應(yīng)的微積分特性階躍響應(yīng)與的關(guān)系。i5=r(5)*sin(w*t+rdir(5))。x=[real(i1),real(i2),real(i3),real(i4),real(i5),real(i6)]。IS33]。Y23=1/R4。w=4。 % 產(chǎn)生x 向量y = sinx。例如 m u = abs 3 + 4i 。矩陣的輸入、輸出、轉(zhuǎn)置、加減、乘和方陣求逆、分塊矩陣的合成與分解等操作都十分方便。)。omega=linspace(0,pi,512)。運(yùn)行結(jié)果為Numerator ploynomial+000+000Denominator ploynomial+000+000+000Ap=As=。disp(‘Numerator ploynomial’)。w=linspace(0,pi,512)。As=25。數(shù)字濾波器從功能上分類(lèi)，可分為低通濾波器、高通濾波器、帶通濾波器、帶阻濾波器。 %給X軸加標(biāo)簽“時(shí)間/s”ylabel(39。.^(0:10)]39。在程序設(shè)計(jì)階段，為了對(duì)程序進(jìn)行調(diào)試或驗(yàn)證算法的正確性，需要一些特性已知的信號(hào)（簡(jiǎn)單的如正弦、方波、三角波等），所以可以由Matlab提供一些波形產(chǎn)生函數(shù)。 %確定自變量n的取值范圍x=cos(*pi*n)。用MATLAB表示一離散序列x[k]時(shí)，可用兩個(gè)向量來(lái)表示。然后在從MATLAB 命令窗口的File 菜單中選擇Open對(duì)話框，則屏幕出現(xiàn)Open對(duì)話框，在Open對(duì)話框中的File Name 框中輸入文件名，或從右邊的directories框中打開(kāi)這個(gè)M文件。用MATLAB語(yǔ)言編寫(xiě)的程序，稱(chēng)為M文件。Stateflow與Simulink和MATLAB緊密集成，可以將Stateflow創(chuàng)建的復(fù)雜控制邏輯有效地結(jié)合到Simulink的模型中[2]。 MATLAB軟件簡(jiǎn)介MATLAB的名稱(chēng)源自Matrix Laboratory，它是一種科學(xué)計(jì)算軟件，專(zhuān)門(mén)以矩陣的形式處理數(shù)據(jù)。實(shí)驗(yàn)中由于不了解電路及器件的工作原理只要稍有問(wèn)題學(xué)生就盲目詢(xún)問(wèn)老師，不利于培養(yǎng)學(xué)生自己解決問(wèn)題的能力。MATLAB作為一種適用于矩陣運(yùn)算及控制和信息處理的科學(xué)計(jì)算軟件，輸入簡(jiǎn)便，運(yùn)算高效和內(nèi)容豐富等特點(diǎn)在大學(xué)教學(xué)和科學(xué)研究中得到了全面的推廣。利用M語(yǔ)言還開(kāi)發(fā)了相應(yīng)的MATLAB專(zhuān)業(yè)工具箱函數(shù)供用戶直接使用。目前RTW支持生成標(biāo)準(zhǔn)的C語(yǔ)言代碼，并且具備了生成其他語(yǔ)言代碼的能力。命令文件對(duì)MATLAB工作空間的變量進(jìn)行操作，而且函數(shù)文件中定義的變量為局部變量，當(dāng)函數(shù)文件執(zhí)行完畢時(shí)，這些變量被清除。在編輯完成后，選擇File菜單中的Save命令可以把這個(gè)編輯過(guò)的M文件報(bào)存下來(lái)[6]。 %確定k的取值范圍x=[(k2)==0]。 余弦序列x[n]=cos()的圖像另外，MATLAB的工具箱中還提供了大量的信號(hào)產(chǎn)生函數(shù)，例如：sawtooth（產(chǎn)生鋸齒波或三角波信號(hào)）、square（產(chǎn)生方波信號(hào)）、chirp（產(chǎn)生調(diào)諧余弦信號(hào)）、gauspuls（產(chǎn)生高斯正弦脈沖信號(hào)）等等，調(diào)用這些函數(shù)可以方便地生成多種復(fù)雜信號(hào)[7]。例如：要產(chǎn)生一頻率為10kHz的周期高斯脈沖信號(hào)，其帶寬為50％。,10e3,)。如果濾波器的輸入、輸出都是離散時(shí)間信號(hào)，則該濾波器沖擊響應(yīng)也必然是離散的，這樣的濾波器定義為數(shù)字濾波器。下面，對(duì)于IIR濾波器和FIR濾波器的設(shè)計(jì)各舉一例以說(shuō)明。 Ws=Ws*Fs。numd=numd/norm。fprintf(‘%.4e\n’,dend)。a=[0,1]。xlabel(39。3 MATLAB在正弦穩(wěn)態(tài)電路分析中的應(yīng)用3. 1 矩陣計(jì)算與線性電路分析 矩陣工具引入電路理論已有半個(gè)多世紀(jì)的歷史。 2, 1, 1。利用函數(shù)ode23 和ode45可進(jìn)行電路瞬態(tài)分析。無(wú)論在理論研究還是實(shí)際應(yīng)用中，對(duì)于正弦穩(wěn)態(tài)電路的分析都是十分重要的，它是變壓器、交流電機(jī)以及電子電路的理論基礎(chǔ)，在實(shí)際應(yīng)用中，許多電氣設(shè)備的設(shè)計(jì)、性能指標(biāo)就是按正弦穩(wěn)態(tài)來(lái)考慮的，因此分析和計(jì)算正弦穩(wěn)態(tài)電路是工程技術(shù)和科學(xué)研究中常常會(huì)碰到的問(wèn)題。XL=w*L。IS22=US*(R2+i*XL)。U3=U(2)。direction=rdir*180/pir=strength*sqrt(2)t=0:pi/10000:。plot(t,i1,t,i2,t,i3,t,i4,t,i5,t,i6)。2)單位序列響應(yīng)的定義，的定義，求解（經(jīng)典法）；若方程右側(cè)是激勵(lì)及其移位序列時(shí)，注意應(yīng)用線性時(shí)不變性質(zhì)求解；階躍響應(yīng)與的關(guān)系。2)拉氏變換的性質(zhì)線性、尺度變換、時(shí)移、頻移；時(shí)域微分（1次、2次）、時(shí)域積分（1次）；時(shí)域卷積定理、初值終值定理。s域代數(shù)方程224。按照 極點(diǎn)的情況進(jìn)行部分分式展開(kāi)224。7)z域框圖k域框圖224。z單位圓外；s虛軸224。 % 繪制f(t)的波形（2）信號(hào)相乘。Nf=length(nf)。 % 計(jì)算f(n)和x(n)的卷積subplot(3,1,1),stem(nf,f,39。subplot(3,1,3),stem(left:right,y,39。,39。 % 輸入序列y1=filter(b,a,fn)。grid onsubplot(1,2,2),stem(n,y2,39。a=[1 2 5 0]。tao=1。 % 計(jì)算傅立葉系數(shù)對(duì)應(yīng)各個(gè)n的值subplot(2,1,1),stem(n,fn,39。filled39。輸入Q=39。gridsubplot(1,2,2),plot(normalizedw,angle(H))。 % 繪制零極點(diǎn)圖subplot(1,3,2),step(b,a)。 % 分子多項(xiàng)式系數(shù)sys=tf(b,a)。 階躍響應(yīng)以上分別包含了時(shí)域，頻域，離散域，復(fù)頻域里信號(hào)與系統(tǒng)的例子，從里面可以充分的體會(huì)到MATLAB這款軟件功能十分強(qiáng)大，用它來(lái)對(duì)信號(hào)與系統(tǒng)進(jìn)行仿真也極其方便簡(jiǎn)單，讓我們對(duì)這門(mén)課的理解也變得簡(jiǎn)單了許多[14]?？删帉?xiě)如下的MATLAB 程序演示信號(hào)疊加噪聲前后的幅值譜分布：t=0::。f=1000*(0:256)/512。)。