TWI415428B

TWI415428B - Channel estimation method and device

Info

Publication number: TWI415428B
Application number: TW098125837A
Authority: TW
Inventors: Yi Ying Liao; Tai Lai Tung
Original assignee: Mstar Semiconductor Inc
Priority date: 2009-07-31
Filing date: 2009-07-31
Publication date: 2013-11-11
Also published as: TW201105069A; US20110026620A1; US8615048B2

Description

通道估算方法與裝置

本發明係有關於通道估算之技術，特別是用於正交分頻多工(OFDM)通訊系統之通道估算方法與裝置。

在無線通訊系統(wireless communication system)中，由於其無線電通道(radio channel)通常有多重路徑衰減(multipath fading)效應，以致於在接收的信號中會有符號間干擾(ISI)之問題。為了消除符號間干擾，一般的作法乃在接收器中設置等化器(equalizer)，而等化器之運作需用到通道脈衝響應(CIR)之資訊，使得通道脈衝響應之估算在行動無線電系統中扮演相當關鍵的角色。

又，OFDM在無線通訊領域中已是一重要通信技術，主要可增加資料傳輸率。例如IEEE 802.11a採用OFDM技術，其資料傳輸率可達54Mbps；而IEEE 802.11b未用OFDM技術，其資料傳輸率僅有11Mbps。故如何有效率地估算OFDM系統之通道脈衝響應以消除符號間干擾，以發揮OFDM高傳輸率之特性，為一重要課題。在OFDM系統中，通道估算，亦即通道脈衝響應估算一般係經由使用傳送器和接收器事前已知的引導符號(pilot symbol)獲得。

請參照圖1，其繪示一習知OFDM通道估算裝置之方塊圖。如圖1所示，該習知OFDM通道估算裝置包含一4096取樣點IFFT運算單元101、一鏡像信號拒斥運算單元102、一4096取樣點FFT運算單元103及一頻域通道響應之平滑化處理單元104。

該4096取樣點IFFT運算單元101係用以對初估頻域通道響應執行4096取樣點反快速傅立葉轉換(IFFT)運算，以產生一時域通道脈衝響應，其中載波之頻率與時間分佈圖請參照圖2，共有4096個時域通道脈衝響應計算值：。該習知OFDM通道共有4096個副載波，其中每8個副載波有一個為引導副載波，用以承載一引導符號，其餘副載波，稱為資料副載波，用以承載資料符號；亦即在該等4096個副載波中，有512個引導副載波用以承載引導符號，有3584個資料副載波用以承載資料符號。在圖3之時間-頻率平面中，繪示了該習知OFDM通道所傳輸之17×4096個符號。該初估頻域通道響應之估算，一般乃在各引導副載波處對其引導符號之頻域傳送值與頻域接收值施行一最小平方差演算而得。亦即該初估頻域通道響應只在頻率位置k對應於引導副載波處有演算值，其他對應於資料副載波處之頻域通道響應值則令其為0。該等512個引導符號在該等4096個副載波中之分佈型式分為偶數型及奇數型二種類型，偶數型為：●○○○○○○○●○○○○○○○●×●×●×●×奇數型為：○○○○●○○○○○○○●○○○×●×●×●×●其中●代表引導符號，其對應之引導副載波頻率位置在偶數型為0，8，16，24，‧‧‧‧‧，4088，在奇數型為4，12，20，28，‧‧‧‧‧，4092；○代表資料符號；×代表重複左側之排列。故該初估頻域通道響應依該等512個引導符號在該等4096個副載波中之分佈型式亦分為偶數型頻域通道響應及奇數型頻域通道響應二種類型。該初估頻域通道響應以偶數型頻域通道響應、奇數型頻域通道響應交錯之方式傳送給該4096取樣點IFFT運算單元101。

該鏡像信號拒斥運算單元102係用以保留該4096個時域通道脈衝響應之前、後各256個時域通道脈衝響應計算值，而將屬於鏡像信號之其他3584個時域通道脈衝響應計算值濾除以產生一時域通道脈衝響應，亦即

