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