【正文】 %只處理周期序列的一個周期的值 Xk=fft(xh,N)。 xn=。 n = 0:1:N1。))。 離散傅里葉變換 DFT N=str2num(get(,39。 end = double(x)39。 xtilde=(xtilde(:))39。 %產(chǎn)生一個長度為 L 的正弦型序列 28 xtilde=x39。 case 6 n=[0:1:N1]。 xtilde=(xtilde(:))39。 %產(chǎn)生一個長度為 L 的實指數(shù)序列 xtilde=x39。 x=.^n.*[(n=0)amp。 stem(xtilde)。*ones(1,5)。 x=(n+1).*(udt(n)udt(nceil(L/2)))+(Ln).*(udt(nceil(L/2))udt(nL))。 stem(xtilde)。*ones(1,5)。 x=(Ln).*(udt(n)udt(nL))+0.*(udt(nL)udt(nN))。 stem(xtilde)。*ones(1,5)。 x=(n+1).*(udt(n)udt(nL))。 stem(n,x)。 x2=[n=L]。 switch val case 1 n=[0:1:N1]。)。 val=get(hObject,39。))。 L=str2num(get(,39。string39。)。 title(39。 angXk=angle([Xk(ceil(N/2+1):N) Xk(1:ceil(N/2+1))])。)。 title(39。 %從 Xk 中分離出其幅度 axes()。 %只處理周期序列的一個周期的值 Xk=fft(xh,N)。 xn=。 %從的 edit text 控件中獲取 a0 的值 n = 0:1:N1。string39。 guidata(hObject, handles)。)。 title(39。 axes()。*ones(1,5)。 x=cos(2*pi/N*n)。)。 title(39。 axes()。*ones(1,5)。 26 x=.^n.*([n=0][n=L])。)。 title(39。 axes()。*ones(1,5)。 x=(n+1).*([n=0][(nceil(L/2))=0])+(Ln).*(([(nceil(L/2))=0])[n=L])。)。 title(39。 axes()。*ones(1,5)。 x=(Ln).*([n=0][n=L])+0.*([n=L][n=N])。)。 title(39。 axes()。*ones(1,5)。 x=(n+1).*([n=0][(nL)=0])。)。 %畫出周期序列 title(39。 axes()。 %將序列 x周 期化，周期為 N xtilde=(xtilde(:))39。 %產(chǎn)生一個長度為 L 的矩形序列 xtilde=x39。)。 %獲得句柄的 value(值 ) 屬性 str = get(hObject, 39。Value39。))。 L=str2num(get(,39。string39。在靜態(tài)界面搭建好之后，編寫各個控件的 callback 函數(shù)，就可以實現(xiàn)界面的功能。 %將經(jīng)過簡化等處理的 x(n)的值顯示在按鈕旁邊的 static text 文本框中， guidata(hObject, handles)。String39。 xbh=char(xbh)。 %通過 iztrans函數(shù)求 的反變換 x(n)， iztrans函數(shù)只能處理由 sym函數(shù)定義的符號變量，它所求出來的 x(n)也是符號變量。 Z=z*(a0*z+a1)/(b0*z^2+b1*z+b2)。)。 %通過 filter 函數(shù)求 的反變換 x(n) stem(n,x,39。 axes()。 zplane(A,B,)。 B=[b0,b1,b2]。))。 b2=str2num(get(,39。string39。))。 b0=str2num(get(,39。string39。))?！?z的正反變換及系統(tǒng)函數(shù)”的靜態(tài)界面如下圖所示， 各個 callback 函數(shù)為： 系統(tǒng)函數(shù) H（ z） global a0 a1 b0 b1 b2 a0=str2num(get(,39。 輸入值后結(jié)果如下圖所示列如單位階躍信號： ① 取樣 ② 傅里葉變換 FT z 的正反變換及系統(tǒng)函數(shù)的設計與實現(xiàn) “ z 的正反變換及系統(tǒng)函數(shù) ” 的靜態(tài)界面的制作方法與前面幾個界面的制作是相似的，同樣是根據(jù)稿紙上設計的草圖，將界面上所用到的 控件從 GUI 設計窗口左側(cè)的工具欄中直接拖拽到對象設計區(qū)，并調(diào)整大小、位置及修改各個控件的string 屬性。 axes()。 %Ff 的共軛 Ff_conj GF=sqrt(Ff*Ff_conj)。)。,39。 Ff=subs(Fw,39。 Fw=simplify(fourier(xn))。 stem(n,yn2)。 t2=(t1)39。 n = [5:: 10]。 len=length(yn)。 取樣 axes()。 %將 x轉(zhuǎn)換為 vpa 格式的字符 = double(y)39。 plot(t,y)。Heaviside(t)39。 y=3*sin(2*pi*t)。 plot(t,y)。Heaviside(t)39。 y=3*2.^t。 plot(t,y)。Heaviside(t)39。 y=[t=0][t=2]。 %將單位階躍信號用 sym 聲明，便于求其傅里葉變換 plot(t,y)。Heaviside(t)39。 y=[t=0]。 axes()。String39。)。 22 各個 callback 函數(shù)為： 取樣信號、矩陣信號、實指數(shù)信號、正弦信號 val = get(hObject,39。在靜態(tài)界面設計好之后，編寫各個按鈕的 callback 函數(shù)，就可以實現(xiàn)界面的功能。 輸入之后如下圖所示： ① 離散時間信號的傅里葉變換（ DTFT） ② 系統(tǒng)頻率響應 信號取樣的設計與實現(xiàn) “信號取樣 ” 的靜態(tài)界面的制作方法與前面的界面有一點不同的是：它內(nèi)置了一張圖片來演示整個信號取樣流程，界面所需要的控件，如按鈕、坐標框、下拉菜單等都是按照圖片指定的位置布局的。相位響應 39。plot(w/pi,angH)。)。title(39。 axes()。 %利用函數(shù) freqz 求系統(tǒng)的頻率響應 magH=abs(H)。 %將序列 x的系數(shù)向量賦給 b a=[a0,a1,a2]。))。 b1=str2num(get(,39。string39。))。 a2=str2num(get(,39。string39。))。 %將序列 x(n)的 DTFT 的相位繪制在坐標 2（ tu2）中 21 系統(tǒng)頻率響應 global a0 a1 b0 a0=str2num(get(,39。相位響應 39。 plot(w/pi,angX)。)。title(39。 axes()。 %求序列 x(n)的離散傅里葉變換（ DTFT） magX=abs(X)。*(exp(j*pi/500)).^(n39。w=(pi/500)*k。 x = 。 %將序列 y1 轉(zhuǎn)換成 double 型的數(shù)據(jù) guidata(hObject, handles)。 x1=f.*sin(2*pi*(n+p))。 x1=f.*2.^(n+p)。 x1=f.*x。 case 2 n=[10:10]。 x=[n==p]。)。 %獲得句柄的 value(值 ) 屬性 str = get(hObject, 39。Value39。))。 p=str2num(get(,39。string39。 %將系統(tǒng)響應 y(n)序列繪制 在坐標 2（ axes2）中，線端是實心圓點 ① 系統(tǒng)響應 ② 系統(tǒng)的單位取樣響應 離散時間信號傅里葉變換及系統(tǒng)頻率響應的設計與實現(xiàn) “離散時間信號傅里葉變換及系統(tǒng)頻率響應 ” 的靜態(tài)界面如下圖所示，在靜態(tài)界面設計好之后，編寫各個按鈕的 callback 函數(shù)，就可以實現(xiàn)界面的功能。fill39。 %通過 filter 函數(shù)求系統(tǒng)的響應 axes()。 x= 。 %將序列 x的系數(shù)向量賦給 A B=[b0 b1]。))。 a2=str2num(get(,39。string39。))。 a0=str2num(get(,39。string39。))。 lisan(handles) % hObject handle to pushbutton3 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global b0 b1 a0 a1 a2 b0=str2num(get(,39。)。 stem(n,h,39。 %令序列 x(n)為單位取樣序列 h=filter(B,A,x)。 n = [10 : 20]。 %將系統(tǒng)響應 y(n)序列繪制在坐標 2（ axes2）中，線端是實心圓點 系統(tǒng)的單位取樣響應 h（ n） global b0 b1 a0 a1 a2 A=[a0 a1 a2]。fill39。 %通過 filter 函數(shù)求系統(tǒng)的響應 axes()。 x= 。 %將序列 x的系數(shù)向量賦給 A B=[b0 b1]。))。 a2=str2num(get(,39。string39。))。 a0=str2num(get(,39。string39。))。 系統(tǒng)響應 y（ n） global b0 b1 a0 a1 a2 b0=str2num(get(,39。%正弦序列 stem(n,y1) end =double(y1)。 %實指數(shù)序列 stem(n,y1) case 4 n=[10:1:20]。 %取樣序列 case 3 n=10:20。(nn2+n4)]。 n4=0。 n0=0。 n2=5。 y1=[zeros(1,10),1,zeros(1,20)]。 axes()。String39。)。 靜態(tài)界面： 各個 callback 函數(shù)如下： 取樣序列 、矩形序列、實指數(shù)、正弦型序列 val=get(hObject,39。 lisan(handles) 輸入值后如下個圖所示： ① 矩形序列 *矩形序列： ② 實指數(shù)序列 *矩形序列： ③ 實指數(shù)序列 *實指數(shù)： 離散時間系統(tǒng)的設計與實現(xiàn) 離散時間系統(tǒng)的靜態(tài)界面如下圖所示，與上一個界面相似，將界面上所用到的控件從 GUI設計窗口左側(cè)的工具欄中直接拖拽到對象設計區(qū)，并調(diào)整大小、位置及修改各個控件的string 屬性，就制作好了離散時間系統(tǒng)的靜態(tài)界面。 stem(n,y2)。 y2=.^n。 stem(n, y2)。 y2=[x1x2]。 x1=[n=0]。 %獲得句柄的 String(內(nèi)容 )屬性 axes()。String39。)。 %從菜單左邊的 edit text 控件中獲取序列的長度 val = get(hObject,39。string39。 stem(n,y)。 n=0:ny1。 y=a*y1。 h=。 a=A。string39。 靜態(tài)界面： 因計算機只能處理有限長序列，因此在本系統(tǒng)中所給出的序列也是有限長的，并且序列的長度可由用戶自定義輸入。 各序列輸入值運行后如下圖所示： ① 取樣序列： ② 階躍序列： ③ 矩形序列： ④ 實指數(shù)序列： ⑤ 正弦型序列： ⑥ 復指數(shù)序列： 16 線性卷積的設計與實現(xiàn) “線性卷積”界面所用到的控