【正文】 x(n)39。ylabel(39。) %畫出復(fù)指數(shù)序列 x的幅度 subplot(3,3,4)。 title(39。)。x(n)39。ylabel(39。)。string39。復(fù)指數(shù)序列定義為 x(n)=exp(an+jwn)39。這里 a 被稱為阻尼系數(shù)， w 為數(shù)字域頻率，單位為弧度 39。復(fù)指數(shù)序列的幅度 |x(n)|=exp(an)，相角 arg[x(n)]=wn39。注意：當 w為π的倍數(shù)時，輸入時π應(yīng)寫成 pi，乘用 *表示，除用 /表示 39。HorizontalAlignment39。left39。FontSize39。 a=str2num(get(,39。))。 subplot(3,3,1)。 title(39。)。x(n)39。ylabel(39。)。 a=str2num(get(,39。))。string39。 p=str2num(get(,39。))。 subplot(3,3,1)。 title(39。)。x(n)/pi39。ylabel(39。)。在靜態(tài)界面搭建好之后，編寫各個控件的 callback函數(shù)，就可以實現(xiàn)界面的功能。各個 callback 函數(shù)如下所示： ① 卷積： A=str2num(get(,39。))。 x=。 y1=conv(x,h)。 ny=length(y1)。 axes()。 ② 矩形序列和實指數(shù)序列 global N N=str2num(get(,39。))。Value39。 %獲得句柄的 value(值 ) 屬性 str = get(hObject, 39。)。 %設(shè)置當前圖形繪制在坐標 1（ axes1）中 switch val case 1 n=[0: N1]。 x2=[n=N]。 17 axes()。 %畫出矩形序列 y1，并且矩形的長度 N是由用戶規(guī)定的 case 2 n=[0: N1]。 axes()。 end ③ 退出 close(gcf)。在靜態(tài)界面設(shè)計好之后，編寫各個按鈕的 callback函數(shù)，就可以實現(xiàn)界面的功能。Value39。 str=get(hObject,39。)。 18 switch val case 1 n=[10:20]。 stem(n,y1) %矩形序列 case 2 n1=10。 n3=20。 n=n1:n3。 y1=[(n=n0+n4)amp。 stem(n,y1)。 y1=power(,n)。 y1=sin(*pi*n) 。 guidata(hObject,handles)。string39。 %從 b0的 edit text控件中獲取 b0的值 b1=str2num(get(,39。))。string39。 a1=str2num(get(,39。))。string39。 A=[a0 a1 a2]。 n = [10 : 20]。 %將用戶所選擇的序列 x(n)賦給 x y=filter(B,A,x)。 stem(n,y,39。)。 B=[b0 b1]。 19 x=[n==0]。 %求系統(tǒng)單位取樣響 應(yīng) axes()。fill39。 function pushbutton3_Callback(hObject, eventdata, handles) close(gcf)。string39。 %從 b0的 edit text控件中獲取 b0的值 b1=str2num(get(,39。))。string39。 a1=str2num(get(,39。))。string39。 A=[a0 a1 a2]。 n = [10 : 20]。 %將用戶所選擇的序列 x(n)賦給 x y=filter(B,A,x)。 stem(n,y,39。)。 靜態(tài)界面： 20 各個 callback 函數(shù)為： 取樣序列、矩形序列、實 指數(shù)序列、正弦序列 global x1 f p f=str2num(get(,39。))。string39。 val = get(hObject,39。)。String39。 %獲得句柄的 String(內(nèi)容 )屬性 switch val case 1 n = [10 : 10]。 %對取樣序列 x進行初始化 x1=f.*x。 x=[n=p][n=p+5]。 case 3 n=[10:10]。 case 4 n=[10:10]。 end = double(x1)39。 %保存句柄 離散時間信號的傅里葉變換（ DTFT） global X n = [10 : 10]。 %將用戶所選擇的序列 x(n)賦給 x k=0:500。 X=x39。*k)。 angX=angle(X)。plot(w/pi,magX)。幅 度響應(yīng) 39。 %將序列 x(n)的 DTFT 的幅度繪制在坐標 1（ tu1）中，離散時間信號的 DTFT 是一個復(fù)數(shù)，且是頻率 的連續(xù)函數(shù) axes()。title(39。)。string39。 %從 a0 的 edit text 控件中獲取 a0 的值 a1=str2num(get(,39。))。string39。 b0=str2num(get(,39。))。string39。 b=[b0 b1]。 [H,w]=freqz(b,a,400)。angH=angle(H)。plot(w/pi,magH)。幅度響應(yīng) 39。 %將系統(tǒng)頻率響應(yīng)的幅度繪制在坐標 1（ tu1）中，系統(tǒng)頻率響應(yīng)是一個復(fù)數(shù)，且是頻率 的連續(xù)函數(shù) axes()。title(39。)。圖片制作完成后，只需將在 GUI 對象設(shè)計區(qū)所布置好的幾個控件的 string 屬性修改完，信號取樣的靜態(tài)界面基本就制作完成了。“信號取樣”的靜態(tài)界面 如下圖所示，界面上展示了信號取樣的詳細流程，可以讓用戶對信號取樣的過程一目了然。Value39。 str = get(hObject, 39。)。 switch val case 1 t = [5:: 10]。 %產(chǎn)生單位階躍信號 x=sym(39。)。 case 2 t=[5::10]。 x=sym(39。)。 case 3 t=[5::10]。 x=sym(39。)。 case 4 t=[5::10]。 x=sym(39。)。 end = vpa(x)39。 guidata(hObject, handles)。 yn=。 t=[1,zeros(1,49)]。 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]。 yn2=yn.*t2。 傅里葉變換 DTFT xn=。 %對 x求傅里葉變換， x為 popupmenu菜單所產(chǎn)生的經(jīng)過 sym 23 函數(shù)聲明了的信號， fourier( )函數(shù)只能處理由 sym 函數(shù)定義的符號變量。2*pi*f39。w39。 Ff_conj=conj(Ff)。 %求傅里葉變換 Ff 的幅度。 ezplot(GF)。在靜態(tài)界面設(shè)計好之后，編寫各個按鈕的 callback 函數(shù)，就可以實現(xiàn)界面的功能。string39。 a1=str2num(get(,39。))。string39。 b1=str2num(get(,39。))。string39。 24 A=[a0,a1]。 n=0:30。 %在坐標 1（ axes1）中畫出 的零極點圖 h=impdt(n)。 x=filter(A,B,h)。fill39。grid on syms z。 x=iztrans(Z)。 xbh=simplify(x)。 set(,39。,xbh)。 輸入之后結(jié)果如下圖 1所示 離散傅里葉變換及其快速算法的 GUI 界面的設(shè)計與實現(xiàn) “離散時傅里葉變換及其快速算法”的靜態(tài)界面如上圖所 2示： 離散傅里葉級數(shù) (DFS)的設(shè) 計與實現(xiàn) “離散傅里葉級數(shù) (DFS)”的靜態(tài)界面如下圖所示。 ① 靜態(tài)界面： ② 輸入之后結(jié)果如下圖所示：列如選擇矩形序列 25 各個 callback 函數(shù)如下： 矩形序列、遞增序列、遞減序列、三角序列、實指數(shù)序列、正弦序列 global n xtilde N=str2num(get(,39。))。string39。 val = get(hObject,39。)。String39。 %獲得句柄的 String(內(nèi)容 )屬性 switch val case 1 x=[ones(1,L),zeros(1,NL)]。*ones(1,4)。 n = [N :1: (3.*N1)]。 stem(n,xtilde)。矩形序列 39。 case 2 n=[0:1:N1]。 %產(chǎn)生一個長度為 L 的遞增序列 xtilde=x39。 %將序列 x周期化，周期為 N xtilde=(xtilde(:))39。 stem(xtilde)。遞增序列 39。 case 3 n=[0:1:N1]。 %產(chǎn)生長度為 L 的遞減序列 xtilde=x39。 xtilde=(xtilde(:))39。 stem(xtilde)。遞減序列 39。 case 4 n=[0:1:N1]。 %產(chǎn)生一個長度為 L 的三角形序列 xtilde=x39。 xtilde=(xtilde(:))39。 stem(xtilde)。三角序列 39。 case 5 n=[0:1:N1]。 %產(chǎn)生一個長度為 L 的實指數(shù)序列 xtilde=x39。 xtilde=(xtilde(:))39。 stem(xtilde)。實指數(shù)序列 39。 case 6 n=[0:1:N1]。 %產(chǎn)生一個長度為 L 的正弦型序列 xtilde=x39。 xtilde=(xtilde(:))39。 stem(xtilde)。正弦型序列 39。 end = double(x)39。 離散傅里葉級數(shù) DFS N=str2num(get(,39。))。 k=[ceil(N/2):ceil(N/2)]。 xh=xn*([n=0][n=N])。 %對周期 序列的一個周期求 DFS,這里用函數(shù) fft( )來求 DFS，其實對于有限長的一個周期序列， DFS 與 DFT 沒有本質(zhì)區(qū)別 magXk=abs([Xk(ceil(N/2+1):N) Xk(1:ceil(N/2+1))])。 stem(k,magXk)。幅度響應(yīng) 39。 axes()。 %從 Xk 中分離出其相位 stem(k,angXk)。相位響應(yīng) 39。 離散傅里葉變換 (DFT)的設(shè)計與實現(xiàn) ① “ 離散傅里葉變換 (DFT)”的靜態(tài)界面如下所示 ② 輸入結(jié)果后，如實指數(shù)序列： 27 各個 callback 函數(shù)為 矩形序列、遞增序列、遞減序列、三角序列、實指數(shù)序列、正弦序列 global N L N=str2num(get(,39。))。string39。 axes()。Value39。 axes()。 x1=[n=0]。 x=x1x2。 %畫出周期序列 case 2 n=[0:1:N1]。 %產(chǎn)生一個長度為 L 的遞增序列 xtilde=x39。 %將序列 x周期化，周期為 N xtilde=(xtilde(:))39。 case 3 n=[0:1:N1]。 %產(chǎn)生長度為 L 的遞減序列 xtilde=x39。 xtilde=(xtilde(:))39。 case 4 n=[0:1:N1]。 %產(chǎn)生一個長度為 L 的三角形序列 xtilde=x39。 xtilde=(xtilde(:))39。 case 5 n=[0:1:N1]。(nN)]。*ones(1,5)。 stem(xtilde)。 x=cos(2*pi/N*n)。*ones(1,5)。 stem(xtilde)。 guidata(hObject, handles)。string39。 %從的 edit text 控件中獲取 a0 的值 axes()。 k=[ceil(N/2):ceil(N/2)]。 xh=xn*(udt(n)udt(nN))。