該4096取樣點FFT運算單元103係用以對該時域通道脈衝響應執行4096取樣點快速傅立葉轉換(FFT)運算，以產生一頻域通道響應,k=0~4095。

該頻域通道響應之平滑化處理單元104係用以依該頻域通道響應之17組資料：執行一算術平均運算，以產生一頻域通道響應，其中

在該習知OFDM通道估算裝置中，該4096取樣點反快速傅立葉轉換(IFFT)運算之取樣值實際上只會在512個取樣點處可能有非零值，其餘皆為零，故其4096之IFFT計算規模顯得過大而缺乏效率，應有改進之空間。因此，本發明乃根據中國移動多媒體廣播(CMMB)OFDM之系統特性，研究降低其IFFT計算規模之方法，以期在一低計算規模下保有相同之通道估算效能，從而提供一低成本之解決方案。

有鑒於此瓶頸，本發明提出一新穎的演算法以完成通道估算，該演算法可用較小規模的IFFT計算、依初估頻域通道響應取樣型式適應性地調整計算公式，而大幅縮小整體計算規模。

本發明之一目的在於提供一低成本、低功耗之通道估算方法，用以適應性地估算通道之脈衝響應，以供接收器執行補償功能。

本發明之另一目的在於提供一低成本、低功耗之通道估算裝置，用以適應性地估算通道之脈衝響應，以供接收器執行補償功能。

本發明又一目的在於提供一具低計算複雜度之通道估算方案，其可依初估頻域通道響應之取樣型式適應性地調整計算公式，以用較小規模之反快速傅立葉轉換(IFFT)計算實現通道脈衝響應之估算。

為達成本發明上述諸目的，一通道估算方法乃被提出，該方法包含對第一數量之初估頻域通道響應值執行具第二數量取樣點之反快速傅立葉轉換運算，以產生第一時域通道脈衝響應，其中該第二數量小於該第一數量；對該第一時域通道脈衝響應執行一時域窗口過濾處理，以產生第二時域通道脈衝響應；對複數組不同時域之該第二時域通道脈衝響應執行一平滑化運算，以產生一平滑時域通道脈衝響應；以及對該平滑時域通道脈衝響應執行具該第一數量取樣點之快速傅立葉轉換運算，以產生一頻域通道響應。

為達成上述諸目的，本發明進一步提出一通道估算裝置，包含一反快速傅立葉轉換運算單元，用以對具第一數量響應值之初估頻域通道響應執行具第二數量取樣點之反快速傅立葉轉換運算，以產生第一時域通道脈衝響應，其中該第二數量小於該第一數量；一窗口過濾處理單元，用以對該第一時域通道脈衝響應執行一時域窗口過濾處理，以產生第二時域通道脈衝響應；一平滑化處理單元，用以對複數組不同時域之該第二時域通道脈衝響應執行一平滑化運算，以產生一平滑時域通道脈衝響應；以及一快速傅立葉轉換運算單元，用以對該平滑時域通道脈衝響應執行具該第一數量取樣點之快速傅立葉轉換運算，以產生一頻域通道響應。

為使　貴審查委員能進一步瞭解本發明之結構、特徵及其目的，茲附以圖式及較佳具體實施例之詳細說明如后。

請參照圖4，本案對中國移動多媒體廣播(CMMB)OFDM系統進行研究後發現，其在4096個副載波中設有保護頻帶(guard band)且位於該保護頻帶內之副載波均不使用，以避免訊號間干擾，故其真正會用到的引導副載波只有192×2個，且其係分佈於不連續之二區段中。此外，CMMB OFDM系統之引導符號分佈可區分成在時域軸上交替出現之偶數型及奇數型兩類。因此，可將初估頻域通道響應之4096取樣點IFFT運算簡化如下：偶數型計算公式：

由於兩組192個引導副載波分別分佈於及處，故只需對及計算。奇數型計算公式：

本發明512個取樣點之IFFT運算，只需針對及等512點實施，其餘3584個點則無須進行IFFT運算，其原因在於非及非間之時域通道脈衝響應為鏡像信號，而該鏡像信號會在進行IFFT運算之後被濾除。

