matlab大作業(yè)Matlab在數(shù)字信號處理及圖像處理中的應(yīng)用

在數(shù)字信號處理及圖像處理中的應(yīng)用>>產(chǎn)生簡單信號程序如下:allall定義的取值范圍y1=sinc(t1;%調(diào)用sinc函數(shù)title('sinc函數(shù)';調(diào)用門函數(shù)title('門函數(shù)';axis（）定義坐標(biāo)軸的范圍y3=tripuls(t1,3;%調(diào)用三角脈沖函數(shù)title三角脈形沖'）axis（）;

（2）求模型的零輸入響應(yīng)+Y=[1.96916.4493]初始條件為程序如下a=[-0.55720];c=[1.969;0]%初始狀態(tài)x010狀態(tài)方程initial(sys,x0

(3）

濾波器設(shè)計

設(shè)計一個數(shù)字低通濾波器，抽樣頻率為1000HZ，濾波器的3db截至頻率為40HZ，阻帶截止頻率為１５０HZ，阻帶最小衰減為６０ｄｂ．繪制該濾波器的頻率響應(yīng)曲線。程序如下：Wp=40/500;%通帶截止頻率Ws=150/500;%阻帶截止頻率[n,Wn]=buttord(Wp,Ws,3,60%求butterworth濾波器的階數(shù)和帶的截止頻率[b,a]=butter(n,Wn;%調(diào)butter函數(shù)設(shè)計模擬濾波器freqz(b,a,512,1000;%頻率響應(yīng)圖title('巴特沃斯濾波器'

２.計一個ｃｈｅｂｙｓｈｅＩ數(shù)字低通濾波器，抽樣頻率為，濾波器的3db截至頻率為，阻帶截止頻率為１５０HZ，阻帶最小衰減為６０ｄｂ．繪制該濾波器的頻率響應(yīng)曲線。程序如下：Ws=150/500;通帶和阻帶的波紋系數(shù)[b,a]=cheby1(n,Rp,Wn;freqz(b,a,512,1000;title('切比雪夫濾波器'(4像處理1.實(shí)現(xiàn)兩圖像相減程序如下：I=imread('pepper.jpg';%用函數(shù)讀入圖片

Z=imsubtract(I,K;%用imsubtract函數(shù)實(shí)現(xiàn)兩圖片相減的功能imshow(Z%顯示相減后的圖片2.實(shí)現(xiàn)輛圖像相加程序如下：I=imread('pepper.jpg';%用函數(shù)讀入圖片K=imadd(I,J,'uint8';%用函數(shù)實(shí)現(xiàn)兩圖片相加的功能imshow(Ktitle('兩圖像相加'3.圖像切割

程序如下：i=imcrop(I,[40,30,50,60];%用imcrop函數(shù)實(shí)現(xiàn)切割圖片的功能subplot(121;imshow(I;title('辣椒原圖';subplot(122;title('切割后的圖';4.圖像旋轉(zhuǎn)程序如下：用imrotate函數(shù)實(shí)現(xiàn)圖片旋轉(zhuǎn)的功能subplot(121;imshow(I;title('辣椒原圖';subplot(122;

title('旋轉(zhuǎn)后的圖';實(shí)驗小結(jié)：通過做這次作業(yè)以后才發(fā)現(xiàn)自己對matlab軟件還不是很熟悉，看到別的同學(xué)做的東西都是很復(fù)雜而且功很強(qiáng)大，在自己做的大作業(yè)卻是相當(dāng)?shù)暮唵?，有時想一下真好慚愧。通過這次大業(yè)，自己又重新熟悉了matlab的一些功能，剛開始的時候本來想用來仿真一個通信系統(tǒng)的，但那軟件不太會用，最后就選擇了直接用matlab來理數(shù)字信號處理中的一些問題，做完以后又發(fā)現(xiàn)做的內(nèi)