取樣速率轉換器及其方法和音頻裝置的制作方法

            文檔序號:7537626閱讀:455來源:國知局
            專利名稱:取樣速率轉換器及其方法和音頻裝置的制作方法
            技術領域
            本申請涉及取樣速率轉換器及其方法以及音頻裝置,所述取樣速率轉換器可以被應用于例如音頻取樣頻率的轉換以及用于圖像的幀的放大或縮小的分辨率的轉換。
            背景技術
            在音頻和圖像的數字信號處理中,經常使用濾波處理。由于線性相位是通過有限個分支而得到的特征,經常利用線性相位FIR(有限脈沖響應)濾波器作為用于濾波處理的濾波器。
            圖1是示出了線性相位FIR濾波器的橫向(transversal)型電路的配置的圖。
            如圖1所示,該線性相位FIR濾波器1具有相對于輸入端TIN串聯連接的且構成移位寄存器的(n-1)個延遲單元2-1到2-n-1,用于將輸入到輸入端TIN的信號和延遲單元2-1到2-n-1的輸出信號乘以濾波系數h(0)到h(n-1)的n個乘法器3-1到3-n,以及用于將n個乘法器3-1到3-n的輸出信號相加并將結果輸出到輸出端TOUT的加法器4。
            作為這種線性相位FIR濾波器的代表性設計方法,例如已知由Parks,T.W.和McClellan、J.H.等人應用于線性相位FIR濾波器的Remez交換算法(參見例如非專利文獻1)。
            Remez交換算法是用于近似的算法,使得加權近似誤差顯示出相對于期望的振幅特性的等波紋形狀。
            作為使用線性相位FIR濾波器的濾波處理的應用,有利用取樣速率轉換的圖像分辨率的轉換以及音頻取樣頻率的轉換。
            例如,在分辨率的轉換中,使用利用內插器的多速率濾波器、分米器(decimeter)和線性相位FIR濾波器作為元件技術(參見例如非專利文獻2)。
            在多速率濾波器中,一般線性相位FIR濾波器被使用與內插器多相分解匹配。內插器和分米器兩者都是循環的時間不變式系統,并且具有與時間不變式系統的特征不同的特性特征。
            由于內插器的循環時間不變性,在圖像分辨率的轉換中會出現稱為“棋盤失真”(chessboard distortion)的格子的失真。
            因此,Harada和Takaie考慮用于避免棋盤失真濾波器零點排列的條件(參見非專利文獻3)。
            沒有棋盤失真的多速率濾波器的傳送函數H(z)通過將由某種方法構出的線性相位FIR濾波器(以下稱為補償器)的傳送函數K(z)與零點傳送函數Z(z)相乘得到,以避免后面的棋盤失真。
            (等式1)H(z)=Z(z)·K(z)(1)(等式2)Z(z)=1+z-1+Z-2+…+Z-(U-1)(2)這里,用于避免棋盤失真的如零點傳送函數Z(z)的事先固定的線性相位FIR濾波器將被稱作預濾波器。
            圖2A到2C示出了通過將通過Remez交換算法設計的補償器與預濾波器和加權近似誤差相乘來避免棋盤失真的多速率濾波器的頻率響應的示例。
            非專利文獻1Parks,T.W.和McClellan,J.H.“ChebyshevApproximation for Nonrecursive Digital Filters with Linear Phase”,IEEE Trans.Circuit Theory,CT-19,2,189-194頁,1972,以及Rabiner,L.R.,Mcclellan,J.H.和Parks,T.W.“FIR Digital Filter DesignTechniques Using Weighted Chebyshev Approximation”,Proc.IEEE,vol.63,4月,595-610頁,1975;非專利文獻2Takaie,Hitoshi,Multi-rate Signal Processing,Shokodo,1997;
            非專利文獻3Harada,Yasuhiro和Takaie,HitoshiMulti-rateFilter not Accompanies by Chessboard Distortion and Zero PointArrangement of Same,Shingaku Giho CAS96-78,1-6頁,1997-01。

            發明內容
            本發明將要解決的問題但是,通過上述方法避免棋盤失真的方法具有如下缺點。
            也就是說,在具有用傳統方法設計的傳送函數H(z)的多速率濾波器中,如圖2C所示,通過Remez交換算法設計的加權近似誤差的等波紋塌陷。
            此外,在用傳統方法設計的多速率濾波器中,如圖2B所示,通帶增益不是常數值,其右端衰減。
            當通過使用這種濾波器進行分辨率的轉換時,圖像的輪廓變得模糊,并影響圖像質量。同樣地,當通過使用這種濾波器進行取樣頻率的轉換時,高頻分量衰減,并且不能進行高精度的音頻復制。
            即使在濾波器系數增加時,通帶處的衰減也不能避免。
            此外,有這樣的情況,即當在線性相位FIR濾波器的設計規格中頻率w為0時,DC增益必須被設定在1。
            但是,在Remez交換算法中,如圖3A到3C所示,不能得到經過任何指定頻率點的振幅特性。
            此外,圖4示出了基于傳統多速率信號處理理論的取樣速率轉換的配置的實例。
            圖4的取樣速率轉換單元10具有向上取樣器11、FIR濾波器12和向下取樣器13。
            此外,在圖4中,“U”和“D”是互為本原的正整數,H(z)指FIR濾波器的傳送函數。此外,向上的箭頭表示用于在信號之間插入(U-1)個零點的向上取樣器,向下的箭頭表示用于以D信號的間隔使信號稀疏的向下取樣器。
            圖4的取樣速率轉換單元10通過向上取樣器11增加信號到Y倍的取樣速率,通過使用FIR濾波器12進行帶限制,最后通過向下取樣器13進行將取樣速率降低到1/D的操作。
            這樣,取樣速率可以被轉換到U/D倍。該FIR濾波器的截止頻率如下。
            (等式3)ωC=π/maxU,D (3)例如,為了基于傳統多速率信號處理理論實現從44.1kHz或直到44.1kHz的取樣速率轉換,如圖4所示,必須準備其中低通濾波器的FIR濾波器2的截止頻率極度嚴格的濾波器。
            一般而言,如下面所示,這種濾波器需要大量的分支,并且極難實現。
            從44.1kHz到48kHz,Up=160,Down=147,并且截止頻率π/160。
            從44.1kHz到32kHz,Up=320,Down=441,并且截止頻率π/441。
            從44.1kHz到24kHz,Up=80,Down=147,并且截止頻率π/147。
            從48kHz到44.1kHz,Up=147,Down=160,并且截止頻率π/160。
            從32kHz到44.1kHz,Up=441,Down =320,并且截止頻率π/441。
            從24kHz到44.1kHz,Up=147,Down =80,并且截止頻率π/147。
            本發明是考慮這種情況做出的,其目的是提供取樣速率轉換器及其方法和音頻裝置,在所述取樣速率轉換器中加權近似誤差的等波紋不塌陷,通帶的增益可以保持為基本為常數值,可以得到通過任何頻率的振幅特性,并且可以實現高精度轉換。
            解決目的的途徑為了實現上述目的,根據本發明第一方面的取樣速率轉換器具有向上取樣器,用于在樣本信號之間插入U-1個零點并增加取樣頻率到U倍;卷積處理單元,包括FIR濾波器,并進行相對于向上取樣器輸出信號的預定卷積處理;以及線性內插塊,用于選擇相對于卷積處理單元的處理結果的兩個樣本點,并根據線性內插找到處于所需位置的值,其中,卷積處理單元的FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            根據本發明第二方面的取樣速率轉換器包括多個卷積處理單元,包括通過多相分解預定的FIR濾波器以及進行輸入樣本信號和被分解成多相的多相濾波器的卷積處理得到的預相濾波器;多個向上取樣器,用于在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍;加法裝置,用于在相加所有信號之后通過調整多個向上取樣器輸出信號的傳播時間產生信號;以及線性內插塊,用于相對于加法裝置得到的信號選擇兩個樣本點以及根據線性內插找到所需位置處的值,其中,FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            根據本發明第三方面的取樣速率轉換器包括卷積處理單元,包括能夠設定通過多相分解預定FIR濾波器而得到的不同濾波器系數的多相濾波器,并進行輸入樣本信號與具有所選系數的多相濾波器的卷積處理;選擇器,用于選擇輸入樣本所需的兩個樣本點,并選擇相應多相濾波器的系數;以及線性內插塊,用于根據線性內插找到處于所需位置的值,其中FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            根據本發明第四方面的取樣速率轉換方法包括第一步驟,在樣本信號之間插入U-1個零點并增加取樣頻率到U倍;第二步驟,進行相對于其取樣頻率被包括FIR濾波器的卷積處理單元乘以U的信號的預定卷積處理,在FIR濾波器中,脈沖響應由有限時間長度表示,脈沖響應變為濾波器系數,并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯;以及第三步驟,相對于處理結果選擇兩個樣本點,并根據線性內插找到所需位置處的值,其中FIT濾波器的濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來計算。
            根據本發明第五方面的取樣速率轉換方法包括第一步驟,進行輸入樣本信號與被多個卷積處理單元分解成多相的多相濾波器的卷積處理,所述多個卷積處理單元包括通過多相分解預定FIR濾波器而得到的多相濾波器;第二步驟,在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍;第三步驟,調整多個信號的傳播時間使取樣頻率被增加到U倍,并產生通過相加所有信號得到的信號;以及第四步驟,相對于第三步驟得到的信號選擇兩個樣本點,并根據線性內插找到所需位置處的值,其中,FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來計算。
            根據本發明第六方面的取樣速率轉換方法包括第一步驟,選擇輸出樣本所需的兩個樣本點,并選擇相應多相濾波器的系數;以及第二步驟,通過卷積處理單元進行樣本信號與具有所選系數的多相濾波器的卷積處理,所述卷積處理單元包括通過多相分解預定FIR濾波器得到的多相濾波器并能夠設定不同的濾波器系數,其中FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來計算。
            本發明的第七方面是一種包括取樣速率轉換器的音頻裝置,其中,取樣速率轉換器包括向上取樣器,用于在樣本信號之間插入U-1個零點并增加取樣頻率到U倍;卷積處理單元,包括FIR濾波器,并進行相對于向上取樣器輸出信號的預定卷積處理;以及線性內插塊,用于相對于卷積處理單元的處理結果選擇兩個樣本點,并根據線性內插找到處于期望位置的值,其中,卷積處理單元的FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            本發明的第八方面是一種包括取樣速率轉換器的音頻裝置,其中,取樣速率轉換器包括多個卷積處理單元,包括通過多相分解預定的FIR濾波器以及進行輸入樣本信號和被分解成多相的多相濾波器的卷積處理得到的預相濾波器;多個向上取樣器,用于在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍;加法裝置,用于在相加所有信號之后通過調整多個向上取樣器輸出信號的傳播時間產生信號;以及線性內插塊,用于相對于加法裝置的信號選擇兩個樣本點以及根據線性內插找到所需位置處的值,其中,FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            本發明的第九方面是一種包括取樣速率轉換器的音頻裝置,其中,取樣速率轉換器包括卷積處理單元,包括能夠設定通過多相分解預定FIR濾波器而得到的不同濾波器系數的多相濾波器,并進行輸入樣本信號與具有所選系數的多相濾波器的卷積處理;選擇器,用于選擇輸入樣本所需的兩個樣本點,并選擇相應多相濾波器的系數,其中FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            根據本發明,例如,設計了FIR濾波器。例如,通過初始設置,線性相位FIR濾波器被設定,帶被設定,預濾波器的系數被設定,將要通過的任何頻率點被輸入,并且初始極值點被設定。
            接下來,用于根據當前極值點和將要通過的頻率點內插振幅特性的內插多項式被產生。
            接下來,新的極值點根據從所產生的內插多項式發現的振幅特性被確定。
            這些重復進行,并且判斷例如極值位置是否近似在期望范圍內。
            然后,從近似的振幅特性中找到濾波器特性。
            在這種方式下,在其中具有系數設定的FIR濾波器中,加權近似誤差變為等波形,并且在通帶處的增益被保持在常數值。
            此外,可以通過指定的頻率點。
            然后,向上取樣器通過在樣本信號之間插入U-1個零點增加取樣頻率到U倍。
            接下來,相對于其取樣頻率被包括如上設計的FIR濾波器的卷積處理單元乘以U的信號進行預定的卷積處理。
            接下來,其選擇相對于卷積處理單元的處理結果的兩個樣本點,并且根據線性內插找到所需位置處的值。
            發明效果根據本發明,具有嚴格截止頻率的取樣速率轉換是可能的。
            此外,可能避免棋盤失真。
            此外,任何預濾波器都可以考慮,并且可能通過任何頻率點。
            此外,可能將處理量抑制到所需的最低限,并且可以實現處理速度的改進。


            圖1是示出了FIR濾波器的橫向型電路的配置的圖。
            圖2A到2C是示出了傳統方法中避免棋盤失真和加權近似誤差的頻率響應的實例的圖。
            圖3A到3C是傳統方法中增益1的頻率響應和附近的放大圖。
            圖4是示出了一般取樣速率轉換器的配置的實例的圖。
            圖5是示出了根據本發明的取樣速率轉換器的第一實施例的配置的圖。
            圖6是概念性地示出了根據本實施例的線性內插塊的線性內插處理的圖。
            圖7是概念性地示出了根據本發明實施例通過線性內插塊的線性內插處理來找到在所需位置處的值的處理的圖。
            圖8A到8D是示出了其中FIR濾波器具有線性相的四種情況的脈沖響應的圖。
            圖9是示出了關于線性相位FIR濾波器的四種情況的Q(ejω)和R的圖。
            圖10是示出了加權切比謝夫(chebyshev)近似的實例的圖。
            圖11是根據本發明考慮了預濾波器的頻率響應的Remez交換算法的流程圖。
            圖12A到12C是說明用于確定加權近似誤差E(ejω)的新極值的方法的圖。
            圖13A到13D是示出了當指定本發明的任何頻率點時的頻率響應的圖及其放大圖。
            圖14是示出了根據本發明的取樣速率轉換器的第二實施例的配置的圖。
            圖15是示出了根據本發明的取樣速率轉換器的第三實施例的配置的圖。
            圖16是示出了根據本發明的取樣速率轉換器的第四實施例的配置的圖。
            圖17是概念性地示出了根據本發明的取樣速率轉換器的第五實施例的配置的圖。
            圖18是更詳細地示出了在圖17中概念性示出的取樣速率轉換器的圖。
            圖19是用于說明線性內插塊中使用的變量的圖。
            圖20A和20B是用于具體說明MasterCount的圖。
            圖21是概念性示出了根據本發明的取樣速率轉換器的第六實施例的圖。
            圖22是更具體地示出了圖21中概念性示出的取樣速率轉換器的圖。
            圖23是說明根據本實施例的取樣速率轉換器的第一安裝方法的流程圖。
            圖24是說明根據本實施例的取樣速率轉換器的第二安裝方法的流程圖。
            圖25是示出了使用根據本發明的取樣速率轉換器的音頻裝置配置的實例的框圖。
            標號說明1--線性相位FIR濾波器,2-1到2-n-1--延遲單元,3-1到3-n--乘法器,4--法器,h(0)到h(n-1)--濾波器系數,TIN--輸入端,TOUT--輸出端,100、100A--樣速率轉換器,101--輸入端,102--LPF,103--向上取樣器,104--卷積處理單元,105--線性內插塊,106--輸出端,200、200A到200E--取樣速率轉換器,201--輸入端,202--LPF,203-1到203-U--卷積處理單元,204-1到204-U--向上取樣器,205-1到205-U-1--延遲單元,206-1到206-U-1--加法器,207--向下取樣器,208--輸出端,209--LPF,210--選擇器,211--輸入端,212--FIR濾波器設計單元,213--第一存儲器,214--LPF卷積處理單元,215--第二存儲器,216--輸入緩沖器,217--卷積處理單元,218--計數器控制單元,219--線性內插處理單元,220--輸出端,230--選擇器,231--輸入端,232--FIR濾波器設計單元,233--第一存儲器,234--輸入緩沖器,235--卷積處理單元,236--計數器控制單元,237--線性內插處理單元,238--LPF,卷積處理單元,239--第二存儲器,240--輸出端,300--音頻裝置,301--輸入端,302--LR分離電路(DSB),303--取樣速率轉換器(SRC),304--衰減器(ATT),305--靜音電路(MUTE),306--輸出端。
            具體實施例方式
            下面,將結合附圖給出本發明優選實施例的詳細描述。
            <第一實施例>
            圖5是示出了根據本發明的取樣速率轉換器的第一實施例的配置的圖。
            在圖5中,“U”和“D”是互為本原的正整數,H(z)表示FIR濾波器的傳送函數。此外,向上的箭頭表示用于在信號之間插入(U-1)個零點的向上取樣器。
            此外,Fsi表示輸入的取樣頻率,FSO表示輸出的取樣頻率。第一實施例示出了其中FSO<Fsi的配置的實例。
            也就是說,如圖5所示,取樣速率轉換器100具有輸入端101、低通濾波器(LPF)102、向上取樣器103、卷積處理單元104、線性內插塊105和輸出端106。
            輸入端101接收具有取樣頻率Fsi的取樣信號x(n)作為輸入。
            LPF 102防止(抑制)從輸入端101輸入的取樣頻率Fsi的取樣信號x(n)發生折疊,因為在輸入的取樣頻率Fsi高于輸出的取樣頻率FSO且如此輸出到向上取樣器103時,會產生重疊部分并會出現折疊。
            向上取樣器103經由LPF 102接收從輸入端101輸入的具有取樣頻率Fsi的樣本數據x(n),插入(內插)U-1個零點,增加取樣頻率Fsi到Up倍(過取樣),并將具有取樣頻率UFsi的樣本信號輸出到卷積處理單元104。
            卷積處理單元104包括基于后面解釋的Remez交換算法設計的FIR濾波器,進行在下面的等式中示出的卷積處理(進行帶限制),并在下一步將處理結果輸出到線性內插塊105。卷積處理單元104通過具有1/Up截止頻率的低通濾波器(傳送頻率H(z))的卷積來內插值。
            y(n′)=Σi=0∞h(i)x(n′-1)---(4)]]>這里,h(n)是FIR濾波器的脈沖響應。卷積的輸出(在向下取樣器之前)是通過向上取樣器插入零點得到的樣本。
            線性內插塊105從卷積處理單元104的如圖6和圖7中所示具有取樣頻率UFsi的輸出信號中選擇兩個樣本點,根據如下等式中的線性內插找到期望位置處的值(圖7),并將其作為樣本信號y(m)從輸出端106輸出。
            (等式5)y(m)=αxIn A+βxIn Bβ=1-α(5)下面,將詳細描述卷積處理單元104的FIR濾波器的設計方法。
            根據本實施例的線性相位FIR濾波器可以等同地使用例如如圖1所示的橫向型電路配置。
            注意,根據近似的振幅發現,如將在下面詳細說明的,濾波器系數h(n)具有如下特性,其擴展Remez交換算法,指定將被通過的頻率點,考慮預濾波器的頻率響應,并切比謝夫近似期望的振幅特性。
            下面,將順序參照

            根據本發明設定線性相位FIR濾波器的系數的具體方法。
            如等式(6)所示,N個分支的線性相位FIR濾波器的傳送函數H(z)是從預濾波器的傳送函數Z(z)和補償器的傳送函數K(z)的成績得出的濾波器。
            (等式6)H(z)=Z(z)·K(z)(6)這里,假設預濾波器和補償器是U個分支和N-(U-1)個分支的線性相位FIR濾波器,并且預濾波器的傳送函數被事先給出。
            此外,在頻率域中通過任何Np數量的頻率點。因此,這里的傳送函數H(z)的濾波器設計就要確定N-(U-1)個分支的補償器的傳送函數K(z),使得被指定的頻率點被通過,并且振幅特性H(ejω)趨近期望的振幅特性D(ejω)。
            分配給具有傳送函數K(z)的補償器的分支數被設定為L=N-(U-1)。
            線性相位FIR濾波器的傳送函數K(z)具有如圖8A到8D所示的線性相位,所以被分成四種情況。
            具體地說,其被分成如圖8A所示的奇數個分支和偶對稱的情況1,如圖8B所示的偶數個分支和偶對稱的情況2,如圖8C所示的奇數個分支和奇對稱的情況3,以及如圖8D所示的偶數個分支和奇對稱的情況4。
            然后,振幅特性函數K(ejω)在情況1中被保持為原樣,在情況2到4中被重寫如下[等式7]情況1Σn-0(L-1)/2a(n)cos(nω)---(7-1)]]>
            情況2Σn-1L/2b(n)cos{(n-1/2)ω}---(7-2)]]>=cos(ω2)Σn=2L/2-1b~(n)cos(nω)]]>情況3Σn=1(L-1)/2c(n)sin(nω)]]>=sin(ω)Σn=0(L-3)/2c~(n)cos(nω)]]>情況4Σn=1L/2d(n)sin{(n-12)ω}]]>=sin(ω2)Σn=0L/2-1d~(n)cos(nω)]]>也就是說,振幅特性系數K(ejω)由固定參數的函數Q(ejω)與包括如圖9所示設計參數的余弦級數P(ejω)的乘積表示。下面,等式(7-1)到等式(7-4)的上限將由R-1+2xNp表示。也就是說,R如圖9中所示來計算。此外,a(n);-b(n);-c(n)和-d(n)在一起被統稱為“p(n)”。
            當期望的振幅特性被定義為D(ejω)并且相對于每一個頻率的權重是W(ejω)時,加權近似誤差被如下定義[等式8]E(ejω)=W(ejω){D(ejω)-H(ejω)}(8)[等式9]H(ejω)=K(ejω)·Z(ejω)=Q(ejω)·P(ejω)·Z(ejω)(9)當將等式(9)代入等式(8)時,結果如下。
            E(ejω)=W^(ejω){D^(ejω)-P(ejω)}---(10)]]>注意, 和 如下。
            W^(ejω)=W(ejω)·Q(ejω)·Z(ejω)---(11)]]>[等式12]D^(ejω)=D(ejω)/Q(ejω)·Z(ejω)---(12)]]>等式(10)表示情況1到情況4的四種情況中線性相位FIR濾波器的加權近似誤差。
            加權切比謝夫近似問題是要確定等式(7-1)到等式(7-4)的這些a(n);-b(n);-c(n)和-d(n),其使等式(8)中指定頻率帶中的|E(ejw)|的最大值最小化。
            下面,將參照具體的實例對此進行說明。
            這里,如下面的描述和圖10所示,定義振幅特性D(ejω)。
            D (ejω)=1(在±δ1誤差范圍內,0<ω<ωp)D(ejω)=0(在±δ2誤差范圍內,ωs<ω<π)(13)注意,當給定R時,δ1和δ2的值不能被自由指定,但是其比率可以被指定。
            W(ejω)在通帶處被設定為常數值W1,而在阻擋頻率下被設定為W2并且被選擇為使得W1δ1=W2δ2成立。例如,其被選擇為使得W1=1且W2=δ1/δ2。這時,下面的反射定理成立。
            <定理>
            第(R-1)級余弦級數P(ejω)成為相對于w的(0,π)部分中目標特性的最佳加權切比謝夫近似所需的充分條件是(1)E(ejω)在(0,π)部分中變為極值至少(R+1)次。此時取得極值所處的頻率被設定為w0<w1<w2<...<wR-1<wR。
            (2)相鄰極值的符號是不同的,所有極值的絕對值相等。也就是說,滿足下面的條件。
            E(ejωi)·E(ejωiii)<0(i=0,1,-,R-1)(E(ejωi)|=|E(ejωiii)|(i=0,1,-,R-1)(14)
            因此,|E(ejwi)|等于該部分中|E(ejw)|的最大值。
            用于獲得最佳切比謝夫近似的技術包括基于反射定理(參見Raber,L.R.,McClellan,J.H.和Parks,T.W.“FIR Digital FilterDesign Techniques Using Weighted Chebyshev Approximation”,Proc.IEEE,vol.63,April,pp.595-610,1975)的Remez交換算法。
            Remez交換算法切比謝夫近似在頻率域內的期望振幅特性,并從近似的振幅特性中找出線性相位FIR濾波器的系數。
            圖11是根據本發明的經過任何頻率點且考慮了預濾波器頻率響應的Remez交換算法的流程圖。
            考慮了具體預濾波器的頻率響應的Remez交換算法如下。
            <步驟0>
            如圖11所示,首先,進行初始設置(F101)。在該初始設置中,線性相位FIR濾波器被設定,帶被設定,預濾波器的系數被設定,將要通過的任何頻率點被輸入,并且初始極值點被設定。
            具體地說,設定項如下。
            ·分支數,·偶對稱或奇對稱的線性相位FIR濾波器,·帶數,·每一個帶的兩個端頻率,·每一個帶的期望振幅值,·相對于每一個帶的加權,·預濾波器系數,·將要通過的頻率點和振幅值(WR+1,D(ejωR+1),i=1,...,Np),·在近似帶中變為極值的頻率w(0)=wk(0)(k=0,...,R)。
            注意,上標(i)表示重復次數。
            <步驟1>
            接下來,用于從當前極值點內插振幅特性的拉格朗日內插多項式被產生(F102)。
            在上面等式(8)中示出的切比謝夫近似的目標函數要變為最小值所需的充分條件由反射定理表示。因此,基于反射定理,下面等式的參數p(n)被發現,使得來自期望振幅特性的加權近似誤差δ(1)在每一個頻率點都是相等的,并且符號相反。
            p(ejω)=Σn=0R-1p(n)cos(nω)---(15)]]>也就是說,等式(9)在頻率點w(i)=wk(i)(k=0,...,R)處的加權近似誤差滿足如下等式[等式16]W^(ejωk(i)){D^(ejωk(i))-P(ejωk(i))}=(-1)kδ(i)(k=0,1,···,R)---(16)]]>下面,為了簡化,省略了下標(i)。在修改等式(16)時,其變成下面的等式。
            P(ejωk)+(-1)kδ/W^(ejωk)=D^(ejωk)]]>Σn=0R-1p(n)cos(nωk)+(-1)kδW^(ejωk)=D^(ejωk)(k=0,1,···,R)---(17)]]>然后,頻率域中將被通過的點的等式作為限制被添加到等式(17)[等式18]P(ejωk)=D^(ejωk)]]>Σn=0R-1p(n)cos(nωk)=D^(ejωk),(k=R+1,···,R+Np)---(18)]]>當用矩陣表示等式(17)和等式(18)時,它們變成下面的等式。

            =D^(ejω0)D^(ejω1)···D^(ejωR-1)D^(ejωR)D^(ejωR+1)···D^(ejωR+Rp)---(19)]]>但是,解該等式需要非常大量的處理,所以將首先分析δ。
            δ=Σj=0RαjD^(e)jω1Σj=0R(-1)jαj/W^(ejω1)---(20)]]>[等式21]αk=Πj≠kR1(xk-xj)---(21)]]>[等式22]xj=cos(ωj)(22)αk是矩陣F中第k行第(R+1)列元素的余項因子。注意,對于 和 使用等式(11)和等式(12)。
            接下來,通過使用該δ來設定下面的等式。
            Ck=D^(ejωk)-(-1)kδW^(ejωk),k=1,...,R---(23)]]>[等式24]Ck=D^(ejωk),(k=R+1,···,R+Np)---(24)]]>為了找到除極值點之外的頻率的振幅特性,這次通過使用極值點和將要通過的頻率點,來使用拉格朗日內插多項式作為用于內插的內插多項式。也就是說,P(ejω)通過進行內插來計算,從而通過使用拉格朗日內插多項式由wk(k=0,...,R+Np)來取得值Ck。
            P(e)jω)=Σk=0R+NpCk(βkx-xk)Σk=0R+Np(βkx-xk)---(25)]]>[等式26]βk=Πj=0j≠kR+Np1(xk-xj)---(26)]]>[等式27]x=cos(ω) (27)該結果符合正被求解的等式(19)。
            <步驟2>
            從于內插多項式中發現的振幅特性中找出新的極值點(F103)。反復判斷是否得到了最佳近似(F104)。
            上述步驟1的結果的每一個極值點wk不總是成為加權誤差函數E(ejω)的極值。有時存在|E(ejw)|變得大于δ(1)的點。因此,新的極值點w(i+1)通過所有點的替換方法確定。
            <所有點的同時替換方法>
            基于下面的等式,在整個近似帶上搜索根據用于內插的極值點算出的加權近似誤差的極值。這被定義為新的極值點w(i+1)=wk(i+1)(k=0,1,...,R),然后程序返回到步驟1的處理。
            E(ejω)=W^(ejω){D^(ejω)-P(ejω)}---(28)]]>假設在極值的位置不再改變時得到最佳近似。這是重復結束的條件。程序然后進行到下一個步驟3的處理。
            圖12A到12C是所有點替換方法的概念圖。
            簡單對此進行說明,圖12A到12C中的黑點表示用于內插的極值點。從該極值點中找到的加權近似誤差E(ejω)對應于實線。
            如圖12A所示,在黑點的極值點處的加權近似誤差的值變為白點,但是實際的極值是由方塊表示的頻率。因此,程序返回到步驟1的處理,同時將由方塊表示的頻率定義為新的極值點。
            此外,如圖12B所示,用于內插的極值點的頻率和實際極值的頻率是偏離的,因此,由方塊表示的頻率被定義為新的極值點,然后程序返回到步驟1的處理。
            然后,如圖12C所示,當用于內插的極值點和實際的加權近似誤差的極值點(白點)相同時,反復結束。
            <步驟3>
            線性相位FIR濾波器的系數從近似的振幅特性中找到(F105)。
            當從最佳近似函數P(ejω)中找到N個分支的脈沖響應h(n)時,是從下面的等式中找,來替代從p(n)中找。
            H(ejω)=P(ejω)·Q(ejω)·Z(ejω) (29)[等式30]情況1h(n)=1N{H(0)}+2Σk=1N-12(-1)kH(2πNk)cos(2πNk(n+12))}---(30)]]>[等式31]情況2
            h(n)=2NΣk=0N2-1(-1)kH(2πN(k+12))sin(2πN(k+12)(n+12))---(31)]]>[等式32]情況3h(n)=-2NΣk=0N-12(-1)kH(2πNk)sin(2πNk(n+12))---(32)]]>[等式33]情況4h(n)=2NΣk=0N2-1(-1)kH(2πN(k+12))cos(2πN(k+12)(n+12))---(33)]]>此外,當找到L=N-(U-1)個分支的補償器的脈沖響應k(n)時,根據如下等式計算。
            K(ejω)=P(ejω)·Q(ejω) (34)[等式35]情況1k(n)=1L{K(0)}+2Σk=1L-12(-1)kK(2πLk)cos(2πLk(n+12))}---(35)]]>[等式36]情況2k(n)=2LΣk=0L2-1(-1)kK(2πL(k+12))sin(2πL(k+12)(n+12))---(36)]]>[等式37]情況3k(n)=-2LΣk=0L-12(-1)kK(2πLk)sin(2πLk(n+12))---(37)]]>[等式38]情況4
            k(n)=2LΣk=0L2-1(-1)kK(2πL(k+12))cos(2πL(k+12)(n+12))---(38)]]>如果預濾波器的傳送函數Z(z)是如下面等式中示出的1,那么與通過任何頻率點的Remez交換算法是一樣的。
            (等式39)Z(z)=1(39)此外,在其中不存在將被通過預濾波器的頻率點的情況下,當Np=0時,其與考慮了預濾波器的頻率響應的Remez交換算法一樣。
            此外,當預濾波器的傳送函數Z(z)如在下面等式中所示為1且不存在將要通過的頻率點時,當Np=0時,其與普通的Remez交換算法是一樣的。
            (等式40]z(z)=1 (40)圖13A到13D是示出了通過擴展Remez交換算法設計的、使得任何頻率點被通過且預濾波器的頻率響應可以按如下規格考慮的低通濾波器的頻率響應的圖。
            注意,在下面的說明中,用于避免棋盤失真的零點將被作為預濾波器處理。
            預濾波器的頻率響應被如下表示。
            Z(ejω)=(1+e-jω+e-2jω+···+e-jω(U-1))]]> 說明將在下面示出。
            <線性相位FIR濾波器>
            ·24個分支·偶對稱
            ·等式(2)的預濾波器U=3(被調節使得DC增益變為U)<設計方法>
            由經過任何頻率點且考慮了預濾波器頻率響應的Remez交換算法設計。

            指定的頻率點

            圖13A是示出了用分貝顯示的頻率響應的圖,圖13B是示出了將值保持為原樣而被顯示的頻率響應的圖。圖13C是增益3鄰近部分的放大圖,圖13D是增益0鄰近部分的放大圖。
            在圖13A到13D中,虛線表示預濾波器的頻率特性和補償器的頻率特性,實曲線表示最后得到的頻率特性(所提出的H(z)),豎實線表示其中H(z)=0必須成立以避免棋盤失真的頻率(零點)黑點表示帶的截點。
            根據圖13A,可以確認通帶的增益保持常數值,并且通過零點,以避免棋盤失真。
            此外,根據圖13C可以確認所指定的頻率點被通過。
            此外,根據圖13C和13D可以確認保持了等波形。
            也就是說,通過Remez交換算法設計的通過任何頻率點且被放大以使得預濾波器的頻率響應可以被考慮的低通濾波器,可以給出良好的頻率響應特性。
            接下來,將說明具有上述配置的取樣速率轉換器的操作。
            從輸入端101輸入的具有取樣頻率FS的樣本數據x(n)被輸入到LPF 102。
            在LPF 102中,從輸入端101輸入的取樣頻率Fsi的樣本信號的折疊出現被避免(抑制),并且信號被輸出到向上取樣器103。
            在向上取樣器103中,(U-1)個零點被插入到信號之間,取樣頻率Fsi被增加到U倍,并且具有取樣頻率UFsi的樣本信號被輸出到卷積處理單元104。
            卷積處理單元104基于等式(3)進行卷積處理,限制樣本信號的帶,并在下一個階段中將其原樣提供給線性內插塊105。
            這時,卷積輸出(在向下取樣器之前)是通過插入零點由向上取樣器得到的樣本。
            然后,線性內插塊105從卷積處理單元104的輸出信號中選擇取樣頻率UFsi的兩個樣本點。然后,根據線性內插找到所需位置處的值。
            據此,具有取樣頻率Fso的樣本信號y(m)被從輸出端106輸出。
            在具有上述功能的取樣速率轉換器100中,卷積處理單元104的FIR濾波器是線性相位FIR濾波器,其中,脈沖響應變為濾波器系數,傳送函數H(z)與預濾波器的傳送函數Z(z)和補償器的傳送函數K(z)相關聯并且與將要通過的頻率點以及預濾波器的頻率響應相關聯,并且濾波器系數是通過使用經過任何頻率點的Remez交換算法且考慮了預濾波器相對于將要通過的頻率點的頻率響應和預濾波器的頻率響應、基于通過進行相對于期望特性的加權近似而得到的補償器的振幅特性被設定的,因此本發明的取樣速率轉換器具有如下優點。
            也就是說,可以避免棋盤失真。此外,可以考慮任何預濾波器,并且可以通過任何頻率點。
            此外,設置是如下組成的向上取樣器103,用于在信號之間插入(U-1)個零點并增加取樣頻率Fsi到U倍,卷積處理單元104,用于通過具有1/up截止頻率的低通濾波器的卷積(傳送函數H(z))內插值,以及線性內插塊105,用于從具有取樣頻率UFsi的卷積處理單元104的輸出信號中選擇兩個樣本點以及根據線性內插找出在所需位置處的值,因此,如圖表3所示,具有如下優點,即其中截止頻率很嚴格的取樣速率轉換器是可能的。
            表3示出了根據本實施例由取樣速率轉換器支持的頻率帶。
            表3.由本發明的取樣速率轉換支持的頻率

            在表3中,輸入的取樣頻率Fsi是8kHz、11.025kHz、12kHz、22.05kHz、24kHz、32kHz、44.1kHz、48kHz和96kHz,輸出的取樣頻率Fso是8kHz、11.025kHz、12kHz、22.05kHz、24kHz、32kHz、44.1kHz和48kHz。
            在表3中,“△”部分是可以由取樣速率轉換器基于一般的多速率信號處理得到的頻率。一般而言,基于多速率信號處理的那些具有更好的衰減特性,并且在通帶等的控制中更容易,但是通過根據本實施例的取樣速率轉換器也是可能實現的。
            此外,根據第一實施例,LPF 102被設置在輸入端101和向上取樣器103之間,因此,當輸入的取樣頻率Fsi高于輸出的取樣頻率Fso時,重疊部分的出現和折疊的出現可以被抑制。
            <第二實施例>
            圖14是示出了根據本發明的取樣速率轉換器的第二實施例的配置的圖。
            第二實施例與上面說明的第一實施例的區別在于LPF 17被設置在線性內插塊105與輸出端106之間,取代了在輸入端101與向上取樣器103之間設置LPF 102,并且配置被制成使得當輸入的取樣頻率Fsi低于輸出的取樣頻率Fso時會出現成像部分,并且避免(抑制)了在原始數據中不存在的頻率部分的出現。
            配置的其余部分與上面說明的第一實施例的相同。
            根據第二實施例,具有如下優點,可以避免(抑制)在原始數據中不存在的頻率分量的出現,并且以與第一實施例相同的方式,可以避免棋盤失真,并且可以考慮可選的預濾波器,可以通過任意頻率點,并且其中截止頻率很嚴格的取樣速率轉換變為可能。
            注意,還可能在輸入端101與向上取樣器103之間提供LPF 102,在線性內插塊105與輸出端106之間提供LPF 107。
            <第三實施例>
            圖15是示出了根據本發明的取樣速率轉換器的第三實施例的配置的圖。
            注意,在圖15中,“U”和“D”是互為本原的正整數,R(z)指多相濾波器的傳送函數。此外,向上的箭頭表示在信號之間插入(U-1)個零點的向上取樣器。
            第三實施例與上面說明的第一實施例的區別在于FIR濾波器被所謂“多相分解”成U個濾波器,并且使用多相配置。
            這里,將首先給出對多相分解的說明。
            <多相分解>
            取樣速率轉換器通過使用第(N-1)級FIR濾波器(傳送函數H(z))限制帶,并內插零點部分的樣本。
            H(z)=Σn=0N-1h(n)z-n---(42)]]>該FIR濾波器的截止頻率ωc變為[等式43]wc=π/maxU,D (43)圖5的配置可以由多相配置如圖15中所示地等同表示。等式(42)的濾波器的傳送函數H(z)和多相濾波器的傳送函數Ri(z)通過如下關系表示。
            H(z)=RU-1(zU)+z-1RU-2(zU)+···+z-(U-1)Ro(zU)]]>=Σn=0U-1Ri(zU)z(U-1-i)---(44)]]>[等式45]Ri(z)=Σn=0N/U-1h(Un-i-1+U)z-n---(45)]]>注意,計算是如下進行的,將整個倍數增加到U,或者,在對NU個分支不充足時,假設存在0系數。
            在圖15中示出的內插的輸入/輸出關系在將多相濾波器的脈沖響應表示為ri(n)時通過卷積和向上取樣器的處理變為如下。
            y(m)=x0(m-(U-1)U),m=kU+(U-1)···xu-2((m-1)U),m=kU+1xU-1(mU),m=kU---(46)]]>注意,k是整數,xi(m)變為如下 ni(m)=Σn=0∞n(n)x(m-n)---(47)]]>根據第三實施例的取樣速率轉換器200包括,如圖15所示,輸入端201,LPF 202,卷積處理單元203-1到203-U,向上取樣器204-1到204-U,延遲單元205-1到205-U-1,加法器206-1到206-U-1,線性內插塊207,和輸出端208。
            在這些部分中,加法裝置是通過延遲單元205-1到205-U-1和加法器206-1到206-U-1構成的。
            輸入端201接收具有取樣頻率FS的樣本信號x(n)作為輸入。
            當輸入的取樣頻率Fsi高于輸出的取樣頻率FSO時,產生重疊部分,并且出現折疊,因此LPF 202防止(抑制)從輸入端201輸入的具有取樣頻率Fsi的樣本信號x(n)出現折疊,并將結果輸出到卷積處理單元203-1到203-U。
            卷積處理單元203-1到203-U包括通過多相分解基于上述Remez交換算法設計的FIR濾波器而得到的多相濾波器,基于上述等式(47)進行輸入樣本信號與多相濾波器之間的卷積處理,并在下一個階段將處理結果輸出到向上取樣器204-1到204-U。
            向上取樣器204-1到204-U接收卷積處理單元203-1到203-U的輸出樣本信號,插入(U-1)個零點,增加取樣頻率Fsi到U倍,并輸出具有取樣頻率UFsi的樣本信號。
            延遲單元205-1通過向上取樣器204-1將具有取樣頻率UFs的樣本信號延遲預定的時間,并將結果輸出到加法器206-1。
            此外,延遲單元205-1到205-U-1將加法器206-1到206-U-2的輸出信號延遲預定的時間,并將結果輸出到加法器206-2到206-U-1。
            加法器206-1將通過顯示取樣器204-2的具有取樣頻率UFsi的樣本信號與被延遲單元205-1延遲的樣本信號相加,并將結果輸出到延遲單元205-2。
            加法器206-2將通過顯示取樣器204-3的具有取樣頻率UFsi的樣本信號與被延遲單元205-2延遲的樣本信號相加,并將結果輸出到延遲單元205-3。
            以相同的方式,加法器206-U-1將通過顯示取樣器204-U的具有取樣頻率UFsi的樣本信號與被延遲單元205-U-1延遲的樣本信號相加,并將結果輸出到線性內插塊207。
            線性內插塊207從加法器206-U-1的輸出信號中選出兩個樣本點,所述輸出信號即通過將具有其取樣頻率被向上取樣器204-1到204-U-1乘以U的樣本信號相加而得到的信號,根據如上面等式(5)中的線性內插找出位于期望位置(圖7)的值,并將結果作為樣本信號y(m)從輸出端208輸出。
            在取樣速率轉換器200中,在折疊的出現在LPF 202被避免(抑制)之后,從輸入端201輸入的具有頻率Fs的樣本信號x(n)被并行輸入到包括多相濾波器的卷積處理單元203-1到203-U。
            卷積處理單元203-1到203-U進行輸入樣本信號和多相濾波器的卷積處理,并將處理結果在下一個階段中提供給向上取樣器204-1到204-U。
            向上取樣器204-1到204-U在卷積處理單元203-1到203-U的輸出樣本信號之間插入(U-1)個零點,增加取樣頻率FsiU-flod,并輸出具有取樣頻率UFsi的樣本信號。
            向上取樣器204-1到204-U的輸出信號被延遲單元205-1到205-U-1和加法器206-1到206-U-1延遲,然后被累積相加并提供給線性內插塊207。
            然后,線性內插塊207從具有取樣頻率UFsi的加法器206-U-1的輸出信號中選擇兩個樣本點。然后,根據線性內插找到在所需位置的值。
            據此,具有取樣頻率Fso的樣本信號y(m)被從輸出端208輸出。
            根據第三實施例,除了上述第一實施例的效果之外,還有如下優點,即處理量可以被降低到所需的最低限,并且可以實現處理速度的改進。
            <第四實施例>
            圖16是示出了根據本發明取樣速率轉換器的第四實施例的配置的圖。
            第四實施例與上述第三實施例的區別在于,LPF 209被設置在線性內插快207與輸出端208之間,而取代了在輸入端201與卷積處理單元203-1到203-U-1之間設置LPF 202,并且配置被做成使得防止(抑制)了在輸入的取樣頻率Fsi低于輸出的取樣頻率Fso時成像部分的產生以及原始數據中不存在的頻率部分的產生。
            配置的其余部分與第三實施例的相同。
            根據第四實施例,除了上述第二實施例的效果之外,還有如下優點,即處理量可以被降低到所需的最低限,并且可以實現處理速度的改進。
            注意,還可能在輸入端201與卷積處理單元203-1到203-U-1之間設置LPF 202,以及在線性內插快207與輸出端208之間提供LPF209。
            <第五實施例>
            圖17是概念性示出了根據本發明的取樣速率轉換器的第五實施例的圖。
            第五實施例與上述第三實施例的區別在于,線性內插所需的樣本的位置可以在線性內插塊中被指定,因此,配置被制成使得將線性內插快看作是選擇器210,選擇相應于選擇器210輸出的樣本而被設定的多相濾波器,并基于上述等式(47)進行多相的卷積處理。
            通過這種方式,可能將線性內插快看作是選擇器210,選擇相應于輸出樣本的多相濾波器,并且只進行等式(47)的多相卷積處理。
            通過這樣做,就不必進行不必要的計算了。
            在該情況下的取樣速率轉換器200B的配置,如圖17所示,基本上包括,具有取樣頻率Fs的樣本信號x(n)被輸入到的輸入端201,用于抑制當輸入的取樣頻率Fsi高于輸出的取樣頻率Fso時重疊部分產生以及折疊的出現的LPF 202,卷積處理單元203(從1到n),用于進行輸入樣本和通過將濾波器分解成多相而得到的多相濾波器的卷積處理(等式(47)),選擇器210,用于選擇相應于輸出樣本的多相濾波器,以及輸出端208,用于將具有被轉換頻率的樣本信號y(m)輸出到期望的取樣頻率。
            圖18是更具體地示出了圖17中概念性示出的取樣速率轉換器200B的圖。
            如圖18所示,該取樣速率轉換器200C包括,具有取樣頻率Fs的樣本信號x(n)被輸入到其中的輸入端211,具有基于圖11的上述算法被設定系數的FIR濾波器設計單元212,第一存儲器213,用于存儲FIR濾波器設計單元212的多相濾波器系數設定,卷積處理單元214,用于進行LPF的卷積處理以抑制當輸入的取樣頻率Fsi高于輸出的取樣頻率Fso時重疊部分的產生以及折疊的出現,用于存儲LPF系數的第二存儲器215,用于存儲卷積處理單元214的輸入數據的輸入緩沖器216,卷積處理單元217,用于基于輸入緩沖器216中存儲的輸入數據以及從第一系數存儲器213中讀出的系數進行輸入數據和分解成多相的多相濾波器的卷積處理,計數器控制單元218,其作用是選擇輸入樣本所需的兩個樣本點并選擇相應的多相濾波器以及找出線性內插的系數,線性內插處理單元219,用于基于計數器控制單元218的線性內插系數從卷積處理單元217的輸出中找到處于所需位置的值,以及輸出端220,用于輸出頻率被轉換成期望取樣頻率的樣本信號。
            計數器控制單元218包含在線性內插塊中,并且包括變量主計數器(MasterCounter)。
            圖19是用于說明線性內插塊中使用的變量的圖。
            如圖19所示,線性內插塊包括如下變量,MasterCount作為用于確定線性內插的操作的主計數器,Count1000作為將被相加到主計數器的小數點或更小的計數器,是將被相加到主計數器的整數部分的相加值,以及CountAmari作為將被相加到Count1000的小數點或更小的誤差的相加值。
            圖20A和20B是具體說明MasterCount的圖。
            MasterCount如圖20A和20B所示地由至少26位組成,其中,從第0位到第15位的16位表示線性內插的系數,從第16位到第22位的7位表示多相系數設定的數量,從第23位到第25位的3位表示將被讀并被取出的樣本數。
            <第六實施例>
            圖21是概念性示出根據本發明的取樣速率轉換器的第六實施例的圖。
            第六實施例與上述第四實施例的差別在配置上,即線性內插所需的樣本的位置可以被指定在線性內插塊中,因此,線性內插塊被認為是選擇器230,相應于選擇器230的樣本輸出而設定的多相濾波器被選擇,并進行基于上述等式(47)的多相卷積處理。
            在這種方式下,通過將線性內插塊看作是選擇器230并選擇相應于輸出樣本的多相濾波器,僅可進行等式(47)的卷積處理。
            通過這樣做,就不必進行不必要的計算。
            如圖21所示,在該情況中的取樣速率轉換器200D的配置基本上包括,具有取樣頻率Fs的樣本信號x(n)被輸入到其中的輸入端201,卷積處理單元203(從1到n),用于進行輸入樣本和被分解成多相的多相濾波器的卷積處理(等式(47)),選擇器230,用于選擇相應于輸出樣本的多相濾波器,LPF 209,用于抑制當輸入的取樣頻率Fsi低于輸出的取樣頻率Fso時成像部分的產生以及原始數據中不存在的頻率部分的出現,以及輸出端208,用于輸出具有被轉換成期望取樣頻率的頻率的樣本信號y(m)。
            圖22是更具體地示出了在圖21中概念性示出的取樣速率轉換器200D的圖。
            如圖18所示,該取樣速率轉換器200E包括,具有取樣頻率Fs的樣本信號x(n)被輸入到其中的輸入端231,具有基于圖11的上述算法被設定系數的FIR濾波器設計單元232,第一存儲器233,用于存儲FIR濾波器設計單元232的多相濾波器系數設定,用于存儲卷積處理單元235的輸入數據的輸入緩沖器234,卷積處理單元235,用于基于輸入緩沖器234中存儲的輸入數據以及從第一系數存儲器233中讀出的系數進行輸入數據和分解成多相的多相濾波器的卷積處理,計數器控制單元236,其作用是選擇輸出樣本所需的兩個樣本點并選擇相應的多相濾波器以及找出線性內插的系數,線性內插處理單元237,用于基于計數器控制單元236的線性內插系數從卷積處理單元235的輸出中找到處于所需位置的值,卷積處理單元238,用于進行LPF的卷積處理以抑制當輸入的取樣頻率Fsi低于輸出的取樣頻率Fso時成像部分的產生以及原始數據中不存在的頻率部分的出現,用于存儲LPF系數的第二存儲器239,以及輸出端240,用于輸出頻率被轉換成期望取樣頻率的樣本信號。
            注意,計數器控制單元236的本質與參照圖19和圖20A和20B說明的那些相同。
            下面,將依次說明根據實施例的取樣速率轉換器的安裝方法和具體實例。
            這里,將依次說明圖5、圖17和圖18的取樣速率轉換器的安裝方法和具體實例,以及圖14、圖21和圖22的取樣速率轉換器的安裝方法和具體實例。
            <取樣速率轉換器(圖5、圖17和圖18)的安裝方法>
            在實現圖5中所示的取樣速率轉換器時,是通過使用圖17中示出的多相配置實現的。在下面說明的實現方法中,通過選擇相對于輸出的所需輸入以及選擇多相濾波器,所需的最少計算是令人滿意的。
            圖23是說明根據本實施例的取樣速率轉換器的安裝方法的流程圖。
            具體的處理如下。注意,為了便于說明,使用類似C語言的表述,但是安裝方法是硬件、軟件或硬件與軟件的混和是不重要的。
            <步驟10>
            如圖23所示,首先,進行初始設置(F201)。在該初始設置中,輸入數據的數量、輸出數據的數量、輸入數據的中心、多相濾波器的分支數以及多相濾波器的系數設定的數量被設定。具體地說,計數器控制單元215的初始值被設定。
            它們具體如下。
            ·輸入數據數InputSample;·分支數Tap;·輸入數據的中心InputOffset=0;·多相濾波器的分支數PolyTap=(Tap+Up-1)/Up;·多相濾波器的系數設定數Coeff Count;以及·線性內插塊找到CountAdd和CountAmari的方法根據輸入的取樣頻率Fsi和輸出的取樣頻率Fso定義。
            只有整數部分*/CountAdd=(int)(Fso/Fsi*2^23);100小數點或更少*/CountAmari=((Fso/Fsi*2^23)-CountAdd)*1000;當頻率被從44.1kHz轉換成48kHz時,具體如下實例轉換44.1kHz→48kHz*/CountAdd=9130457;CountAmari=687;Count1000的初始值被如下確定。
            Count1000=500;/*for rounding*/<步驟11>
            接下來,低通FIR濾波器通過已經參照例如圖11說明的Remez交換算法設計(F202)。
            這里,將省略其詳細描述。
            <步驟12>
            接下來,準備多相濾波器(F203)。即,根據FIR濾波器的傳送函數H(z),通過使用上述等式(42)的發展等式找到多相濾波器的傳送函數Ri(z)。
            如在下面的等式中,找到每一個多相的標準化系數RiNormal。
            RiNormal=Σk=0PolyTap-1Coef(k)---(48)]]>然后,進行卷積的系數倒置。也就是說,為了進行上述等式(47)的卷積,多相濾波器的系數被重新排列成相反的順序。
            第一多相的系數設定使得分支的中心位置成為第一輸入*/有必要選擇系數設定。*/CoeffCount=(PolyTap*/Up/2)%Up;設定MasterCount的初始值。*/包括MasterCount=H(z)的中間值的多相的系數設定數CoeffCount)<<16注意,<<16表示16位的改變。
            <步驟13>
            接下來,確定多相濾波器和輸入的中心(F204)。
            發現線性內插的系數LinearCord以及多相系數設定數CoefCount如下/*系數設定和線性內插系數根據MasterCount找到*/CoefCount=((MasterCount&0x007fffff)>>16)&0x0000007f;LinearCoef=((MasterCount&0x007fffff)&0x0000ffff);<步驟14>
            進行通過低通濾波器LPF的濾波處理(F209)具體地,在卷積如等式(49)進行并且接下來如等式(50)被舍入之后,其如等式(51)中被標準化,然后如在等式(52)中進行剪輯。
            卷積
            Yconv=Σk=0LPFTap-1LPF(k)×x(InputOffset)+k-LPFTap/2---(49)]]>[等式50]舍入結束ytmp=Yconv+(LPFNormal/2)Yconv≥0Yconv-(LPFNormal/2)Yconv<0---(50)]]>[等式51]標準化YNormal=Ytmp/LPFNormal(51)[等式52]剪輯Src(n)=CLIP_MAXYNormal≥CLIP_MAXCLIP_MINYNormal≤CLIP_MINYNormalelse---(52)]]>在16位的情況下,CLIP MAX=32767和CLIP MIN=32768成立。
            在24位的情況下,CLIP MAX=8388607和CLIP MIN=8388608成立。
            在有符號8位的情況下,CLIP MAX=127和CLIP MIN=128成立。
            在無符號8位的情況下,CLIP MAX=255和CLIP MIN=0成立。
            <步驟15>
            接下來,進行上述等式(47)的卷積(F205)。
            找到線性內插塊中需要的兩個樣本點。在卷積計算中,多相系數設定以及輸入數據的數量通過使用表4中左邊樣本InA和右邊樣本InB的下列值來計算。


            具體地,在如等式(53)進行卷積并且結果如等式(54)中舍入結束之后,結果被如等式(55)中標準化,然后如等式(56)中剪輯。
            卷積Inconv=Σk=0PolyTap-1CoefCount(k)xSrc(InputOffset+(k-PolyTap/2))---(53)]]>[等式54]舍入結束Intap=Inconv+(RNormal/2)xconv≥0Inconv-(RNormal/2)xconv<0---(54)]]>[等式55]標準化InNormal=InTap/RNormal(55)[等式56]剪輯

            在16位的情況下,CLIP MAX=32767和CLIP MIN=32768成立。
            在24位的情況下,CLIP MAX=8388607和CLIP MIN=8388608成立。
            在有符號8位的情況下,CLIP MAX=127和CLIP MIN=128成立。
            在無符號8位的情況下,CLIP MAX=255和CLIP MIN=0成立。
            <步驟16>
            接下來,進行線性內插(F206)。
            如下進行一個實例的線性內插。
            PCM的負最小值24位*/#定義SRC PCM 24BIT MINUS-8388608加上PCM的最大值24位*/#定義SRC PCM 24BIT PLUS 8388607線性Alpha、Beta系數的最大值到128位*/#定義SRC LINEAR MAX COEF 65535線性內插的16位的舍入*/#定義SRC LINEAR NORMAL 32768線性內插的16位的改變量*/#定義SRC LINEAR SHIFT 16進行線性內插的1實例alpha,/*16位輸入*/IN A,/24位輸入*/IN B,/*24位輸入*/OUT Y/*24位輸出*/lldata=alpha*In B+(SRC線性MAX COEF-alpha)*/INA;/*舍入*/if(lldata>=0)lldata+=(SRC線性標準化);else lldata-=(SRC線性標準化);/*標準化*/Out Ya=(lldata>>SRC線性位移);剪輯*/
            Y(n)=(Out Ya<SRC PCM 24BIT MINUS)?SRC PCM 24BITMINUS((Out Ya>SRC PCM 24BIT PLUS)?SRC PCM 24BITPLUSOut Ya);<步驟17>
            更新計數器值和樣本中心(F208)。
            將被讀的MasterCount和樣本數Readsample被如下更新。
            將計數值相加到控制主計數器*/Mas terCount+=CountAdd;小數點的余數或更少*/1000次*/Count1000+=CountAmari;if(Count1000>=1000){Mas terCount+=1;Count1000-=1000;}InputOffset+=((MasterCount)>>23)&0x00000007;MasterCount=MasterCount&0x007fffff;<步驟18>
            這里,進行是否滿足結束條件的判斷(F207)。當條件滿足時,結束處理。當不滿足時,程序返回到F208和F209的處理。
            接下來,將說明圖14、圖21和圖22的取樣速率轉換器的安裝方法和具體實例。
            <取樣速率轉換器(圖14、圖21和圖22)的安裝方法>
            在實現圖14中所示的取樣速率轉換器時,是通過使用圖21中示出的多相配置實現的。在下面說明的實現方法中,通過選擇相對于輸出的所需輸入以及選擇多相濾波器,所需的最少計算是令人滿意的。
            圖24是說明根據本實施例的取樣速率轉換器的安裝方法的流程圖。
            具體的處理如下。注意,為了便于說明,使用類似C語言的表述,但是安裝方法是硬件、軟件或硬件與軟件的混和是不重要的。
            圖24的處理與圖23的相同,只是圖23處理中F209的處理在F206處理之后進行。
            步驟中處理的內容與參照圖23說明的內容相同,因為這里省略了重疊的部分。
            <第七實施例>
            作為第七實施例,將說明每一個上述取樣速率轉換器的音頻裝置。
            圖25的框圖示出了使用根據本發明的取樣速率轉換器的音頻裝置的配置的實例。
            音頻裝置300包括輸入端301,LR分離電路(DSB)302,取樣速率轉換器(SRC)303,衰減器(ATT)304,靜音電路(MUTE)305,以及輸出端306。
            關于取樣速率轉換器(SRC)303,使用如上述第一到第六實施例中參照圖5、圖14、圖17、圖18、圖21和圖22說明的、取樣速率轉換器100、100A、200和200A。
            在這種情況下,關于與上述說明不同的地方,進行如下處理。
            進行結束點處理。也就是說,當進行第一濾波器計算時,缺少Tap/2個樣本。0被補充用于Tap/2個樣本。
            此外,進行先前數據的保持處理。關于音頻特有的處理,有必要在應用濾波器之后相對于寬度如下復制輸入數據用于下一個濾波器計算。
            Src(i)=Src(Width+i),i=0,…,Tap-1 (57)在該音頻裝置300中,例如,從輸入端301輸入的48kHz的PCM數字信號在分離電路302被進行預定的分離,然后取樣速率在取樣速率轉換器303被轉換,并且結果作為例如44.1kHz的信號被輸出。
            然后,在衰減器304進行衰減處理,并且PCM數字音頻信號經由靜音電路305從輸出端306輸出。
            根據音頻裝置300,因為其具有其中避免了棋盤失真、可以考慮任何預濾波器、可以通過任何頻率點、并且處理量可以被降低到所需最低限的取樣速率轉換器,所以具有如下優點,噪聲耐量被改善,線性增益的偏差可以避免,并且可以實現處理速度的改進。
            <第八實施例>
            作為第八實施例,將說明利用取樣速率轉換器的圖像處理裝置。
            在使用取樣速率轉換用于使用線性相位FIR濾波器的濾波處理應用的圖像分辨率轉換中,使用由利用內插器、分米器和線性相位FIR濾波器的多速濾波器組成作為元件技術。
            在這種情況下,進行結束點處理。關于兩端的處理,進行使得零被首先補充、然后輸入被反射(映射)數據、最后保持兩端象素的處理。
            此外,進行初始化。具體地,在將處理進行到圖像末端時,CoefCount=0被初始化。
            根據本圖像處理裝置,因為其具有其中避免了棋盤失真、可以考慮任何預濾波器、可以通過任何頻率點、并且可能將處理量降低到所需最低限的取樣速率轉換器,所以具有如下優點,即圖像不顯示晶格狀態,可以避免線性增益的偏差,并且可以實現處理速度的改進。
            工業實用性根據本發明,可以得到通過任何頻率的振幅特性,并且可以不考慮截止頻率來實現高精度轉換,因此本發明可以應用于例如音頻取樣頻率的轉換以及例如圖像幀的放大或縮小的分辨率轉換。
            權利要求
            1.一種取樣速率轉換器,包括向上取樣器,用于在樣本信號之間插入U-1個零點并增加取樣頻率到U倍;卷積處理單元,包括FIR濾波器,并進行相對于向上取樣器輸出信號的預定卷積處理;以及線性內插塊,用于相對于卷積處理單元的處理結果選擇兩個樣本點,并根據線性內插找到處于所需位置的值,其中,卷積處理單元的FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似來設定。
            2.如權利要求1所述的取樣速率轉換器,其中,濾波器系數基于通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似而得到的補償器的振幅特性而被設定。
            3.如權利要求1所述的取樣速率轉換器,其中,加權近似使用考慮了預濾波器的頻率響應的Remex交換算法相對于期望特性來進行。
            4.如權利要求1所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時防止出現重疊部分和出現折疊的低通濾波器。
            5.如權利要求1所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和以及非原始頻率部分出現的低通濾波器。
            6.一種取樣速率轉換器,包括向上取樣器,用于在樣本信號之間插入U-1個零點并增加取樣頻率到U倍;卷積處理單元,包括FIR濾波器,并進行相對于向上取樣器輸出信號的預定卷積處理;以及線性內插塊,用于相對于卷積處理單元的處理結果選擇兩個樣本點,并根據線性內插找到處于所需位置的值,其中,卷積處理單元的FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數的FIR濾波器,并且濾波器系數通過使用增加了限制條件以通過任何頻率點的算法進行相對于期望特性的加權近似來設定。
            7.如權利要求6所述的取樣速率轉換器,其中,加權近似是使用通過任何頻率點的Remez交換算法相對于期望特性進行的。
            8.如權利要求6所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時防止出現重疊部分和出現折疊的低通濾波器。
            9.如權利要求6所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和以及非原始頻率部分出現的低通濾波器。
            10.一種取樣速率轉換器,包括向上取樣器,用于在樣本信號之間插入U-1個零點并增加取樣頻率到U倍;卷積處理單元,包括FIR濾波器,并進行相對于向上取樣器輸出信號的預定卷積處理;以及線性內插塊,用于相對于卷積處理單元的處理結果選擇兩個樣本點,并根據線性內插找到處于期望位置的值,其中,卷積處理單元的FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率點和預濾波器的頻率響應相關的期望特性的加權近似來設定。
            11.如權利要求10所述的取樣速率轉換器,其中,濾波器系數基于通過進行相對于與將要通過的頻率點和預濾波器的頻率響應相關的期望特性的加權近似而得到的補償器的振幅特性而被設定。
            12.如權利要求10所述的取樣速率轉換器,其中,加權近似使用通過任何頻率點且考慮了預濾波器的頻率響應的Remex交換算法相對于期望特性來進行。
            13.如權利要求10所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時防止出現重疊部分和出現折疊的低通濾波器。
            14.如權利要求10所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和以及非原始頻率部分出現的低通濾波器。
            15.一種取樣速率轉換器,包括多個卷積處理單元,包括通過多相分解預定的FIR濾波器以及進行輸入樣本信號和被分解成多相的多相濾波器的卷積處理得到的預相濾波器,多個向上取樣器,用于在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍,加法裝置,用于在相加所有信號之后通過調整多個向上取樣器輸出信號的傳播時間產生信號,以及線性內插塊,用于相對于加法裝置得到的信號選擇兩個樣本點以及根據線性內插找到所需位置處的值,其中,FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似來設定。
            16.如權利要求15所述的取樣速率轉換器,其中,濾波器系數基于通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似而得到的補償器的振幅特性而被設定。
            17.如權利要求15所述的取樣速率轉換器,其中,加權近似使用考慮了預濾波器的頻率響應的Remex交換算法相對于期望特性來進行。
            18.如權利要求15所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時防止出現重疊部分和出現折疊的低通濾波器。
            19.如權利要求15所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和非原始頻率部分出現的低通濾波器。
            20.一種取樣速率轉換器,包括多個卷積處理單元,包括通過多相分解預定的FIR濾波器以及進行輸入樣本信號和被分解成多相的多相濾波器的卷積處理得到的預相濾波器,多個向上取樣器,用于在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍,加法裝置,用于在相加所有信號之后通過調整多個向上取樣器輸出信號的傳播時間產生信號,以及線性內插塊,用于相對于加法裝置得到的信號選擇兩個樣本點以及根據線性內插找到所需位置處的值,其中,FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數的FIR濾波器,并且濾波器系數是通過使用增加了限制條件以通過任何頻率點的算法進行相對于期望特性的加權近似來設定的。
            21.如權利要求20所述的取樣速率轉換器,其中,加權近似是使用通過任何頻率點的Remez交換算法相對于期望特性進行的。
            22.如權利要求20所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時防止出現重疊部分和出現折疊的低通濾波器。
            23.如權利要求20所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和以及非原始頻率部分出現的低通濾波器。
            24.一種取樣速率轉換器,包括多個卷積處理單元,包括通過多相分解預定的FIR濾波器以及進行輸入樣本信號和被分解成多相的多相濾波器的卷積處理得到的預相濾波器,多個向上取樣器,用于在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍,加法裝置,用于在相加所有信號之后通過調整多個向上取樣器輸出信號的傳播時間產生信號,以及線性內插塊,用于相對于加法裝置得到的信號選擇兩個樣本點以及根據線性內插找到所需位置處的值,其中,FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率點和預濾波器的頻率響應相關的期望特性的加權近似來設定。
            25.如權利要求24所述的取樣速率轉換器,其中,濾波器系數基于通過進行相對于與將要通過的頻率點和預濾波器的頻率響應相關的期望特性的加權近似而得到的補償器的振幅特性而被設定。
            26.如權利要求24所述的取樣速率轉換器,其中,加權近似使用通過任何頻率點且考慮了預濾波器的頻率響應的Remex交換算法相對于期望特性來進行。
            27.如權利要求24所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時避免出現重疊部分和出現折疊的低通濾波器。
            28.如權利要求24所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和非原始頻率部分出現的低通濾波器。
            29.一種取樣速率轉換器,包括卷積處理單元,包括能夠設定通過多相分解預定FIR濾波器而得到的不同濾波器系數的多相濾波器,并進行輸入樣本信號與具有所選系數的多相濾波器的卷積處理,選擇器,用于選擇輸出樣本所需的兩個樣本點,并選擇相應多相濾波器的系數,以及線性內插塊,用于根據線性內插找到處于所需位置的值,其中FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似來設定。
            30.如權利要求29所述的取樣速率轉換器,其中,濾波器系數基于通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似而得到的補償器的振幅特性而被設定。
            31.如權利要求29所述的取樣速率轉換器,其中,加權近似使用考慮了預濾波器的頻率響應的Remex交換算法相對于期望特性來進行。
            32.如權利要求29所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時防止出現重疊部分和出現折疊的低通濾波器。
            33.如權利要求29所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和非原始頻率部分出現的低通濾波器。
            34.如權利要求29所述的取樣速率轉換器,其中,選擇器包括計數器,通過該計數器至少找到線性內插系數、多個多相系數設定以及多個輸入樣本。
            35.一種取樣速率轉換器,包括卷積處理單元,包括能夠設定通過多相分解預定FIR濾波器而得到的不同濾波器系數的多相濾波器,并進行輸入樣本信號與具有所選系數的多相濾波器的卷積處理,選擇器,用于選擇輸出樣本所需的兩個樣本點,并選擇相應多相濾波器的系數,以及線性內插塊,用于根據線性內插找到處于所需位置的值,其中FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數的FIR濾波器,并且濾波器系數通過使用增加了限制條件以通過任何頻率點的算法進行相對于期望特性的加權近似來設定。
            36.如權利要求35所述的取樣速率轉換器,其中,濾波器系數是基于補償器的振幅特性而被設定的,所述振幅特性通過使用增加了限制條件以通過任何頻率點的算法進行相對于期望特性的加權近似得到。
            37.如權利要求35所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時防止出現重疊部分和出現折疊的低通濾波器。
            38.如權利要求35所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和非原始頻率部分出現的低通濾波器。
            39.如權利要求35所述的取樣速率轉換器,其中,選擇器包括計數器,通過該計數器至少找到線性內插系數、多個多相系數設定以及多個輸入樣本。
            40.一種取樣速率轉換器,包括卷積處理單元,包括能夠設定通過多相分解預定FIR濾波器而得到的不同濾波器系數的多相濾波器,并進行輸入樣本信號與具有所選系數的多相濾波器的卷積處理,選擇器,用于選擇輸出樣本所需的兩個樣本點,并選擇相應多相濾波器的系數,以及線性內插塊,用于根據線性內插找到處于所需位置的值,其中,上述FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率點和預濾波器的頻率響應相關的期望特性的加權近似來設定。
            41.如權利要求40所述的取樣速率轉換器,其中,濾波器系數基于通過進行相對于與將要通過的頻率點和預濾波器的頻率響應相關的期望特性的加權近似而得到的補償器的振幅特性而被設定。
            42.如權利要求40所述的取樣速率轉換器,其中,加權近似使用通過任何頻率點且考慮了預濾波器的頻率響應的Remex交換算法相對于期望特性來進行。
            43.如權利要求40所述的取樣速率轉換器,還包括當輸入的取樣頻率低于輸出的取樣頻率時防止出現重疊部分和出現折疊的低通濾波器。
            44.如權利要求40所述的取樣速率轉換器,還包括當輸入的取樣頻率高于輸出的取樣頻率時防止成像部分出現和以及非原始頻率部分出現的低通濾波器。
            45.如權利要求40所述的取樣速率轉換器,其中,選擇器包括計數器,通過該計數器至少找到線性內插系數、多個多相系數設定以及多個輸入樣本。
            46.一種取樣速率轉換方法,包括第一步驟,在樣本信號之間插入U-1個零點并增加取樣頻率到U倍,第二步驟,進行相對于其取樣頻率被包括FIR濾波器的卷積處理單元乘以U的信號的預定卷積處理,在所述FIR濾波器中,脈沖響應由有限時間長度表示,脈沖響應變為濾波器系數,并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯,以及第三步驟,相對于處理結果選擇兩個樣本點,并根據線性內插找到所需位置處的值,其中FIT濾波器的濾波器系數通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似來計算。
            47.一種取樣速率轉換方法,包括第一步驟,在樣本信號之間插入U-1個零點并增加取樣頻率到U倍,第二步驟,進行相對于其取樣頻率被包括FIR濾波器的卷積處理單元乘以U的信號的預定卷積處理,在所述FIR濾波器中,脈沖響應由有限時間長度表示,脈沖響應變為濾波器系數,以及第三步驟,相對于處理結果選擇兩個樣本點,并根據線性內插找到所需位置處的值,其中濾波器系數通過使用增加了限制條件以通過任何頻率點的算法進行相對于期望特性的加權近似來計算。
            48.一種取樣速率轉換方法,包括第一步驟,在樣本信號之間插入U-1個零點并增加取樣頻率到U倍,第二步驟,進行相對于其取樣頻率被包括FIR濾波器的卷積處理單元乘以U的信號的預定卷積處理,在所述FIR濾波器中,脈沖響應由有限時間長度表示,脈沖響應變為濾波器系數,并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯,以及第三步驟,相對于處理結果選擇兩個樣本點,并根據線性內插找到所需位置處的值,其中所述FIR濾波器的濾波器系數通過進行相對于與將要通過的頻率點和預濾波器的頻率響應相關的期望特性的加權近似來計算。
            49.一種取樣速率轉換方法,包括第一步驟,進行輸入樣本信號與被多個卷積處理單元分解成多相的多相濾波器的卷積處理,所述多個卷積處理單元包括通過多相分解預定FIR濾波器而得到的多相濾波器,第二步驟,在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍,第三步驟,調整多個信號的傳播時間使取樣頻率被增加到U倍,并產生通過相加所有信號得到的信號,以及第四步驟,相對于第三步驟得到的信號選擇兩個樣本點,并根據線性內插找到所需位置處的值,其中,FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似來計算。
            50.一種取樣速率轉換方法,包括第一步驟,進行輸入樣本信號與被多個卷積處理單元分解成多相的多相濾波器的卷積處理,所述多個卷積處理單元包括通過多相分解預定FIR濾波器而得到的多相濾波器,第二步驟,在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍,第三步驟,調整多個信號的傳播時間使取樣頻率被增加到U倍,并產生通過相加所有信號得到的信號,以及第四步驟,相對于第三步驟得到的信號選擇兩個樣本點,并根據線性內插找到所需位置處的值,其中,FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數的FIR濾波器,并且濾波器系數通過使用增加了限制條件以通過任何頻率點的算法進行相對于期望特性的加權近似來計算。
            51.一種取樣速率轉換方法,包括第一步驟,進行輸入樣本信號與被多個卷積處理單元分解成多相的多相濾波器的卷積處理,所述多個卷積處理單元包括通過多相分解預定FIR濾波器而得到的多相濾波器,第二步驟,在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍,第三步驟,調整多個信號的傳播時間使取樣頻率被增加到U倍,并產生通過相加所有信號得到的信號,以及第四步驟,相對于第三步驟得到的信號選擇兩個樣本點,并根據線性內插找到所需位置處的值,其中,上述FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率點和預濾波器的頻率響應相關的期望特性的加權近似來計算。
            52.一種取樣速率轉換方法,包括第一步驟,選擇輸出樣本所需的兩個樣本點,并選擇相應多相濾波器的系數,以及第二步驟,通過卷積處理單元進行樣本信號與具有所選系數的多相濾波器的卷積處理,所述卷積處理單元包括通過多相分解預定FIR濾波器得到的多相濾波器并能夠設定不同的濾波器系數,其中上述FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與預濾波器的頻率響應相關的期望特性的加權近似來計算。
            53.一種取樣速率轉換方法,包括第一步驟,選擇輸出樣本所需的兩個樣本點,并選擇相應多相濾波器的系數,以及第二步驟,通過卷積處理單元進行樣本信號與具有所選系數的多相濾波器的卷積處理,所述卷積處理單元包括通過多相分解預定FIR濾波器得到的多相濾波器并能夠設定不同的濾波器系數,其中上述FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數的FIR濾波器,并且濾波器系數通過使用增加了限制條件以通過任何頻率點的算法進行相對于期望特性的加權近似來計算。
            54.一種取樣速率轉換方法,包括第一步驟,選擇輸出樣本所需的兩個樣本點,并選擇相應多相濾波器的系數,以及第二步驟,通過卷積處理單元進行樣本信號與具有所選系數的多相濾波器的卷積處理,所述卷積處理單元包括通過多相分解預定FIR濾波器得到的多相濾波器并能夠設定不同的濾波器系數,其中上述FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和預濾波器的頻率響應相關的期望特性的加權近似來計算。
            55.一種包括取樣速率轉換器的音頻裝置,其中,取樣速率轉換器包括向上取樣器,用于在樣本信號之間插入U-1個零點并增加取樣頻率到U倍;卷積處理單元,包括FIR濾波器,并進行相對于向上取樣器輸出信號的預定卷積處理;以及線性內插塊,用于相對于卷積處理單元的處理結果選擇兩個樣本點,并根據線性內插找到處于期望位置的值,其中,卷積處理單元的FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            56.一種包括取樣速率轉換器的音頻裝置,其中,取樣速率轉換器包括多個卷積處理單元,包括通過多相分解預定的FIR濾波器以及進行輸入樣本信號和被分解成多相的多相濾波器的卷積處理得到的預相濾波器,多個向上取樣器,用于在相應卷積處理單元的輸出信號之間插入U-1個零點并增加取樣頻率到U倍,加法裝置,用于在相加所有信號之后通過調整多個向上取樣器輸出信號的傳播時間產生信號,以及線性內插塊,用于相對于加法裝置得到的信號選擇兩個樣本點以及根據線性內插找到所需位置處的值,其中,上述FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率點和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            57.一種包括取樣速率轉換器的音頻裝置,其中,取樣速率轉換器包括卷積處理單元,包括能夠設定通過多相分解預定FIR濾波器而得到的不同濾波器系數的多相濾波器,并進行輸入樣本信號與具有所選系數的多相濾波器的卷積處理,選擇器,用于選擇輸出樣本所需的兩個樣本點,并選擇相應多相濾波器的系數,其中上述FIR濾波器是其中脈沖響應由有限時間長度表示、脈沖響應變為濾波器系數、并且傳送函數H(z)與預濾波器的傳送函數Z(z)相關聯的FIR濾波器,并且濾波器系數通過進行相對于與將要通過的頻率點和/或預濾波器的頻率響應相關的期望特性的加權近似來設定。
            全文摘要
            提供了一種取樣速率轉換裝置,其能夠得到通過任意頻率的振幅特性從而實現高精確度的轉換而不依賴截止頻率。所述取樣速率轉換裝置包括向上取樣器(103),用于在樣本之間插入(U-1)個零點并增加取樣頻率F
            文檔編號H03H17/06GK1765051SQ200580000109
            公開日2006年4月26日 申請日期2005年1月18日 優先權日2004年1月30日
            發明者茂木幸彥, 西崎譽 申請人:索尼株式會社
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品