依照以上之公式推導，本案提出一種通道估算方法，其可大幅縮小計算規模及電路成本。請參照圖5，其繪示本發明通道估算一實施例之流程圖。如圖5所示，該實施例包含以下步驟：對第一數值個響應值之初估頻域通道響應執行第二數值取樣點之IFFT運算，以產生第一時域通道脈衝響應(步驟a)；對該第一時域通道脈衝響應執行一時域窗口(time-domain windowing)過濾處理，以產生第二時域通道脈衝響應(步驟b)；對複數組不同時域之該第二時域通道脈衝響應執行一平滑化運算，以產生一平滑時域通道脈衝響應(步驟c)；以及對該平滑時域通道脈衝響應執行第一數值取樣點之FFT運算，以產生一頻域通道響應(步驟d)。以下為各步驟之詳細說明。

在步驟a，對第一數值個響應值之初估頻域通道響應執行第二數值取樣點之IFFT運算，以產生第一時域通道脈衝響應。由上述公式可知，使用512取樣點之IFFT即可達成與使用4096點IFFT相同之效果，故當該第一數值為4096，第二數值可為512。且該512取樣點之IFFT運算包含該偶數型計算公式及該奇數型計算公式，其中該IFFT運算之偶數型計算公式在k=8m+2及k=8m+2559，其中m=0~191等384個取樣點處有非零之取樣值，在k之其餘取樣點處之取樣值則設為零；該IFFT運算之奇數型計算公式在k=8m+6及k=8m+2563，其中m=0~191等384個取樣點處有非零之取樣值，在k之其餘取樣點處之取樣值則設為零。其中該第二數值取樣點之IFFT運算包含：，其所需之512個取樣值分別為：

(奇數型)，亦即取樣值只在m介於0~191間有非零值，在192至511間則為零。本發明512個取樣點之IFFT運算，只需針對及等512點實施，其餘3584個點則無須進行IFFT運算，其原因在於非及非間之時域通道脈衝響應為鏡像信號，而該鏡像信號會在進行IFFT運算之後被濾除。

在步驟b，對該第一時域通道脈衝響應執行一時域窗口(time-domain windowing)過濾處理，以產生第二時域通道脈衝響應。該時域窗口過濾處理係用以濾除雜訊。請參照圖6，該窗口過濾處理根據一臨界數值及一窗口長度，過濾該第一時域通道脈衝響應，以濾除該窗口長度以外及低於該臨界數值之時域通道脈衝響應，而產生該第二時域通道脈衝響應。其中該窗口長度係為回音延遲(echo delay)之一函數，與信號之多重路徑有關；而該臨界數值係為窗口內各時域通道脈衝響應之加權平均值。

在步驟c，對複數組不同時域之該第二時域通道脈衝響應執行一平滑化運算，以產生一平滑時域通道脈衝響應。在本實施例中，使用17組時域連續之第二時域通道脈衝響應，其包括一待估測時間之通道脈衝響應及其時間軸前後各8組之通道脈衝響應，而該平滑化運算可為一平均運算，或一加權運算。亦即步驟c係對該第二時域通道脈衝響應之17組資料進行該平滑化運算，以分別計算出及等512點之響應值，而等3584點的響應值則補零，從而產生4096點平滑時域通道脈衝響應。更詳細地說，當本發明採用平均運算時，其運算如下：

當本發明採用加權運算時，則會考慮時間軸上通道變化，以決定加權數。當通道在時間軸上變化劇烈時，本發明會調高在待估測時間之通道脈衝響應之權重，而降低在待估測時間前、後時間之通道脈衝響應之權重，使離待測時間越遠之通道脈衝響應的權重降越低，以依此加權平均而產生在及等512點有響應值，而在等3584點之響應值則為零之該平滑時域通道脈衝響應。

在步驟d中，對該平滑時域通道脈衝響應執行第一數值取樣點之FFT運算，以產生一頻域通道響應。其中該第一數值為4096，與初估頻域通道響應具有相同之頻域通道響應數目，亦即步驟d在k=0~4095之所有副載波頻率處皆產生頻域通道響應值。

