李亞文

（商洛學(xué)院 物理與電子信息工程系，陜西商洛 726000）

MATLAB的離散系統(tǒng)的數(shù)學(xué)理論很早已經(jīng)形成，直到20世紀(jì)90年代計(jì)算機(jī)應(yīng)用和發(fā)展，才使得其得到了廣泛的應(yīng)用[1]。離散系統(tǒng)的研究存在眾多科學(xué)領(lǐng)域，比如：信號(hào)處理與通訊、圖像處理、信號(hào)檢測(cè)、地質(zhì)勘探、道路檢測(cè)等[2-3]，前期通過(guò)對(duì)離散信號(hào)特性、離散系統(tǒng)的特性進(jìn)行數(shù)學(xué)演算分析，可以判斷系統(tǒng)是否是物理可實(shí)現(xiàn)的系統(tǒng)，并能進(jìn)一步優(yōu)化系統(tǒng)特性[4]。本文摒棄了常規(guī)使用數(shù)學(xué)數(shù)值計(jì)算或者羅斯準(zhǔn)則求解零極點(diǎn)，進(jìn)而分析系統(tǒng)特性；利用時(shí)域分析和頻域分析法研究離散系統(tǒng)響應(yīng)，借助MATLAB軟件通過(guò)分析系統(tǒng)函數(shù)特征,間接研究離散系統(tǒng)的因果穩(wěn)定性幅頻特性和相頻特性，為研究離散系統(tǒng)中的數(shù)字信號(hào)傳輸和離散信號(hào)檢測(cè)提供前提基礎(chǔ)。

1 離散系統(tǒng)

系統(tǒng)分為連續(xù)系統(tǒng)和離散系統(tǒng)，信號(hào)在系統(tǒng)中傳輸主要研究的物理量是：激勵(lì)、系統(tǒng)傳輸特性和系統(tǒng)響應(yīng)。其中，當(dāng)系統(tǒng)各個(gè)物理量隨時(shí)間變化的規(guī)律不能用連續(xù)函數(shù)描述時(shí)，而只在離散的瞬間給出數(shù)值，這種系統(tǒng)稱為離散系統(tǒng)。

如圖1為信號(hào)通過(guò)系統(tǒng)傳輸、處理的框圖。f(·)是系統(tǒng)的輸入（激勵(lì)），y(·)是系統(tǒng)的輸出（響應(yīng)），h(·)是系統(tǒng)特性的一種描述?！啊ぁ笔切盘?hào)的自變量，可以是連續(xù)變量t，也可以是離散變量n。

圖1 系統(tǒng)分析框圖

由于離散系統(tǒng)精度高、自適應(yīng)性靈活、抗干擾能力強(qiáng)、可靠性好、集成化程度高且成本低[2]，因此其在通信、交通、航天航空、生物醫(yī)學(xué)、遙感等領(lǐng)域得到了廣泛的應(yīng)用。

1.1 LTI離散系統(tǒng)時(shí)域響應(yīng)分析

LTI的離散系統(tǒng)滿足可分解、線性以及非時(shí)變特性。線性時(shí)不變連續(xù)系統(tǒng)是由常系數(shù)微分方程描述的，而線性時(shí)不變離散系統(tǒng)是由常系數(shù)差分方程描述的。在差分方程中構(gòu)成方程的各項(xiàng)包含有未知離散變量的 y(n)，以及 y(n+1)，y(n+2)，…，y(n-1)，y(n-2)，…。

圖2 離散系統(tǒng)響應(yīng)框圖

N階LTI離散系統(tǒng)的數(shù)學(xué)模型是常系數(shù)N階線性差分方程一般形式[1]

a0y(n)+a1y(n-1)+…+aNy(n-N)=b0x(n)+b1x(n-1)+…+bmx(n-M)或者表示為

分析N階LTI離散系統(tǒng)時(shí)域響應(yīng)就是求解N階線性差分方程的過(guò)程。時(shí)域分析法是將離散系統(tǒng)的響應(yīng)劃分為零輸入響應(yīng)和零狀態(tài)響應(yīng)[3]，即y(k)=yzi(k)+yzs(k)，其中零輸入響應(yīng)是齊次差分方程的解，零狀態(tài)響應(yīng)可由卷積積分的方法求解,

1.2 LTI離散系統(tǒng)頻域分析

直接求解離散時(shí)間信號(hào)通過(guò)系統(tǒng)的響應(yīng)較為復(fù)雜，因此，引入頻域分析去求解離散系統(tǒng)的時(shí)域響應(yīng)，也稱為Z域響應(yīng)。Z變換的定義[4]可由抽樣信號(hào)的拉氏變換引出，表示為

