[發明專利]可配置變換長度DFT的頻域補償方法及裝置

無效

• 申請號: 200910048198.X
• 申請日: 2009-03-25
• 公開（公告）號: CN101521649A
• 公開（公告）日: 2009-09-02
• 發明（設計）人: 呂正德
• 申請（專利權）人: 呂正德
• 主分類號: H04L27/26
• 分類號: H04L27/26;G06F17/14
• 代理公司: 上海漢聲知識產權代理有限公司
• 代理人: 胡 晶

1.一種可配置變換長度DFT的頻域補償方法，適用于補償因子的幅度可以用式A(k)＝P(x₁)·S(x₂)+Q(x₁)逼近的情況，其中S(x₂)是反映前處理單元的基本幅度頻率特性的倒數的函數，P(x₁)和Q(x₁)為多項式，變量x₁和x₂和頻率值k成線性關系，其特征在于，包括以下步驟：

(1)建立并存儲包括補償因子擬合系數查找表在內的補償信息，所述補償因子擬合系數查找表保存變換長度N和與之對應的擬合系數，所述補償信息中還包括擬合曲線公式所需自變量x₁和x₂的計算方法；

(2)每一次頻域補償計算時，先根據變換長度N查找補償因子擬合系數查找表獲得擬合系數，再計算出擬合自變量，后分別通過曲線擬合公式和直線擬合公式計算出頻率點k處的頻域補償因子的幅度A(k)和補償因子的相位φ(k)，其中，

其中k′為需要進行頻域補償的頻率點的序號；

(3)根據頻率點k處頻域補償因子的幅度和相位，將復數的頻域補償因子轉化為笛卡爾坐標表示形式，然后根據共軛對稱關系得到最后所需要的頻率點k′處的補償因子；

(4)將得到的頻域補償因子與可配置FFT快速傅立葉變換處理器輸出的對應頻率點序號為k′處的數據進行復數相乘，完成頻域補償過程。

2.如權利要求1所述的方法，其特征在于，步驟(1)建立并存儲補償信息進一步包括：建立并保存S(x₂)＝1/sinc²(cx₂)函數查找表，其中c是曲線擬合的系數，根據其自變量x＝cx₂的取值范圍，將自變量在這個范圍內進行若干等分，得到對應的1/sinc²x函數值，(x，1/sinc²x)的關系便可以構成一個查找表，其中由x可計算出查找表的對應的地址，1/sinc²x對應查找表的內容。?

3.如權利要求1所述的方法，其特征在于，步驟(1)建立并存儲補償信息進一步包括：建立并保存cosx函數查找表，將余弦函數cosx的自變量x在[0，π/2]范圍內進行若干等分，得到對應的cosx函數值，(x，cosx)的關系便可以構成一個查找表，其中由x可計算出查找表的對應的地址，cosx對應查找表的內容。

4.如權利要求1所述的方法，其特征在于，還包括：

對于每一種變換長度N，只需要計算頻率序號為0，1，…，?的頻域補償因子，而?…，N-1范圍內的頻域補償因子利用共軛對稱特性C(k)＝C^*(N-k)由前述的補償因子得到。

5.如權利要求1所述的方法，其特征在于，步驟(2)進一步包括：對于頻率序號k在0，1，…，?范圍內的?個頻域補償因子的幅度，采用曲線擬合公式A(k)＝P(x₁)·S(x₂)+Q(x₁)計算得到；曲線擬合公式中P(x₁)和Q(x₁)是以x₁為自變量的多項式，S(x₂)是以x₂為自變量并且是與可配置變換長度DFT離散傅立葉變換的前處理單元的幅頻特性有關的函數。

6.如權利要求1所述的方法，其特征在于，步驟(2)進一步包括：

采用A(k)＝P(x₁)·S(x₂)+Q(x₁)計算補償因子幅度，其中P(x₁)和Q(x₁)使用一次多項式，即P(x₁)＝a+bx₁和Q(x₁)＝d+ex₁；對于可配置變換長度DFT離散傅立葉變換前處理單元使用線性插值的情況下，S(x₂)＝1/sinc²(cx₂)；

其中a，b，c，d，e為計算補償因子幅度所用到的擬合系數。