與習知相比，本案使用512取樣點之IFFT即可達成與使用4096點IFFT相同之效果，且其可由該偶數型計算公式及該奇數型計算公式得到證明，因此本案可使該IFFT運算、該時域窗口過濾處理及該平滑化運算所需之計算量，大為降低。

請參照圖7，其繪示本發明通道估算裝置一實施例之方塊圖。如圖7所示，該實施例包含一512取樣點IFFT運算單元701、一窗口過濾處理單元702、一時域通道脈衝響應之平滑化處理單元703及一4096取樣點FFT運算單元704。

該512取樣點IFFT運算單元701係用以對第一數值個響應值之初估頻域響應，依該偶數型公式或該奇數型公式執行512取樣點IFFT運算，以產生第一時域通道脈衝響應。其中該第一數值為4096，該IFFT運算之偶數型計算公式只在k=8m+2及k=8m+2559,m=0~191等384個取樣點處有非零之取樣值(引導符號所在處)，在k之其餘取樣點處其取樣值則為零；該IFFT運算之奇數型計算公式只在k=8m+6及k=8m+2563,m=0~191等384個取樣點處可能有非零之取樣值，在k之其餘取樣點處其取樣值則為零。其中該具第二數值取樣點之IFFT運算包含：，其所需之512個取樣值分別為：

(奇數型)，亦即取樣值只在m介於0~191間有非零值，在192至511間則為零。本發明512個取樣點之IFFT運算，只需針對及等512點實施，其餘3584個點無須進行IFFT運算，因為在非及非間之時域通道脈衝響應為鏡像信號，而該鏡像信號會在進行IFFT運算之後被濾除。

該窗口過濾處理單元702係用以濾除雜訊。該窗口過濾處理單元702根據一臨界數值及一窗口長度，過濾該第一時域通道脈衝響應，以濾除窗口長度以外及低於該臨界數值之時域通道脈衝響應，而產生該第二時域通道脈衝響應。其中該窗口長度為回音延遲(echo delay)之一函數，與信號之多重路徑有關；而該臨界數值為窗口內各時域通道脈衝響應值之加權平均值。

該時域通道脈衝響應之平滑化處理單元703係用以對複數組不同時域之該第二時域通道脈衝響應執行一平滑化運算，以產生一平滑時域通道脈衝響應。在本實施例中，使用17組時域連續之第二時域通道脈衝響應，其係對應於待估測時間之通道脈衝響應及其時間軸前後各8組之通道脈衝響應，而該平滑化運算可為一平均運算或一加權運算，以在及等512點分別產生平滑響應值，而其餘等3584點的響應值補零，從而產生4096點平滑時域通道脈衝響應。更詳細地說，當採用平均運算時，其運算如下：

當採用加權運算時，本發明考慮時間軸上通道變化，以決定加權數。當通道在時間軸上變化劇烈時，待估測時間之通道脈衝響應之權重會調高，待估測時間前、後時間之通道脈衝響應之權重會降低，離待測時間越遠之通道脈衝響應的權重會越低，以依此加權平均而在及等512點產生平滑響應值，而其餘等3584點響應值則補零，以產生平滑時域通道脈衝響應。

該4096取樣點FFT運算單元704係用以對該平滑時域通道脈衝響應執行4096取樣點FFT運算，以產生一頻域通道響應。其中該頻域通道響應在k=0~4095之所有副載波頻率處皆有頻域通道響應值。

請參照圖8，其繪示圖7該512取樣點IFFT運算單元701一實施例之方塊圖。如圖8所示，該實施例包含二192取樣點緩衝器801及802、二512取樣點IFFT計算器803及804、二乘法器805及806及一加法器807。

該192取樣點緩衝器801及802係用以緩衝該初估頻域通道響應。當該初估頻域通道響應之取樣點分佈為偶數型時，該192取樣點緩衝器801及802分別用以緩衝二不連續區域之初估頻域通道響應；當該初估頻域通道響應之取樣點分佈為奇數型時，該192取樣點緩衝器801及802則分別用以緩衝二不連續區域之初估頻域通道響應。

