陳紹榮，朱行濤，徐 舜，沈建國

（陸軍工程大學通信士官學校，重慶 400035）

0 引言

在實時信號處理中，通常用FFT的算法來完成線性卷和運算，因此需要對長序列進行分段，其分段方法有無重疊等長度分段和重疊等長度分段，對應的兩種計算線性卷和的方法，分別是重疊相加法和重疊保留法。著作[1-3]用圖解的方式導出了重疊相加法和重疊保留法，既不嚴謹又不便于理解。本文基于著作[4]，利用矩形窗函數(shù)對序列的分段進行了描述，給出了一種導出重疊相加法和重疊保留法的簡易方法，揭示了兩種方法的聯(lián)系和區(qū)別。解決了兩種方法在實時信號處理中的信號拼接問題，給出了重疊保留法在快速計算自相關(guān)函數(shù)方面的應用。

1 圓周卷和代替線性卷和的條件

設(shè)矩形窗函數(shù)為：

式中，ε(n)為單位階躍序列。

利用矩形窗函數(shù)，起點為n1（n1為整數(shù)），長度為N1的有限長序列x1(n)可表示為：

同理，起點為n2（n2為整數(shù)），長度為N2的有限長序列x2(n)可表示為：

有限長序列x1(n)和x2(n)的線性卷和yl(n)表示為：

考慮到式（2）、式（3）及式（4），則有：

考慮到式（5）及式（6），則有：

由式（7）可知，線性卷和yl(n)的序列長度為：

結(jié)論1：

兩個有限長序列的線性卷和yl(n)是一個有限長序列，該有限長序列始于兩個有限長序列的始點之和，止于兩個有限長序列的終點之和，其序列長度為兩個有限長序列的長度相加再減1。

利用矩形窗函數(shù)，考慮到式（7）及式（8），則線性卷和yl(n)可以表示為：

設(shè)周期為N的周期沖激序列為：

在式（2）及式（3）描述的有限長序列x1(n)及x2(n)之后，分別補上N-N1及N-N2個零值點，再延拓成周期為N的周期序列，即：

由兩個周期序列的周期卷和的定義，并考慮到周期序列在一周期內(nèi)的求和與起始位無關(guān)及式（11），則兩個周期序列的周期卷和可以表示為：

結(jié)論2：

兩個周期序列作周期卷和等價于一個周期序列任一周期內(nèi)的序列和另一個周期序作線性卷和，也等價于對兩個周期序列任一周期內(nèi)的序列的線性卷和作同周期的周期延拓。

將兩個有限長序列x1(n)及x2(n)之后，先分別補上N-N1及N-N2個零值點，得到了兩個長度為N的序列，再分別作周期為N的周期延拓，得到了周期為N的周期序列的周期卷和取主值，即取區(qū)間[0,N-1]對應的值，定義為N點序列x1(n)及x2(n)的N點圓周卷和y(n)。

考慮到兩個N點序列的N點圓周卷和的定義及式（13），則N點序列x1(n)及x2(n)的N點圓周卷和y(n)可以表示為：

結(jié)論3：

兩個N點序列作N點圓周卷和，等價于兩個N點序列的線性卷和作N點的周期延拓，再取其主值序列。

考慮到式（15）及式（9），若滿足條件：

則有：

特別地

若滿足條件：

則有：

若滿足條件：

則有：

結(jié)論4：

若序列x1(n)=x1(n)RN1(n)，x2(n)=x2(n)RN2(n)，當N=N1+N2-1 時，則滿足式（21）的條件，那么可用序列x1(n)和x2(n)的N點圓周卷和y(n)來代替序列x1(n)和x2(n)的線性卷和yl(n)。

若序列x1(n)=x1(n)RN1(n)，x2(n)=x2(n)RN2(n)，當 式（21）的條件得不到滿足時，由式（14）可知，將序列x1(n)和x2(n)的線性卷和yl(n)作周期延拓，即yl(n)*δN(n)=yl(n)RN1+N2-1(n)*δN(n)，那么，每周內(nèi)將出現(xiàn)重疊。雖然增加點數(shù)N可避免序列x1(n)和x2(n)的N點圓周卷和y(n)出現(xiàn)重疊，但是使用FFT來計算N點圓周卷和y(n)時，將會增加運算量。

由于序列x(n)通過單位沖激響應為h(n)的LSI離散時間系統(tǒng)時，系統(tǒng)的零狀態(tài)響應yzs(n)可以表示成yzs(n)=x(n)*h(n)，有時序列x(n)可能會很長，如果將序列x(n)存儲完畢后，再做線性卷和將產(chǎn)生兩個問題：一是要求計算機的存儲量過大，二是要等待序列x(n)的輸入時間過長，即使h(n)是有限長的，也不能實現(xiàn)信號的“實時處理”。為了解決這一問題，如果h(n)的長度為N2，選取長度N1，保證M為最小的正整數(shù)時，滿足條件N=N1+N2-1=2M，將長度為L的輸入序列x(n)按等長度N1進行分段，則第i段的序列可以表示為xi(n)=x(n)RN1(n-iN1)，i=0,1,2,…,，其中表示大于或等于L/N1的正整數(shù)。將序列xi(n)及h(n)做線性卷和運算得到序列yi(n)，然后將各段序列yi(n)按一定的規(guī)則相加，即可得到完整的零狀態(tài)響應yzs(n)，這樣做有可能實現(xiàn)信號的實時處理。

