在訊號處理中,常需從帶有噪聲的訊號中提取原始訊號。一般用濾波處理方法濾除或削弱噪聲干擾以及其他不需要的訊號。對於疊加性組合訊號,可用線性濾波器將它們分離開。對於實用中常見的非疊加性 組合訊號(如乘積性訊號和褶積性訊號),靠線性濾波器分離或處理這些訊號分量往往是無效的,這時應採用非線性濾波,即要用同態濾波處理系統進行訊號處理。在輸入輸出運算相同的情況下,同態系統可分為相乘訊號的同態濾波處理和褶積訊號的同態濾波處理兩種。在許多實際問題中,訊號為兩個或多個分量的乘積(如在有衰落的傳輸通道中,衰落效應可看作一個緩變分量和傳輸訊號相乘)。對這類相乘訊號,如用線性系統來分離訊號各成分或單獨地改善某一訊號成分往往是無效的。但利用相乘訊號的同態濾波處理,就可以取得較好的濾波效果。在多徑或混響環境中進行通訊、定位或記錄,產生失真的效果可以看成是干擾與所需訊號的褶積,對這類訊號可用褶積訊號的同態濾波處理。在語音、影象、雷達、聲吶、地震勘探以及生物醫學工程等領域中,同態訊號處理獲得廣泛的應用。
同態系統服從廣義迭加原理
用符號 T[·]表示系統變換,用“囗”表示輸入訊號的組成分量的廣義迭加,用“:“表示輸入訊號向量與標量C之間的一種廣義乘法運算,用“○”和“”分別表示輸出訊號的組成分量的廣義相加和標乘的運算,於是這類同態系統可以用三個子系統D□相級聯組成(圖1)。其中第一個子系統D□具有如下性質: D□【x1(n)□x2(n)】=D□【x1(n)】+D□【x2(n)】 =憫1(n)+憫2(n) D□【C:x1(n)】=CD□【x1(n)】=C憫1(n) 可見,子系統D□服從廣義迭加原理,其輸入運算是“囗”,輸出運算是“+”。子系統D□ 的作用是將由訊號分量按運演算法則“囗”組合起來的訊號變換成按常規線性相加的訊號分量D□【x1(n) 】和D□【x2(n)】。因此,D□是一個同態系統。 系統L是一個普通的線性系統,於是有 L【憫1(n)+憫2(n)】=L【憫1(n)】+L【憫2(n)】 =尳1(n)+尳2(n) L【C憫1(n)】=CL【憫1(n)】=C尳1(n) 最後一個子系統 D劸是將系統L的輸出按運算規則“○”進行逆運算變換為系統輸出,有 D劸【尳1(n)+尳2(n)】=D劸【尳1(n)】+D劸【尳2(n)】 =y1(n)+y2(n) D劸【C尳(n)】=CD劸【尳(n)】=C尳(n) 由於子系統D□是按運算“囗”和“:“確定的一種特徵系統,因此稱為對運算“囗”的特徵系統。同樣,D0是對運算“○“的特徵系統。顯然,具有相同運算規則的第一子系統和第三子系統的所有同態系統,只是其線性系統部分有所不同。換言之,特徵系統確定之後,剩下的就是線性濾波問題了。例如,要從訊號x(n)=x1(n)□x2(n)中恢復出有用訊號x1(n),首先就要找出一個特徵系統D□,它能把x(n)變為 D□【x(n)】=D□【x1(n)】+D□【x2(n)】 =憫1(n)+憫2(n) 然後,適當選擇和設計線性系統L,只讓D□【x(n)】中的憫1(n)=D□【n1(n)】分量透過,理想的情況是 尳(n)=憫1(n)=D□【x1(n)】 然後,取D0=D□,得到逆特徵系統的輸出為 y(n)=D劸【尳(n)】=D劸{D□【x1(n)】}=x1(n) 因此,為了分離x1(n)和x2(n),必須用一個線性濾波器L來完全分離x1(n)和x2(n)。理想的分離取決於運演算法則“囗”以及輸入訊號分量x1(n)和x2(n)的性質。 同態濾波處理系統 在輸入和輸出運算相同的情況下,同態系統可分為相乘訊號的同態濾波處理和褶積訊號的同態濾波處理兩種。
相乘訊號的同態濾波處理
應用
在許多實際問題中,訊號為兩個或多個分量的乘積。例如,在有衰落的傳輸通道中,衰落效應可以看作是一個緩變分量和傳輸訊號相乘。又如,調幅訊號可表示為載頻訊號與包絡函式的乘積,在接收機內需要分離載波和包絡。在這一類相乘訊號中,用線性系統來分離訊號各成分或單獨地改善某一訊號成分往往是無效的。但利用相乘訊號的同態濾波處理,就可以取得較好的濾波效果。
假設信
,且對所有n值x1(n)>0和x2(n)>0,則有 log【x1(n)·x2(n)】=log{x1(n)+log【x2(n)】} 但是,輸入訊號x(n)不一定是正的,而且往往是複數訊號。這時,就要用到複對數函式,於是輸入和輸出均為乘法的同態系統,如圖2所示,其中序列x(n)、憫(n)、尳(n)及y(n)一般均為複數。 令x(n)=|x(n)|exp{jarg【x(n)】}表示一個複數序列,則x(n)的複對數為 log【x(n)】=log|x(n)|+jarg【x(n)】複對數log【x(n)】 的逆是復指數,即 y(n)=exp{log【x(n)】} =exp【log|x(n)|】·exp{jarg【x(n)】} 複對數的虛部arg【x(n)】加上2π的任意整數倍,並不改變上式的結果。因此,如果不另加限制,複對數並不是唯一性變換,會出現多值性問題。為此,必須選擇能消除模糊的arg【x(n)】。但是系統log【x(n)】服從廣義迭加原理,x(n)=x1(n)x2(n)必須存在下列關係 log|x(n)|=log|x1(n)|+log|x2(n)| arg【x(n)】=arg【x1(n)】+arg【x2(n)】 只要保證上述相角關係成立,就能消除arg【x(n)】中的模糊。 褶積訊號的同態濾波處理 在多徑或混響環境中進行通訊、定位或記錄,產生失真的效果可以看成是干擾與所需訊號的褶積。在語音訊號處理中,經常要分離激勵源與聲道衝激響應,至少在一段短時間內可以認為語音波形是由兩者的褶積構成的。地震記錄資料是地震子波與含有岩層結構資訊的反射係數序列的褶積組合。離散褶積組合訊號可以表示為 x(n)=x1(n)*x2(n) 利用同態濾波處理可以解褶積(圖3)。特徵系統D*具有如下性質:D*【x1(n)*x2(n)】=D*【x1(n)】+D*【x2(n)】 =憫1(n)+憫2(n) D*【C:x(n)】=CD*【x(n)】=C憫(n) 系統L是一個線性系統,D*-1是D*的逆系統。 褶積的Z變換是 x(Z)=x1(Z)x2(Z) 即Z變換運算也可以看作是一個系統,它的輸入運算為褶積,輸出運算是乘法的同態變換。因此,利用Z變換就能把褶積組合變換成乘法組合,從而可以利用上述乘法同態濾波系統來處理。但是,通常函式x(Z)是個複數,故不得不採用複對數運算。
在訊號處理中,常需從帶有噪聲的訊號中提取原始訊號。一般用濾波處理方法濾除或削弱噪聲干擾以及其他不需要的訊號。對於疊加性組合訊號,可用線性濾波器將它們分離開。對於實用中常見的非疊加性 組合訊號(如乘積性訊號和褶積性訊號),靠線性濾波器分離或處理這些訊號分量往往是無效的,這時應採用非線性濾波,即要用同態濾波處理系統進行訊號處理。在輸入輸出運算相同的情況下,同態系統可分為相乘訊號的同態濾波處理和褶積訊號的同態濾波處理兩種。在許多實際問題中,訊號為兩個或多個分量的乘積(如在有衰落的傳輸通道中,衰落效應可看作一個緩變分量和傳輸訊號相乘)。對這類相乘訊號,如用線性系統來分離訊號各成分或單獨地改善某一訊號成分往往是無效的。但利用相乘訊號的同態濾波處理,就可以取得較好的濾波效果。在多徑或混響環境中進行通訊、定位或記錄,產生失真的效果可以看成是干擾與所需訊號的褶積,對這類訊號可用褶積訊號的同態濾波處理。在語音、影象、雷達、聲吶、地震勘探以及生物醫學工程等領域中,同態訊號處理獲得廣泛的應用。
同態系統服從廣義迭加原理
用符號 T[·]表示系統變換,用“囗”表示輸入訊號的組成分量的廣義迭加,用“:“表示輸入訊號向量與標量C之間的一種廣義乘法運算,用“○”和“”分別表示輸出訊號的組成分量的廣義相加和標乘的運算,於是這類同態系統可以用三個子系統D□相級聯組成(圖1)。其中第一個子系統D□具有如下性質: D□【x1(n)□x2(n)】=D□【x1(n)】+D□【x2(n)】 =憫1(n)+憫2(n) D□【C:x1(n)】=CD□【x1(n)】=C憫1(n) 可見,子系統D□服從廣義迭加原理,其輸入運算是“囗”,輸出運算是“+”。子系統D□ 的作用是將由訊號分量按運演算法則“囗”組合起來的訊號變換成按常規線性相加的訊號分量D□【x1(n) 】和D□【x2(n)】。因此,D□是一個同態系統。 系統L是一個普通的線性系統,於是有 L【憫1(n)+憫2(n)】=L【憫1(n)】+L【憫2(n)】 =尳1(n)+尳2(n) L【C憫1(n)】=CL【憫1(n)】=C尳1(n) 最後一個子系統 D劸是將系統L的輸出按運算規則“○”進行逆運算變換為系統輸出,有 D劸【尳1(n)+尳2(n)】=D劸【尳1(n)】+D劸【尳2(n)】 =y1(n)+y2(n) D劸【C尳(n)】=CD劸【尳(n)】=C尳(n) 由於子系統D□是按運算“囗”和“:“確定的一種特徵系統,因此稱為對運算“囗”的特徵系統。同樣,D0是對運算“○“的特徵系統。顯然,具有相同運算規則的第一子系統和第三子系統的所有同態系統,只是其線性系統部分有所不同。換言之,特徵系統確定之後,剩下的就是線性濾波問題了。例如,要從訊號x(n)=x1(n)□x2(n)中恢復出有用訊號x1(n),首先就要找出一個特徵系統D□,它能把x(n)變為 D□【x(n)】=D□【x1(n)】+D□【x2(n)】 =憫1(n)+憫2(n) 然後,適當選擇和設計線性系統L,只讓D□【x(n)】中的憫1(n)=D□【n1(n)】分量透過,理想的情況是 尳(n)=憫1(n)=D□【x1(n)】 然後,取D0=D□,得到逆特徵系統的輸出為 y(n)=D劸【尳(n)】=D劸{D□【x1(n)】}=x1(n) 因此,為了分離x1(n)和x2(n),必須用一個線性濾波器L來完全分離x1(n)和x2(n)。理想的分離取決於運演算法則“囗”以及輸入訊號分量x1(n)和x2(n)的性質。 同態濾波處理系統 在輸入和輸出運算相同的情況下,同態系統可分為相乘訊號的同態濾波處理和褶積訊號的同態濾波處理兩種。
相乘訊號的同態濾波處理
應用
在許多實際問題中,訊號為兩個或多個分量的乘積。例如,在有衰落的傳輸通道中,衰落效應可以看作是一個緩變分量和傳輸訊號相乘。又如,調幅訊號可表示為載頻訊號與包絡函式的乘積,在接收機內需要分離載波和包絡。在這一類相乘訊號中,用線性系統來分離訊號各成分或單獨地改善某一訊號成分往往是無效的。但利用相乘訊號的同態濾波處理,就可以取得較好的濾波效果。
假設信
,且對所有n值x1(n)>0和x2(n)>0,則有 log【x1(n)·x2(n)】=log{x1(n)+log【x2(n)】} 但是,輸入訊號x(n)不一定是正的,而且往往是複數訊號。這時,就要用到複對數函式,於是輸入和輸出均為乘法的同態系統,如圖2所示,其中序列x(n)、憫(n)、尳(n)及y(n)一般均為複數。 令x(n)=|x(n)|exp{jarg【x(n)】}表示一個複數序列,則x(n)的複對數為 log【x(n)】=log|x(n)|+jarg【x(n)】複對數log【x(n)】 的逆是復指數,即 y(n)=exp{log【x(n)】} =exp【log|x(n)|】·exp{jarg【x(n)】} 複對數的虛部arg【x(n)】加上2π的任意整數倍,並不改變上式的結果。因此,如果不另加限制,複對數並不是唯一性變換,會出現多值性問題。為此,必須選擇能消除模糊的arg【x(n)】。但是系統log【x(n)】服從廣義迭加原理,x(n)=x1(n)x2(n)必須存在下列關係 log|x(n)|=log|x1(n)|+log|x2(n)| arg【x(n)】=arg【x1(n)】+arg【x2(n)】 只要保證上述相角關係成立,就能消除arg【x(n)】中的模糊。 褶積訊號的同態濾波處理 在多徑或混響環境中進行通訊、定位或記錄,產生失真的效果可以看成是干擾與所需訊號的褶積。在語音訊號處理中,經常要分離激勵源與聲道衝激響應,至少在一段短時間內可以認為語音波形是由兩者的褶積構成的。地震記錄資料是地震子波與含有岩層結構資訊的反射係數序列的褶積組合。離散褶積組合訊號可以表示為 x(n)=x1(n)*x2(n) 利用同態濾波處理可以解褶積(圖3)。特徵系統D*具有如下性質:D*【x1(n)*x2(n)】=D*【x1(n)】+D*【x2(n)】 =憫1(n)+憫2(n) D*【C:x(n)】=CD*【x(n)】=C憫(n) 系統L是一個線性系統,D*-1是D*的逆系統。 褶積的Z變換是 x(Z)=x1(Z)x2(Z) 即Z變換運算也可以看作是一個系統,它的輸入運算為褶積,輸出運算是乘法的同態變換。因此,利用Z變換就能把褶積組合變換成乘法組合,從而可以利用上述乘法同態濾波系統來處理。但是,通常函式x(Z)是個複數,故不得不採用複對數運算。