該512取樣點IFFT計算器803係耦接至該192取樣點緩衝器801，用以執行512取樣點之第一IFFT計算以產生一第一IFFT數值，其中或。

該512取樣點IFFT計算器804係耦接至該192取樣點緩衝器802，用以執行512取樣點之第二IFFT計算以產生一第二IFFT數值，其中或。

該乘法器805係用以依一第一相位變數e ^j ² ^π ⁿ ⁽² ⁺ ⁴ ⁱ ⁾ ^/ ⁴ ⁰ ⁹ ⁶ 及該第一IFFT數值之乘積，以產生一第一部分時域通道脈衝響應。其中該第一相位變數e ^j ² ^π ⁿ ⁽²⁺ ⁴ ⁱ ⁾ ^/ ⁴⁰⁹⁶ 分為i=0,1兩種情況；i=0代表該初估頻域通道響應其取樣點分佈屬於該偶數型，i=1代表該初估頻域通道響應其取樣點分佈屬於該奇數型。

該乘法器806係用以依一第二相位變數e ^j ^2π ⁿ ⁽ ²⁵⁵⁹⁺⁴ ⁱ ⁾ ^/4096 及該第二IFFT數值之乘積，產生一第二部分時域通道脈衝響應。其中該第二相位變數e ^j ^2π ⁿ ⁽ ²⁵⁵⁹⁺⁴ ⁱ ⁾ ^/4096 分為i=0,1兩種情況；i=0代表該初估頻域通道響應其取樣點分佈屬於該偶數型，i=1代表該初估頻域通道響應其取樣點分佈屬於該奇數型。

該加法器807係用以結合該第一部分時域通道脈衝響應及該第二部分時域通道脈衝響應以產生該第一時域通道脈衝響應。

請參照圖9，其繪示圖7該512取樣點IFFT運算單元701另一實施例之方塊圖。如圖9所示，該實施例包含一192取樣點緩衝器901、一512取樣點IFFT計算器902、一乘法器903、一選擇開關904、一時域通道脈衝響應緩衝器905及一加法器906。

該192取樣點緩衝器901係用以緩衝該初估頻域通道響應。當該初估頻域通道響應之取樣點分佈為偶數型時，該192取樣點緩衝器901乃用以依序緩衝及；當該初估頻域通道響應之取樣點分佈為奇數型時，該192取樣點緩衝器901則用以依序緩衝及。本案只需在下一批192個初估頻域通道響應值送至該192取樣點緩衝器901前完成IFFT計算，而使得該192取樣點緩衝器901可從容緩衝每一批192個初估頻域通道響應值。

該512取樣點IFFT計算器902係耦接至該192取樣點緩衝器901，用以執行512取樣點之IFFT計算以產生一IFFT數值，其中或。

該乘法器903係用以依一相位變數e ^j ^2π ⁿ ⁽ ²⁺⁴ ⁱ ⁾ ^/4096 (或e ^j ² ^π ⁿ ⁽ ²⁵⁵⁹⁺⁴ ⁱ ⁾ ^/4096 )及該IFFT數值之乘積，產生一第一部分時域通道脈衝響應(或一第二部分時域通道脈衝響應)。

該選擇開關904係用以自一第一相位變數及一第二相位變數中擇一作為該相位變數。其中該第一相位變數為e ^j ² ^π ⁿ ⁽ ² ⁺⁴ ⁱ ⁾ ^/4096 ,i=0,1；i=0代表該初估頻域通道響應其取樣點分佈屬於該偶數型，i=1代表該初估頻域通道響應其取樣點分佈屬於該奇數型。其中該第二相位變數為e ^j ² ^π ⁿ ⁽ ²⁵⁵⁹⁺⁴ ⁱ ⁾ ^/ ⁴⁰⁹⁶ ,i=0,1；i=0代表該初估頻域通道響應其取樣點分佈屬於該偶數型，i=1代表該初估頻域通道響應其取樣點分佈屬於該奇數型。

該時域通道脈衝響應緩衝器905係耦接至該乘法器903，用以緩衝一第一部份時域通道脈衝響應及一第二部份時域通道脈衝響應。

