【正文】 由于我的學術水平有限，所寫論文難免有不足之處，懇請各位老師和學友批評和指正。致謝在這次畢業(yè)設計中，我首先要感謝學校給我們提供的機會，其次我要感謝我的老師蔡艷艷老師，在畢業(yè)設計上給予我的指導、提供給我的支持和幫助，這是我能順利完成這次論文的主要原因，更重要的是老師幫我解決了許多理論上的難題，讓我能把系統(tǒng)做得更加完善。身為電子信息工程專業(yè)的學生，論文是我們必備的一步。 總 結通過這次畢業(yè)論文，我認識到了MATLAB功能非常的強大，使得我們在使用的時候直接調用這些庫函數(shù)并賦予實際參數(shù)就能解決實際問題，具有極高的變成效率。如果仿真參數(shù)設置為：階躍時間為1 s，初始值為0，終止值為1，數(shù)設置為：階躍時間為2 s，初始值為0，終止值為10，[16]。在數(shù)字圖像處理中常用的邊緣檢測算子有Sobel 算子、Robert 算子、Prewitt 算子、LOG 算子、Canny 算子，在MATLAB 中只需要將語句BW＝edge （I，’sobel’） 中的sobel 修改為roberts，prewitt，log 和canny 算子就可以了，這樣就可以對同一幅圖像進行不同的邊緣檢測，各個算子檢測的結果以圖形的形式顯示，程序如下：I＝imread（‘’）；subplot（2，3，1），imshow（I）title（‘原始圖像’）BW1＝edge（I， ‘sobel’）；BW2＝edge（I， ‘roberts’）；BW3＝edge（I， ‘prewitt’）；BW4＝edge（I， ‘log’）；subplot（2，3，2），imshow（BW1）title（‘sobel 算子’）subplot（2，3，3），imshow（BW2）title（‘roberts 算子’）subplot（2，3，4），imshow（BW3）title（‘prewitt 算子’）subplot（2，3，5），imshow（BW4）title（‘log 算子’）subplot（2，3，6），imshow（BW5）title（‘canny 算子’）。 title(39。 title(39。 title(39。 title(39。Ps=abs(Ys)。示例1：通過傅里葉變換找出被隨機噪聲淹沒的信號頻率。 % 系統(tǒng)函數(shù)分子多項式系數(shù)step(b,a)。 % 電容C=470pFb=[L*C 0 1]。 % 系統(tǒng)函數(shù)分母多項式系數(shù)sys=tf(b,a)。幅頻特性曲線39。)。n39。 周期矩形脈沖程序如下：clear allsyms t n T tao A T=4。grid on。零狀態(tài)響應39。 % 差分方程的系數(shù)n=0:15。 % 定義f(n)fz=ztrans(fn,n,z) % 對f(n)進行z變換運行結果如下：fz =(zcos(OMEGA))*z/(z^22*z*cos(OMEGA)+1)（5）設方程 ，試求零狀態(tài)響應 程序如下： yzs=dsolve(39。)。 % 確定卷積序列的起點right=nf(Nf)+nh(Nh)。)。 % 定義符號變量tf=cos(18*pi*t)+cos(20*pi*t)。8)S域與z域的關系s左半平面224。求；系統(tǒng)函數(shù)轉化為差分方程。3)ZT的性質線性、移位性質（單邊右移）、z域尺度、k域卷積定理、k域反轉、部分和、初值終值定理（因果序列）。5)s域框圖時域框圖224。4)系統(tǒng)的頻率響應周期信號輸入，可用傅立葉級數(shù)法 。3)卷積積分定義：激勵、零狀態(tài)響應、沖激響應之間關系。i6=r(6)*sin(w*t+rdir(6))。y=[imag(i1),imag(i2),imag(i3),imag(i4),imag(i5),imag(i6)]。U=Y\I。Y32=Y23。L=4。 % 求y 向量 ph = plotx ,y。 % 求復數(shù)模 f uj = angle3 + 4i 。例如 A = 1, 1, 1。grid。mag=freqz(h,[1],omega)。 例一數(shù)字濾波器的增益響應例2：試用Kaiser窗設計滿足下列指標的FIR高通濾波器Ωp=，Ωs=，δs=解：設計滿足上述指標數(shù)字濾波器的MATLAB程序如下：%Kaiser窗設計FIR高通濾波器Rs=。fprintf(‘%.4e\n’,numd)。h=freqz(numd,dend,w)。Fs=1。從濾波器的網絡結構或者從單位脈沖響應分類，可分為IIR濾波器（即無限長單位沖擊響應濾波器）和FIR濾波器（即有限長單位沖擊響應濾波器）。幅值39。 y=pulstran(t,d,39。例如在上節(jié)中提到的方波發(fā)生器square（）、鋸齒波和三角波發(fā)生器sawtooth（）等[8]。 stem(x) %作圖xlabel(‘n’)。其中一個向量表示自變量k的取值范圍，另一個向量表示序列x[k]的值，在命令窗口直接輸入表示兩個向量的命令語句即可。在M文件所在的目錄，再從File Name 下面的列表框中選中這個文件，然后按OK按鈕即打開這個M文件。M文件有兩類：命令文件和函數(shù)文件。在MATLAB產品族中，自動化的代碼生成工具主要有RealTime Workshop（RTW）和Stateflow Coder，這兩種代碼生成工具可以直接將Simulink的模型框圖和Stateflow的狀態(tài)圖轉換成高效優(yōu)化的程序代碼。MATLAB將高性能的數(shù)值計算和可視化集成在一起，并提供了大量的內置函數(shù)，從而被廣泛地應用于科學計算、控制系統(tǒng)、信息處理等領域的分析、仿真和設計工作，而且利用MATLAB產品的開放式結構，可以非常容易地對MATLAB的功能進行擴充，從而在不斷深化對問題認識的同時，不斷完善MATLAB產品以提高產品自身的競爭能力。由于電力電子器件自身的非關聯(lián)性，給電力電子電路和系統(tǒng)的分析了一定的復雜性和困難，一般常用波形分析和分段性化處理的方法來研究電力電子字電路?，F(xiàn)代計算機仿真技術為電力電子電路和系統(tǒng)的分析提供了嶄新的方法，可以使復雜的電力電子系統(tǒng)的分析和設計變得更加容易和有效。MATLAB是MATLAB產品家族的基礎，它提供了基本的數(shù)學算法，例如矩陣運算、數(shù)值分析算法，MATLAB集成了2D和3D圖形功能，以完成相應數(shù)值可視化的工作，并且提供了一種交互式的高級編程語言——M語言，利用M語言可以通過編寫腳本或者函數(shù)文件實現(xiàn)用戶自己的算法[1]。利用RTW生成的代碼簡潔、可靠、易讀。兩者區(qū)別在于：命令文件沒有輸入?yún)?shù)，也不返回輸出參數(shù)；而函數(shù)文件可以輸入?yún)?shù)，也可以返回輸出參數(shù)。在M文件窗口可以對打開的M文件進行編輯修改。例如利用MATLAB表示單位脈沖序列δ[k2]在4≤k≤4范圍內各點的取值，可在命令窗口輸入如下語句：k=[4:4]。 %在X軸添加標簽：n。此外，對于采樣處理，MATLAB也提供了一些簡單的處理函數(shù)，例如函數(shù)resample（）可用于改變信號的采樣率、函數(shù)decimat（）用于經低通濾波后信號的下采樣以及函數(shù)interp（）用于經低通濾波后信號的上采樣等等。gauspuls39。) %給Y軸加標簽“幅值” MATLAB在數(shù)字濾波器中的應用濾波器是指用來對輸入信號進行濾波的硬件或軟件。對于每種濾波器的設計，方法也不盡相同，在此不予贅述，可查閱相關參考文獻。 %抽樣頻率（Hz）%確定模擬BW指標Wp=WP*Fs。%幅度歸一化DF的幅度響應norm=max(abs(h))。disp(‘Denominator ploynomial’)。f=[,]。plot(omega/pi,20*log10(abs(mag)))。 5, 4, 3。 % 求復數(shù)幅角3. 2 微分方程求解 MATLAB提供了常微分方程初值問題的數(shù)值解法。 %繪制正弦曲線 簡單的正弦穩(wěn)態(tài)分析計算 正弦穩(wěn)態(tài)電路圖介紹正弦激勵的動態(tài)電路中，若各電壓、。C=3。Y33=1/R1+1/R41/(i*XC)。U2=U(1)。[rdir strength]=cart2pol(x,y)。figure 。（2）離散系統(tǒng)的時域分析1)離散系統(tǒng)的響應差分方程的迭代法求解；差分方程的經典法求解：齊次解+特解（代入初始條件求系數(shù)）；全響應=零輸入響應+ 零狀態(tài)響應。（4）連續(xù)系統(tǒng)的S域分析1)單邊拉普拉斯變換的定義及ROC ROC： 幾個結論。s域框圖（零狀態(tài)）224。4)逆z變換的求解長除法、部分分式展開法、留數(shù)定理法；重點：部分分式展開法；步驟：224。求解零狀態(tài)響應。z單位圓內；s右半平面224。 % 計算符號函數(shù)f(t)=cos(18*pi*t)+cos(20*pi*t)ezplot(f,[0 pi])。 % 加注波形標題（3）求卷積和：若，到，求程序如下：nf=5:30。 % 確定卷積序列的終點y=conv(f,h)。 % 繪制x(n)的圖形axis([0 40 0 ])。D2y+5*Dy+6*y=2*exp(t)39。 % 序列的個數(shù)fn=.^n。)。 差分方程相應響應（7）部分分式展開：程序如下： b=[2 1]。A=1。)。subplot(2,1,2),stem(n,abs(fn),39。 周期矩形脈沖頻譜與幅度譜二階低通濾波器特性為：即：和令和1時，分別求幅頻特性和相頻特性程序如下： Q=inpu