將長序列分成短序列做線性卷和有兩種方法，一是重疊相加法，二是重疊保留法。

2 重疊相加法

設(shè)序列x(n)=x(n)RL(n)，h(n)=h(n)RN2(n)，L和N2為正整數(shù)，并且L＞＞N2。

選取長度N1，保證M為最小正整數(shù)時，滿足條件N=N1+N2-1=2M，將序列x(n)按等長度N1分段，則有：

式中：

考慮到式（23）及式（24），則有：

由 于N=N1+N2-1=2M，h(n)=h(n)RN2(n)，xi(n+iN1)=x(n+iN1)RN1(n)，由結(jié)論4 可知，可用序列x(n+iN1)RN1(n)和h(n)RN2(n)的N點圓周卷和來代替線性卷和。那么，式（25）中的yi(n)可表示為：

考慮到式（26），則可以利用FFT來計算LSI離散時間系統(tǒng)的零狀態(tài)響應yzs(n)，其實現(xiàn)步驟為：

（1）計算N 點FFT

（2）計算N 點FFT

（3）相乘

（4）計算N 點IFFT

（5）將各段yi(n)按式（25）相加，即按式（25）完成LSI 離散時間系統(tǒng)的零狀態(tài)響應信號的拼接。

討論：

由結(jié)論1 可知，第i段xi(n)和h(n)RN2(n)的線性卷和yi(n-iN1)=xi(n)*h(n)RN2(n)和h(n)RN2(n)中的變量n的取值范圍為：

由結(jié)論1 可知，第i-1 段xi-1(n)和h(n)RN2(n)的線性卷和y(i-1)[n-(i-1)N1]=xi-1(n)*h(n)RN2(n)中的變量n的取值范圍為：

由式（31）及式（32）可知，第i-1 段的線性卷和y(i-1)[n-(i-1)N1]的末端與第i段的線性卷和yi(n-iN1)的前端將出現(xiàn)重疊，其重疊的長度為：

結(jié)論5：

由式（33）可知，按式（25）的重疊相加法來計算LSI 離散時間系統(tǒng)的零狀態(tài)響應時，相鄰段的線性卷和將出現(xiàn)N2-1 點的重疊相加運算，故得名重疊相加法。

考慮到h(n)=h(n)RN2(n)，則有：

于是：

式中，i=1,2,…,。

由式（35）可知，yzs(iN1)由xi(n)的首位數(shù)據(jù)xi(iN1)=x(iN1)及首位之前的N2-1 個點的輸入數(shù)據(jù)x(m)(m=iN1-N2+1,iN1-N2+2,…,iN1-1)

共同確定。

結(jié)論6：

用對輸入x(n)等長度N1分段的線性卷和來計算LSI 離散時間系統(tǒng)的零狀態(tài)響應時，出現(xiàn)重疊相加運算的本質(zhì)是第i段首位之前的N2-1 點的輸入數(shù)據(jù)x(m)，對第i段前N2-1 點的線性卷和有影響。

3 重疊保留法

