【正文】 ：自然科學版2010（1）.[16] ：自然科學版2010.[17]鐘麟、：國防工業(yè)出版社，2004.。寫作畢業(yè)論文是一次再系統(tǒng)學習的過程，畢業(yè)論文的完成，同樣也意味著新的學習生活的開始。最后再一次感謝所有在設計中曾經幫助過我的良師益友和同學。最后，我要感謝幫助過我的同學，他們也為我解決了不少我不太明白的理論上的難題。感謝我的同學和朋友，在論文的撰寫和排版等過程中提供給我熱情的幫助。感謝這篇論文所涉及到的各位學者。相信畢業(yè)論文答辯過后，我的理論水平和實踐動手能力都會有一個質的飛躍，也為今后的工作和學習做了鋪墊。因此我們在論文前必須做好充分的準備，認真查找詳細的資料，為我們論文的成功打下堅實的基礎。做畢業(yè)論文雖然很累，但卻是苦中有樂，讓人很充實。畢業(yè)論文恰恰給我提供了一個應用自己所學知識去驗證總結一些東西。能夠順利完成該畢業(yè)論文，令我覺得受益匪淺，也是大學期間收獲最大的一次實踐。另外，我們也可以嘗試使用其他計算機仿真軟件來編寫簡單實用的信號與系統(tǒng)實驗系統(tǒng)。不僅了解了MATLAB的強大的仿真功能，而且復習鞏固了之前學的電子電路技術的知識。我也熟悉了MATLAB的工作環(huán)境，可以很熟練的對MATLAB進行常規(guī)的操作，快速進行程序編輯和仿真。從上到下是解調波形、原始波形、調制波形。示例2：。MATLAB 在通信領域得到越來越廣泛的應用，運用MATLAB 通信仿真工具箱可以對通信系統(tǒng)中的調制與解調，信號的頻譜分析以及FIR 濾波器的設計等進行仿真。 仿真結構圖 仿真曲線圖 仿真曲線圖從上面例子可以看出，采用MATLAB 對自動控制原理中的內容進行仿真是很方便的，不用搭建任何電子物理器件就可以進行仿真，同時仿真結果以圖形的形式顯示。在SIMULINK 仿真界面上作仿真結構圖（）時域仿真，仿真參數(shù)可以進行設置修改，和現(xiàn)實用仿真儀器進行參數(shù)調整是一樣的，這樣就可以得到在不同情況下的仿真曲線了。此仿真環(huán)境一方面是利用SIMULINK 提供的輸入信號（信號源模塊）對結構圖所描述的系統(tǒng)施加激勵；另一方面是利用SIMULINK 提供的輸出裝置（輸出口模塊）獲得系統(tǒng)的輸出響應，即數(shù)據(jù)或時間響應曲線，成為圖形化、模塊化方式的控制系統(tǒng)仿真。中的參數(shù)值。 檢測示例圖 MATLAB 在自動控制原理中的應用MATLAB 語言開發(fā)了控制系統(tǒng)工具箱，即CONTROL 工具箱，并進一步設計了控制系統(tǒng)結構圖程序設計工具箱，即SIMULINK 工具箱，使得MATLAB 語言成為控制系統(tǒng)仿真分析的必要工具。下面以邊緣檢測為例來說明MATLAB 在數(shù)字圖像處理課程中的應用。所以在數(shù)字圖像處理課程中引進MATLAB軟件。 幅值譜分布上面的程序在MATLAB ，其中可以看到淹沒在時域噪聲中明顯的信號峰值，這就是頻域分析的功效。 含噪聲信號的幅值譜39。)。xlabel(39。subplot(2,2,4)。 不含噪聲信號的幅值譜39。)。xlabel(39。subplot(2,2,3)。 含噪聲信號39。)。xlabel(39。subplot(2,2,2)。 不含噪聲信號39。)。xlabel(39。subplot(2,2,1)。 Pn=abs(Yn)。Yn=fft(sn,512)。sn=signal+noise。signal=sin(2*pi*50*t)+cos(2*pi*100*t)。設采樣頻率f=1000Hz，信號頻率分布在50 ～ 100Hz，然后疊加隨機噪聲。通過上述語句，可見MATLAB 用于信號與系統(tǒng)的實驗仿真中是很方便的[15]。二維傅立葉變換是一個很復雜的數(shù)學運算，如果采用C 或Visual C＋＋語言編寫此變換的仿真程序，需要對傅立葉變換有一個深入的理解，而且需要將數(shù)學公式按步驟編寫程序，需要十幾條語句，但在MATLAB 中，就不需要把精力放在編寫程序上，只需要一條語句就可以實現(xiàn)傅立葉變換的仿真，語法格式如下：B＝fft2（I） （）B＝fft2（I，m，n） （）其中B＝fft2（I）返回數(shù)據(jù)I 的二維fft 變換矩陣，輸入數(shù)據(jù)I和輸出數(shù)據(jù)B 維數(shù)相同。5 MATLAB 在電子信息課程中的應用 MATLAB 在信號與系統(tǒng)中的應用信號與系統(tǒng)是電子信息類專業(yè)的一門重要的專業(yè)基礎課。 % 階躍響應s(t)。 % 系統(tǒng)函數(shù)分母多項式系數(shù)a=[1 ]。 % 繪制對數(shù)幅頻特性和對數(shù)相頻特性曲線[13]。 % 傳遞函數(shù) H(s)subplot(1,2,1),pzmap(sys)。 % 分母多項式系數(shù)a=[L*C R*C 1]。 % 電感L=50uHC=470*(10^12)。 零極點與相關響應圖，已知， 帶阻二階系統(tǒng)a. 畫出零極點圖b. 畫出幅頻特性和相頻特性(對數(shù))其中：，系統(tǒng)函數(shù)： （中心頻率）即： 和程序如下：R=50。 % 階躍響應s(t)subplot(1,3,3),impulse(b,a)。 % 傳遞函數(shù) H(s)subplot(1,3,1),pzmap(sys)。 % 系統(tǒng)函數(shù)分子多項式系數(shù)a=[1 4 5]。)。 % 繪制相頻特性曲線title(39。)。 % 繪制幅頻特性曲線title(39。H=1./(1normalizedw.^2+j*normalizedw/Q)。)。 周期矩形脈沖頻譜與幅度譜二階低通濾波器特性為：即：和令和1時，分別求幅頻特性和相頻特性程序如下： Q=input(39。axis([20 20 0 ])。周期矩形脈沖的幅度譜39。)。subplot(2,1,2),stem(n,abs(fn),39。周期矩形脈沖的頻譜39。 % 繪制頻譜line([20 20],[0 0])。filled39。)。 % 給定頻譜的整數(shù)自變量,eps代表0fn=subs(fn,n,39。 % 計算傅立葉系數(shù)fn=simple(fn)。f=A*exp(j*n*2*pi/T*t)。A=1。 信號調制波形如圖410所示周期矩形脈沖，試求其幅度譜。 % 定義符號變量tf=(2+2*sin(4*pi*t))*cos(50*pi*t) % 計算符號函數(shù)f(t)=(2+2*sin(4*pi*t))*cos(50*pi*t)ezplot(f,[0 pi])。[r p k]=residue(b,a)運行結果如下： r = + p = + 0 k = [] MATLAB在信號調制、濾波及求解中的應用由于信號與系統(tǒng)比較復雜，因此對于許多問題的求解都非常困難，然而MATLAB將會對問題的解答帶來很大的方便，下面會以具體事例給出展示[12]。 差分方程相應響應（7）部分分式展開：程序如下： b=[2 1]。)。title(39。filled39。)。title(39。filled39。 % 零狀態(tài)響應y2=impz(b,a,16)。 % 序列的個數(shù)fn=.^n。a=[1 1 ]。y(0)=0,Dy(0)=039。D2y+5*Dy+6*y=2*exp(t)39。 卷積和（4）求z變換 程序如下： syms n z OMEGA % 定義符號變量fn=cos(n*OMEGA)。 % 繪制y(n)的圖形axis([0 40 0 5])。filled39。 % 繪制x(n)的圖形axis([0 40 0 ])。filled39。 % 繪制f(n)的圖形axis([0 40 0 1])。filled39。 % 確定卷積序列的終點y=conv(f,h)。 % 確定h(n)序列值left=nf(1)+nh(1)。Nh=length(nh)。 % 確定f(n)的序號向量和區(qū)間長度f=.^(nf5)。 % 加注波形標題（3）求卷積和：若，到，求程序如下：nf=5:30。sinc(t)*cos(20*pi*t)39。 % 計算函數(shù)f(t)=sinc(t)*cos(20*pi*t)plot(t,f)。 信號相乘結果程序如下：t=5::5。 % 計算符號函數(shù)f(t)=cos(18*pi*t)+cos(20*pi*t)ezplot(f,[0 pi])。 信號相加結果程序如下：syms t。通過MATLAB可以對信號與系統(tǒng)進行各種處理，具體如下。z單位圓且，s與z之間是多對一的映射關系。z單位圓內；s右半平面224。若給定初始狀態(tài)，可由系統(tǒng)函數(shù)得齊次差分方程；進一步求得。響應的象函數(shù)224。z域框圖（零狀態(tài)）224。求解零狀態(tài)響應。系統(tǒng)函數(shù)224。5)利用ZT求解差分方程（零輸入響應、零狀態(tài)響應、全響應）差分方程利用單邊ZT的移位性質得到z域代數(shù)方程，整理成