JP2007202157A - サンプリングレート変換 - Google Patents

サンプリングレート変換 Download PDF

Info

Publication number
JP2007202157A
JP2007202157A JP2007014341A JP2007014341A JP2007202157A JP 2007202157 A JP2007202157 A JP 2007202157A JP 2007014341 A JP2007014341 A JP 2007014341A JP 2007014341 A JP2007014341 A JP 2007014341A JP 2007202157 A JP2007202157 A JP 2007202157A
Authority
JP
Japan
Prior art keywords
signal
sampling rate
time
upsampled
sampling
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2007014341A
Other languages
English (en)
Other versions
JP5541832B2 (ja
Inventor
Gerhard Uwe Schmidt
ウベ シュミット ゲルハルト
Mohamed Krini
クリニ モハメド
Martin Roessler
ロースラー マルティン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harman Becker Automotive Systems GmbH
Original Assignee
Harman Becker Automotive Systems GmbH
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harman Becker Automotive Systems GmbH filed Critical Harman Becker Automotive Systems GmbH
Publication of JP2007202157A publication Critical patent/JP2007202157A/ja
Application granted granted Critical
Publication of JP5541832B2 publication Critical patent/JP5541832B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/06Non-recursive filters
    • H03H17/0621Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing
    • H03H17/0635Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing characterized by the ratio between the input-sampling and output-delivery frequencies
    • H03H17/0642Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing characterized by the ratio between the input-sampling and output-delivery frequencies the ratio being arbitrary or irrational
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/06Non-recursive filters
    • H03H17/0621Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing
    • H03H17/0635Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing characterized by the ratio between the input-sampling and output-delivery frequencies
    • H03H17/0685Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing characterized by the ratio between the input-sampling and output-delivery frequencies the ratio being rational

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Transmission Systems Not Characterized By The Medium Used For Transmission (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

【課題】第1のサンプリングレートの入力信号を処理するための方法を提供すること。
【解決手段】本発明は、例えば、第1のサンプリングレートの入力信号を処理するための本発明の方法であって、入力信号をアップサンプリングして第2の信号を得るステップ;第2信号を低域通過フィルタリングして第3の信号を得るステップ;第2の信号を、時間遅延し、その後、時間遅延された第2の信号を低域通過フィルタリングして第4の信号を得るステップ;第1の重み因子によって、第3の信号を重み付けして第5の信号を得るステップ;第1の重み因子とは異なる第2の重み因子によって、第4の信号を重み付けして第6の信号を得るステップ;および、第1のサンプリングレートとは異なる第2のサンプリングレートの出力信号を得るために、第5の信号および第6の信号を加算するステップ;を包含する、方法を提供する。
【選択図】図1

Description

本発明は、帯域制限された信号の処理に関し、特に、第1のサンプリングレートにてサンプリングされた電気信号を、第1のサンプリング信号とは異なる第2のサンプリングレートを有する電気信号に変換することに関する。より詳細には、本発明は、一時的に変化する目的のサンプリングレートの場合における、サンプリングレート変換に関する。
サンプリングレートの変換は、信号処理における主問題である。2つの電子デバイスが電気信号の交換によって通信する場合、通常、一つのデバイスのサンプリングレートが、もう一方のデバイスのサンプリングレートとは異なるという問題が生じる。例えば、第1のデバイスは、第1のサンプリングレートを用いて信号を検出し得、これらの信号を処理し得る。処理後、これらの信号は、第2のデバイスに出力される。第2のデバイスは、第1のデバイスからくる入力信号を処理し、第1のデバイスのサンプリングレートよりも高いサンプリングレートにてそれらを出力することが意図され得る。この場合、一部のアップサンプリングは第2のデバイスによって実行される必要がある。そのような通信の一例は、Bluetooth(登録商標)技術を介した、セル式電話とハンズフリーセットとの間における信号の交換である。
電気信号の一部のアップサンプリングのニーズに関する別の例は、電話との関連において、帯域制限された信号の拡張によって与えられる。通常の電話の帯域制限されたオーディオ(特に、音声信号)は、300Hz〜3.4kHzのみの帯域幅を示す。より低い周波数およびより高い周波数を有する生の信号の一部を取り除くことによって、音声品質における低下(特に、了解度の低下)が生じるゆえに、制限された帯域幅をより広い帯域幅に拡張することがさらに所望される。補間によってサンプリングレートを増加すると、通常は、拡張されたより不自然でない信号を生成するために、帯域制限された信号を分析する前に、一部の処理を実行する必要がある。
所定のサンプリングレート(すなわち、離散的な信号を生成するために、連続した信号から単位時間ごとに取得された個々のサンプルの数)の増加は、離散的な時間サンプル間における信号の補間を必要とする。周知の補間方法は、Lagrange補間を含み、それは、n+1個の所定のデータポイントを通過するn次の多項式を見出す古典的手法である。三次スプラインとして知られている別の手法は、データポイントのうちの一つにおいて所定の傾きを達成するように、2つのポイントを介する3次の多項式にフィッティングする。
特に、一時的に変化する目的のサンプリングレートを用いたサンプリングレート変換に関連して、以下の2つの方法が過去に用いられてきた。第1の方法に従うと、極めてオーバーサンプリングされた信号が入力信号から生成される。それぞれの目的の離散的サンプルに対して、最も近似している人工的に生成された時間サンプルが選択される。そのオーバーサンプリングレートの半分は、最大誤差を表す。明らかであるように、オーバーサンプリングレートがより高く選択されるほど、それに対応する誤差は低くなる。しかしながら、この強引な方法は、高度な計算およびメモリ容量を必要とする。
第2の方法は、連続的で完全積分可能な信号x(t)の個々のサンプルx(n)の擬連続補間を示し、ここで、tは秒単位の時間であり、nは整数の範囲である。本方法に従い、シンク(sine cardinal)関数が、所望のサンプリング時間Tにおいて、信号値を生成するために用いられる。
Figure 2007202157
これらの式において、Tinはオリジナル(入力)サンプリング時間(期間)を表し、finはサンプリングレート(fin=1/Tin)を表す。
目的のサンプリングレートがオリジナルのサンプリングレートよりも小さい場合、低域通過のカットオフは、エイリアシングを避けるために、新たなより低いサンプリングレートの半分以下に置かれなければならない。
しかしながら、シンク補間は、目的のサンプリングレートの恣意的な選択に対して、比較的高品質を保証するが、それはまた、多量の計算時間を必要とする。シンク関数の詳細な補間に対して、多量のメモリリソースが必要とされる。さらに、実際には、再生信号の質は、式1における有限和のニーズによって制限される。
従って、上記の欠点を克服し、かつ従来技術に比べてコンピュータリソースに対してより小さな要求でサンプリングレートを変換する方法を提供するための課題が本発明に内在する。
上記の課題は、項目1に従った方法によって解決される。第1のサンプリングレートの入力信号を処理するための本発明の方法は、アップサンプリングされた信号を得るために該入力信号をアップサンプリングするステップ;第1のフィルタリングされたアップサンプリングされた信号を得るために、該アップサンプリングされた信号を低域通過フィルタリングするステップ;第2のフィルタリングされたアップサンプリングされた信号を得るために、該アップサンプリングされた信号を、時間遅延し、その後、該時間遅延されたアップサンプリングされた信号を低域通過フィルタリングするステップ;第1の重み付けされた信号を得るために、第1の重み因子によって、第1のフィルタリングされたアップサンプリングされた信号を重み付けするステップ;第2の重み付けされた信号を得るために、該第1の重み因子とは異なる第2の重み因子によって、第2のフィルタリングされたアップサンプリングされた信号を重み付けするステップ;および、第1のサンプリングレートとは異なる第2のサンプリングレートの第1の出力信号を得るために、該第1の重み付けされた信号および該第2の重み付けられた信号を加算するステップ;を包含する。
第1のサンプリングレートの入力信号は、所定のサンプリングレート(すなわち、離散的なデジタル信号を生成するために連続的な信号から取得される時間あたりの個々のサンプル数)を用いて類似信号からサンプリングされるデジタル信号であり得る。入力信号は、第1の出力信号よりも、多くのまたは少ない、1つの帯域間隔に対するサンプリングインスタント、従って、個々のサンプルを含み得る。
本発明の方法は、比較的小さいコンピュータリソース(特に、メモリ)のみを必要とする。従って、本方法は、ハードウェアインプリメンテーションに対して非常に実用的である。要するに、本方法は、単に、低次の1つの畳み込みと、少ないさらなる演算を含む。
アップサンプリングされた信号と、時間遅延されたアップサンプリングされた信号との両方の低域通過フィルタリングは、好ましくは、有限インパルス応答(FIR)フィルタによって実行され、それによって、そのようにしなれば生じ得るエイリアシング成分を削除する。
該時間遅延されたアップサンプリングされた信号の該低域通過フィルタリングは、アップサンプリングされた信号の低域通過フィルタリングに対して同じであるが信号の時間遅延に従って遅延されているインパルス応答を使用して、実行され得る。有利なことに、時間遅延するステップは、一サンプリングインスタント(一つの離散的なサンプル)の時間遅延という結果となる。この場合、時間遅延されたアップサンプリングされた信号の低域通過フィルタリングは、遅延されていない信号をフィルタリングするために使用されるインパルス応答であるが一サンプリングインスタントだけ遅延ているインパルス応答によって実行される。
第1の重み因子および第2の重み因子は、時間に依存し得る。第2の重み因子の値は、1から該第1の重み因子の値を引いたものによって与えられ得る。従って、y(n)およびy(n−1)がそれぞれ、第1のフィルタリングされたアップサンプリングされた信号、および第2のフィルタリングされたアップサンプリングされた信号を示し、a(n)が第1の重み因子を示しているとすると、第1の出力信号v(n)は、v(n)=a(n)y(n)+(1−a(n))y(n−1)によって得られる。
アップサンプリングされた信号が、該第1のサンプリングレートの入力信号のそれぞれのサンプリングインスタントに対する個々のサンプル間に、L−1個のヌルを挿入すること(すなわち、mod(n,L)≠0の場合、ヌルを挿入することであって、nは複数の整数の範囲に亘り、かつ離散的な時間インデックス(すなわち、サンプリングインスタント)を表す、こと)であって、Lは、第1の信号をアップサンプリングするのに用いる因子である、ことと、
Figure 2007202157
に従った低域通過FIRフィルタhTP,iのインパルス応答を用いた該入力信号の拡張畳み込みとによって得られた場合、該低域通過フィルタリングは、非常に効果的に実行され得る。
式3においては、パラメータNは、低域通過フィルタの長さ(フィルタ係数の数)を示す。
この場合において、第2のサンプリングレートの第1の出力信号は、
Figure 2007202157
に従った該低域通過フィルタリングのために使用される時間依存のフィルタ係数hTP(n)を用いて得られ得る。
ここで、上付きインデックスTは転置演算を示し、
Figure 2007202157
である。
通常、入力信号の一部のアップサンプリングが望まれた場合、アップサンプリング後の第1の出力信号は、最終的に望まれるサンプリングレートより上のサンプリングレートを示す。従って、第3のサンプリングレートであって、該第2のサンプリングレートとは異なり、かつ通常第1のサンプリングレートとも異なる、第3のサンプリングレートの第2の(宛先)出力信号(v(n))を得るために、該第2のサンプリングレートの該第1の出力信号(v(n))がダウンサンプリングされ得る。
この場合において、該第1の重み因子は、
Figure 2007202157
によって得られ得、
が、
Figure 2007202157
として定義され、δはクロネッカーのデルタ関数であり、Lは、該第2の(宛先)出力信号の該第3のサンプリングレートfoutへと、該第1の(オリジナル)サンプリングレートfinの該入力(オリジナル)信号をアップサンプリングするのに用いる因子である。記号
Figure 2007202157
は、それぞれ、それよりも大きなもののうち最も近い整数、それよりも小さなもののうちの最も近い整数に丸めることを示す。
本発明の方法の上記局面は、入力信号が変換される時間依存の第3のサンプリングレートに対して有用になるように修正もされ得る。第3のサンプリングレートが時間に依存する場合、第1のサンプリングレートの入力信号は第1のバッファ(特に、第1のリングバッファ)に格納され得、および/または第3のサンプリングレートの第2の出力信号は第2のバッファ(特に、第2のリングバッファ)に格納され得る。それぞれの書取および読取ポインタをモニタリングすることによって、複数のバッファのうちの1つが空になるか、またはオーバーフロー状態に近いかが決定され得る。これらの場合においては、第3のサンプリングレートは、入力信号の誤りのない処理を保証するために、補正され得る。
時間依存の第3のサンプリングレートは、
Figure 2007202157
に従って調整され得る時間依存の補正項△fout(n)とともに一定の名目上の値を含み得る。ここで、△は、一サンプリングインスタントから次のサンプリングインスタントへのサンプリングレートの最大の変化である。
式9に与えられる補正項△fout(n)を利用して、第3のサンプリングレートの第2の出力信号v(n)は、ダウンサンプリング因子
Figure 2007202157
を用いて
Figure 2007202157
によって得られ得る。
時間依存の第3の(宛先)サンプリングレートの場合におけるサンプリングレート変換に有用な本発明の方法の実施形態は、第1の(オリジナル)サンプリングレートが時間に依存している場合にも、利用され得る。
本発明は、信号処理手段であって、本発明の方法の上述された例のうちの1つに従った処理ステップを実行することによって、第1のサンプリングレートとは異なる第2のサンプリングレートの第1の出力信号を、第1のサンプリングレートの入力信号から生成するように構成されている信号処理手段、および/または、第3のサンプリングレートであって、該第1のサンプリングレートおよび/または該第2のサンプリングレートとは異なる、第3のサンプリングレートの第2の出力信号を、該第1のサンプリングレートの入力信号から生成するように構成されている信号処理手段も提供する。
更に、電子デバイスであって、上記の信号処理手段を含んでおり、他の電子デバイスに信号を送信し、および/または他の電子デバイスから信号を受信するように構成されている、電子デバイスが提供される。電子デバイスの例としては、セル式電話、ヘッドセット、およびハンズフリーセットが挙げられる。
更に、本発明は、コンピュータプログラム製品であって、上述されたようなサンプリングレート変換のための本発明の方法の実施形態のステップを実行するためのコンピュータ実行可能命令を有している1つ以上のコンピュータ可読媒体を備えている、コンピュータプログラム製品を提供する。
本発明の更なる特徴および利点は、図面を参照して説明される。以下の説明においては、本発明の好ましい実施形態を示すことを意図する添付の図に対して参照が行われる。そのような実施形態が、添付の特許請求の範囲によって規定される本発明の範囲全体を表さないことが理解される。
本発明は、さらに、以下の手段を提供する。
(項目1)
第1のサンプリングレートの入力信号(x(n))を処理する方法であって、
アップサンプリングされた信号(x(n))を得るために該入力信号(x(n))をアップサンプリングするステップ(1)と、
第1のフィルタリングされたアップサンプリングされた信号(y(n))を得るために該アップサンプリングされた信号(X(n))を低域通過フィルタリングするステップ(2)と、
第2のフィルタリングされたアップサンプリングされた信号(y(n−1))を得るために、該アップサンプリングされた信号(x(n))を、特に、一サンプリングインスタントだけ時間遅延し(3)、その後、該時間遅延されたアップサンプリングされた信号(x(n−1))を低域通過フィルタリングするステップ(4)と、
第1の重み付けされた信号を得るために、第1の重み因子によって、第1のフィルタリングされたアップサンプリングされた信号(y(n))を重み付けするステップと、
第2の重み付けされた信号を得るために、該第1の重み因子とは異なる第2の重み因子によって、第2のフィルタリングされたアップサンプリングされた信号(y(n−1))を重み付けするステップと、
該第1のサンプリングレートとは異なる第2のサンプリングレートの第1の出力信号(v(n))を得るために、該第1の重み付けされた信号および該第2の重み付けられた信号を加算するステップと
を包含する、方法。
(項目2)
上記アップサンプリングされた信号および上記時間遅延されたアップサンプリングされた信号の両方(x(n),x(n−1))の上記低域通過フィルタリングするステップ(2,4)が、有限インパルス応答フィルタを用いて実行される、項目1に記載の方法。
(項目3)
上記時間遅延されたアップサンプリングされた信号(x(n−1))が、時間遅延されていない上記アップサンプリングされた信号(x(n))に対して一サンプリングインスタントだけ遅延されており、該時間遅延されたアップサンプリングされた信号(x(n−1))の上記低域通過フィルタリングするステップ(4)は、時間遅延されていない該アップサンプリングされた信号(x(n))の上記低域通過フィルタリングするステップ(2)に対して使用されるものと同じであるが一サンプリングインスタントだけ遅延されているインパルス応答を使用して実行される、項目1または2に記載の方法。
(項目4)
上記第1の重み因子および上記第2の重み因子が、時間に依存するか一定かであり、該第2の重み因子の値が、1から該第1の重み因子の値を引いたものによって与えられる、項目1〜3のいずれか一項に記載の方法。
(項目5)
上記アップサンプリングされた信号(x(n))が、上記第1のサンプリングレートの入力信号の個々のサンプル間にL−1個のヌルを挿入することによって得られ、Lは、該入力信号(x(n))をアップサンプリングするのに用いる因子であり、上記第1のフィルタリングされたアップサンプリングされた信号(y(n))が、
Figure 2007202157
に従った有限インパルス応答フィルタのインパルス応答(hTP)を用いた該入力信号(x(n))の拡張畳み込みによって得られ、Lは、該入力信号(x(n))をアップサンプリングするのに用いる因子であり、Nは、上記低域通過フィルタのフィルタ長を示す、項目2〜4のいずれか一項に記載の方法。
(項目6)
上記第2のサンプリングレートの上記第1の出力信号(v(n))が、
Figure 2007202157
を用いた
Figure 2007202157
に従った上記低域通過フィルタリングのために使用される時間依存のフィルタ係数(hTP)を用いて得られ、a(n)は、上記第1の重み因子であり、上記上付きインデックスTは、ベクトルまたはマトリックスの転置を示す、項目5に記載の方法。
(項目7)
上記第2のサンプリングレートとは異なる第3のサンプリングレートの第2の出力信号(v(n))を得るために、該第2のサンプリングレートの上記第1の出力信号(v(n))がダウンサンプリングされる(6)、項目1〜6のいずれか一項に記載の方法。
(項目8)
上記第1の重み因子が、
Figure 2007202157
によって与えられ、
が、
Figure 2007202157
として規定され、δはクロネッカーのデルタ関数であり、Lは、上記第1のサンプリングレートfinの上記入力信号(x(n))をアップサンプリングするのに用いる因子であり、foutは上記第2の出力信号(v(n))の該第3のサンプリングレートであり、記号
Figure 2007202157
は、それぞれ、それよりも大きなもののうち最も近い整数、それよりも小さなもののうちの最も近い整数に丸めることを示す、項目7に記載の方法。
(項目9)
上記第3のサンプリングレートが時間に依存し、上記入力信号(x(n))が第1のバッファに格納され、および/または上記第2の出力信号(v(n))が第2のバッファに格納される、項目7または8に記載の方法。
(項目10)
上記第3のサンプリングレートが、
Figure 2007202157
によって与えられる時間依存の寄与と一定の寄与とを含み、上記第2の出力信号(v(n))が、
Figure 2007202157
を用いて
Figure 2007202157
によって得られる、項目9に記載の方法。
(項目11)
上記第1のサンプリングレートが時間に依存する、項目8〜10のいずれか一項に記載の方法。
(項目12)
信号処理手段であって、項目1〜6のいずれか一項に記載の処理ステップを実行することによって、上記第1のサンプリングレートとは異なる第2のサンプリングレートの第1の出力信号(v(n))を、第1のサンプリングレートの入力信号(x(n))から生成するように構成されている、信号処理手段。
(項目13)
項目12に記載の信号処理手段であって、項目7〜11のいずれか一項に記載の処理ステップを実行することによって、第3のサンプリングレートであって、上記第1のサンプリングレートおよび/または上記第2のサンプリングレートとは異なる、第3のサンプリングレートの第2の出力信号(v(n))を、該第1のサンプリングレートの入力信号(x(n))から生成するように構成されている、項目12に記載の信号処理手段。
(項目14)
特に、セル式電話またはヘッドセットまたはハンズフリーセットである電子デバイスであって、項目12または13に記載の信号処理手段を備えている、電子デバイス。
(項目15)
コンピュータプログラム製品であって、項目1〜11のうちのいずれか一項に記載の方法のステップを実行するためのコンピュータ実行可能命令を有している1つ以上のコンピュータ可読媒体を備えている、コンピュータプログラム製品。
(摘要)
本発明は、第1のサンプリングレートの入力信号を処理するための本発明の方法であって、アップサンプリングされた信号を得るために入力信号をアップサンプリングするステップ;第1のフィルタリングされたアップサンプリングされた信号を得るために、アップサンプリングされた信号を低域通過フィルタリングするステップ;第2のフィルタリングされたアップサンプリングされた信号を得るために、アップサンプリングされた信号を、時間遅延し、その後、時間遅延されたアップサンプリングされた信号を低域通過フィルタリングするステップ;第1の重み付けされた信号を得るために、第1の重み因子によって、第1のフィルタリングされたアップサンプリングされた信号を重み付けするステップ;第2の重み付けされた信号を得るために、第1の重み因子とは異なる第2の重み因子によって、第2のフィルタリングされたアップサンプリングされた信号を重み付けするステップ;および、第1のサンプリングレートとは異なる第2のサンプリングレートの出力信号を得るために、第1の重み付けされた信号および第2の重み付けられた信号を加算するステップ;を包含する、方法に関する。
本発明により、上記の欠点を克服し、かつ従来技術に比べてコンピュータリソースに対してより小さな要求でサンプリングレートを変換する方法が提供され得る。
図1は、オリジナルの(入力)サンプリングレートfinの入力信号x(n)から、目的の(出力)サンプリングレートfoutの出力信号v(n)(ここでnは離散的な時間インデックスを意味する)を得るために、サンプリングレート変換の本発明の方法の一例を示す。まず、入力信号x(n)は、Lの因子によって、サンプリングレートfup>finにまでアップサンプリングされ(参照番号1)、ここでLは整数であり、すなわち、fup=Lfinである。アップサンプリングされた信号x(n)は、オリジナルの信号のサンプル間において、L−1個のヌルを挿入することによって得られる。すなわち、mod(n,L)=0の場合、x(n)=x(n/L)であり(modはモジュロ関数)、それ以外の場合は、x(n)=0。
他方、x(n)は、一つの処理ブランチにおいて、係数hTP,iを有する低域通過フィルタのインパルス応答を用いた畳み込み(参照番号2)の対象となる。
Figure 2007202157
他方、x(n)は、別の処理ブランチにおいて、一つのサンプリングインスタントによって、時間遅延され(参照番号3)、その後、その時間遅延された信号は、第1の処理ブランチに用いられるのと同じであるが一つのサンプリングインスタントによって遅延されたフィルタベクトルを有するインパルス応答を用いた畳み込み(参照番号4)の対象となる。エイリアシングを防ぐために、FIRフィルタは低域通過フィルタリングプロセスに用いられる。
さらに、x(n)は、ヌルを挿入することによって生成されているため、そのフィルタリングされた非遅延のアップサンプリングされた信号y(n)は、以下の式によって得られうる。
Figure 2007202157
(n)はゼロになることのない値に対してのみ用いられるので、
Figure 2007202157
である。
ここで、記号
Figure 2007202157
は、それぞれ、それよりも大きなもののうちの最も近い整数、それよりも小さなもののうちの最も近い整数に丸めることを意味する。
制御手段5は、第1の処理ブランチに対して時間依存因子a(n)、および第2の処理ブランチに対して時間依存因子1−a(n)を提供し、y(n)およびy(n−1)の重み付け合計によって、入力信号のうちの一つよりもはるかに高いサンプリングレートの出力信号v(n)を得る:
(n)=a(n)y(n)+(1−a(n))y(n−1)
である。
重み付けプロセスおよび畳み込みプロセスは、
Figure 2007202157
の信号ベクトルと、
Figure 2007202157
という用いられる低域通過フィルタの時間依存フィルタベクトルとを用いて、
(n)=x(n) TP(n)
によって実行され得る。
基本的には、処理手段5によるa(n)の適切な選択によって、それぞれの任意のサンプリングインスタントに対する信号値(すなわち、信号の個々のサンプルが取られる離散的な時間点)は、少なくとも近似的に計算され得る。
出力信号v(n)のサンプリングレートfupは、通常、所望のサンプリングレートfoutよりも高く、したがって、v(n)はダウンサンプリングされなければならない。特に、a(n)およびダウンサンプリング因子R(n)の適切な選択によって、オリジナルの信号の各サンプリングインスタントは、以下のように、所望のサンプリングレートの一つのサンプリングインスタントに正確に対応することが確保され得る
Figure 2007202157
ここで、クロネッカーのデルタ関数δを用い、f
Figure 2007202157
として定義されている。
したがって、入力信号x(n)よりも高いが、v(n)のサンプリングレートよりも低いサンプリングレートを有する所望の出力信号v(n)は、結果として、
Figure 2007202157
として得られる。
図2は、整数mによって示されるような11,025kHzのサンプリングレートの目的の信号(下の行)への、整数nによって示されるような8kHzのオリジナルのサンプリングレートを有する信号(上の行)のサンプリングレート変換の一例を示す。8個のサンプリングインスタント(n=0〜7)が、比較可能な帯域間隔において、10個のサンプリングインスタント(m=0〜9)に対応するオリジナルの信号x(n)に対して示されている。オリジナルの信号は、重み付け因子a(m)および1−a(m)(第2の行)と、時間依存のダウンサンプリングR(m)とを用いて変換される。アップサンプリングに対する因子は、理解しやすくするためにL=1に設定されており、第2の行のより良い明瞭さのために、重み付け因子は、互いに僅かにシフトされている。
図3は、サンプリングレート変換の本発明の方法のための別の例を示し、目的のサンプリングレートは時間に依存する。図3に示されるように、オリジナルの(入力)信号
Figure 2007202157
および目的の(出力)信号v(n)は、リングバッファにおいてバッファされる。あるいは、オリジナルの信号と目的の信号とのうちの一つのみがバッファに入力される。書き込み/読み出しポインタをモニタリングすることによって、バッファのうちの一つが空になるか、またはオーバーフローのリスクを負うか否かが決定され得る。バッファのうちの一つの空状態を避けるために、またはオーバーフロー状態を避けるために、目的のサンプリングレートの補正は以下のように実行され得る。目的のレートfoutは、定数部分
Figure 2007202157
および時変(time−varying)補正
Figure 2007202157
に分割され、ここで、
Figure 2007202157
であり、Δは、一つのサンプリングインスタントから次のサンプリングインスタントへのサンプリングレートの最大変化であり、例えば、10−7〜0.01Hzの範囲内であり得る。
補正項であるΔfout(n)を含んでいる時間依存の目的のサンプリングレートを用いて、図1に関連して記載された計算は、目的の信号を得るために実行される。その結果、
Figure 2007202157
のダウンサンプリング因子を用いて、
Figure 2007202157
が得られる。
さらに、デジタル入力信号x(n)が、連続した類似の信号からサンプリングされ得る時変入力サンプリングレートに対して、上述の信号処理がまた、実行され得る。
以上のように、本発明の好ましい実施形態を用いて本発明を例示してきたが、本発明は、この実施形態に限定して解釈されるべきものではない。本発明は、特許請求の範囲によってのみその範囲が解釈されるべきであることが理解される。当業者は、本発明の具体的な好ましい実施形態の記載から、本発明の記載および技術常識に基づいて等価な範囲を実施することができることが理解される。
図1は、アップサンプリングされたオリジナルの信号と、時間遅延されたアップサンプリングされたオリジナルの信号の低域通過フィルタリング、および、低域通過フィルタリングされた信号の重み付けを含む、本発明に従ったサンプリングレート変換の方法の一例の原理を示すブロック図である。 図2は、オリジナルの信号、時間依存の重み付け、時間依存のサンプリングレート・デシメーション、および目的の信号を示す。 図3は、時間依存の目的のサンプリングレートの場合に対して、本発明に従ったサンプリングレート変換の一例を示す。この例は、オリジナルの信号と目的の信号とのバッファリングを用いる。

Claims (15)

  1. 第1のサンプリングレートの入力信号(x(n))を処理する方法であって、
    アップサンプリングされた信号(x(n))を得るために該入力信号(x(n))をアップサンプリングするステップ(1)と、
    第1のフィルタリングされたアップサンプリングされた信号(y(n))を得るために該アップサンプリングされた信号(X(n))を低域通過フィルタリングするステップ(2)と、
    第2のフィルタリングされたアップサンプリングされた信号(y(n−1))を得るために、該アップサンプリングされた信号(x(n))を、特に、一サンプリングインスタントだけ時間遅延し(3)、その後、該時間遅延されたアップサンプリングされた信号(x(n−1))を低域通過フィルタリングするステップ(4)と、
    第1の重み付けされた信号を得るために、第1の重み因子によって、第1のフィルタリングされたアップサンプリングされた信号(y(n))を重み付けするステップと、
    第2の重み付けされた信号を得るために、該第1の重み因子とは異なる第2の重み因子によって、第2のフィルタリングされたアップサンプリングされた信号(y(n−1))を重み付けするステップと、
    該第1のサンプリングレートとは異なる第2のサンプリングレートの第1の出力信号(v(n))を得るために、該第1の重み付けされた信号および該第2の重み付けられた信号を加算するステップと
    を包含する、方法。
  2. 前記アップサンプリングされた信号および前記時間遅延されたアップサンプリングされた信号の両方(x(n),x(n−1))の前記低域通過フィルタリングするステップ(2,4)が、有限インパルス応答フィルタを用いて実行される、請求項1に記載の方法。
  3. 前記時間遅延されたアップサンプリングされた信号(x(n−1))が、時間遅延されていない前記アップサンプリングされた信号(x(n))に対して一サンプリングインスタントだけ遅延されており、該時間遅延されたアップサンプリングされた信号(x(n−1))の前記低域通過フィルタリングするステップ(4)は、時間遅延されていない該アップサンプリングされた信号(x(n))の前記低域通過フィルタリングするステップ(2)に対して使用されるものと同じであるが一サンプリングインスタントだけ遅延されているインパルス応答を使用して実行される、請求項1または2に記載の方法。
  4. 前記第1の重み因子および前記第2の重み因子が、時間に依存するか一定かであり、該第2の重み因子の値が、1から該第1の重み因子の値を引いたものによって与えられる、請求項1〜3のいずれか一項に記載の方法。
  5. 前記アップサンプリングされた信号(x(n))が、前記第1のサンプリングレートの入力信号の個々のサンプル間にL−1個のヌルを挿入することによって得られ、Lは、該入力信号(x(n))をアップサンプリングするのに用いる因子であり、前記第1のフィルタリングされたアップサンプリングされた信号(y(n))が、
    Figure 2007202157
    に従った有限インパルス応答フィルタのインパルス応答(hTP)を用いた該入力信号(x(n))の拡張畳み込みによって得られ、Lは、該入力信号(x(n))をアップサンプリングするのに用いる因子であり、Nは、前記低域通過フィルタのフィルタ長を示す、請求項2〜4のいずれか一項に記載の方法。
  6. 前記第2のサンプリングレートの前記第1の出力信号(v(n))が、
    Figure 2007202157
    を用いた
    Figure 2007202157
    に従った前記低域通過フィルタリングのために使用される時間依存のフィルタ係数(hTP)を用いて得られ、a(n)は、前記第1の重み因子であり、前記上付きインデックスTは、ベクトルまたはマトリックスの転置を示す、請求項5に記載の方法。
  7. 前記第2のサンプリングレートとは異なる第3のサンプリングレートの第2の出力信号(v(n))を得るために、該第2のサンプリングレートの前記第1の出力信号(v(n))がダウンサンプリングされる(6)、請求項1〜6のいずれか一項に記載の方法。
  8. 前記第1の重み因子が、
    Figure 2007202157
    によって与えられ、
    が、
    Figure 2007202157
    として規定され、δはクロネッカーのデルタ関数であり、Lは、前記第1のサンプリングレートfinの前記入力信号(x(n))をアップサンプリングするのに用いる因子であり、foutは前記第2の出力信号(v(n))の該第3のサンプリングレートであり、記号
    Figure 2007202157
    は、それぞれ、それよりも大きなもののうち最も近い整数、それよりも小さなもののうちの最も近い整数に丸めることを示す、請求項7に記載の方法。
  9. 前記第3のサンプリングレートが時間に依存し、前記入力信号(x(n))が第1のバッファに格納され、および/または前記第2の出力信号(v(n))が第2のバッファに格納される、請求項7または8に記載の方法。
  10. 前記第3のサンプリングレートが、
    Figure 2007202157
    によって与えられる時間依存の寄与と一定の寄与とを含み、前記第2の出力信号(v(n))が、
    Figure 2007202157
    を用いて
    Figure 2007202157
    によって得られる、請求項9に記載の方法。
  11. 前記第1のサンプリングレートが時間に依存する、請求項8〜10のいずれか一項に記載の方法。
  12. 信号処理手段であって、請求項1〜6のいずれか一項に記載の処理ステップを実行することによって、前記第1のサンプリングレートとは異なる第2のサンプリングレートの第1の出力信号(v(n))を、第1のサンプリングレートの入力信号(x(n))から生成するように構成されている、信号処理手段。
  13. 請求項12に記載の信号処理手段であって、請求項7〜11のいずれか一項に記載の処理ステップを実行することによって、第3のサンプリングレートであって、前記第1のサンプリングレートおよび/または前記第2のサンプリングレートとは異なる、第3のサンプリングレートの第2の出力信号(v(n))を、該第1のサンプリングレートの入力信号(x(n))から生成するように構成されている、請求項12に記載の信号処理手段。
  14. 特に、セル式電話またはヘッドセットまたはハンズフリーセットである電子デバイスであって、請求項12または13に記載の信号処理手段を備えている、電子デバイス。
  15. コンピュータプログラム製品であって、請求項1〜11のうちのいずれか一項に記載の方法のステップを実行するためのコンピュータ実行可能命令を有している1つ以上のコンピュータ可読媒体を備えている、コンピュータプログラム製品。
JP2007014341A 2006-01-24 2007-01-24 サンプリングレート変換 Active JP5541832B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP06001452A EP1811663B1 (en) 2006-01-24 2006-01-24 Sampling rate conversion
EP06001452.9 2006-01-24

Publications (2)

Publication Number Publication Date
JP2007202157A true JP2007202157A (ja) 2007-08-09
JP5541832B2 JP5541832B2 (ja) 2014-07-09

Family

ID=36579680

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007014341A Active JP5541832B2 (ja) 2006-01-24 2007-01-24 サンプリングレート変換

Country Status (4)

Country Link
US (1) US7554465B2 (ja)
EP (1) EP1811663B1 (ja)
JP (1) JP5541832B2 (ja)
AT (1) ATE524875T1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9621130B2 (en) * 2011-10-27 2017-04-11 Avago Technologies General Ip (Singapore) Pte. Ltd. Configurable generic filter hardware block and methods
US8947993B2 (en) * 2012-07-24 2015-02-03 Telefonaktiebolaget L M Ericsson (Publ) Methods of transmitting using filtering in the time domain and related devices
WO2016108648A1 (ko) * 2014-12-30 2016-07-07 주식회사 쏠리드 분산 안테나 시스템의 노드 유닛 및 신호 처리 방법
US10341042B2 (en) * 2015-07-27 2019-07-02 Qualcomm Incorporated Narrowband operation with reduced sampling rate
EP3716479A1 (en) 2019-03-26 2020-09-30 Bang & Olufsen A/S A method for sampling rate conversion
US20230015697A1 (en) * 2021-07-13 2023-01-19 Citrix Systems, Inc. Application programming interface (api) authorization

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6477327A (en) * 1987-09-18 1989-03-23 Toshiba Corp Sample rate converting circuit
JP2000031947A (ja) * 1998-07-09 2000-01-28 Sony Corp 標本化周波数変換装置およびこれを具備した電子機器
JP2004282612A (ja) * 2003-03-18 2004-10-07 Sony Corp データ処理装置およびその方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6427157B1 (en) * 1998-07-31 2002-07-30 Texas Instruments Incorporated Fir filter structure with time- varying coefficients and filtering method for digital data scaling
US6711214B1 (en) * 2000-04-07 2004-03-23 Adc Broadband Wireless Group, Inc. Reduced bandwidth transmitter method and apparatus
DE602004017120D1 (de) * 2004-02-24 2008-11-27 Accent S P A Verfahren zur Implementierung eines fraktionalen Abtastratenumsetzers und entsprechende Umsetzerarchitektur

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6477327A (en) * 1987-09-18 1989-03-23 Toshiba Corp Sample rate converting circuit
JP2000031947A (ja) * 1998-07-09 2000-01-28 Sony Corp 標本化周波数変換装置およびこれを具備した電子機器
JP2004282612A (ja) * 2003-03-18 2004-10-07 Sony Corp データ処理装置およびその方法

Also Published As

Publication number Publication date
JP5541832B2 (ja) 2014-07-09
EP1811663B1 (en) 2011-09-14
US20080151984A1 (en) 2008-06-26
EP1811663A1 (en) 2007-07-25
US7554465B2 (en) 2009-06-30
ATE524875T1 (de) 2011-09-15

Similar Documents

Publication Publication Date Title
EP1639703B1 (en) Rational sample rate conversion
JP5541832B2 (ja) サンプリングレート変換
JP5638787B2 (ja) サブバンド信号処理
US8849883B2 (en) Method and apparatus for adaptive control of the decimation ratio in asynchronous sample rate converters
JP2006197599A (ja) 非同期サンプルレートコンバータ
JP2005287041A (ja) サンプリングレート変換システム、サンプリングレート変換方法及びサンプリングレート変換装置
KR101102410B1 (ko) 샘플링 레이트 컨버터, 변환 방법 및 컴퓨터 판독 가능한 기록 매체치
JP2002271174A (ja) サンプリングレート変換装置及び方法
JP2006345508A (ja) デジタル信号のサンプリング周波数を変換するための方法および装置
US20080315940A1 (en) Mixing device for plural digital data having different sampling rates
CN109976660B (zh) 基于线性插值的任意信号采样率重构方法和数据采样系统
EP3704796B1 (en) Low delay decimator and interpolator filters
JP4449007B2 (ja) サンプリング周波数変換装置
US8949303B2 (en) Filter
JPH07235861A (ja) 加重平均を使用したサンプリング周波数変換方法
JP2000165664A (ja) 画像の解像度変換装置及び方法
JP4364598B2 (ja) フィルタ処理装置,フィルタ処理方法及びそのプログラム
JP2007067578A (ja) オーディオ信号処理装置及び方法
JP2008219560A (ja) デシメーションフィルタ
JP7333423B2 (ja) サンプリングレート変換器、サンプリングレート変換方法及び通信システム
CN111585543A (zh) 一种Farrow结构实现音频采样率转换的方法
JP2008089791A (ja) オーディオ信号処理装置
JP5774475B2 (ja) インパルス応答の抽出合成によるsrcの装置及び方法
JP2008109279A (ja) オーディオ信号処理装置及びオーディオ信号処理方法
JP4276258B2 (ja) 速度とフレキシビリティを最大にするためにハードウエアとソフトウエアを用いてサンプルレートコンバータを実装するシステム及び方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20091217

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20101209

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20110819

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120605

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20120904

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20120907

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20120928

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20121003

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20121030

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20121102

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20121203

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130625

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130814

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20131120

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140217

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140415

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140502

R150 Certificate of patent or registration of utility model

Ref document number: 5541832

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250