由于用對x(n)等長度N1分段的線性卷和來計算LSI 離散時間系統(tǒng)的零狀態(tài)響應時，出現(xiàn)了重疊相加運算，為了避免重疊相加運算，可以將第i段首位之前的N2-1 點輸入數(shù)據(jù)x(n)保留下來與第i段的序列xi(n) 依次銜接起來，構(gòu)成一個長度為N=N1+N2-1 的 新 序 列pi(n)(i=1,2,…,，考慮到x(n)=x(n)RL(n)，只須在x0(n)=x(n)RN1(n)段的前端補上N2-1 個零值點，便可構(gòu)成新序列p0(n)，這樣各段pi(n)均為N1+N2-1 點的序列。即：

式中，i=1,2,…,。

考慮到式（36）及式（37），則長度為N1+N2-1 的序列pi(n)和LSI 離散時間系統(tǒng)單位沖激響應h(n)的線性卷和yli(n)(i=0,1,2,…,，可以表示為：

下面確定n在什么范圍內(nèi)取值時，才能使yli(n)=yzs(n)。

對比式（38）和式（34）可知，為使線性卷和yli(n)=yzs(n)，則要求：

式中，m為整數(shù)。

顯然，若滿足條件：

則式（39）恒成立。

即，若滿足條件：

則式（39）恒成立。

由式（38）和式（34）可知，mmin=n-(N2-1)，mmax=n。

于是式（41）可以表示為：

解得：

由式（43）可知，其區(qū)間長度為：

考慮到式（43）及式（44），則有：

考慮到式（45），則有：

分析表明，將各段pi(n)和LSI 離散時間系統(tǒng)單位沖激響應h(n)的線性卷和yli(n)，按式（46）逐段銜接（或拼接）起來，就得到了LSI 離散時間系統(tǒng)的零狀態(tài)響應yzs(n)。

下面的任務是考察是否可用圓周卷和來替代線性卷和yli(n)。

考慮到h(n)=h(n)RN2(n)，式（36）及式（37），并注意到N=N1+N2-1=2M，由結(jié)論1 可知，pi(n)和LSI 離散時間系統(tǒng)單位沖激響應h(n)的線性卷和yli(n)(i=0,1,2,…,，可以表示為：

于是：

即：

考慮到式（49），由結(jié)論3 可知，pi(n+iN1)和h(n)的N點圓周卷和yi(n+iN1)可以表示為：

由式（50）可得：

考慮到式（51），則式（46）可以寫成：

考慮到式（50），則可以利用FFT來計算LSI離散時間系統(tǒng)的零狀態(tài)響應yzs(n)，其實現(xiàn)步驟為：

（1）計算N 點FFT

（2）基于式（36）及式（37），計算N 點FFT

（3）相乘

（4）計算N 點IFFT

（5）將各段yi(n+iN1)按式（52）逐段銜接（或拼接）起來，就得到了LSI 離散時間系統(tǒng)的零狀態(tài)響應yzs(n)。

由式（36）及式（37）可知，將第i段首位之前的N2-1 點輸入數(shù)據(jù)x(n)保留下來，與第i段xi(n)依次銜接起來構(gòu)成了pi(n)，從而使pi-1(n)的末端與pi(n)的前端N2-1 點的數(shù)據(jù)相同，即pi-1(n)與pi(n)之間有N2-1 點的數(shù)據(jù)重疊，故得名重疊保留法。由式（50）可知，pi(n+iN1)和h(n)的N點圓周卷和yi(n+iN1)末端出現(xiàn)了N2-1 點的重疊數(shù)據(jù)，由式（51）可知，將pi(n+iN1)和h(n)的N點圓周卷和yi(n+iN1)末端N2-1 點的重疊數(shù)據(jù)舍去，才保證了pi(n+iN1)和h(n)的圓周卷和等于線性卷和，故又得名重疊舍去法。

4 重疊保留法在快速計算自相關(guān)方面的應用

設(shè)序列x(n)=x(n)RL(n)，先將x(n)按長度N2分段，則第j段可以表示為：

式中，j=0,1,2,…,。

于是：

現(xiàn)在，構(gòu)成各段pi(n)均為N=N1+N2-1=2M點的序列，即：

式中，i=1,2,…,。

令：

考慮到式（57），則式（61）可以表示為：

式（62）表明，hj(n)是一個長度為N2點的序列，并且滿足hj(n)=hj(n)RN2(n)。

考慮到序列自相關(guān)函數(shù)的定義及式（58），則序列x(n)的自相關(guān)函數(shù)Rx(n)可以表示為：

令：

則有：

考慮到式（61），則式（65）可以寫成：

式中，i=0,1,2,…,-1。

由于式（66）與式（47）同形，考慮到式（52），則有：

式 中，Rpixj[n+iN1-jN2-(N2-1)]為pi(n+iN1) 和hj(n)的N點圓周卷和。

考慮到式（66），由結(jié)論3 可知：

考慮到式（67），則有：

將式（69）代入式（63），可得序列x(n)的自相關(guān)函數(shù)Rx(n)，即：

考慮到式（68），則可以利用FFT 來計算序列x(n)的自相關(guān)函數(shù)Rx(n)，其實現(xiàn)步驟為：

（1）基于式（59）及式（60），計算N 點FFT

（2）基于式（61），計算N 點FFT

（3）相乘

（4）計算N 點IFFT

（5）將各段Rpixj[n+iN1-jN2-(N2-1)]，按式（70）相加，即完成自相關(guān)函數(shù)的拼接，就得到了自相關(guān)函數(shù)Rx(n)。

5 結(jié)語

本文利用矩形窗函數(shù)對輸入序列的等長度N1分段進行了描述，分析了用對輸入x(n)等長度N1分段的線性卷和來計算LSI 離散時間系統(tǒng)的零狀態(tài)響應時，出現(xiàn)重疊相加運算的原因。不僅給出了一種導出重疊相加法和重疊保留法的簡易方法，而且揭示了重疊相加法和重疊保留法的聯(lián)系和區(qū)別。解決了兩種方法在實時信號處理中的信號拼接問題，最后給出了重疊保留法在快速計算自相關(guān)函數(shù)方面的應用。