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