在Z域中分析N階離散系統(tǒng)的時(shí)域響應(yīng)，步驟如下：

2）對(duì)N階的離散差分方程兩邊同時(shí)取Z變換：

其中y(l)是初始條件。

3）代入系統(tǒng)的初始條件y(l)，進(jìn)而可得系統(tǒng)響應(yīng) Y(z)。

4）對(duì)系統(tǒng)響應(yīng)Y(z)進(jìn)行逆Z變換，最終可得離散系統(tǒng)的時(shí)域響應(yīng)y(n)。

2 離散系統(tǒng)特性分析

2.1 系統(tǒng)函數(shù)

當(dāng)外加的激勵(lì)信號(hào)x(n)是因果序列時(shí)，并且系統(tǒng)初始條件為零（y(l)=0,-N≤l≤-1），因此，系統(tǒng)的響應(yīng)就是零狀態(tài)響應(yīng)，即：Y(z)=Yzs(z)。

稱H(z)為系統(tǒng)（傳輸）函數(shù)，定義LTI離散系統(tǒng)輸出Z變換與輸入Z變換之比為系統(tǒng)函數(shù)[5]

根據(jù)系統(tǒng)函數(shù)對(duì)應(yīng)的單位沖擊響應(yīng)，將離散系統(tǒng)劃分為：有限脈沖響應(yīng)FIR（finite impulse response）和無(wú)線脈沖響應(yīng)IIR[6-7]（infinite impulse response）。

2.2 離散系統(tǒng)特性

H(z)與h(n)是一對(duì)Z變換對(duì)，所以只要知道H(z)在z平面上的零、極點(diǎn)分布情況,從而得到系統(tǒng)的脈沖響應(yīng)h(n)的變化規(guī)律。系統(tǒng)函數(shù)不但能反應(yīng)系統(tǒng)的內(nèi)部結(jié)構(gòu)，并且能較好的反應(yīng)系統(tǒng)的因果穩(wěn)定性[7]。

1）因果性

由因果系統(tǒng)的時(shí)域條件n＜0，h(n)=0及H(z)的定義，可知因果系統(tǒng)的H(z)只有z的負(fù)冪項(xiàng)，其收斂區(qū)為RH-＜|z|≤∞。因此收斂區(qū)包含無(wú)窮時(shí)，必為因果系統(tǒng)。

2）穩(wěn)定性

3 離散系統(tǒng)分析

由系統(tǒng)的電路模型，根據(jù)電路關(guān)系中電流KCL和電壓KVL關(guān)系得到輸入外加激勵(lì)和系統(tǒng)響應(yīng)的數(shù)學(xué)關(guān)系，進(jìn)而求解系統(tǒng)響應(yīng)；引入頻域分析，研究離散系統(tǒng)的系統(tǒng)傳輸函數(shù)分析零極點(diǎn)分布進(jìn)而分析離散系統(tǒng)因果穩(wěn)定性，系統(tǒng)響應(yīng)的幅頻特性和相頻特性[8-9]。舉例分析一個(gè)離散系統(tǒng)。

某離散系統(tǒng)的數(shù)學(xué)模型

通過(guò)分析系統(tǒng)函數(shù)，則可研究離散系統(tǒng)特性，由離散系統(tǒng)的數(shù)學(xué)模型，對(duì)差分方程兩邊同時(shí)取Z變換，得到系統(tǒng)傳輸函數(shù)

繪制零極點(diǎn)圖，matlab部分關(guān)鍵語(yǔ)句

B=[0.30.40.20.30.4];

A=[1-1.21.2-0.60.4];

[r,p,k]=residue(B,A)

zplane(B,A)

p1=abs(p)

a1=angle(p/pi)

零極點(diǎn)圖分布如圖3所示。

圖3 零極點(diǎn)圖分布

分析零極點(diǎn)圖圖3發(fā)現(xiàn)，極點(diǎn)都位于單位圓內(nèi)，因此該離散系統(tǒng)是個(gè)因果穩(wěn)定的系統(tǒng)[10]。

