TWI876412B - Compensated composite filter and signal processing method thereof - Google Patents
Compensated composite filter and signal processing method thereof Download PDFInfo
- Publication number
- TWI876412B TWI876412B TW112124987A TW112124987A TWI876412B TW I876412 B TWI876412 B TW I876412B TW 112124987 A TW112124987 A TW 112124987A TW 112124987 A TW112124987 A TW 112124987A TW I876412 B TWI876412 B TW I876412B
- Authority
- TW
- Taiwan
- Prior art keywords
- filter
- group delay
- prototype
- transfer function
- shaping filters
- Prior art date
Links
Images
Landscapes
- Networks Using Active Elements (AREA)
Abstract
Description
本揭示內容係關於一種補償式複合濾波器及其訊號處理方法,特別是關於一種具有線性相位響應的補償式複合濾波器及其訊號處理方法。The present disclosure relates to a compensating composite filter and a signal processing method thereof, and in particular to a compensating composite filter with linear phase response and a signal processing method thereof.
配置具有完美頻率響應(即平坦通帶、尖銳過渡帶、高抑制阻帶及線性相位響應)的濾波器始終是任何數位訊號處理器(Digital Signal Processor,以下簡稱DSP)從業者的最終目標。於習知技術中,具有線性相位響應的高性能濾波可透過有限脈衝響應(Finite Impulse Response,以下簡稱FIR)濾波器實現。Configuring a filter with perfect frequency response (i.e., flat passband, sharp transition band, highly suppressed stopband, and linear phase response) is always the ultimate goal of any digital signal processor (DSP) practitioner. In the prior art, high-performance filtering with linear phase response can be achieved through finite impulse response (FIR) filters.
然而,FIR濾波器含有許多的加法器和乘法器,導致DSP的系統複雜度提高。與FIR濾波器相比,無限脈衝響應(Infinite Impulse Response,IIR)濾波器透過自身的遞歸結構可基於低系統複雜度來實現高性能濾波,但也由於遞歸結構,使得IIR濾波器無法提供線性相位響應,並限制IIR濾波器的適用性。由此可知,目前市場上缺乏一種具有低複雜度且線性相位響應的濾波器及其訊號處理方法,故相關業者均在尋求其解決之道。However, the FIR filter contains many adders and multipliers, which increases the system complexity of the DSP. Compared with FIR filters, infinite impulse response (IIR) filters can achieve high-performance filtering based on low system complexity through their own recursive structure. However, due to the recursive structure, IIR filters cannot provide linear phase response and limit the applicability of IIR filters. It can be seen from this that there is currently a lack of a filter and signal processing method with low complexity and linear phase response on the market, so relevant industries are looking for solutions.
因此,本揭示內容之目的在於提供一種補償式複合濾波器及其訊號處理方法。補償式複合濾波器不僅具有低複雜度的電路結構,且還能在通帶上實現線性相位響應。補償式複合濾波器之訊號處理方法可完整地保留輸入訊號的波形,並同時過濾干擾訊號。Therefore, the purpose of the present disclosure is to provide a compensating composite filter and a signal processing method thereof. The compensating composite filter not only has a low-complexity circuit structure, but also can achieve a linear phase response in the passband. The signal processing method of the compensating composite filter can completely retain the waveform of the input signal and filter the interference signal at the same time.
依據本揭示內容的結構態樣的一實施方式提供一種補償式複合濾波器,其包含一原型濾波器以及一補償濾波器。原型濾波器用以將一輸入訊號乘以一原型轉移函數而產生一原型濾波輸出訊號,且原型濾波器於一通帶中具有一第一群延遲。補償濾波器串聯連接原型濾波器,且用以將原型濾波輸出訊號乘以一補償轉移函數而產生一輸出訊號。補償濾波器包含級聯連接的複數整形濾波器,並於通帶中具有一第二群延遲。第一群延遲相加第二群延遲而產生一複合群延遲,複合群延遲為一常數。According to an implementation method of the structural aspect of the present disclosure, a compensated composite filter is provided, which includes a prototype filter and a compensation filter. The prototype filter is used to multiply an input signal by a prototype transfer function to generate a prototype filter output signal, and the prototype filter has a first group delay in a passband. The compensation filter is connected in series with the prototype filter, and is used to multiply the prototype filter output signal by a compensation transfer function to generate an output signal. The compensation filter includes cascade-connected complex shaping filters and has a second group delay in the passband. The first group delay is added to the second group delay to produce a composite group delay, which is a constant.
藉此,本揭示內容的補償式複合濾波器透過將原型濾波器串聯連接含有級聯連接的多個整形濾波器的補償濾波器,並利用補償濾波器的第二群延遲補償原型濾波器的第一群延遲,達到線性相位響應且同時具有低複雜度。Thus, the compensated composite filter of the present disclosure achieves linear phase response and low complexity by connecting the prototype filter in series with the compensating filter including a plurality of cascade-connected shaping filters and using the second group delay of the compensating filter to compensate for the first group delay of the prototype filter.
前述實施方式之其他實施例如下:前述原型轉移函數與補償轉移函數相乘而產生一複合轉移函數。原型轉移函數表示為 ,補償轉移函數表示為 ,複合轉移函數表示為 ,且複合轉移函數符合下式: ;其中 為一變數, 為此些整形濾波器的一數量, 為此些整形濾波器之一者的一整形轉移函數, 為正整數且介於1至 之間。 Other embodiments of the above-mentioned implementation are as follows: The above-mentioned prototype transfer function and the compensation transfer function are multiplied to generate a composite transfer function. The prototype transfer function is represented as , the compensation transfer function is expressed as , the composite transfer function is expressed as , and the composite transfer function conforms to the following formula: ;in is a variable, is the number of these shaping filters, is an integer transfer function of one of these shaping filters, Is a positive integer between 1 and between.
前述實施方式之其他實施例如下:前述第一群延遲表示為 ,第二群延遲表示為 ,複合群延遲表示為 ,且複合群延遲符合下式: ;其中 為一角頻率, 為此些整形濾波器的一數量, 為此些整形濾波器之一者的一群延遲, 為正整數且介於1至 之間。 Other embodiments of the above-mentioned implementation are as follows: The above-mentioned first group delay is expressed as , the second group delay is expressed as , the composite group delay is expressed as , and the complex group delay meets the following formula: ;in is the angular frequency, is the number of these shaping filters, is the group delay of one of these shaping filters, Is a positive integer between 1 and between.
前述實施方式之其他實施例如下:前述此些整形濾波器之此者的群延遲符合下式: ;其中, 為此些整形濾波器之此者的一濾波器參數, 為此些整形濾波器之此者的另一濾波器參數;其中此些整形濾波器的複數個群延遲彼此相異。 Other embodiments of the aforementioned embodiment are as follows: The group delay of one of the aforementioned shaping filters conforms to the following formula: ;in, is a filter parameter for one of these shaping filters, is another filter parameter of one of these shaping filters; wherein a plurality of group delays of these shaping filters are different from each other.
前述實施方式之其他實施例如下:前述原型濾波器為一巴特沃斯濾波器(Butterworth filter)、一橢圓濾波器(Elliptic filter)、一無限脈衝響應濾波器(Infinite Impulse Response filter,IIR filter)或一切比雪夫濾波器(Chebyshev filter),且各整形濾波器為一全通濾波器(All-pass filter)。Other embodiments of the aforementioned implementation method are as follows: the aforementioned prototype filter is a Butterworth filter, an elliptic filter, an infinite impulse response filter (IIR filter) or a Chebyshev filter, and each shaping filter is an all-pass filter.
依據本揭示內容的方法態樣的一實施方式提供一種補償式複合濾波器之訊號處理方法,其包含以下步驟:一第一濾波步驟以及一第二濾波步驟。第一濾波步驟係驅動一原型濾波器接收一輸入訊號,並將輸入訊號乘以一原型轉移函數而產生一原型濾波輸出訊號,且原型濾波器於一通帶中具有一第一群延遲。第二濾波步驟係驅動一補償濾波器接收原型濾波輸出訊號,並將原型濾波輸出訊號乘以一補償轉移函數而產生一輸出訊號,且補償濾波器包含級聯連接的複數整形濾波器,並於通帶中具有一第二群延遲。第一群延遲相加第二群延遲而產生一複合群延遲,複合群延遲為一常數。According to an embodiment of the method aspect of the present disclosure, a signal processing method of a compensating composite filter is provided, which includes the following steps: a first filtering step and a second filtering step. The first filtering step is to drive a prototype filter to receive an input signal and multiply the input signal by a prototype transfer function to generate a prototype filter output signal, and the prototype filter has a first group delay in a passband. The second filtering step is to drive a compensation filter to receive the prototype filter output signal and multiply the prototype filter output signal by a compensation transfer function to generate an output signal, and the compensation filter includes cascade-connected complex shaping filters and has a second group delay in the passband. The first group delay is added to the second group delay to generate a complex group delay, and the complex group delay is a constant.
藉此,本揭示內容的補償式複合濾波器之訊號處理方法根據原型轉移函數與補償轉移函數轉換輸入訊號為輸出訊號,進而可有效地抑制通帶外的干擾訊號,並以高保真度保留通帶中的訊號。Thus, the signal processing method of the compensated composite filter disclosed in the present invention converts the input signal into the output signal according to the prototype transfer function and the compensated transfer function, thereby effectively suppressing the interference signal outside the passband and retaining the signal in the passband with high fidelity.
前述實施方式之其他實施例如下:前述原型轉移函數與補償轉移函數相乘而產生一複合轉移函數。原型轉移函數表示為 ,補償轉移函數表示為 ,複合轉移函數表示為 ,且複合轉移函數符合下式: ;其中 為一變數, 為此些整形濾波器的一數量, 為此些整形濾波器之一者的一整形轉移函數, 為正整數且介於1至 之間。 Other embodiments of the above-mentioned implementation are as follows: The above-mentioned prototype transfer function and the compensation transfer function are multiplied to generate a composite transfer function. The prototype transfer function is represented as , the compensation transfer function is expressed as , the composite transfer function is expressed as , and the composite transfer function conforms to the following formula: ;in is a variable, is the number of these shaping filters, is an integer transfer function of one of these shaping filters, Is a positive integer between 1 and between.
前述實施方式之其他實施例如下:前述第一群延遲表示為 ,第二群延遲表示為 ,複合群延遲表示為 ,且複合群延遲符合下式: ;其中 為一角頻率, 為此些整形濾波器的一數量, 為此些整形濾波器之一者的一群延遲, 為正整數且介於1至 之間。 Other embodiments of the above-mentioned implementation are as follows: The above-mentioned first group delay is expressed as , the second group delay is expressed as , the composite group delay is expressed as , and the complex group delay meets the following formula: ;in is the angular frequency, is the number of these shaping filters, is the group delay of one of these shaping filters, Is a positive integer between 1 and between.
前述實施方式之其他實施例如下:前述此些整形濾波器之此者的群延遲符合下式: ;其中, 為此些整形濾波器之此者的一濾波器參數, 為此些整形濾波器之此者的另一濾波器參數;其中此些整形濾波器的複數個群延遲彼此相異。 Other embodiments of the aforementioned embodiment are as follows: The group delay of one of the aforementioned shaping filters conforms to the following formula: ;in, is a filter parameter for one of these shaping filters, is another filter parameter of one of these shaping filters; wherein a plurality of group delays of these shaping filters are different from each other.
前述實施方式之其他實施例如下:前述原型濾波器為一巴特沃斯濾波器(Butterworth filter)、一橢圓濾波器(Elliptic filter)、一無限脈衝響應濾波器(Infinite Impulse Response filter,IIR filter)或一切比雪夫濾波器(Chebyshev filter),且各整形濾波器為一全通濾波器(All-pass filter)。Other embodiments of the aforementioned implementation method are as follows: the aforementioned prototype filter is a Butterworth filter, an elliptic filter, an infinite impulse response filter (IIR filter) or a Chebyshev filter, and each shaping filter is an all-pass filter.
以下將參照圖式說明本揭示內容之複數個實施例。為明確說明起見,許多實務上的細節將在以下敘述中一併說明。然而,應瞭解到,這些實務上的細節不應用以限制本揭示內容。也就是說,在本揭示內容部分實施例中,這些實務上的細節是非必要的。此外,為簡化圖式起見,一些習知慣用的結構與元件在圖式中將以簡單示意的方式繪示之;並且重複之元件將可能使用相同的編號表示之。The following will describe multiple embodiments of the present disclosure with reference to the drawings. For the sake of clarity, many practical details will be described together in the following description. However, it should be understood that these practical details should not be used to limit the present disclosure. In other words, in some embodiments of the present disclosure, these practical details are not necessary. In addition, in order to simplify the drawings, some commonly known structures and components will be depicted in a simple schematic manner in the drawings; and repeated components may be represented by the same number.
此外,本文中當某一元件(或單元或模組等)「連接/連結」於另一元件,可指所述元件是直接連接/連結於另一元件,亦可指某一元件是間接連接/連結於另一元件,意即,有其他元件介於所述元件及另一元件之間。而當有明示某一元件是「直接連接/連結」於另一元件時,才表示沒有其他元件介於所述元件及另一元件之間。而第一、第二、第三等用語只是用來描述不同元件,而對元件本身並無限制,因此,第一元件亦可改稱為第二元件。且本文中之元件/單元/電路之組合非此領域中之一般周知、常規或習知之組合,不能以元件/單元/電路本身是否為習知,來判定其組合關係是否容易被技術領域中之通常知識者輕易完成。In addition, in this article, when a certain component (or unit or module, etc.) is "connected/linked" to another component, it may refer to that the component is directly connected/linked to another component, or it may refer to that a certain component is indirectly connected/linked to another component, that is, there are other components between the component and the other component. When it is clearly stated that a certain component is "directly connected/linked" to another component, it means that there are no other components between the component and the other component. The terms first, second, third, etc. are only used to describe different components, and there is no restriction on the components themselves. Therefore, the first component can also be renamed as the second component. Moreover, the combination of components/units/circuits in this article is not a generally known, conventional or familiar combination in this field. Whether the components/units/circuits themselves are known cannot be used to determine whether their combination relationship is easy to be completed by ordinary knowledgeable people in the technical field.
請參閱第1圖,第1圖係繪示依照本揭示內容的第一實施方式的補償式複合濾波器100的方塊示意圖。如第1圖所示,補償式複合濾波器100包含原型濾波器200以及補償濾波器300。原型濾波器200用以將輸入訊號S
in乘以一原型轉移函數而產生原型濾波輸出訊號S
pro,且原型濾波器200於一通帶(Passband)中具有一第一群延遲。補償濾波器300串聯連接原型濾波器200,且用以將原型濾波輸出訊號S
pro乘以一補償轉移函數而產生一輸出訊號S
out。補償濾波器300包含級聯連接的複數整形濾波器310、320、330,並於通帶中具有一第二群延遲。原型濾波器200的第一群延遲相加補償濾波器300的第二群延遲而產生一複合群延遲,且複合群延遲為一常數。
Please refer to FIG. 1, which is a block diagram of a compensated
詳細地說,當輸入訊號S
in輸入至補償式複合濾波器100時,補償式複合濾波器100根據原型濾波器200的原型轉移函數與補償濾波器300的補償轉移函數對輸入訊號S
in進行濾波,並會經歷一個恆定的複合群延遲後產生輸出訊號S
out,且輸出訊號S
out可完整地保留與輸入訊號S
in相同的波形。通帶以外的其他波段受到原型濾波器200抑制,且由於複合群延遲於通帶中為常數,因此補償式複合濾波器100的相位響應(Phase response)於通帶中可呈現線性。
Specifically, when the input signal S in is input to the compensated
原型轉移函數與補償轉移函數相乘而產生一複合轉移函數。原型轉移函數表示為 ,補償轉移函數表示為 ,複合轉移函數表示為 ,且複合轉移函數可符合下式子(1): (1)。 The prototype transfer function is multiplied by the compensation transfer function to generate a composite transfer function. The prototype transfer function is expressed as , the compensation transfer function is expressed as , the composite transfer function is expressed as , and the composite transfer function can satisfy the following formula (1): (1).
於式子(1)中,
為一變數(即Z轉換中的一複數),
為整形濾波器310、320、330的一數量,
為整形濾波器310、320、330之一者的一整形轉移函數,
為正整數且介於1至
之間。
In formula (1), is a variable (i.e. a complex number in the Z transform), is the number of
原型濾波器200的第一群延遲表示為
,補償濾波器300的第二群延遲表示為
,複合群延遲表示為
,且複合群延遲可符合下式子(2):
(2)。
The first group delay of the
於式子(2)中,
為一角頻率,
為整形濾波器310、320、330的數量,
為整形濾波器310、320、330之一者的一群延遲,
為正整數且介於1至
之間。由式子(1)與式子(2)可知,本揭示內容的補償式複合濾波器100透過將原型濾波器200串聯連接含有級聯連接的多個整形濾波器310、320、330的補償濾波器300,並利用補償濾波器300的第二群延遲補償原型濾波器200的第一群延遲,達到線性相位響應且同時具有低複雜度。以下段落將配合圖式以詳細說明本揭示內容的補償式複合濾波器100的運作機制與配置規範。
In formula (2), is the angular frequency, is the number of
請一併參閱第1圖、第2圖及第3圖,其中第2圖係繪示第1圖之補償式複合濾波器100的整形濾波器310、320、330的群延遲比較圖;以及第3圖係繪示第1圖之補償式複合濾波器100與原型濾波器200的群延遲比較圖。如第2圖所示,整形濾波器310、320、330的複數個群延遲彼此相異,即整形濾波器310、320、330各自的相位響應彼此相異。Please refer to FIG. 1, FIG. 2 and FIG. 3 together, wherein FIG. 2 is a group delay comparison diagram of the
整形濾波器310、320、330之一者的群延遲可符合下式子(3):
(3)。
The group delay of one of the
於式子(3)中,
為整形濾波器310、320、330之此者的一濾波器參數,
為整形濾波器310、320、330之此者的另一濾波器參數。整形濾波器310、320、330的濾波器參數
之預設值可分別為
、
及
,其中參數
係對應於群延遲曲線的峰值,參數
控制曲線的形狀。由於
,使得整形濾波器310、320、330趨於穩定,因此整形濾波器310、320、330在所有頻段中的群延遲可為正值。濾波器參數
的選擇能提供整形補償式複合濾波器100的群延遲之自由度。在其他實施例中,本揭示內容的補償濾波器亦可使用更多的整形濾波器(即整形濾波器的數量大於4),但勢必會增加補償式複合濾波器內部的乘法器數量及加法器數量。
In formula (3), is a filter parameter of one of the
具體而言,原型濾波器200可為一巴特沃斯濾波器(Butterworth filter)、一橢圓濾波器(Elliptic filter)、一無限脈衝響應濾波器(Infinite Impulse Response filter,IIR filter)或一切比雪夫濾波器(Chebyshev filter),且各整形濾波器310、320、330為一全通濾波器(All-pass filter)。具有原型轉移函數
的原型濾波器200需滿足振幅響應的配置規範。具有補償轉移函數
的補償濾波器300主要用以重構補償式複合濾波器100的相位響應。
Specifically, the
儘管習知的FIR濾波器在整個頻帶上可具有線性相位響應,但對於頻寬受限的訊號是不必要的,因此補償式複合濾波器100的相位響應僅在通帶PB(繪示於第3圖)中呈線性,這是因為其他波段已被原型濾波器200抑制。Although conventional FIR filters may have a linear phase response over the entire frequency band, this is not necessary for bandwidth-limited signals. Therefore, the phase response of the compensated
如第3圖所示,通帶PB介於0至通帶邊緣頻率
之間。用以將原型濾波器200的第一群延遲D
1補償至複合群延遲D
C的一群延遲裕度M
GD可由補償濾波器300的第二群延遲D
2所提供,且第二群延遲D
2可由整形濾波器310的群延遲D
21、整形濾波器320的群延遲D
22及整形濾波器330的群延遲D
23相互疊加所構成。補償式複合濾波器100在通帶PB中的複合群延遲D
C之一平均值表示為
,且可符合下式子(4):
(4)。
As shown in Figure 3, the passband PB is between 0 and the passband edge frequency. A group delay margin MGD for compensating the first group delay D1 of the
此外,通帶PB上的複合群延遲D C之一平坦度定義為複合轉移函數 的一均方根(Root Mean Square,RMS)群延遲誤差,並可表示為 且符合下式子(5): (5)。 In addition, the flatness of the complex group delay DC over the passband PB is defined as the complex transfer function The root mean square (RMS) group delay error can be expressed as And it meets the following formula (5): (5).
均方根群延遲誤差
可作為補償式複合濾波器100的配置規範。補償式複合濾波器100的配置規範可包含規範(i)通帶邊緣頻率
(rad./sample);規範(ii)阻帶邊緣頻率
(rad./sample);規範(iii)通帶漣波峰對峰值小於0.05 dB;規範(iv)阻帶抑制不小於73 dB;規範(v)通帶PB上的均方根群延遲誤差
小於0.2;規範(vi)補償式複合濾波器100的階數不大於20;及規範(vii)整形濾波器310、320、330的數量大於等於1且小於等於3,即
。
RMS Group Delay Error It can be used as the configuration specification of the compensated
請一併參閱第1-3、4及5圖,其中第4圖係繪示本揭示內容的原型濾波器200與習知的FIR濾波器的振幅響應比較圖;以及第5圖係繪示本揭示內容的原型濾波器200與習知的FIR濾波器的群延遲比較圖。規範(i)-(iv)主要作為原型濾波器200的配置規範。如第4圖所示,習知的FIR濾波器(即最小二乘法有限脈衝響應濾波器Least-Squares FIR filter,以下簡稱最小二乘法FIR濾波器)、最小多項式無限脈衝響應濾波器(Least
pth-norm IIR filter,以下簡稱最小多項式IIR濾波器)及切比雪夫II型濾波器均符合上述規範(i)-(iv),但本揭示內容為了降低系統複雜度,故選擇階數較少的最小多項式IIR濾波器(8階)與切比雪夫II型濾波器(12階)作為原型濾波器200。如第5圖所示,與切比雪夫II型濾波器相比,最小多項式IIR濾波器具有較大的一群延遲裕度,但其可由整形濾波器310、320、330填充。
Please refer to Figures 1-3, 4 and 5 together, wherein Figure 4 is a graph showing a comparison of the amplitude response of the
請一併參閱第1-5、6、7及8圖,其中第6圖係繪示本揭示內容的原型濾波器200的均方根群延遲誤差的曲線示意圖;第7圖係繪示本揭示內容的補償式複合濾波器100與習知的FIR濾波器的相位響應比較圖;以及第8圖係繪示第1圖之補償式複合濾波器100、原型濾波器200及整形濾波器310、320、330的群延遲的曲線示意圖。如第6圖所示,與最小多項式IIR濾波器相比,選擇切比雪夫II型濾波器作為原型濾波器200可達到最小的均方根群延遲誤差,這是因為切比雪夫II型濾波器所需補償的群延遲裕度小於最小多項式IIR濾波器的群延遲裕度,因此補償式複合濾波器100的優選配置如後。選擇切比雪夫II型濾波器作為原型濾波器200,整形濾波器310、320、330的數量為3(
),且整形濾波器310、320、330的濾波器參數
可分別為
、
及
。原型濾波器200的均方根群延遲誤差為0.1196。基於上述優選配置的參數,第7圖的補償式複合濾波器100的相位響應不僅可如同最小二乘法FIR濾波器的相位響應在通帶PB中呈現線性,還能具有低複雜度。第8圖的補償式複合濾波器100的複合群延遲D
C可為58。
Please refer to Figures 1-5, 6, 7 and 8 together, wherein Figure 6 is a curve diagram showing the root mean square group delay error of the
請接續參閱表一,表一係列示濾波器(i)-(iv)各自的複雜度(例如乘法器數量及加法器數量),其中濾波器(i)-(iv)分別為補償式複合濾波器100、最小二乘法FIR濾波器、等值漣波FIR濾波器以及採用濾波器設計工具(例如MATLAB)中的窗口設計方法(Window design method)的FIR濾波器。
表一
如表一所示,本揭示內容的補償式複合濾波器100相比其他三個FIR濾波器具有較少的乘法器數量與加法器數量,因此補償式複合濾波器100在通帶PB中不僅可實現線性相位響應,且同時具有低複雜度的電路結構。As shown in Table 1, the compensated
請一併參閱第1、9A、9B、10A、10B及10C圖,其中第9A圖係繪示本揭示內容的輸入訊號S
in的波形示意圖;第9B圖係繪示本揭示內容的輸入訊號S
in經過離散時間傅立葉轉換(Discrete-Time Fourier Transform,DTFT)後的頻譜示意圖;第10A圖係繪示第9A圖之輸入訊號S
in經由切比雪夫濾波器過濾後所產生的輸出訊號S
1的波形示意圖;第10B圖係繪示第9A圖之輸入訊號S
in經由補償式複合濾波器100過濾後所產生的輸出訊號S
out的波形示意圖;以及第10C圖係繪示第9A圖之輸入訊號S
in經由最小二乘法FIR濾波器過濾後所產生的輸出訊號S
2的波形示意圖。
Please refer to Figures 1, 9A, 9B, 10A, 10B and 10C together, wherein Figure 9A is a waveform diagram of the input signal S in of the present disclosure; Figure 9B is a spectrum diagram of the input signal S in of the present disclosure after discrete-time Fourier transform (DTFT); Figure 10A is a waveform diagram of the output signal S 1 generated after the input signal S in of Figure 9A is filtered by a Chebyshev filter; Figure 10B is a waveform diagram of the output signal S out generated after the input signal S in of Figure 9A is filtered by a compensated
如第10A圖所示,雜訊雖已從輸出訊號S
1過濾,但輸出訊號S
1的波形明顯地產生多處跳動點,並相異於輸入訊號S
in。如第10B圖與第10C圖所示,具有線性相位響應的補償式複合濾波器100與最小二乘法FIR濾波器均以高保真度保留輸入訊號S
in的波形,並同時過濾干擾訊號。由於最小二乘法FIR濾波器的階數為200並具有較大的群延遲,因此與最小二乘法FIR濾波器相比,本揭示內容的補償式複合濾波器100在低複雜度的配置下仍可具有較低的通帶內群延遲。
As shown in FIG. 10A , although the noise has been filtered from the output signal S 1 , the waveform of the output signal S 1 obviously produces multiple jump points and is different from the input signal S in . As shown in FIG. 10B and FIG. 10C , the compensated
請一併參閱第1圖與第11圖,其中第11圖係繪示第1圖之補償式複合濾波器100的零極點分佈圖。由於整形濾波器310、320、330的數量為3,因此整形濾波器310、320、330共有六個(2×3)個零極點對P
1、P
2、P
3、P
4、P
5、P
6。如第11圖所示,零極點對P
1、P
2、P
3、P
4、P
5、P
6散佈在頻帶
上,且其對應的六個極點均位於單位圓內,因此補償式複合濾波器100具有高穩定性。
Please refer to FIG. 1 and FIG. 11 together, wherein FIG. 11 is a diagram showing the zero-pole distribution of the compensated
請一併參閱第1-8及12圖,其中第12圖係繪示依照本揭示內容的第二實施方式的補償式複合濾波器之訊號處理方法400的流程示意圖。如第12圖所示,補償式複合濾波器之訊號處理方法400可應用於第一實施方式的補償式複合濾波器100,且包含以下步驟:第一濾波步驟S02以及第二濾波步驟S04。Please refer to Figures 1-8 and 12 together, wherein Figure 12 is a flow chart showing a
第一濾波步驟S02係驅動原型濾波器200接收輸入訊號S
in,並將輸入訊號S
in乘以原型轉移函數而產生原型濾波輸出訊號S
pro,且原型濾波器200於通帶PB中具有第一群延遲D
1。第二濾波步驟S04係驅動補償濾波器300接收原型濾波輸出訊號S
pro,並將原型濾波輸出訊號S
pro乘以補償轉移函數而產生輸出訊號S
out。補償濾波器300包含級聯連接的整形濾波器310、320、330,並於通帶PB中具有第二群延遲D
2。於第二濾波步驟S04中,第二群延遲D
2疊加第一群延遲D
1以形成複合群延遲D
C,且複合群延遲D
C於通帶PB中為常數。因此,本揭示內容的補償式複合濾波器之訊號處理方法400根據原型轉移函數與補償轉移函數轉換輸入訊號S
in為輸出訊號S
out,達到有效地抑制通帶PB外的干擾訊號,且以高保真度保留通帶PB中的訊號,藉以令輸出訊號S
out的波形與輸入訊號S
in相同。
The first filtering step S02 is to drive the
綜上所述,本揭示內容提供的補償式複合濾波器100及補償式複合濾波器之訊號處理方法400具有下列優點:其一,於數位訊號處理器中,補償式複合濾波器可替換習知的FIR濾波器,進而降低系統複雜度。其二,在低複雜度的電路結構下,補償式複合濾波器仍可實現線性相位響應。其三,補償式複合濾波器具有高穩定性。In summary, the compensated
雖然本發明已以實施方式揭露如上,然其並非用以限定本發明,任何熟習此技藝者,在不脫離本發明之精神和範圍內,當可作各種之更動與潤飾,因此本發明之保護範圍當視後附之申請專利範圍所界定者為準。Although the present invention has been disclosed in the above embodiments, it is not intended to limit the present invention. Anyone skilled in the art can make various changes and modifications without departing from the spirit and scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the scope defined in the attached patent application.
100:補償式複合濾波器
200:原型濾波器
300:補償濾波器
310,320,330:整形濾波器
400:補償式複合濾波器之訊號處理方法
D
1:第一群延遲
D
2:第二群延遲
D
21,D
22,D
23:群延遲
D
C:複合群延遲
M
GD:群延遲裕度
PB:通帶
S02:第一濾波步驟
S04:第二濾波步驟
S
in:輸入訊號
S
pro:原型濾波輸出訊號
S
out,S
1,S
2:輸出訊號
:通帶邊緣頻率
100: Compensated composite filter 200: Prototype filter 300:
第1圖係繪示依照本揭示內容的第一實施方式的補償式複合濾波器的方塊示意圖; 第2圖係繪示第1圖之補償式複合濾波器的整形濾波器的群延遲比較圖; 第3圖係繪示第1圖之補償式複合濾波器與原型濾波器的群延遲比較圖; 第4圖係繪示本揭示內容的原型濾波器與習知的有限脈衝響應濾波器的振幅響應比較圖; 第5圖係繪示本揭示內容的原型濾波器與習知的有限脈衝響應濾波器的群延遲比較圖; 第6圖係繪示本揭示內容的原型濾波器的均方根群延遲誤差的曲線示意圖; 第7圖係繪示本揭示內容的補償式複合濾波器與習知的有限脈衝響應濾波器的相位響應比較圖; 第8圖係繪示第1圖之補償式複合濾波器、原型濾波器及整形濾波器的群延遲的曲線示意圖; 第9A圖係繪示本揭示內容的輸入訊號的波形示意圖; 第9B圖係繪示本揭示內容的輸入訊號經過離散時間傅立葉轉換後的頻譜示意圖; 第10A圖係繪示第9A圖之輸入訊號經由切比雪夫濾波器過濾後所產生的輸出訊號的波形示意圖; 第10B圖係繪示第9A圖之輸入訊號經由補償式複合濾波器過濾後所產生的輸出訊號的波形示意圖; 第10C圖係繪示第9A圖之輸入訊號經由最小二乘法有限脈衝響應濾波器過濾後所產生的輸出訊號的波形示意圖; 第11圖係繪示第1圖之補償式複合濾波器的零極點分佈圖;以及 第12圖係繪示依照本揭示內容的第二實施方式的補償式複合濾波器之訊號處理方法的流程示意圖。 FIG. 1 is a block diagram of a compensated composite filter according to the first embodiment of the present disclosure; FIG. 2 is a group delay comparison diagram of a shaping filter of the compensated composite filter of FIG. 1; FIG. 3 is a group delay comparison diagram of the compensated composite filter of FIG. 1 and a prototype filter; FIG. 4 is an amplitude response comparison diagram of the prototype filter of the present disclosure and a known finite pulse impulse response filter; FIG. 5 is a group delay comparison diagram of the prototype filter of the present disclosure and a known finite pulse impulse response filter; Figure 6 is a curve diagram showing the root mean square group delay error of the prototype filter of the present disclosure; Figure 7 is a phase response comparison diagram showing the compensated composite filter of the present disclosure and the known finite pulse impulse response filter; Figure 8 is a curve diagram showing the group delay of the compensated composite filter, the prototype filter and the shaping filter of Figure 1; Figure 9A is a waveform diagram showing the input signal of the present disclosure; Figure 9B is a spectrum diagram showing the input signal of the present disclosure after discrete time Fourier transformation; FIG. 10A is a waveform diagram of the output signal generated after the input signal of FIG. 9A is filtered by the Chebyshev filter; FIG. 10B is a waveform diagram of the output signal generated after the input signal of FIG. 9A is filtered by the compensated composite filter; FIG. 10C is a waveform diagram of the output signal generated after the input signal of FIG. 9A is filtered by the least squares finite pulse impulse response filter; FIG. 11 is a zero-pole distribution diagram of the compensated composite filter of FIG. 1; and FIG. 12 is a flow diagram of the signal processing method of the compensated composite filter according to the second embodiment of the present disclosure.
100:補償式複合濾波器 100: Compensatory composite filter
200:原型濾波器 200: Prototype filter
300:補償濾波器 300: Compensation filter
310,320,330:整形濾波器 310,320,330: Shaping filter
Sin:輸入訊號 S in : input signal
Spro:原型濾波輸出訊號 S pro : prototype filter output signal
Sout:輸出訊號 S out : output signal
Claims (8)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW112124987A TWI876412B (en) | 2023-07-04 | 2023-07-04 | Compensated composite filter and signal processing method thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW112124987A TWI876412B (en) | 2023-07-04 | 2023-07-04 | Compensated composite filter and signal processing method thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TW202504247A TW202504247A (en) | 2025-01-16 |
| TWI876412B true TWI876412B (en) | 2025-03-11 |
Family
ID=95152640
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW112124987A TWI876412B (en) | 2023-07-04 | 2023-07-04 | Compensated composite filter and signal processing method thereof |
Country Status (1)
| Country | Link |
|---|---|
| TW (1) | TWI876412B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4310803A (en) * | 1978-08-23 | 1982-01-12 | Hiroshi Kurihara | Bandpass filter circuit |
| US4395688A (en) * | 1981-08-11 | 1983-07-26 | Harris Corporation | Linear phase filter with self-equalized group delay |
| US5627899A (en) * | 1990-12-11 | 1997-05-06 | Craven; Peter G. | Compensating filters |
| US20050119772A1 (en) * | 2003-11-28 | 2005-06-02 | Yuji Yamada | Digital signal processing apparatus, method thereof and headphone apparatus |
| US9197462B2 (en) * | 2008-10-28 | 2015-11-24 | Freescale Semiconductor Inc. | Single amplifier filter for constant group delay in radio frequency transmitters |
-
2023
- 2023-07-04 TW TW112124987A patent/TWI876412B/en active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4310803A (en) * | 1978-08-23 | 1982-01-12 | Hiroshi Kurihara | Bandpass filter circuit |
| US4395688A (en) * | 1981-08-11 | 1983-07-26 | Harris Corporation | Linear phase filter with self-equalized group delay |
| US5627899A (en) * | 1990-12-11 | 1997-05-06 | Craven; Peter G. | Compensating filters |
| US20050119772A1 (en) * | 2003-11-28 | 2005-06-02 | Yuji Yamada | Digital signal processing apparatus, method thereof and headphone apparatus |
| US9197462B2 (en) * | 2008-10-28 | 2015-11-24 | Freescale Semiconductor Inc. | Single amplifier filter for constant group delay in radio frequency transmitters |
Also Published As
| Publication number | Publication date |
|---|---|
| TW202504247A (en) | 2025-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Podder et al. | Design and implementation of Butterworth, Chebyshev-I and elliptic filter for speech signal analysis | |
| CN107563032A (en) | Computerization wave filter design system | |
| Darak et al. | Efficient implementation of reconfigurable warped digital filters with variable low-pass, high-pass, bandpass, and bandstop responses | |
| JP3223188B2 (en) | How to reduce digital filter bank power consumption by reducing the number of products. | |
| CN112737544A (en) | Surface acoustic wave filter | |
| TWI876412B (en) | Compensated composite filter and signal processing method thereof | |
| Getu | Digital IIR filter design using bilinear transformation in MATLAB | |
| Kaur et al. | Digital filteration of ECG signals for removal of baseline drift | |
| WO2025256334A1 (en) | Filter and multiplexer | |
| WO2004008637A1 (en) | Digital filter designing method, digital filter designing program, digital filter | |
| Roy et al. | A new design strategy of sharp cut-off fir filter with powers-of-two coefficients | |
| Sutthikarn et al. | Tunable bandpass/bandstop digital filters basedon 1 st-order allpass network instead of unit delay | |
| Lehto et al. | Synthesis of wideband linear-phase FIR filters with a piecewise-polynomial-sinusoidal impulse response | |
| Barcellos et al. | Design of FIR filters combining the frequency-response masking and the WLS-Chebyshev approaches | |
| WO2023148540A1 (en) | A recursive fir digital filter | |
| KR20060119891A (en) | Synthesis filter and passband ripple cancellation method | |
| TW202504248A (en) | Cascading digital filter and signal processing method thereof | |
| JP4214391B2 (en) | How to design a digital filter | |
| Jovanovic-Dolecek et al. | Multiplier-free FIR filter design based on IFIR structure and rounding | |
| JPH11298202A (en) | Band pass filter and double multiplication circuit | |
| JPH1051269A (en) | Low-pass filter | |
| Jovanovic-Dolecek et al. | Design of FIR lowpass filters using stepped triangular approximation | |
| Holland | Signal Conditioning with Almost No Multiplications via Filter Sharpening | |
| Webb et al. | A new approach to designing computationally efficient interpolated FIR filters | |
| CN110518894B (en) | High-speed low-complexity second-order all-pass filter |