內容簡介
本教程中的全部例題程式和演示程式都可以免費從網上下載。本教程把《數位訊號處理》課程定位為講述“用數字方法處理(任何)信號”的技術,因此,對模擬信號的計算機處理給予了充分的篇幅,而這正是目前許多數位訊號處理教材的弱點。
本書適合作為電子信息專業的大學本科學生教材,書中採用了大小兩種字號混排,目的是適應不同專業和不同要求的本科學生。對於學過數位訊號處理課程,但沒有用MATLAB來解決有關問題的研究生和工程技術人員,也有很好的參考價值。
出版背景
本書第1版於2004年11月出版以後,得到了很好的反應。因為這本書與其他傳統的教材不同,它是一本用MATLAB全面闡述數位訊號處理教材,給出了所有概念闡述的形象演示和分析設計的計算機算法,而且提供了例題和演示的程式集,書中所有的理論和公式都可以通過執行程式得到驗證,理論聯繫實際,學完了會快速地用計算機解決問題。心中有底。讀者反映這本書的缺點主要是對本科生而言,內容有些偏多。
在修訂出版第2版的時候,開始我們想把書改淺一些,薄一點,為此又專門徵求了一些曾仔細閱讀使用過本書的教授的意見。現把他們的意見摘出幾條:“我覺得您已出版的‘釋義’雖然厚一些,但寫得易懂。若從頭寫一本更薄的,恐怕篇幅小,難以達到“釋義”那本書的效果。”,“您的書是國內用MATLAB闡述DSP的第一本好書。最好不重新另寫一本。完全可以在原來基礎上進一步完善。”……經過探討,大家覺得本書在對物理和數學概念的闡述上毫不遜於傳統的教材,例如“關於模擬頻率和數字頻率的定義”,“關於角頻率的正負”,“對DFT的序號k與實際頻率f的關係”等等,講得很細緻,很便於本科生接受。其實,作者收到過一些本科生讚揚的來信,他們的經驗是,不要怕書厚,不要怕MATLAB,一個一個問題看下去,這本書講得很清楚,而且馬上可以用程式驗證,很好懂。
目前,大學本科用這本書作教材的還不太多,第一個原因是很多大學生(甚至部分老師)MATLAB基礎差;見到程式就頭痛了,不敢用MATLAB教學;但是從國內外發展趨勢來看,要搞信號處理,不用計算機是不行的,用MATLAB的學校只會愈來愈多。令我們高興的是:隨著時間推移,本書的大量程式已為多本其他教材所引用。我們不能自己放棄或減弱這個特色。第二個原因是書中確有少數偏深的內容,這些內容雖然超越了本科大綱,但有利於已經畢業的讀者在設計中更好地運用MATLAB,也有利於大學生在後續課程或畢業設計中套用。考慮到第1版的讀者群並不限於本科生,很多是畢業了的工程師或研究生,這些內容仍是重要的,不能刪除。為了兼顧大學生閱讀方便,在本版中把它們用小號字印刷。
這樣,第2版的修訂量就比較小,主要包括了以下三方面:(1)改正了第1版中的一些印刷錯誤和圖文標註;(2)由於MATLAB軟體的升級,對相關的部分按新版本作了修訂;(3)儘量減少大號字的篇幅,把它們改為小號字,使本科生能抓住要點,閱讀方便。修訂後程式集內容基本不變,新程式集名稱改為dsk04n,修改中儘量考慮了新舊版本的兼容問題。
北京理工大學的王大倫教授曾仔細地閱讀過本書第1版,對書中的大小錯誤都作了批註,並反饋給作者,對第2版的修訂幫助很大,在此特表謝意。
MATLAB已經可以把信號處理的概念設計變成了非常方便的工作。近若干年來,Mathworks公司致力於擴展MATLAB的功能,目的是用基於模型的流程,使用統一的軟體工具,把概念設計變成工程設計。這就更提高了把MATLAB作為設計工具的重要性。我們在本版中加了第10章,對此作了一個引導性的介紹。北京邁斯沃克公司(Mathworks China)的徐正高高級工程師寫了該章中的10.4和10.5兩節。
作者簡介
陳懷琛教授,1934年生於上海,1953年畢業於軍事電信工程學院,留校任教至今。期間曾到美國賓夕法尼亞大學系統工程系做訪問學者,他先後在西安電子科技大學機械系,自動控制系,和電子工程系講授過十多門課程,並曾任副校長。主管科研和研究生教育。他還曾兼任中國自動化學會理事,中國電子學會教育學會副理事長,陝西人民對外友協副會長,陝西自動化學會副理事長等職務。
陳懷琛教授1994年卸任後,繼續任教並致力於推動大學課程和教學的計算機化。目標是使教師及學生都用計算機取代計算器來解決各課程的問題。出版 著作多本,並多次應邀在各研究所和大學做系列講座。
第1版序
隨著電子技術迅速地向數位化發展,《數位訊號處理》越來越成為廣大理工科,特別是IT領域的學生和技術人員的必修內容,先前也出版過很多有關這方面的書籍。但國內這門課程本身的現代化進程,卻並不令人滿意。這裡的“現代化”指的是利用最新的計算技術和網路手段進行教學。我編寫這本書的目的,就是想突出利用MATLAB和Internet來進行數位訊號處理的教學。
信號處理涉及大量的運算,可以說離開了計算機及相應的軟體,就不可能解決任何稍微複雜的實際套用問題。早在20世紀90年代中期,MATLAB就已成為國際公認的信號處理的標準軟體和開發平台。從1996年後,美國新出版的信號處理教材就沒有一本是不用MATLAB的(見參考文獻[1,3,7,8,9])。而國內把MATLAB貫徹始終的信號處理教材卻還很少,在教學中真正用MATLAB的就更少,兩者相比落後了近八年。本書的一個意圖就是使之成為一本自始至終使用MATLAB闡述問題和進行計算的數位訊號處理教材,更好地與國際接軌。
關於《數位訊號處理》課程的內涵,存在著兩種不同看法:一種看法把“數位訊號”限定為“處理”的對象;另一種看法則把“數位技術”作為一種方法,研究如何用它來“處理信號”。從第一種看法出發,“數位訊號處理”課程只講離散信號。經驗證明,這種劃分會造成學生不會用數位技術處理連續信號,形成學科內容的“空白區”。而從第二種看法出發,“數位訊號處理”課程應該同時考慮連續和離散信號,重點是如何用計算機進行分析和計算。我贊同後一種看法。因此,在本書多處加強了離散信號和連續信號等價性的討論,並詳細介紹了連續信號頻譜的數值計算問題。
中國教材一般將快速傅立葉變換列為一章,篇幅約30~40頁。美國教材則通常列為一節,篇幅一般只有10頁左右。這是什麼原因?我認為在目前所有的科學計算軟體和數位訊號處理晶片中,FFT已經是標準的模組,一般都沒有必要自己編寫FFT程式,所以重點應該是FFT的用法而不是算法。國外教材對這部分內容經歷了“薄—厚—薄”的歷史過程,我國則慢了一拍。所以本教材也精簡了這方面的內容。
在把數位訊號處理概念與MATLAB語言相結合的方面,國外教材有多種模式。有的以講概念為主,MATLAB放在作業和實驗中(見參考文獻[9]);有的以MATLAB為主,概念則建議讀者另選教材(見參考文獻[1]);有的兩者都講得不少,但書就變得很厚(見參考文獻[8]):本書則大體採納了參考文獻[3]的模式,以系統地講解信號處理的原理為主,用MATLAB程式的核心語句幫助解釋概念。要做到這點,學生必須在修本課程以前已經學過並用過MATLAB。當然,這並不是說,沒有學過MATLAB就不能用這本書。我曾在西安電子科技大學和桂林電子工業學院兩次使用MATLAB進行數位訊號雙語教學,教材是英文書(見參考文獻[1]),課堂90%用英語講,同學中多數沒有學過MATLAB。課還是正常進行,完成得還可以。但這些學生都已通過英語四級,因此屬於班級中程度較高者,而且多加一個學分。儘管如此,在教學中有些同學仍提出“究竟著重概念還是編程”這樣的問題,說明他們對於將語言和概念結合起來仍感到困難。所以,經驗告訴我們:從大面積課堂教育而言,最好不要把MATLAB放到高年級和信號處理混在一起學。
我一貫的指導思想是:學生應在大學一年級或二年級學習MATLAB,這對各門課程都有好處。信號處理課是大學三年級開設的。此時學生應該已經能用矩陣建模並用MATLAB語言來思考和概括數學問題了。有了這樣的基礎,才能對本書有較深的體會。建議未學過MATLAB的讀者,至少用10小時左右先學習其入門並實際做一些套用。為此,我們新製作了一套“MATLAB入門”(兩張)講課光碟。該組光碟以MATLAB6.x為背景(也介紹一點MATLAB 7.0),可在計算機上實時播送。本課程的老師統一購買了這套光碟給大家播放,可以不用費勁去備MATLAB課。與之配套的教材是參考文獻[6]或參考文獻[16]。其中,參考文獻[16]是為大學一、二年級同學設計的,在套用方面結合的是數學、物理、力學、電學等;參考文獻[6]是為大學二、三年級同學設計的,結合了電路、信號和系統、信號處理和控制原理四門課程。這套光碟由本人主講,西安電子科技大學網路學院製作。有關發行購買事項,請各位教師注意下載程式中另附的通知。
本書利用電子工業出版社的博文視點公司網站向讀者免費提供完整的子程式集、例題程式集和示教程式集(即附錄B,C,D中的全部程式)共180餘個,程式集名為dsk04,可免費下載。
單擊資源下載,再按書名查找即可。下載前要按提問填寫下載人的必要信息,在此網頁上,還將提供作者其他有關書籍程式集和光碟獲取的信息。由於讀者可以免費獲得全部的程式,從而也可得到程式產生的圖形。我們把書中的程式都用小字排版,圖形也縮排了,以降低篇幅和書價。希望讀者理解。
為了使教師能夠較方便地利用本書進行教學,我們準備了一套教師用光碟。其中包括兩部分內容:一是為講授本書用的PowerPoint幻燈片,約900張;二是本書部分習題的題解,因為有些教師認為若有題解流入學生中,出題的效果將打折扣。為了兼顧對題解的不同看法,我們只給了三分之一的題解。預計這張教師光碟將於本書出版後三個月內完成。由“西安電子科技大學電工電子教學基地”進行著作權管理。我們將在下載程式中附上有關教師用光碟的發行辦法和申購表,請讀者注意。
為了適用不同類型的專業和學生對數位訊號處理的不同要求,本書在排版上分了大小兩種字型。大字是達到基本要求應掌握的內容,小字排版的部分是可以跳過不看的,不會影響對本課程基本概念的掌握。整個第9章都屬於任選內容,但沒有用小字,因為這個會引起誤解,同時使讀者看書時可以舒服一些。每章的最後都給出了複習思考題和習題,以利於讀者自學複習。
為了闡述數位訊號處理中的重要概念,本書採用了一些其他教材中未採用的術語和方法。例如,離散時間傅立葉變換、採樣序列的等價序列、符幅特性、在無窮遠鄰域的收斂性、離散最小二乘設計法等。書中也有一些作者本人的研究工作內容。比如,信號流圖的解析形式及MATLAB解、信號流圖與狀態空間法的轉換和MATLAB解、係數向量量化的子程式以及一些有圖形、動畫和聲音表現的演示程式等。其目的是引起大家討論,更好地解放思想,用創新的態度來對待這門課程的教學。
這本書採用了參考文獻[4,3,1,5,8]中不少的內容和例題,用到了這些作者的寶貴經驗,也借用了一些圖形、例題和習題,乘此機會向這些教材的作者們表示謝意。在2003年8月的“第十一屆信號處理學術年會”上,我曾把編寫本書的指導思想,與許多信號處理界的專家(如袁保宗、程佩青、張賢達等教授)交換過意見,得到了他們的支持。本書初稿形成後,又曾廣泛徵求了本校的丁玉美、趙樹傑、廖桂生、高西全、史林等長期從事數位訊號處理教學的教授們對初稿的意見,特別是丁玉美教授仔細閱讀了稿件,給作者很大的幫助。西安電子科技大學的段寶岩校長還特別資助了作者的學術交流經費,保錚院士對本書的出版也給予了關心。此外Mathworks公司對本書的出版也給予了軟體新版本的支持,電子工業出版社的博文視點公司的各位編輯對本書一路綠燈,提供了出版和網頁方面的一切方便,使本書能及早問世。對大家的這些幫助,作者表示衷心的感謝。
由於作者水平和編著時間的限制,書中肯定有很多不足之處,歡迎廣大讀者指正,我很希望得到讀者對本書的反饋信息。歡迎讀者將教學實踐中積累的經驗告訴本人,並希望大家不斷改進和充實網上的演示程式!
陳懷琛2004年8月31日
於西安電子科技大學
圖書目錄
第1章 信號和信號處理 1
1.1 信號的分類 1
1.2 信號處理 5
1.2.1 信號處理的舉例 5
1.2.2 信號處理的一般任務 6
1.2.3 數位訊號處理 7
1.3 信號採樣和復原實例 8
1.4 本書概貌的說明 10
複習思考題 12
習題 12
第2章 時域中的離散信號和系統 13
2.1 模擬信號與其採樣序列的關係 13
2.2 基本序列及序列的基本運算 14
2.2.1 序列的表示方法 14
2.2.2 常用的典型序列 15
2.2.3 序列的運算和變換 17
2.3 序列的奇偶分解、卷積和相關 21
2.3.1 周期序列 21
2.3.2 實對稱序列 21
2.3.3 複數域的偶序列和奇序列 22
2.3.4 序列的卷積 23
2.3.5 序列的相關 25
2.4 時域離散系統 26
2.4.1 線性系統 27
2.4.2 線性時不變離散系統 27
2.4.3 LTI離散系統的數學描述——線性常係數差分方程 28
2.4.4 用MATLAB的遞推解法 29
2.4.5 用時移運算元z1表示的離散系統 29
2.5 脈衝回響——用時間序列表征的系統特性 30
2.5.1 脈衝回響的特性 30
2.5.2 用解析卷積方法求系統輸出 31
2.5.3 用向量矩陣乘法進行卷積計算 32
2.5.4 用脈衝回響判斷系統的因果性 34
2.5.5 用脈衝回響判斷系統的穩定性 34
2.5.6 用脈衝回響分類數字濾波器 35
2.6 序列運算的MATLAB函式列表 36
複習思考題 37
習題 38
第3章 離散序列的兩種傅立葉變換 40
3.1 序列信號的周期性和頻率 40
3.1.1 序列信號的周期性 40
3.1.2 序列信號的頻率 41
3.1.3 連續信號的頻譜 43
3.1.4 有限頻寬序列 44
3.2 離散時間傅立葉變換(DTFT) 45
3.2.1 離散時間傅立葉變換(DTFT)的定義 45
3.2.2 DTFT的性質 46
3.2.3 實數奇偶序列的DTFT 48
3.2.4 複數奇偶序列的DTFT 48
3.2.5 DTFT的解析計算示例 49
3.2.6 DTFT的數值計算示例 51
3.2.7 用MATLAB計算DTFT 52
3.2.8 MATLAB計算舉例 54
3.3 離散系統的頻率特性 57
3.3.1 定義 57
3.3.2 用MATLAB計算系統頻率特性 59
3.4 周期序列的頻譜——離散傅立葉級數 60
3.4.1 從一個計算實例出發 60
3.4.2 周期序列的傅立葉級數 62
3.4.3 用MATLAB計算離散傅立葉級數 63
3.5 離散傅立葉變換(DFT) 64
3.5.1 離散傅立葉變換(DFT)的定義 64
3.5.2 用MATLAB計算DFT和IDFT 66
3.5.3 兩種離散傅立葉變換(DFT與DTFT)的關係 66
3.6 離散傅立葉變換(DFT)的特性 68
3.6.1 隱含周期性 68
3.6.2 線性性質 69
3.6.3 循環移位性質 69
3.6.4 循環摺疊性質 71
3.6.5 DFT的循環對稱性 72
3.6.6 循環卷積特性 75
3.6.7 帕瑟瓦爾定理 77
3.7 用DFT計算線性卷積 78
3.7.1 循環卷積與線性卷積 78
3.7.2 分段卷積 81
3.7.3 分段卷積的MATLAB函式 82
複習思考題 84
習題 84
第4章 信號頻譜的高效計算 88
4.1 各種傅立葉變換及其相互關係 88
4.1.1 傅立葉變換的主要類型 88
4.1.2 DFT與其他變換的直接關係 90
4.2 快速傅立葉變換(FFT) 91
4.2.1 時域抽取(DIT)基二FFT算法 91
4.2.2 用MATLAB程式描述FFT算法 95
4.2.3 基四和其他基的FFT算法 96
4.2.4 其他FFT算法 98
4.3 用FFT計算離散時間序列的頻譜 99
4.3.1 有限長離散時間序列的頻譜計算 99
4.3.2 用補零方法由FFT求DTFT 100
4.3.3 無限長序列的頻譜計算 101
4.4 連續信號採樣中的頻譜變換 103
4.4.1 時域採樣定理(奈奎斯特定理) 103
4.4.2 由離散序列重構模擬信號 106
4.4.3 理想預濾波作用的定量分析 109
4.5 連續信號的頻譜計算 110
4.5.1 非周期連續信號的頻譜計算 110
4.5.2 連續周期信號的頻譜計算 116
4.5.3 循環計算中對應頻點的確定方法 117
4.6 從頻譜計算序列信號(IDFT) 119
4.6.1 頻率域採樣定理 119
4.6.2 從頻譜計算離散時間序列 120
4.6.3 從頻譜計算連續時間信號 123
4.7 用FFT計算能量 125
4.8 小結 126
複習思考題 127
習題 128
第5章 變換域中的離散時間系統 130
5.1 z變換 130
5.1.1 z變換的定義與收斂性 130
5.1.2 z變換的重要特性和定理 133
5.1.3 z變換的計算實例 135
5.2 z反變換和差分方程的解 136
5.2.1 用留數定理求z反變換 136
5.2.2 用部分分式法求z反變換 137
5.2.3 冪級數法(長除法)求z反變換 140
5.2.4 用z變換解差分方程 142
5.3 z域中對系統的描述 146
5.3.1 系統函式(傳遞函式) 146
5.3.2 正冪系統函式和零極增益形式 149
5.3.3 零極點分布與系統的因果性和穩定性 151
5.4 z平面上的譜分析 153
5.4.1 零極點與頻率特性的關係 153
5.4.2 非單位圓周上的頻譜分析 154
5.4.3 z平面螺旋線上的譜分析(Chirp-z變換) 155
5.5 理想濾波器 158
5.5.1 理想的因果低通濾波器的頻率特性 158
5.5.2 線性相位條件對脈衝回響對稱性的要求 160
5.5.3 考慮因果性和對稱性對理想濾波器的影響 161
5.6 低階數字濾波器 162
5.6.1 一階數字濾波器 162
5.6.2 二階數字濾波器 164
5.6.3 梳狀濾波器的零極點分布 165
5.7 全通濾波器 166
複習思考題 167
習題 168
第6章 數字濾波器結構 171
6.1 概述 171
6.2 濾波器的信號流圖結構及解法 172
6.2.1 信號流圖 172
6.2.2 信號流圖方程的計算機解法 173
6.2.3 信號流圖與數字系統的網路結構圖 175
6.3 FIR基本網路結構 175
6.3.1 概述 175
6.3.2 直接形式 176
6.3.3 級聯形式 176
6.3.4 線性相位形式 177
6.3.5 頻率樣本結構 178
6.4 IIR基本網路結構 181
6.4.1 直接形式 182
6.4.2 級聯形式 182
6.4.3 並聯形式 184
6.5 一些其他類型的濾波器結構 189
6.5.1 全零點格型(Lattice)濾波器 189
6.5.2 全極點格型濾波器 191
6.5.3 用信號流圖與MATLAB結合求解格型濾波器 193
6.6 狀態變數分析法 195
6.6.1 狀態方程和輸出方程 195
6.6.2 狀態方程和其他結構參數的轉換 196
6.6.3 信號流圖方程與狀態方程的關係 197
6.6.4 例題:正餘弦序列產生器 197
6.7 各種結構形式的比較 200
6.7.1 從計算的複雜性做比較 200
6.7.2 從減小係數量化影響做比較 200
複習思考題 200
習題 201
第7章 FIR濾波器設計 204
7.1 數字濾波器設計的基本要求 204
7.1.1 數字濾波器的設計基本步驟 204
7.1.2 數字濾波器的主要指標 204
7.1.3 為什麼以低通濾波器為典型 207
7.2 線性相位濾波器的特性 207
7.2.1 線性相位和符幅特性(Amplitude Response) 207
7.2.2 符幅特性的計算舉例 211
7.2.3 計算符幅特性的MATLAB程式 212
7.2.4 線性相位濾波器的零極點分布 214
7.3 設計FIR濾波器的窗函式法 215
7.3.1 對脈衝回響進行截斷的分析方法 215
7.3.2 矩形窗截斷的定量分析 216
7.3.3 幾種其他窗函式的介紹 219
7.3.4 MATLAB提供的窗函式 220
7.3.5 用窗函式設計FIR濾波器的步驟 222
7.4 利用頻率樣本法設計FIR濾波器 228
7.4.1 預期頻率特性的設定方法 228
7.4.2 頻率樣本法的設計過程 229
7.4.3 MATLAB中頻率樣本法的設計函式 232
7.5 FIR濾波器的最優設計法 233
7.5.1 設計FIR濾波器的離散最小二乘法 234
7.5.2 最小最大波動法設計濾波器 238
7.5.3 MATLAB中最優等波動濾波器的設計函式 241
7.6 FIR濾波器設計的一些深入問題 244
7.6.1 關於四種線性相位濾波器類型的套用範圍 244
7.6.2 數字微分器的設計 244
7.6.3 Hilbert變換器的設計 247
複習思考題 249
習題 250
第8章 IIR濾波器設計 253
8.1 關於模擬濾波器的預備知識 253
8.1.1 模擬濾波器的指標 253
8.1.2 以模擬濾波器為基礎的設計方法 255
8.2 模擬濾波器原型的設計 256
8.2.1 巴特沃斯低通濾波器 256
8.2.2 切比雪夫I型低通濾波器 261
8.2.3 切比雪夫II型低通濾波器 266
8.2.4 橢圓低通濾波器 268
8.2.5 濾波器原型的相位回響 270
8.3 由模擬濾波器變換為數字濾波器 270
8.3.1 脈衝回響不變法 270
8.3.2 雙線性變換法 276
8.4 模擬濾波器的頻帶變換 284
8.4.1 低通到低通的頻帶變換 285
8.4.2 低通到高通的頻帶變換 285
8.4.3 低通到帶通的頻帶變換 286
8.4.4 低通到帶阻的頻帶變換 287
8.4.5 信號處理工具箱中的頻帶變換函式 288
8.5 數字濾波器的頻帶變換 291
8.6 MATLAB中直接設計IIR數字濾波器的函式 296
8.7 IIR和FIR數字濾波器的比較 297
8.8 MATLAB中的濾波器設計工具(fdatool) 298
複習思考題 300
習題 301
第9章 綜合套用 303
9.1 信號處理工具(sptool)的介紹 303
9.2 語音數位化量化噪聲的改善 306
9.3 係數量化和運算量化的影響 310
9.3.1 數的浮點和定點表示方法 310
9.3.2 係數量化對FIR濾波器的影響舉例 311
9.3.3 係數量化對IIR濾波器的影響舉例 312
9.3.4 運算量化和溢出問題 313
9.3.5 濾波器設計工具中的量化處理 314
9.3.6 定點(Fix-point)工具箱 315
9.4 數位訊號處理在雙音撥號系統中的套用 316
9.4.1 雙音多頻撥號 316
9.4.2 Goertzel算法 316
9.4.3 檢測DTMF信號的DFT參數選擇 318
9.4.4 生成和檢測DTMF信號的演示程式 319
9.5 正餘弦信號的頻譜分析 320
9.6 音樂信號處理 326
9.6.1 時域處理方法 326
9.6.2 頻域處理方法 327
9.7 變採樣率數字濾波 329
9.7.1 信號的整數倍內插 329
9.7.2 信號的整數倍抽取 331
9.7.3 分數倍變採樣率和相應的MATLAB函式 332
9.7.4 整數倍抽取和內插的套用舉例 333
9.8 稀疏天線陣列設計 333
複習思考題 335
習題 336
第10章 信號處理從概念設計到晶片實現 338
10.1 嵌入式系統及其開發流程 338
10.2 系統仿真(Simulink)工具簡介 339
10.2.1 概述 339
10.2.2 模組庫及模型的建立 340
10.2.3 仿真的運行 341
10.3 MATLAB與嵌入式代碼的銜接 342
10.4 基於模型的開發流程 343
10.5 嵌入式代碼開發的一個實例 345
10.6 數位訊號處理的未來 349
附錄A 信號處理工具箱函式集 351
附錄B 本書自編子程式索引表 357
附錄C 例題程式索引表 358
附錄D 圖形生成和演示程式索引表 362
參考文獻 363