可以看出，通過(guò)判斷系統(tǒng)極點(diǎn)是否位于單位圓內(nèi)來(lái)判斷系統(tǒng)的穩(wěn)定性，對(duì)于一個(gè)復(fù)雜的系統(tǒng)來(lái)說(shuō)，求極點(diǎn)并不容易，以往總是使用羅斯（Routh）準(zhǔn)則[4]判斷連續(xù)系統(tǒng)的穩(wěn)定性，利用勞斯（Jury）準(zhǔn)則[4]判斷離散系統(tǒng)的穩(wěn)定性。利用matlab軟件分析系統(tǒng)穩(wěn)定性，只需在零極點(diǎn)圖上觀察即可，判斷極點(diǎn)是否位于s的右半平面，或是否有極點(diǎn)在Z平面的單位圓外。

為了分析系統(tǒng)時(shí)域響應(yīng)，因此研究系統(tǒng)的幅頻特性和相頻特性，如圖4所示。

繪制幅頻相頻特性，matlab關(guān)鍵語(yǔ)句

B=[0.30.40.20.30.4];

A=[1-1.21.2-0.60.4];

[H,w]=freqz(b,a,1000,’whole’);

Subplot(2,1,1); plot(w/pi,abs(H));

Ylabel(‘|H|’);title(‘幅頻特性’)

Subplot(2,1,2); plot(w/pi,angle(H));

Ylabel(‘a(chǎn)ng|H|’);title(‘相頻特性’)

Xlabel（‘omega/pi’）;

圖4 幅頻特性和相頻特性

通過(guò)這個(gè)離散系統(tǒng)例子可以看出，借助matlab分析離散系統(tǒng)響應(yīng)和系統(tǒng)特性比較方便，而且以圖形形式顯示，較為直觀,比傳統(tǒng)的數(shù)值計(jì)算分析系統(tǒng)特性更為精確。

4 結(jié)語(yǔ)

離散系統(tǒng)研究成為通信系統(tǒng)和信號(hào)處理方面研究的核心內(nèi)容，既使是連續(xù)系統(tǒng)也是將其以適當(dāng)?shù)牟介L(zhǎng)進(jìn)行時(shí)間離散后，在計(jì)算機(jī)存儲(chǔ)精度下的離散系統(tǒng)，因此實(shí)際的通信和信號(hào)處理都是離散信號(hào)和離散系統(tǒng)的研究。本文主要分析了求解離散系統(tǒng)時(shí)域響應(yīng)的兩種方法，即時(shí)域分析法和頻域分析法，通過(guò)分析離散系統(tǒng)系統(tǒng)的特性。借助matlab分析離散系統(tǒng)使過(guò)程變得簡(jiǎn)便，并且圖形展示直觀形象，這種方法準(zhǔn)確率高，簡(jiǎn)單易于實(shí)現(xiàn)，可操作性強(qiáng)，為研究數(shù)字通信系統(tǒng)的研究和離散信號(hào)分析提供前提條件。

[1]張德豐.matlab在數(shù)字信號(hào)處理中的應(yīng)用[M].北京:清華大學(xué)出版社,2010:13-15.

[2]鄭君里,應(yīng)啟珩,楊為理.信號(hào)與系統(tǒng)[M].北京:高等教育出版社,2000:35-46.

[3]陳懷琛.MATLAB及其在理工課程中的應(yīng)用指南[M].西安:西安電子科技大學(xué)出版社,2000:9-22.

[4]AV奧本海姆,AS威斯基.信號(hào)與系統(tǒng)[M].劉樹(shù)棠,譯.西安:西安交通大學(xué)出版社,1998:75-84.

[5]師五喜,霍 偉,吳宏鑫.一類未知非線性離散系統(tǒng)的直接自適應(yīng)模糊預(yù)測(cè)控制[J].自動(dòng)化學(xué)報(bào),2004,30(5):664-670.

[6]于 馳,王宏偉,宮明龍.多輸入多輸出離散系統(tǒng)的輸出反饋滑模控制[J].東北大學(xué)學(xué)報(bào):自然科學(xué)版,2012,33(2):89-97.

[7]WANG Y Z,ZHANG H G,YANG J.Robust adaptive control for interval time-delay systems[J].Journal of Control Theroy and Application,2006,4:367-371.

[8]肖衛(wèi)初,陳偉宏.MATLAB在IIR濾波器設(shè)計(jì)與仿真中的應(yīng)用[J].湖南工程學(xué)院學(xué)報(bào):自然科學(xué)版,2005,15(1):64-71.

[9]李吉祥,武俊峰,樊麗穎.離散系統(tǒng)降階控制器的設(shè)計(jì)[J].電機(jī)與控制學(xué)報(bào),2011,15(6):78-83.

[10]黨宏社.信號(hào)系統(tǒng)實(shí)驗(yàn)[M].MATLAB版.西安:西安電子科技大學(xué)出版社,2007:98-115.