該加法器906係用以結合該第一部份時域通道脈衝響應及該第二部份時域通道脈衝響應，以產生該第一時域通道脈衝響應。

請參照圖10，其繪示圖7該時域通道脈衝響應之平滑化處理單元703一實施例之方塊圖。如圖10所示，該實施例包含一時域通道脈衝響應儲存單元1001及一平滑計算器1002。

該時域通道脈衝響應儲存單元1001係用以儲存該第二時域通道脈衝響應之17組資料，包含待估測時間之通道脈衝響應及其時間軸前後各8組之通道脈衝響應，由於本案僅使用512點IFFT，因此，該時域通道脈衝響應儲存單元1001只需17*512點的儲存空間，可大幅減少儲存單元的需求。

該平滑計算器1002係耦接至該時域通道脈衝響應儲存單元1001，用以對該第二時域通道脈衝響應之17組資料執行平滑化運算。該平滑化運算可採一平均運算，亦可為一加權運算，用以對該第二時域通道脈衝響應之17組資料進行運算，以分別產生及等512點之響應值，而其餘等3584點的響應值則補零，從而產生4096點平滑時域通道脈衝響應。更詳細地說，當採用平均運算時，執行以下運算：

以產生該平滑時域通道脈衝響應。當採用加權運算時，考慮時間軸上通道變化，以決定加權數。當通道在時間軸上變化劇烈時，待估測時間之通道脈衝響應之權重會調高，前後時間之通道脈衝響應之權重會降低，離待測時間通道越遠之通道脈衝響應的權重會越低，以在及等512點產生平滑響應值，其餘等3584點之響應值則補零，從而產生平滑時域通道脈衝響應。

是故，經由上述本案較佳實施例之實施，即可提供一具低計算複雜度、低成本、低功耗之通道估算方案，從而改進習知技術之缺點。

本案所揭示者，乃較佳實施例，舉凡局部之變更或修飾而源於本案之技術思想而為熟習該項技藝之人所易於推知者，例如引導副載波之排列方式、取樣點數之變更、緩衝器尺寸之變更、平均組數之變更、平滑化處理方式之變更等，俱不脫本案之專利權範疇。

綜上所陳，本案無論就目的、手段與功效，在在顯示其迥異於習知之技術特徵，且其首先發明合於實用，亦在在符合發明之專利要件，懇請　貴審查委員明察，並祈早日賜予專利，俾嘉惠社會，實感德便。

101．．．4096取樣點IFFT運算單元

102．．．鏡像信號拒斥運算單元

103、704．．．4096取樣點FFT運算單元

104．．．頻域通道響應之平滑化處理單元

701．．．512取樣點IFFT運算單元

702．．．窗口過濾處理單元

703．．．時域通道脈衝響應之平滑化處理單元

801、802、901．．．192取樣點緩衝器

803、804、902．．．512取樣點IFFT計算器

805、806、903．．．乘法器

807、906．．．加法器

904．．．選擇開關

905．．．時域通道脈衝響應緩衝器

1001．．．時域通道脈衝響應儲存單元

1002．．．平滑計算器

圖1為一示意圖，其繪示一習知OFDM通道估算裝置之方塊圖。

圖2為一示意圖，其繪示圖1之IFFT運算結果時域通道脈衝響應圖。

圖3為一示意圖，其繪示習知技術之CMMB OFDM初估頻域通道響應於時間-頻率平面之引導取樣分佈圖。

圖4為一示意圖，其繪示本發明之CMMB OFDM初估頻域通道響應於時間-頻率平面之引導取樣分佈圖。

圖5為一示意圖，其繪示本發明通道估算一實施例之流程圖。

圖6為一示意圖，其繪示本發明之時域窗口過濾處理過程。

圖7為一示意圖，其繪示本發明通道估算裝置一實施例之方塊圖。

圖8為一示意圖，其繪示圖7其512取樣點IFFT運算單元一實施例之方塊圖。

圖9為一示意圖，其繪示圖7其512取樣點IFFT運算單元另一實施例之方塊圖。

圖10為一示意圖，其繪示圖7其時域通道脈衝響應之平滑化處理單元一實施例之方塊圖。

