【正文】 )。N=str2num(get(,39。string39。))。 %將幅度響應(yīng)繪制在坐標(biāo)6（axes6）內(nèi)漢明窗設(shè)計(jì)global f p N Rp Asf=str2num(get(,39。幅度響應(yīng)39。 plot(w/pi,mag)。)。title(39。 %將對(duì)數(shù)幅度響應(yīng)繪制在坐標(biāo)4（axes4）內(nèi)axes()。對(duì)數(shù)幅度響應(yīng)39。 plot(w/pi,db)。)。title(39。 %將哈明窗函數(shù)繪制在坐標(biāo)2（axes2）內(nèi)axes()。矩形窗39。 stem(n,w_ham)。)。title(39。 %實(shí)際的最小阻帶衰減axes()。 Rp=(min(db(1:1:wp/delta_w+1)))。 %FIR濾波器的沖激響應(yīng)[db,mag,pha,grd,w]=myfreqz(h,[1])。 %理想低通的沖激響應(yīng)w_ham=(rectwin(N))39。wc=(ws+wp)/2。 ws=p*pi。))。As=str2num(get(,39。string39。))。N=str2num(get(,39。string39。))。其靜態(tài)界面和輸入之后如下圖所示列如當(dāng)為矩形窗設(shè)計(jì)時(shí)候：（其他窗也可行，就不一一枚舉） 各個(gè)callback函數(shù)為矩形窗設(shè)計(jì)global f p N Rp Asf=str2num(get(,39。untitledA(handles) 窗函數(shù)法設(shè)計(jì)FIR低通濾波器的設(shè)計(jì)與分析“窗函數(shù)法設(shè)計(jì)FIR低通數(shù)字濾波器”的界面列出了六種窗函數(shù)，用戶在確定各參數(shù)之后，可以點(diǎn)擊相應(yīng)的窗函數(shù)按鈕進(jìn)行設(shè)計(jì)。)。title(39。axes()。相位響應(yīng)39。plot(w/pi,pha)。對(duì)數(shù)幅度39。plot(w/pi,db)。幅度響應(yīng)39。plot(w/pi,mag)。 %bilinear()是MATLAB內(nèi)部提供的實(shí)現(xiàn)雙線性變換法設(shè)計(jì)的函數(shù)[db,mag,pha,grd,w]=myfreqz(b,a)。)。[b1,a1]=cheby1(N,Rp,Wn,39。s39。 %將數(shù)字濾波器的指標(biāo)利用雙線性變換法轉(zhuǎn)換成模擬域的指標(biāo)omegas=(2/T)*tan(Ws/2)。Fs=1/T。Ws=p*pi。))。Rs=str2num(get(,39。string39。))。 %將從 的edit text控件中獲取的值賦給f,下面三條語句的功能是類似的p=str2num(get(,39。string39。)。title(39。 %畫出相位響應(yīng)axes()。相位響應(yīng)39。plot(w/pi,pha)。)。title(39。 %畫出幅度響應(yīng)axes()。幅度響應(yīng)39。plot(w/pi,mag)。 %imp_invr ()是自寫的沖激響應(yīng)不變法設(shè)計(jì)IIR數(shù)字濾波器的函數(shù)[db,mag,pha,grd,w]=myfreqz(b,a)。)。[b1,a1]=cheby1(N,Rp,Wn,39。s39。Ws=ws/T。 T=1。wp=f*pi。string39。))。 Rp=str2num(get(,39。string39。))。 切比雪夫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。群延遲39。plot(w/pi,grd)。)。title(39。)axes()。title(39。)axes()。title(39。axes()。[b,a]=bilinear(b1,a1,)。s39。)。[N,Wn]=buttord(Wp,Ws,Rp,Rs,39。omegap=(2/T)*tan(Wp/2)。T=1。Wp=f*pi。string39。))。 Rp=str2num(get(,39。string39。))。 %畫出群延遲雙線性變換法設(shè)計(jì)global f p Rp Rsf=str2num(get(,39。群延遲39。plot(w/pi,grd)。)。title(39。 %畫出對(duì)數(shù)幅度axes()。對(duì)數(shù)幅度39。plot(w/pi,db)。)。title(39。 %計(jì)算數(shù)字濾波器頻率響應(yīng)axes()。[b,a]=impinvar(b1,a1,T)。s39。)。 %將數(shù)字濾波器的指標(biāo)轉(zhuǎn)換成模擬域的指標(biāo)[N,Wn]=buttord(Wp,Ws,Rp,Rs,39。Wp=wp/T。 ws=p*pi。))。Rs=str2num(get(,39。string39。))。 %將從 的edit text控件中獲取的值賦給f,下面三條語句的功能是類似的p=str2num(get(,39。string39。stem(k,phaX)。axes()。magX=abs(X)。 %對(duì)離散時(shí)間信號(hào)xn 進(jìn)行FFT變換，信號(hào)譜分析實(shí)際上是FFT分析k=0:N1。N=2.^ceil(log(M)/log(2))。xn=。string39。))。譜分析global Nglobal T FT=str2num(get(,39。end %將連續(xù)時(shí)間信號(hào)轉(zhuǎn)換成double型數(shù)據(jù)=double(x)39。 x=r.*udt(r)。n=0:1:ceil((N1)/10)*10。 %連續(xù)斜變信號(hào)plot(t,y)。 t = [10:: ceil((N1)/10)*10]。 %取樣后的階躍序列stem(n,x)。r=n*T。axes()。y=udt(t)。switch valcase 1axes()。Value39。 %對(duì)1/(T*F)取整N=2.^ceil(log(M)/log(2))。))。F=str2num(get(,39。string39。在本系統(tǒng)中，所處理的離散時(shí)間信號(hào)時(shí)從連續(xù)時(shí)間信號(hào)中取樣得到的，取樣周期可由用戶定義，的頻率分辨率也可以由用戶定義，連續(xù)時(shí)間信號(hào)可以從幾個(gè)典型的連續(xù)時(shí)間信號(hào)中選擇。 信號(hào)譜分析的設(shè)計(jì)與實(shí)現(xiàn)信號(hào)譜分析是快速傅里葉變換（FFT）的一個(gè)重要應(yīng)用。相位響應(yīng)39。 %從Xk中分離出其相位stem(k,angXk)。axes()。幅度響應(yīng)39。 %從Xk中分離出其幅度stem(k,magXk)。 %只處理周期序列的一個(gè)周期的值Xk=fft(xh,N)。 xn=。 n = 0:1:N1。))?？焖俑道锶~變換FFTN=str2num(get(,39。退出close(gcf)。fill39。 y=ifft(fft(x1,N).*fft(x2,N))。axes()。x1=。string39。)。 str = get(hObject, 39。Value39。 guidata(hObject,handles)。 %產(chǎn)生一個(gè)長度為L的三角形序列stem(n,y1)。 stem(n,y1) case 5 %三角形序列n=[0:1:20]。 stem(n,y1) case 4 %遞增序列 n=(0:1:20)。 stem(n,y1) case 3 %正弦序列 n=[0::N1]。 stem(n,y1) case 2 %實(shí)指數(shù)序列 n=0:N1。 %y1=[ones(1,5),zeros(1,N4)]。axes()。String39。)。val=get(hObject,39。string39。循環(huán)卷積的具體的實(shí)現(xiàn)方法見“循環(huán)卷積”按鈕的callback函數(shù)。N實(shí)際是兩個(gè)序列的實(shí)際長度，當(dāng)兩個(gè)序列的長度小于N時(shí)，須補(bǔ)0延長至N；當(dāng)兩序列的長度大于N時(shí)，系統(tǒng)會(huì)報(bào)錯(cuò)，因?yàn)檫@是不符合邏輯的，所以用戶在輸入LLN時(shí)，一定要注意（LL2）必須小于或等于N。循環(huán)卷積要求進(jìn)行卷積的兩個(gè)序列必須長度相同，卷積后的結(jié)果的長度與進(jìn)行卷積的序列的長度也是相同的。退出close(gcf)。 %從Xk中分離出其相位stem(k,angXk)。 axes()。 %對(duì)周期序列的一個(gè)周期求DFS,這里用函數(shù)fft( )來求DFS，其實(shí)對(duì)于有限長的一個(gè)周期序列，DFS與DFT沒有本質(zhì)區(qū)別magXk=abs([Xk(ceil(N/2+1):N) Xk(1:ceil(N/2+1))])。 xh=xn*(udt(n)udt(nN))。k=[ceil(N/2):ceil(N/2)]。 %從的edit text控件中獲取a0的值axes()。string39。guidata(hObject, handles)。stem(xtilde)。*ones(1,5)。x=cos(2*pi/N*n)。stem(xtilde)。*ones(1,5)。(nN)]。case 5n=[0:1:N1]。xtilde=(xtilde(:))39。 %產(chǎn)生一個(gè)長度為L的三角形序列xtilde=x39。case 4n=[0:1:N1]。 xtilde=(xtilde(:))39。 %產(chǎn)生長度為L的遞減序列xtilde=x39。 case 3n=[0:1:N1]。 %將序列x周期化，周期為Nxtilde=(xtilde(:))39。 %產(chǎn)生一個(gè)長度為L的遞增序列xtilde=x39。 %畫出周期序列case 2n=[0:1:N1]。x=x1x2。x1=[n=0]。axes()。Value39。axes()。string39。))。 離散傅里葉變換(DFT)的設(shè)計(jì)與實(shí)現(xiàn)①“離散傅里葉變換(DFT)”的靜態(tài)界面如下所示 ②輸入結(jié)果后，如實(shí)指數(shù)序列： 各個(gè)callback函數(shù)為矩形序列、遞增序列、遞減序列、三角序列、實(shí)指數(shù)序列、正弦序列g(shù)lobal N LN=str2num(get(,39。相位響應(yīng)39。 %從Xk中分離出其相位stem(k,angXk)。axes()。幅度響應(yīng)39。 stem(k,magXk)。 %對(duì)周期序列的一個(gè)周期求DFS,這里用函數(shù)fft( )來求DFS，其實(shí)對(duì)于有限長的一個(gè)周期序列，DFS與DFT沒有本質(zhì)區(qū)別magXk=abs([Xk(ceil(N/2+1):N) Xk(1:ceil(N/2+1))])。 xh=xn*([n=0][n=N])。k=[ceil(N/2):ceil(N/2)]。))。離散傅里葉級(jí)數(shù)DFSN=str2num(get(,39。end = double(x)39。正弦型序列39。stem(xtilde)。xtilde=(xtilde(:))39。 %產(chǎn)生一個(gè)長度為L的正弦型序列 xtilde=x39。case 6n=[0:1:N1]。實(shí)指數(shù)序列39。stem(xtilde)。xtilde=(xtilde(:))39。 %產(chǎn)生一個(gè)長度為L的實(shí)指數(shù)序列 xtilde=x39。case 5n=[0:1:N1]。三角序列39。stem(xtilde)。xtilde=(xtilde(:))39。 %產(chǎn)生一個(gè)長度為L的三角形序列xtilde=x39。case 4n=[0:1:N1]。遞減序列39。stem(xtilde)。 xtilde=(xtilde(:))39。 %產(chǎn)生長度為L的遞減序列xtilde=x39。case 3n=[0:1:N1]。遞增序列39。stem(xtilde)。 %將序列x周期化，周期為Nxtilde=(xtilde(:))39。 %產(chǎn)生一個(gè)長度為L的遞增序列xtilde=x39。case 2n=[0:1:N1]。矩形序列39。stem(n,xtilde)。 n = [N :1: (3.*N1)]。*ones(1,4)。 %獲得句柄的String(內(nèi)容)屬性switch valcase 1 x=[ones(1,L),zeros(1,NL)]。String39。)。val = get(hObject,39。string39。))。①靜態(tài)界面： ② 輸入之后結(jié)果如下圖所示：列如選擇矩形序列各個(gè)callback函數(shù)如下：矩形序列、遞增序列、遞減序列、三角序列、實(shí)指數(shù)序列、正弦序列g(shù)lobal n xtildeN=str2num(get(,39。輸入之后結(jié)果如下圖1所示 “離散時(shí)傅里葉變換及其快速算法”的靜態(tài)界面如上圖所2示： 離散傅里葉級(jí)數(shù)(DFS)的設(shè)計(jì)與實(shí)現(xiàn)“離散傅里葉級(jí)數(shù)(DFS)”的靜態(tài)界面如下圖所示。,xbh)。 set(,39。xbh=simplify(x)。x=iztrans(Z)。grid onsyms z。fill39。