用MatLab進(jìn)行數(shù)字信號(hào)處理

附錄B 用MatLab進(jìn)行數(shù)字信號(hào)處理
一、 Matlab簡(jiǎn)介
在科學(xué)研究和工程應(yīng)用中，往往要進(jìn)行大量的數(shù)學(xué)計(jì)算，這些運(yùn)算一般來(lái)說(shuō)難以用手工精確和快捷地進(jìn)行，而要借助計(jì)算機(jī)編制相應(yīng)的程序做近似計(jì)算。目前流行用Basic、Fortran和c語(yǔ)言編制計(jì)算程序, 既需要對(duì)有關(guān)算法有深刻的了解，還需要熟練地掌握所用語(yǔ)言的語(yǔ)法及編程技巧。在本講義的附錄一中，提供了一部分算法的C代碼，但這些功能代碼段僅僅是在數(shù)字信號(hào)處理過(guò)程中一個(gè)很小的部分。一方面實(shí)驗(yàn)時(shí)間有限，另一方面，各個(gè)專(zhuān)業(yè)的同學(xué)對(duì)編程有不同程度的要求。往往實(shí)驗(yàn)的大部分時(shí)間用來(lái)調(diào)試程序，真正用于理解課程內(nèi)容的時(shí)間相對(duì)被壓縮。
為克服上述困難，美國(guó)Mathwork公司于1967年推出了"Matrix Laboratory"（縮寫(xiě)為Matlab）軟件包，并不斷更新和擴(kuò)充，是一種功能強(qiáng)、效率高便于進(jìn)行科學(xué)和工程計(jì)算的交互式軟件包。其中包括：一般數(shù)值分析、矩陣運(yùn)算、數(shù)字信號(hào)處理、建模和系統(tǒng)控制和優(yōu)化等應(yīng)用程序，并集應(yīng)用程序和圖形于一便于使用的集成環(huán)境中。在此環(huán)境下所解問(wèn)題的Matlab語(yǔ)言表述形式和其數(shù)學(xué)表達(dá)形式相同，不需要按傳統(tǒng)的方法編程。
雖然Matlab作為一種新的計(jì)算機(jī)語(yǔ)言，但由于使用Matlab編程運(yùn)算與人進(jìn)行科學(xué)計(jì)算的思路和表達(dá)方式完全一致，所以不象學(xué)習(xí)其它高級(jí)語(yǔ)言--如Basic、Fortran和C等那樣難于掌握，用Matlab編寫(xiě)程序猶如在演算紙上排列出公式與求解問(wèn)題。Matlab語(yǔ)言具有編程效率高，易學(xué)易懂；使用方便；繪圖方便等一系列優(yōu)點(diǎn)。
二、 MatLab基本繪圖
（一）繪圖函數(shù)簡(jiǎn)介
有關(guān)MatLab運(yùn)算、函數(shù)、邏輯控制等可以參考一些Matlab的書(shū)籍、或參閱多媒體通信實(shí)驗(yàn)室主頁(yè)中有關(guān)MatLab的部分，在本講義中不再贅述。這里介紹以下MatLab下二維圖形的繪制。下表是MatLab中不同繪圖函數(shù)的名稱(chēng)和功能簡(jiǎn)介：

（二）Stem函數(shù)繪圖
各種不同的繪圖函數(shù)分別適用于不同的場(chǎng)合，使用“stem”繪制針狀圖最簡(jiǎn)單，從附錄中提供的MatLab原代碼可以看出，只需要將需要繪制的數(shù)據(jù)存放在一個(gè)數(shù)組中，然后將這個(gè)數(shù)組作為參數(shù)傳遞給“stem”函數(shù)就可以得到輸出圖形。例如，下面的代碼可以繪制正弦函數(shù)的圖形：
n=1:50; %定義序列的長(zhǎng)度是50
x=zeros(1,50); %注意：MATLAB中數(shù)組下標(biāo)從1開(kāi)始
x(1)=1; %沖擊函數(shù)
stem(x); %繪制函數(shù)圖形
得到的函數(shù)圖形如下圖所示：