704．．．4096取樣點FFT運算單元

701．．．512取樣點IFFT運算單元

702．．．窗口過濾處理單元

703．．．時域通道脈衝響應之平滑化處理單元

Claims

一種通道估算方法，用於中國移動多媒體廣播(CMMB)系統，該方法包含以下步驟：對第一數量之初估頻域通道響應值執行具第二數量取樣點之反快速傅立葉轉換運算，以產生第一時域通道脈衝響應，其中該第一數量大於該第二數量；對該第一時域通道脈衝響應執行一時域窗口過濾處理，以產生第二時域通道脈衝響應；對複數組不同時域之該第二時域通道脈衝響應執行一平滑化運算，以產生一平滑時域通道脈衝響應；以及對該平滑時域通道脈衝響應執行具該第一數量取樣點之快速傅立葉轉換運算，以產生一頻域通道響應；其中該第一數量為4096，該第二數量為512，該第一時域通道脈衝響應係由512個時域點之時域響應值組成，且所述具第二數量取樣點之反快速傅立葉轉換運算在所述512個時域點以外之時域點所產生之時域通道響應為鏡像信號。
如申請專利範圍第1項之方法，其中該窗口過濾處理係依一臨界數值及一窗口長度，對該第一時域通道脈衝響應進行處理，濾除窗口長度外及低於該臨界數值者，其中該窗口長度係回音延遲(echo delay)之一函數，而該臨界數值為窗口內該第一時域通道脈衝響各數值之一加權平均值。
如申請專利範圍第1項之方法，其中該平滑化運算為一平均運算及一加權運算其中之一。
一種通道估算裝置，用於中國移動多媒體廣播(CMMB)系統，該裝置包含：一反快速傅立葉轉換運算單元，用以對具第一數量響應值之初估頻域通道響應執行具第二數量取樣點之反快速傅立葉轉換運算，以產生第一時域通道脈衝響應，其中該第一數量大於該第二數量；一窗口過濾處理單元，用以對該第一時域通道脈衝響應執行一時域窗口過濾處理，以產生第二時域通道脈衝響應；一平滑化處理單元，用以對複數組不同時域之該第二時域通道脈衝響應執行一平滑化運算，以產生一平滑時域通道脈衝響應；以及一快速傅立葉轉換運算單元，用以對該平滑時域通道脈衝響應執行具該第一數量取樣點之快速傅立葉轉換運算，以產生一頻域通道響應；其中該第一數量為4096，該第二數量為512，該第一時域通道脈衝響應係由512個時域點之時域響應值組成，且所述具第二數量取樣點之反快速傅立葉轉換運算在所述512個時域點以外之時域點所產生之時域通道響應為鏡像信號。
如申請專利範圍第4項之裝置，其中該窗口過濾處理係依據一臨界數值及一窗口長度，對該第一時域通道脈衝響應進行過濾處理，濾除窗口長度外及低於該臨界數值者，其中該窗口長度係回音延遲(echo delay)之一函數，而該臨界數值為窗口內該第一時域通道脈衝響各數值之一加權平均值。
如申請專利範圍第4項之裝置，其中該平滑化運算為一平均運算及一加權運算其中之一。
如申請專利範圍第4項之裝置，其中該反快速傅立葉轉換運算單元包含：一緩衝器，其具有一第一緩衝部及一第二緩衝部，用以緩衝該初估頻域通道響應；一第一反快速傅立葉轉換計算器，耦接至該第一緩衝部，用以執行具該第二數量取樣點之第一反快速傅立葉轉換計算，以產生一第一反快速傅立葉轉換數值；一第一乘法器，用以將一第一相位變數及該第一反快速傅立葉轉換數值相乘，以產生一第一部份時域通道脈衝響應；一第二反快速傅立葉轉換計算器，耦接至該第二緩衝部，用以執行具該第二數量取樣點之第二反快速傅立葉轉換計算，以產生一第二反快速傅立葉轉換數值；一第二乘法器，用以將一第二相位變數及該第二反快速傅立葉轉換數值相乘，以產生一第二部份時域通道脈衝響應；以及一加法器，用以結合該第一部份時域通道脈衝響應及該第二部份時域通道脈衝響應，以產生該第一時域通道脈衝響應。
如申請專利範圍第7項之裝置，其中該第一緩衝部及該第二緩衝部各緩衝有該初估頻域通道響應之192個取樣點響應值。
如申請專利範圍第8項之裝置，其中該第一相位變數為e ^{j
2πn
(2+4i
)/4096} ,i=0,1；當i=0代表該初估頻域通道響應其取樣點分佈屬於該偶數型，當i=1代表該初估頻域通道響應其取樣點分佈屬於該奇數型。
如申請專利範圍第8項之裝置，其中該第二相位變數為e ^{j
2πn
(2559+4i
)/4096} ,i=0,1；當i=0代表該初估頻域通道響應其取樣點分佈屬於該偶數型，當i=1代表該初估頻域通道響應其取樣點分佈屬於該奇數型。
如申請專利範圍第4項之裝置，其中該反快速傅立葉轉換運算單元包含：一緩衝器，用以緩衝該初估頻域通道響應；一反快速傅立葉轉換計算器，耦接至該緩衝器，用以執行具該第二數量取樣點之反快速傅立葉轉換計算，以產生一反快速傅立葉轉換數值；一選擇開關，用以依據該通道之位置，自一第一相位變數及一第二相位變數中擇一作為一輸出相位變數；一乘法器，用以將該輸出相位變數及該反快速傅立葉轉換數值相乘，產生一暫時時域通道脈衝響應；一時域通道脈衝響應緩衝器，耦接至該乘法器，用以緩衝該暫時時域通道脈衝響應，其中該暫時時域通道脈衝響應包含與第一相位變數有關之一第一部份時域通道脈衝響應及與第二相位變數有關之一第二部份時域通道脈衝響；以及一加法器，用以結合該第一部份時域通道脈衝響應及該第二部份時域通道脈衝響應，以產生該第一時域通道脈衝響應。
如申請專利範圍第11項之裝置，其中該緩衝器緩衝有該初估頻域通道響應之192個取樣點響應值。
如申請專利範圍第11項之裝置，其中該第一相位變數為e ^{j
2πn
(2+4i
)/4096} ,i=0,1；當i=0代表該初估頻域通道響應其取樣點分佈屬於該偶數型，當i=1代表該初估頻域通道響應其取樣點分佈屬於該奇數型。
如申請專利範圍第11項之裝置，其中該第二相位變數為e ^{j
2πn
(2559+4i
)/4096} ,i=0,1；當i=0代表該初估頻域通道響應其取樣點分佈屬於該偶數型，當i=1代表該初估頻域通道響應其取樣點分佈屬於該奇數型。
如申請專利範圍第4項之裝置，其中該平滑化處理單元包含有：一時域通道脈衝響應儲存單元，用以儲存該複數組不同時域之該第二時域通道脈衝響應；以及一平滑計算器，耦接至該時域通道脈衝響應儲存單元，用以對該複數組不同時域之該第二時域通道脈衝響應執行該平滑化運算，以產生該平滑時域通道脈衝響應。
一種通道估算方法，用於中國移動多媒體廣播(CMMB)系統，該方法包含以下步驟：對初估頻域通道響應執行一反快速傅立葉轉換運算，以產生第一時域通道脈衝響應，其中該初估頻域通道響應具有第一數量之響應值，該反快速傅立葉轉換運算具有第二數量之取樣點，且該第一數量大於該第二數量；對該第一時域通道脈衝響應執行一時域窗口過濾處理，以產生第二時域通道脈衝響應；對複數組不同時域之該第二時域通道脈衝響應執行一平滑化運算，以產生一平滑時域通道脈衝響應；以及對該平滑時域通道脈衝響應執行一具該第一數量取樣點之快速傅立葉轉換運算，以產生一頻域通道響應；其中該第一數量為4096，該第二數量為512，該第一時域通道脈衝響應係由512個時域點之時域響應值組成，且所述具第二數量取樣點之反快速傅立葉轉換運算在所述512個時域點以外之時域點所產生之時域通道響應為鏡像信號。