數字信號課程設計報-資料下載頁

　　

【正文】 text控件中獲取a0的值axes()。 n = 0:1:N1。k=[ceil(N/2):ceil(N/2)]。 xn=。 xh=xn*([n=0][(nN)=0])。 %只處理周期序列的一個周期的值Xk=fft(xh,N)。 %對周期序列的一個周期求DFS,這里用函數fft( )來求DFS，其實對于有限長的一個周期序列，DFS與DFT沒有本質區別magXk=abs([Xk(ceil(N/2+1):N) Xk(1:ceil(N/2+1))])。 %從Xk中分離出其幅度stem(k,magXk)。 title(39。幅度響應39。)。axes()。angXk=angle([Xk(ceil(N/2+1):N) Xk(1:ceil(N/2+1))])。 %從Xk中分離出其相位stem(k,angXk)。title(39。相位響應39。)。 信號譜分析的設計與實現信號譜分析是快速傅里葉變換（FFT）的一個重要應用。在利用FFT對模擬信號進行譜分析時，應將模擬信號離散化以得到離散時間信號，同時要考慮譜分析中參數的選擇。在本系統中，所處理的離散時間信號時從連續時間信號中取樣得到的，取樣周期可由用戶定義，的頻率分辨率也可以由用戶定義，連續時間信號可以從幾個典型的連續時間信號中選擇。①靜態界面 ②輸入值后如下圖所示（單位斜變信號）各個callback函數為單位階躍信號、單位斜變信號global T FT=str2num(get(,39。string39。))。F=str2num(get(,39。string39。))。M=ceil(1/(T*F))。 %對1/(T*F)取整N=2.^ceil(log(M)/log(2))。 %確定一個記錄長度中的取樣數val=get(hObject,39。Value39。)。switch valcase 1axes()。 t = [10:: ceil((N1)/10)*10]。y=udt(t)。 %連續階躍信號plot(t,y)。axes()。n=0:1:(N1)。r=n*T。 %對坐標n取樣x=udt(r)。 %取樣后的階躍序列stem(n,x)。case 2axes()。 t = [10:: ceil((N1)/10)*10]。y=t.*udt(t)。 %連續斜變信號plot(t,y)。 axes()。n=0:1:ceil((N1)/10)*10。r=n*T。 x=r.*udt(r)。 %取樣后的斜變序列stem(n,x)。end %將連續時間信號轉換成double型數據=double(x)39。guidata(hObject, handles)。譜分析global Nglobal T FT=str2num(get(,39。string39。))。F=str2num(get(,39。string39。))。xn=。M=ceil(1/(T*F))。N=2.^ceil(log(M)/log(2))。 X=fft(xn,N)。 %對離散時間信號xn 進行FFT變換，信號譜分析實際上是FFT分析k=0:N1。axes()。magX=abs(X)。stem(k,magX)。axes()。phaX=angle(X)。stem(k,phaX)。 %對周期序列的一個周期求DFS,這里用函數fft( )來求DFS，其實對于有限長的一個周期序列，DFS與DFT沒有本質區別①“數字濾波器的設計”的主界面如下所示 ②“巴特沃斯低通濾波器”的靜態界面 巴特沃斯低通濾波器的設計與實現①沖擊響應不變法設計 ② 雙線性變換法設計 各個callback函數為沖擊響應不變法global f p Rp Asf=str2num(get(,39。string39。))。 %將從 的edit text控件中獲取的值賦給f,下面三條語句的功能是類似的p=str2num(get(,39。string39。))。 Rp=str2num(get(,39。string39。))。Rs=str2num(get(,39。string39。))。wp=f*pi。 ws=p*pi。 T=1。Wp=wp/T。Ws=ws/T。 %將數字濾波器的指標轉換成模擬域的指標[N,Wn]=buttord(Wp,Ws,Rp,Rs,39。s39。)。[b1,a1]=butter(N,Wn,39。s39。)。[b,a]=impinvar(b1,a1,T)。 %imp_invr ()是自寫的沖激響應不變法設計IIR數字濾波器的函數[db,mag,pha,grd,w]=myfreqz(b,a)。 %計算數字濾波器頻率響應axes()。plot(w/pi,mag)。title(39。幅度響應39。)。 %畫出幅度響應axes()。plot(w/pi,db)。title(39。對數幅度39。)。 %畫出對數幅度axes()。plot(w/pi,pha)。title(39。相位響應39。)。 %畫出相位響應axes()。plot(w/pi,grd)。title(39。群延遲39。)。 %畫出群延遲雙線性變換法設計global f p Rp Rsf=str2num(get(,39。string39。))。 %將從 的edit text控件中獲取的值賦給f,下面三條語句的功能是類似的p=str2num(get(,39。string39。))。 Rp=str2num(get(,39。string39。))。Rs=str2num(get(,39。string39。))。Wp=f*pi。Ws=p*pi。T=1。Fs=1/T。omegap=(2/T)*tan(Wp/2)。 %將數字濾波器的指標利用雙線性變換法轉換成模擬域的指標omegas=(2/T)*tan(Ws/2)。[N,Wn]=buttord(Wp,Ws,Rp,Rs,39。s39。)。[b1,a1]=butter(N,Wn,39。s39。)。[b,a]=bilinear(b1,a1,)。 %bilinear()是MATLAB內部提供的實現雙線性變換法設計的函數[db,mag,pha,grd,w]=myfreqz(b,a)。axes()。plot(w/pi,mag)。title(39。幅度響應39。)axes()。plot(w/pi,db)。title(39。對數幅度39。)axes()。plot(w/pi,pha)。title(39。相位響應39。)。axes()。plot(w/pi,grd)。title(39。群延遲39。)。 切比雪夫I型低通濾波器的設計與實現①沖擊響應不變法設計 ② 雙線性變換法設計 各個callback函數為沖擊響應不變法設計global f p Rp Asf=str2num(get(,39。string39。))。 %將從 的edit text控件中獲取的值賦給f,下面三條語句的功能是類似的p=str2num(get(,39。string39。))。 Rp=str2num(get(,39。string39。))。Rs=str2num(get(,39。string39。))。wp=f*pi。 ws=p*pi。 T=1。Wp=wp/T。Ws=ws/T。 %將數字濾波器的指標轉換成模擬域的指標[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,39。s39。)。[b1,a1]=cheby1(N,Rp,Wn,39。s39。)。[b,a]=impinvar(b1,a1,T)。 %imp_invr ()是自寫的沖激響應不變法設計IIR數字濾波器的函數[db,mag,pha,grd,w]=myfreqz(b,a)。 %計算數字濾波器頻率響應axes()。plot(w/pi,mag)。title(39。幅度響應39。)。 %畫出幅度響應axes()。plot(w/pi,db)。title(39。對數幅度39。)。 %畫出對數幅度axes()。plot(w/pi,pha)。title(39。相位響應39。)。 %畫出相位響應axes()。plot(w/pi,grd)。title(39。群延遲39。)。 %畫出群延遲雙線性變換法設計global f p Rp Rsf=str2num(get(,39。string39。))。 %將從 的edit text控件中獲取的值賦給f,下面三條語句的功能是類似的p=str2num(get(,39。string39。))。 Rp=str2num(get(,39。string39。))。Rs=str2num(get(,39。string39。))。Wp=f*pi。Ws=p*pi。T=1。Fs=1/T。omegap=(2/T)*tan(Wp/2)。 %將數字濾波器的指標利用雙線性變換法轉換成模擬域的指標omegas=(2/T)*tan(Ws/2)。[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,39。s39。)。[b1,a1]=cheby1(N,Rp,Wn,39。s39。)。[b,a]=bilinear(b1,a1,)。 %bilinear()是MATLAB內部提供的實現雙線性變換法設計的函數[db,mag,pha,grd,w]=myfreqz(b,a)。axes()。plot(w/pi,mag)。title(39。幅度響應39。)axes()。plot(w/pi,db)。title(39。對數幅度39。)axes()。plot(w/pi,pha)。title(39。相位響應39。)。axes()。plot(w/pi,grd)。title(39。群延遲39。)。3退出close(gcf)。untitledA(handles) 窗函數法設計FIR低通濾波器的設計與分析“窗函數法設計FIR低通數字濾波器”的界面列出了六種窗函數，用戶在確定各參數之后，可以點擊相應的窗函數按鈕進行設計。 窗函數法設計FIR數字濾波器是從時域出發，把理想的濾波器的單位取樣響應用合適的窗函數截短成為有限長的，并使逼近理想的，以實現所設計的濾波器的頻率響應逼近于理想濾波器的頻率響應。其靜態界面和輸入之后如下圖所示列如當為矩形窗設計時候：（其他窗也可行，就不一一枚舉） 各個callback函數為矩形窗設計global f p N Rp Asf=str2num(get(,39。string39。))。p=str2num(get(,39。string39。))。N=str2num(get(,39。string39。))。Rp=str2num(get(,39。string39。))。As=str2num(get(,39。string39。))。wp=f*pi。 ws=p*pi。 n=[0:1:N1]。wc=(ws+wp)/2。 %理想低通的截止頻率hd = ideallp(wc,N)。 %理想低通的沖激響應w_ham=(rectwin(N))39。 %哈明窗函數h=hd.*w_ham。 %FIR濾波器的沖激響應[db,mag,pha,grd,w]=myfreqz(h,[1])。delta_w=2*pi/1000。 Rp=(min(db(1:1:wp/delta_w+1)))。 %實際的通帶衰減As=round(max(db(ws/delta_w+1:1:501)))。 %實際的最小阻帶衰減axes()。 stem(n,hd)。title(39。理想沖激響應39。)。 %將理想沖激響應繪制在坐標1（axes1）內axes()。 stem(n,w_ham)。title(39。矩形窗39。)。 %將哈明窗函數繪制在坐標2（axes2）內axes()。 stem(n,h)。title(39。實際沖激響應39。)。 %將實際沖激響應繪制在坐標3（axes3）內axes()。 plot(w/pi,db)。title(39。對數幅度響應39。)。 %將對數幅度響應繪制在坐標4（axes4）內axes()。 plot(w/pi,pha)。title(39。相位響應39。)。 %將相位響應繪制在坐標5（axes5）內axes()。 plot(w/pi,mag)。title(39。幅度響應39。)。 %將幅度響應繪制在坐標6（axes6）內漢明窗設計global f p N Rp Asf=str2num(get(,39。string39。))。p=str2num(get(,39。string39。))。N=str2num(get(,39。string39。)