【正文】 。filled39。),title(39。實部39。)。 subplot(2,2,3), stem(k,Xi,39。filled39。),title(39。虛部39。)。 subplot(2,2,2), stem(k,Xa,39。filled39。),title(39。模39。)。 subplot(2,2,4),stem(k,Xn,39。filled39。),title（’相角’）。輸入dfzsu(0,20,pi/4) ，結果如圖： 離散時間序列卷積及MATLAB實現(xiàn)例6：已知某LTI 離散系統(tǒng)，其單位響應h(k) = e(k) e(k 4)，求該系統(tǒng)在激勵為f (k) = e(k) e(k 3)時的零狀態(tài)響應y(k)，并繪出其時域波形圖。解：由LTI 變離散系統(tǒng)的分析可得y(k) = h(k) * f (k)因此，我們可調用dconv()函數來解決此問題，相應的MATLAB 命令為：f1=ones(1,4)。k1=0:3。f2=ones(1,3)。k2=0:2。[f,k]=dconv(f1,f2,k1,k2)運行結果為：f =1 2 3 3 2 1k =0 1 2 3 4 5 時域波形圖例7：已知兩連續(xù)時間信號如下圖所示，試用MATLAB 求f（t）= f1（t） * f2（ t），并繪出f (t)的時域波形圖。解: 我們可以調用前述的函數sconv( )來解決此問題，即首先設定取樣時間間隔p，并對連續(xù)信號f1（t）和f2（t）的非零值區(qū)間以時間間隔p 進行抽樣，產生離散序列f1和f2，然后構造離散序列f1 和f2 所對應的時間向量k1 和k2，最后再調用sconv 函數即可求出f1（t） * f2（ t）的數值近似，并繪出其時域波形圖。實現(xiàn)上述過程的 MATLAB 命令如下：p=。k1=0:p:2。f1=*k1。k2=k1。f2=f1。[f,k]=sconv(f1,f2,k1,k2,p)上述命令繪制的波形如圖所示。圖中分別給出了取樣時間間隔p= 和p=。可見，當取樣時間p 足夠小時，函數sconv( )的計算結果就是連續(xù)時間卷積f（t）=f1（t）*f2（ t）的較好的數值近似 例8：已知描述某連續(xù)系統(tǒng)的微分方程為：2 y162。162。(t ) + y162。(t) + 8 y(t) = f (t )試用 MATLAB 繪出該系統(tǒng)的沖激響應和階躍響應的波形。解：直接調用函數impulse()和step()即可解決此問題，對應的MATLAB 命令如下：b=[1]。a=[2 1 8]。subplot(1,2,1)impulse(b,a)subplot(1,2,2)step(b,a）上述命令繪制的系統(tǒng)沖激響應和階躍響應的波形如圖所示。例9：已知描述某離散系統(tǒng)的差分方程如下：2 y(k) 2y(k 1) + y(k 2) = f (k ) + 3 f (k 1) + 2 f (k 2)試用 MATLAB 繪出該系統(tǒng)0~50 時間范圍內單位響應的波形。解：首先用向量a 和b 表示出該離散系統(tǒng)，然后再調用impz()函數即可解決此問題。實現(xiàn)上述過程的MATLAB 命令如下：a=[2 2 1]b=[1 3 2]。impz(b,a)繪制的系統(tǒng)單位響應波形圖如圖所示例10：已知描述離散系統(tǒng)的差分方程為：y(k ) y(k 1) + y(k 2) = f (k ) + f (k 1)且知該系統(tǒng)輸入序列為f (k) = （）^k*ε(k)試用 MATLAB 實現(xiàn)下列分析過程：（1） 畫出輸入序列的時域波形；（2） 求出系統(tǒng)零狀態(tài)響應在0~20 區(qū)間的樣值；（3） 畫出系統(tǒng)的零狀態(tài)響應波形圖。解：我們可以調用filter()函數來解決此問題。實現(xiàn)這一過程的MATLAB 命令如下：a=[1 ]。b=[1 1]。t=0:20。x=(1/2).^t。y=filter(b,a,x)。subplot(2,1,1)。stem(t,x)。title(‘輸入序列’)。subplot(2,1,2)。stem(t,y)。title(‘響應序列’) 小露一手：x=cos(5/16*pi*n)。X_16=fft(x,16)。X_32=fft(x,32)。subplot(211)。hold on。plot(n/16,abs(X_16),39。r39。)。plot(L/32,abs(X_32))。hold off。w=linspace(0,pi,512)。h=freqz(x,1,w)。subplot(212)。plot(w/pi,abs(h))。課程設計體會在短短兩周的課程設計里，我學到了MATLAB的基礎應用知識，和簡單的信號與系統(tǒng)分析及MATLAB實現(xiàn)。在以前的課程中，對信號總覺得非常抽象。尤其是卷積和在傅里葉變換、拉普拉斯變換，S域模型等，都只是通過一系列的數學推導和運算來說明其性質和規(guī)律，但自己并不能很好的理解，由于沒有其他物理規(guī)律只是那么直觀，所以只能靠記憶來學習。往往還會混淆，這非常不利于知識的延伸和升華。通過這次課程設計，在對MATLAB的應用基礎熟悉之后，按照書本上的例題，自己動手實際操作，通過計算機的模擬實踐，我們能清楚直觀地看到以前覺得抽象的數學符號運算變?yōu)榫唧w的圖形，使知識變得簡單易懂。特別有助于我們學習和記憶。但是，由于時間較短，我們不可能學會相關的編程，只是在一步一步驗證以前學過的理論。雖然這些看起來很簡單，但是對于我們來說，MATLAB是一個全新的東西，正如大一時的VisualBasic，一時間我們很難去把握它。所以在這些驗證過程中，往往還是出錯，很多時候就是因為一個標點符號。這在一定程度上告訴我們，做學問一定要嚴謹，不能有絲毫的差錯。在看到有關書籍和網絡上介紹有關MATLAB的更加強大的功能之后，更覺得我們學到的真是微乎其微。在驚嘆這個軟件強大之余，也萌生了要找時間好好學習使用這個非常實用的軟件。指導教師評語及設計成績 評 語 課程設計成績： 指導教師： 日期： 年 月 日 參考文獻 [1] 吳大正, 楊林耀, 張永瑞. 信號與線性系統(tǒng)分析（第三版）. 高等教育出版社. [2] 鄭君里，應啟珩，楊為理. 信號與系統(tǒng). 高等教育出版社. 北京. 2000[3] 吳新余，周井泉，沈元隆. 信號與系統(tǒng)——時域、頻域分析及MATLAB 軟件的應用. 電子工業(yè)出版社. 北京. 1999[4] 閻鴻森，王新鳳，田惠生. 信號與線性系統(tǒng). 西安交通大學出版社. 西安. 1999[5] 程衛(wèi)國，馮峰，姚東等. 應用指南. 人民郵電出版社. 北京. 1999[6] 樓順天，李博菡. 基于MATLAB 的系統(tǒng)分析與設計——信號處理. 西安電子科技大學出版社. 西安. 1999[7] 陳懷琛，王朝英，高西全等譯. 數字信號處理及其MATLAB 實現(xiàn). 電子工業(yè)出版社. 北京. 1998[8] 劉樹棠譯. 信號與系統(tǒng)計算機練習——利用MATLAB. 西安交通大學出版社. [9] 高俊斌. MATLAB 語言與程序設計. 華中理工大學出版社. 武漢. 1998[10] 張志涌等. 版本. 北京航空航天大學出版. 北京. 200028