JPH05145428A - デジタル信号の符号化−復号化方法およびその装置 - Google Patents
デジタル信号の符号化−復号化方法およびその装置Info
- Publication number
- JPH05145428A JPH05145428A JP4107053A JP10705392A JPH05145428A JP H05145428 A JPH05145428 A JP H05145428A JP 4107053 A JP4107053 A JP 4107053A JP 10705392 A JP10705392 A JP 10705392A JP H05145428 A JPH05145428 A JP H05145428A
- Authority
- JP
- Japan
- Prior art keywords
- supplied
- input terminal
- multiplication
- calculating
- digital signal
- 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
Links
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
- H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
- H03M7/3053—Block-companding PCM systems
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/147—Discrete orthonormal transforms, e.g. discrete cosine transform, discrete sine transform, and variations therefrom, e.g. modified discrete cosine transform, integer transforms approximating the discrete cosine transform
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Discrete Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Complex Calculations (AREA)
Abstract
手順を簡易化でき、リアルタイム処理に適したコストお
よび信頼性を得る。 【構成】 連続する複数のディジタルサンプルを重み付
けする窓付け操作を行った後、Cの式で偶数階数の変換
係数Y2kを計算するために、サンプルの複数のブロッ
クについて修正ディスクリートコサイン変換を実行す
る。次に、k=0,・・・,N/4に対してDに示す可
逆複素変換の形として変換係数Y2kを表し、Eに示す
補助計算方程式に基づいて可逆複素変換を計算して、複
数のディジタルサンプルで表されるディジタル信号を符
号化する。
Description
ジタル信号のビットレートの減少、または記憶に要求さ
れる容量を減少するための上記デジタル信号の符号化−
復号化方法およびその装置に関する。
−復号化技術の1つは、周波数変換コーディングであ
る。このタイプの変換器では、デジタル信号は、初期
に、時間−周波数変換が必要であり、符号化され、かつ
伝送される周波数変換の係数である。同一の伝送される
信号の品質に対して、伝送速度における援助は、デコー
ディングにおいて実行される逆変換とともに、コンコー
ディングに用いられる変換の選択に依存する。したがっ
て、時間領域別名キャンセル(time domain aliasingca
ncellation;TDAC,すなわち、時間領域における除
去歪)に基づく修正ディスクリートコサイン変換(MD
CT)は、特に、ディスクリートフーリエ変換(DF
T)またはディスクリートコサイン変換(DCT)のよ
うな従来の変換と比較して効果的である。これらのタイ
プの処理に関係する動作の比較のために、参考文献とし
て、J.P.プリンセン(Princen)およびA.B.ブ
ラッドリ(Bradley)による”時間領域別名キャンセルに
基づく分析/合成フィルタバンク設計(Analysis/synth
esis filter bank design based on time domain alias
ingcancellation)”,IEEE Transacti
ons on ASSP,vol.34,pp.153
−1161,10月,1986年がある。
較すると、この処理および符号化の技術はなおさら強く
みえる。さらに、上述した符号化技術は、比較的、低圧
縮比(すなわち、データ速度が高速度ビットに代り12
8kbit/sに減少させられる高品質のデジタルオー
ディオ信号)にとっては良くないようであるが、その使
用は、さらなる高圧縮比においては避けられないように
みえる。しかしながら、そのような状況の下では、上記
圧縮動作は、複雑になる。リアルタイムで動作すべくM
DCTを用いたコーディングアルゴリズムを実現するこ
と、すなわち、コストおよび信頼性を満足する条件下で
実行することは、前述した変換に対する高速コーディン
グおよび算出方法を確立するための本質であり、準備さ
れたコーディングプロセッサを実現するために、取付け
および実行だけでなく、上述したプロセッサのリアルタ
イム動作を容易にするためである。
られた上記方法において、変換係数は、 (1)… y(m,k)=Σxm(n)・h(N-1-n)・cos(2π・(k+1/2)(n+N0)/N) n=0,...,N-1 によって定義される。なお、k=0,.....,N−
1である。この式において、Nは変換ブロックにおける
サンプルx(n)の大きさまたは数である。h(n)は
空間または時間重付けによって得られるブロック重付け
窓である。また、mは変換ブロックの数であり、xm
(n)は考慮中のm番目のブロックにおけるn番目のサ
ンプル信号を示す。位相シフトN0は原信号の完全な再
構成のために必要とされる。
め、周波数領域における独立係数y(m,k)の数はN
/2に等しい。原オーディオ信号の再構成または合成に
必要とされる逆変換MDCT-1は、元のシーケンスを提
出できない。周波数領域における不完全なサンプリング
は、「アライジング(aliasing)」として知られる時間領
域における歪を生じさせる。それにもかかわらず、合成
は、連続し重複したブロック変換の助けによって表すこ
とができる。ブロックmoの信号合成のために、再構成
されたオーディオ信号xrmo(n)のn番目のサンプ
ルは、 xrmo(n)=f(n+N/2)・xmo-1(n+N/2)+f(n)・xmo(n) で書かれ、n=0,...,N/2−1となる。
o(n)は変換係数y(mo,k)の逆変換である。すなわち、xmo
(n)=MDCT-1y(mo,k)であり、 (2)… Xm(n)=(1/N)Σy(m,k)・cos(2π・(k+1/2)(n+N0)/N) n=0,...,N-1 となる。原信号は正確に再構成される。すなわち、分解
窓h(n)、合成窓f(n)および位相シフト項N0に
基づく条件を負わすことによって、xrmo(n)=xmo(n)と
なる。
ックまでのアライジング項は同一の単位および反対の符
号を有し、ゆえに、取消される。分解、合成窓h(n)およ
びf(n)が、同一長、すなわち同一で、N/2の回りで対
称であるならば、以下の式を満足する。 (3)… f(n+N/2)・h(n)-f(n)・h(N/2+n)=0 (4)… f2(n+N/2)+f(n)=2 なお、n=0,...,N/2−1である。 (5)… no=N/4+1/2 例えば、次式の分解および合成窓は、50%インターブ
ロックを有し、かつ、上記条件を満足する。 (6)… h(n)=f(n)=√2・sin(π・(n+1/2)/N) なお、n=0,...,N−1である。
れるべき係数当たりの有効ビットレートを減少させるこ
となく、周波数領域において上述した高性能な窓を用い
ることができるということにある。上記窓が組込まれた
MDCTは、信号成分の完全な分離を確実なものとし、
ゆえに、よりよいエネルギ集中を提供する。スペクトラ
ム表現における上記改善点は、広帯域における主要な減
少のため、知覚現象を考慮することを容易にする。この
スペクトラム分割は、例えば、矩形の窓が用いられる場
合ではない。
(1)式に定義されたMDCTは、次式に書かれる上記
MDCTを用いて、同一サイズの高速フーリエ変換(F
FT)によって計算される。 (7)… y(m,k)= Re[ej2・pai・(k+1/2)no/N・Σxm(n)・h(N-1-n)・ej2・pai・(k+1/2)n] pai=π,n=0,...,N-1 ここで、Re[...]は、括弧の間の複雑な表現の実際の部
分を明確に示すと、 (8)… y(m,k)=Re[ej2 (k+1/2)no/N・DEF(xm(n)・h(N-1-n)・ej2 n/N] となる。
的な応用は、次の演算操作数を要求するということを認
める。x(n)・ej n/Nの演算、すなわち実数と虚数とのN
個の乗算と、長さNのFFTと、長さNの2つの複雑な
ベクトルの項毎の乗算の実数部の演算とが要求される。
「Implementation of split Radix FFT algorithms f
or complex, real andreal symmetric data」,ASSPに基
づくIEEE発行,34巻,285-295頁,1986年4月号,な
どに述べられているN=2pおよびもっとも効率的な演
算アルゴリズムのために、(8)式の算出は、窓処理に
必要とされる空間−時間重付け演算を無視して、2N+N・
(p-3+4)+2N個の乗算および3N(p-1)+4+N個の加算、すな
わち、合計、N(4p−1)+8個の演算を必要とす
る。変換ブロックのサイズが大きいとき、すなわち、デ
ジタルオーディオ信号を符号化−デコーディングするた
め、N=1024であるとき、上記演算の数は、リアル
タイムにおける使用のために、符号化−複号化法のMD
CTの使用を不利な立場に置くことができる。
もので、本願発明のデジタルサンプルx(n)によって
表されるデジタル信号(x)を符号化−複号化するため
の方法および装置の目的は、上述した欠点を解決するこ
とにある。すなわち、本願発明の他の目的は、符号化に
必要とされる演算数を実際上、減少させる、デジタル信
号を符号化−複号化するための方法および装置を提供す
ることにある。また、本願発明の他の目的は、処理手順
を簡易化する理由で、リアルタイム処理に適合したコス
トおよび信頼性の条件の下で、上記処理手順を、提供さ
れたプロセッサ内に設けられることを可能とした、デジ
タル信号の符号化−複号化方法およびその装置を提供す
ることにある。
るために、請求項1記載の発明では、h(N−1−n)
が決定された重み付け法則を表すときの以下に示す式に
よって表されるブロック内のそれぞれのサンプルの位置
nの関数として、連続するN個のサンプルからなる複数
のブロックにひとまとめにされた複数のディジタルサン
プルx(n)を空間的あるいは、時間的に重み付けする
窓付け操作を行い、 yn=x(n)・h(N−1−n) YN−k=−Yk−1としたときのk∈〔0,・・・,
N/2−1〕に対して、位相シフト項が前記変換係数の
階数kの関数として分布された数式で表される偶数階数
の変換係数Y2kを計算するために、サンプルの複数の
ブロックについて修正ディスクリートコサイン変換を実
行し、 y’n=x2n・h2n、y’’n=(xN−2n−
1)・(hN−2nー1)、W4N=cos(2π/4
N)+jsin(2π/4N)としたときのk=0,・
・・,N/4に対して可逆複素変換の形として前記計算
された偶数階数の変換係数Y2kを表し、 n=0,・・・,N/4−1に対してzn=(y2n−y
N/2−1−2n)+j(yN-1−2n+yN+2+2n)
であるときの補助計算方程式に基づいて前記可逆複素変
換を計算して、前記複数のディジタルサンプルx(n)
で表されるディジタル信号(x)を符号化することを特
徴とする。
ディジタル信号の符号化方法において、前記可逆複素変
換は、 n=0,・・・,N/4−1に対して前記znを計算す
る第1のステップと、 n=0,・・・,N/4−1に対してz’n=zn・W
Nnを計算する第2のステップと、サイズN/4の逆高
速フーリエ変換FFT-1によって以下に示す式で表され
る対応する係数Z’kを計算する第3のステップと、 Z’k=FFT-1・z’n k=0,・・・,N/4−1に対して以下に示す方程式
を計算する第4のステップと、 Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k Zkの虚数部Im〔Zk〕=Y2kと、Zkの実数部R
e〔Zk〕=(Y2k+N/4)とからなる偶数階数の
変換係数を連続的に決定する第5のステップと、前記デ
ィジタル信号(x)を構成する前記複数のディジタルサ
ンプルx(n)の代わりに、前記偶数階数の変換係数Y
2kおよび(Y2k+N/4)を符号化された信号の形
で伝送する第6のステップとからなる連続するステップ
により計算されることを特徴とする。
ジタル信号の符号化方法において、前記補助計算方程式
は、Zk=(-1)k+1・W8-1・WNk・Σ(zn・W4N(4n+1)・WN/4nk)で
表され(n=0〜N/4−1)、前記第2のステップ
と、前記第4のステップとが交換され、W4N(4n+1)に
よる虚数乗法が前記計算の最初あるいは、最後に設けら
れていることを特徴とする。
係数の実数部(YN/4+2k)および虚数部Y2kを
交換し、Re〔Zk〕およびIm〔Zk〕の形式のフォ
ワードMDCTの複素数部を計算することによって前記
交換された変換係数を符号化処理し、前記複素数部の実
数部および虚数部を交換し、前記複数のディジタルサン
プルx(n)の再構成をするためにオーバーラップ計算
を実行して、請求項1記載のディジタル信号の符号化方
法によって符号化された前記複数のディジタルサンプル
x(n)で表されるディジタル信号(x)を復号化する
ことを特徴とする。
ディジタル信号の復号化方法において、前記フォワード
MDCTの複素数部を計算することによる符号化処理と
同一の前記符号化処理は、 n=0,・・・,N/4−1に対してz’n=zn・W
Nnを計算する第2のステップと、サイズN/4の逆高
速フーリエ変換FFT-1によって以下に示す式で表され
る係数Z’kを計算する第3のステップと、 Z’k=FFT-1・z’n k=0,・・・,N/4−1に対して以下に示す方程式
を計算する第4のステップと Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k からなることを特徴とする。
ディジタル信号の復号化方法において、可逆複素変換を
形成する一連の前記第2ないし第4のステップととも
に、前記第2ないし第4のステップによって構成された
逆変換は、逆MDCTのコアを構成し、前記逆MDCT
は、逆入力信号を生成するために、以下に示す式により
入力信号の実数部と虚数部とを交換する第7のステップ
と、 z=a+jb→z*=b+ja 出力信号Y*kを生成するために、逆入力信号のフォー
ワードMDCTを計算する第8のステップと、逆出力信
号Ykを生成するために、出力信号Y*kの実数部と虚
数部とを交換する第9のステップとからなる対応するフ
ォーワードMDCTによって計算されることを特徴とす
る。
サンプルからなる複数のブロックにひとまとめにされた
複数のディジタルサンプルx(n)の連続するN/2個
のサンプルを記憶する第1の記憶手段と、前記複数のデ
ィジタルサンプルx(n)の連続するN/2個のサンプ
ルをN/2遅延して記憶する第2の記憶手段と、 h(N−1−n)が決定された重み付け法則を表すとき
の以下に示す式によって表される複数のサンプルynの
空間的重み付けを実行することからなる前記複数のディ
ジタルサンプルx(n)の窓付け操作を行い、 yn=x(n)・h(N−1−n) YN−k=−Yk−1としたときのk∈〔0,・・・,
N/2−1〕に対して、位相シフト項が前記変換係数の
階数kの関数として分布された偶数階数の変換係数Y2
kを計算するために、サンプルの複数のブロックについ
て修正ディスクリートコサイン変換を実行し、 n=0,・・・,N/4−1に対してzn=(y2n−y
N/2−1−2n)+j(yN-1−2n+yN+2+2n)
であるときの補助計算方程式を用いた可逆複素変換(Y
2k+N/4+jY2k)から偶数階数の変換係数Y2
kを得る第1の計算手段とを具備し、前記複数のディジ
タルサンプルx(n)で表されるディジタル信号(x)
を符号化することを特徴とする。
ディジタル信号の符号化装置において、前記第1の計算
手段は、前記znを計算する第1の複素数計算手段と、 n=0,・・・,N/4−1に対してz’n=zn・W
Nnを計算する第2の複素数計算手段と、(Z’k=F
FT-1・z’n)を計算する第1の逆高速フーリエ変換
計算手段と、 k=0,・・・,N/4−1に対して以下に示す方程式
を計算し、 Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k Zkの虚数部Im〔Zk〕=Y2kと、Zkの実数部R
e〔Zk〕=(Y2k+N/4)とからなる偶数階数の
変換係数を連続的に決定する第3の複素数計算手段とが
縦続接続されていることを特徴とする。
ディジタル信号の符号化装置において、前記第1の複素
数計算手段は、前記第1および第2の記憶手段に接続さ
れ、重み付けされたサンプル値y2nが前記第1の記憶
手段から第1の入力端に直接供給され、重み付けされた
サンプル値(yN/2−1−2n)が(−1)を乗算す
る乗算器を介して第2の入力端に供給される第1の加算
手段と、重み付けされたサンプル値(yN/2+2n)
が前記第1の記憶手段から第1の入力端に供給され、重
み付けされたサンプル値(yN−1−2n)が前記第2
の記憶手段から第2の入力端に供給される第2の加算手
段と、それぞれ記憶容量N/4を有し、前記第1および
第2の加算手段のそれぞれの出力端にそれぞれ接続さ
れ、zn=(y2n−yN/2−1−2n)+j(yN-1−2
n+yN+2+2n)の実数部および虚数部をそれぞれ
記憶する第1および第2の補助記憶手段とを具備するこ
とを特徴とする。
のディジタル信号の符号化装置において、前記第2の複
素数計算手段は、定数(WNn=cos(2πn/N)
+jsin(2πn/N))を記憶する第1の定数記憶
手段と、(zn=a+jb)の実数部aが第1の入力端
に供給され、前記定数WNnの実数部(cos(2πn
/N))が第2の入力端に供給され、積(a・cos
(2πn/N)を出力する第1の乗算手段と、(zn=
a+jb)の実数部aが第1の入力端に供給され、前記
定数WNnの虚数部(sin(2πn/N))が第2の
入力端に供給され、積(a・sin(2πn/N))を
出力する第2の乗算手段と、(zn=a+jb)の虚数
部bが第1の入力端に供給され、前記WNnの実数部
(cos(2πn/N))が第2の入力端に供給され、
積(b・cos(2πn/N))を出力する第3の乗算
手段と、(zn=a+jb)の虚数部bが第1の入力端
に供給され、前記WNnの虚数部(sin(2πn/
N))が第2の入力端に供給され、積(b・sin(2
πn/N))を出力する第4の乗算手段と、前記積(a
・cos(2πn/N))が前記第1の乗算手段から第
1の入力端に供給され、前記積(b・sin(2πn/
N))が前記第4の乗算手段から(−1)を乗算する乗
算器を介して第2の入力端に供給され、a’=a・co
s(2πn/N)−b・sin(2πn/N)、b’=
a・sin(2πn/N)+b・cos(2πn/
N))としたときのa’を出力する第3の加算手段と、
前記積(a・sin(2πn/N))が前記第2の乗算
手段から第1の入力端に供給され、前記積(b・cos
(2πn/N))が前記第3の乗算手段から第2の入力
端に供給され、前記(z’n=a’+jb’)の虚数部
b’を出力する第4の加算手段と、前記z’nの実数部
a’および虚数部b’をそれぞれ記憶する第3および第
4の補助記憶手段とを具備することを特徴とする。
のディジタル信号の符号化装置において、前記第3の複
素数計算手段は、定数W8-1およびW4N4k+1の実数部
および虚数部を記憶する第2の定数記憶手段と、式Z’
kの実数部が第1の入力端に供給され、前記定数WNn
およびW4N4k+1の実数部が前記第2の定数記憶手段か
ら第2の入力端に供給され、対応する第1の積を出力す
る第5の乗算手段と、式Z’kの実数部が第1の入力端
に供給され、前記定数WNnおよびW4N4k+1の虚数部
が前記第2の定数記憶手段から第2の入力端に供給さ
れ、対応する第2の積を出力する第6の乗算手段と、式
Z’kの虚数部が第1の入力端に供給され、前記定数W
NnおよびW4N4k+1の実数部が前記第2の定数記憶手
段から第2の入力端に供給され、対応する第3の積を出
力する第7の乗算手段と、式Z’kの虚数部が第1の入
力端に供給され、前記定数WNnおよびW4N4k+1の虚
数部が前記第2の定数記憶手段から第2の入力端に供給
され、対応する第4の積を出力する第8の乗算手段と、
前記第1の積が前記第5の乗算手段から第1の入力端に
供給され、前記第4の積が前記第8の乗算手段から(−
1)を乗算する乗算器を介して第2の入力端に供給さ
れ、式Zkの実数部Re〔Zk〕=YN/4+2kを出
力する第5の加算手段と、前記第2の積が前記第2の乗
算手段から第1の入力端に供給され、前記第3の積が前
記第7の乗算手段から第2の入力端に供給され、前記式
Zkの虚数部Im〔Zk〕=Y2kを出力する第6の加
算手段とを具備することを特徴とする。
換係数の実数部Re〔Zk〕=YN/4+2kおよび虚
数部Im〔Zk〕=Y2kを記憶する変換係数記憶手段
と、該変換係数記憶手段に記憶された前記伝送された変
換係数を交換して読み出す読出手段と、前記z’nを計
算する第4の複素数計算手段と、逆高速フーリエ変換
(FFT-1・z’n)を計算する第2の逆高速フーリエ
変換計算手段と、以下に示す方程式を計算する方程式計
算手段とが縦続接続され、フォワードMDCTの複素数
部を計算する第2の計算手段と Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k を具備し、連続するN個のサンプルからなる複数のブロ
ックにひとまとめにされ、請求項7記載のディジタル信
号の符号化装置によって符号化された複数のディジタル
サンプルx(n)で表されるディジタル信号(x)を復
号化することを特徴とする。
載のディジタル信号の復号化装置において、前記第2の
計算手段は、請求項7ないし請求項11のいずれかに記
載の符号化装置の複数のモジュールにそれぞれ等しい複
数のモジュールから構成されていることを特徴とする。
((2π・(4n+1))/4N)、s=sin((2
π・(4n+1)/N)としたときの偶数のサンプルx
2nおよび奇数のサンプル(xN/2−1−2n)に基
づいて、以下に示す方程式を虚数乗法によって計算する
第3の計算手段と、 z’n=((x2n・c−xN/2−1−2n・s) +jzo-N/2・(x2n・s+xN/2−2n−1・c))・(c+js ) 式Z’(N/8−k)=FFT-1・z’nを計算するた
めに、前記式z’nに対して逆高速フーリエ変換FFT
-1を実行する逆高速フーリエ変換計算手段と、変換の奇
数係数Y2kおよび偶数係数(YN/2−2K−1)を
計算するために、虚数乗法を実行する第4の計算手段と
を具備し、連続するN個のサンプルからなる複数のブロ
ックにひとまとめにされた複数のディジタルサンプルx
(n)で表されるディジタル信号(x)を符号化および
復号化することを特徴とする。
載のディジタル信号の符号化および復号化装置におい
て、前記第3の計算手段は、偶数サンプル(xN/2−
1/2n)および奇数サンプルx2Nをそれぞれ記憶す
る第3および第4の記憶手段と、前記値cおよび前記値
sをそれぞれ記憶する第5および第6の記憶手段と、前
記偶数サンプルx2Nが第1の入力端に供給され、前記
値cが(−1)を乗算する乗算器を介して第2の入力端
に供給される第9の乗算手段と、前記奇数サンプル(x
N/2−1/2n)が第1の入力端に供給され、前記値
cが直接第2の入力端に供給される第10の乗算手段
と、前記偶数サンプルx2Nが第1の入力端に供給さ
れ、前記値sが前記第6の記憶手段から直接第2の入力
端に供給される第11の乗算手段と、前記奇数サンプル
(xN/2−1/2n)が第1の入力端に供給され、前
記値sが前記第6の記憶手段から(−1)を乗算する乗
算器を介して第2の入力端に供給される第12の乗算手
段と、前記第9の乗算手段から出力される積と前記第1
1の乗算手段から出力される積とが供給される第7の加
算手段と、前記第10の乗算手段から出力される積と前
記第12の乗算手段から出力される積とが供給される第
8の加算手段と、前記第7の加算手段から出力される加
算結果が直接第1の入力端に供給され、前記加算結果が
ディレイラインおよびインバータを経て第2の入力端に
供給される第9の加算手段と、前記第8の加算手段から
出力される加算結果が直接第1の入力端に供給され、前
記加算結果がディレイラインを経て第2の入力端に供給
される第10の加算手段と、前記奇数サンプル(xN/
2−1/2n)がインバータを経て第1の入力端に供給
され、前記第9の加算手段から出力される加算結果が直
接第2の入力端に供給され、前記式z’nの実数部を出
力する第11の加算手段と、前記偶数サンプルx2Nが
インバータを経て第1の入力端に供給され、前記第10
の加算手段から出力される加算結果が直接第2の入力端
に供給され、前記式z’nの虚数部を出力する第12の
加算手段と、前記式z’nの実数部および虚数部をそれ
ぞれ記憶するRAMタイプの第7および第8の記憶手段
とを具備することを特徴とする。
載のディジタル信号の符号化および復号化装置におい
て、前記第3の計算手段は、前記値cおよび前記値sを
それぞれ記憶する第9および第10の記憶手段と、前記
式z’nに対する前記逆高速フーリエ変換FFT-1・
z’nの演算結果の実数部Re〔Z’N/8−k〕が第
1の入力端に供給され、前記値cが第2の入力端に供給
される第13の乗算手段と、前記式z’nに対する前記
逆高速フーリエ変換FFT-1・z’nの演算結果の虚数
部Im〔Z’N/8−k〕が第1の入力端に供給され、
前記値cが第2の入力端に供給される第14の乗算手段
と、前記式z’nに対する前記逆高速フーリエ変換FF
T-1・z’nの演算結果の実数部Re〔Z’N/8+
k〕が第1の入力端に供給され、前記値cが第2の入力
端に供給される第15の乗算手段と、前記式z’nに対
する前記逆高速フーリエ変換FFT-1・z’nの演算結
果の虚数部Im〔Z’N/8+k〕が第1の入力端に供
給され、前記値cが第2の入力端に供給される第16の
乗算手段と、前記式z’nに対する前記逆高速フーリエ
変換FFT-1・z’nの演算結果の実数部Re〔Z’N
/8−k〕が第1の入力端に供給され、前記値sが第2
の入力端に供給される第17の乗算手段と、前記式z’
nに対する前記逆高速フーリエ変換FFT-1・z’nの
演算結果の虚数部Im〔Z’N/8−k〕が第1の入力
端に供給され、前記値sが第2の入力端に供給される第
18の乗算手段と、前記式z’nに対する前記逆高速フ
ーリエ変換FFT-1・z’nの演算結果の実数部Re
〔Z’N/8+k〕が第1の入力端に供給され、前記値
sが第2の入力端に供給される第19の乗算手段と、前
記式z’nに対する前記逆高速フーリエ変換FFT-1・
z’nの演算結果の虚数部Im〔Z’N/8+k〕が第
1の入力端に供給され、前記値sが第2の入力端に供給
される第20の乗算手段と、前記第13の乗算手段の乗
算結果が第1の入力端に供給され、前記第14の乗算手
段の乗算結果が第2の入力端に供給される第13の加算
手段と、前記第17の乗算手段の乗算結果がインバータ
を経て第1の入力端に供給され、前記第14の乗算手段
の乗算結果が第2の入力端に供給される第14の加算手
段と、前記第15の乗算手段の乗算結果が第1の入力端
に供給され、前記第20の乗算手段の乗算結果がインバ
ータを経て第2の入力端に供給される第15の加算手段
と、前記第19の乗算手段の乗算結果が第1の入力端に
供給され、前記第16の乗算手段の乗算結果が第2の入
力端に供給される第16の加算手段と、前記第13およ
び第14の加算手段のそれぞれの加算結果が変換の奇数
係数(YN/2−2k−1)として記憶される第9の記
憶手段と、前記第15および第16の加算手段のそれぞ
れの加算結果が変換の偶数係数Y2kとして記憶される
第10の記憶手段とを具備することを特徴とする。
−n)が決定された重み付け法則を表すときの以下に示
す式によって表されるブロック内のそれぞれのサンプル
の位置nの関数として、連続するN個のサンプルからな
る複数のブロックにひとまとめにされた複数のディジタ
ルサンプルx(n)を空間的あるいは、時間的に重み付
けする窓付け操作を行い、 yn=x(n)・h(N−1−n) YN−k=−Yk−1としたときのk∈〔0,・・・,
N/2−1〕に対して、位相シフト項が前記変換係数の
階数kの関数として分布された偶数階数の変換係数Y2
kを計算するために、サンプルの複数のブロックについ
て修正ディスクリートコサイン変換を実行し、y’n=
x2n・h2n、y’’n=(xN−2n−1)・(h
N−2nー1)、W4N=cos(2π/4N)+js
in(2π/4N)としたときのk=0,・・・,N/
4に対して可逆複素変換の形として前記計算された偶数
階数の変換係数Y2kを表し、n=0,・・・,N/4
−1に対してzn=(y2n−yN/2−1−2n)+j
(yN-1−2n+yN+2+2n)であるときの補助計算
方程式に基づいて前記可逆複素変換を計算して、前記複
数のディジタルサンプルx(n)で表されるディジタル
信号(x)を符号化する。
ば、伝送された変換係数の実数部(YN/4+2k)お
よび虚数部Y2kを交換し、Re〔Zk〕およびIm
〔Zk〕の形式のフォワードMDCTの複素数部を計算
することによって前記交換された変換係数を符号化処理
し、前記複素数部の実数部および虚数部を交換し、前記
複数のディジタルサンプルx(n)の再構成をするため
にオーバーラップ計算を実行して、複数のブロックにひ
とまとめにされた複数のディジタルサンプルx(n)に
よって表されるデジタル信号xを復号化する。
よれば、連続するN個のサンプルからなる複数のブロッ
クにひとまとめにされた複数のディジタルサンプルx
(n)の連続するN/2個のサンプルを第1の記憶手段
に記憶し、前記複数のディジタルサンプルx(n)の連
続するN/2個のサンプルをN/2遅延して第2の記憶
手段に記憶する。また、第1の計算手段により、h(N
−1−n)が決定された重み付け法則を表すときの以下
に示す式によって表される複数のサンプルynの空間的
重み付けを実行することからなる前記複数のディジタル
サンプルx(n)の窓付け操作を行い、 yn=x(n)・h(N−1−n) YN−k=−Yk−1としたときのk∈〔0,・・・,
N/2−1〕に対して、位相シフト項が前記変換係数の
階数kの関数として分布された偶数階数の変換係数Y2
kを計算するために、サンプルの複数のブロックについ
て修正ディスクリートコサイン変換を実行し、n=0,
・・・,N/4−1に対してzn=(y2n−yN/2−1
−2n)+j(yN-1−2n+yN+2+2n)であると
きの補助計算方程式を用いた可逆複素変換(Y2k+N
/4+jY2k)から偶数階数の変換係数Y2kを得
て、前記複数のディジタルサンプルx(n)で表される
ディジタル信号(x)を符号化する。
送された変換係数の実数部Re〔Zk〕=YN/4+2
kおよび虚数部Im〔Zk〕=Y2kを変換係数記憶手
段に記憶し、該変換係数記憶手段に記憶された前記伝送
された変換係数を交換して読出手段により読み出す。ま
た、前記z’nを第4の複素数計算手段によって計算
し、逆高速フーリエ変換(FFT-1・z’n)を第2の
逆高速フーリエ変換計算手段によって計算し、以下に示
す方程式を方程式計算手段によって計算するとともに、
フォワードMDCTの複素数部を第2の計算手段によっ
て計算して、 Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k 連続するN個のサンプルからなる複数のブロックにひと
まとめにされ、請求項7記載のディジタル信号の符号化
装置によって符号化された複数のディジタルサンプルx
(n)で表されるディジタル信号(x)を復号化する。
ば、c=cos((2π・(4n+1))/4N)、s
=sin((2π・(4n+1)/N)としたときの偶
数のサンプルx2nおよび奇数のサンプル(xN/2−
1−2n)に基づいて、以下に示す方程式を虚数乗法に
よって第3の計算手段により計算する。 z’n=((x2n・c−xN/2−1−2n・s) +jzo-N/2・(x2n・s+xN/2−2n−1・c))・(c+js ) 式Z’(N/8−k)=FFT-1・z’nを計算するた
めに、前記式z’nに対して逆高速フーリエ変換FFT
-1を逆高速フーリエ変換計算手段によって実行し、変換
の奇数係数Y2kおよび偶数係数(YN/2−2K−
1)を計算するために、虚数乗法を第4の計算手段によ
って実行して、連続するN個のサンプルからなる複数の
ブロックにひとまとめにされた複数のディジタルサンプ
ルx(n)で表されるディジタル信号(x)を符号化お
よび復号化する。
されるデジタル信号(s)を符号化するための本願発明
の方法を説明する。本願発明の方法は、N個の連続した
サンプルの複数のブロックに一体に集められた複数のサ
ンプルに適用される。上記方法は、ブロック内の各サン
プルの位置nの関数のようなサンプルの空間または時間
重み付けよりなる窓処理において顕著なものである。上
記ブロックは、y(n)=x(n)・h(N-1-n)であり、h(N-1-n)は
確定された重み付け法を表し、修正ディスクリートコサ
イン変換(MDCT)は、次式の変換の偶数階数係数を
計算するために、サンプルブロックレベルで実行され
る。
いる位相シフト項YN-k=-Yk-1において、 (9)… Y2K=Σyn・cos(2π/(4N)・(2n+1)・(4k+1)+(4k+1)・π/4 n=0,...,N-1 計算された偶数階数変換係数Y2Kは逆複素数変換の形
式に表される。 (10)… Y2k+N/4+jY2K= (-1)k+1・W8-1・Σ((y'n-y''n+N/4)+j(y''n+y'n+N/4))・W4N(4k+1)(4n+1) n=0,...,N/4-1 なお、k=0,...,N/4−1であり、N=2に対
して、 y’n=x2n・h2n y’’n=xN−2n−1・hN−2N−1 W4N=cos(2π/4N)+j・sin(2π/4N) W8=W4Nである。
れたデジタルサンプルx(n)によって表されるデジタ
ルサンプル(x)を符号化するための装置は、N/2個
の連続したサンプルを記憶する第1および第2の入力メ
モリを備えることにおいて顕著な特徴を有する。第2の
入力メモリは、N/2分、遅延されるサンプルを格納す
るために用いられる。上記サンプルに関する窓操作を下
記に示す。この法則は、h(n-1-n)が確定された重み付け
法を表すところのサンプルyn=x(n)・h(N−1
−n)をサンプリングする空間を有する。演算回路は、
下記の形式の偶数階数変換係数を計算するために、修正
ディスクリートコサイン変換MDCTを実行可能として
いる。
ト項YN-k=-Yk-1において、k∈[0,...,n/2−
1]であり、逆複素数変換Y2k+N/4+jY2kを
用いている。上記逆複素数変換は次式から得ている。 Zk=(-1)k+1・W8-1・W4N4k+1・Σ(zn・WNn)・WN/4nk n=0,...,N/4-1 ここで、n=0,...,N/4−1に対して zn=(y2n-yN/2-1-2n)+j(yN-1-2n+yN/2+2n) である。
ルx(n)によって表されるデジタル信号xを符号化す
ることに関して述べられている。サンプルは、N個の連
続したサンプルのブロックに一体に集められている。従
来の方法では、この発明のエンコディング方法は、図1
に示すように、まず、注目しているブロック内のサンプ
ル位置nの関数として、サンプルの空間または時間重み
付けによる窓処理の実行から成る。図1において、符号
Aが付けられた第1のステップは、注目しているN個の
サンプルのブロックのサンプルの入力の処理から成り、
実際上、特別な条件下の注目しているサンプルを格納す
ることに一致する入力処理が以下に述べられている。
を、yn=x(n)・h(N-1-n)で書かれる重み付けが行なわれた
サンプルへ変換するための上述した窓処理から成るステ
ップBによって続けられる。上記式において、h(N-1-n)
は確定された重み付け法を表す。当然、上記重み付け法
は、上記条件を満足し、かつ、上記重み付け法は上述し
た(6)式において満足される形式を有する。図1の処
理Cによって表されるように、この発明の符号化法は、
Y2k=Σyn・cos(2π/(4N)・(2n+1)・(4k+1)+(4k+1)・π/4な
る式の変換の偶数階数係数を計算するための、サンプル
のブロック、特に、ステップBから得られた重み付けさ
れたサンプルに基づく、修正ディスクリートコサイン変
換MDCTの実行から成る。
階数変換は、k∈[0,...,N/2-1]に対して算出される。な
お、非対象なためYN-k=-YK-1を満足する偶数項のみを計
算する必要がある。 k∈[0,...,N/2-1] 位相シフト項が、注目している変換係数の位置kの関数
として分散されていることが上記(2)式に観測される
であろう。本願発明の符号化方法の格段の利点におい
て、上述した偶数階数変換係数は、図1のステップDに
逆複素数変換の形式で表される。
れる。この式において、次の表現は以下に与えられる値
を有する。 y’n=x2n・h2n y’’n=xN−2n−1・hN−2n−1 W4N=cos(2π/4n)+j・sin(2π/4N)
すことに気づくべきである。また、図1に示されるよう
に、上記逆複素数変換は、下記の補助計算式を用いて計
算される。 Zk=(-1)k+1・W8-1・W4N4k+1・Σ(zn・WNn)・WN/4nk n=0,...,N/4-1 上記式において、zn項は、n=0,...,N/4−
1に対して同じく算出される上記式とともに、zn=(y2n-
yN/2-1-2n)+j(yN-1-2n+yN/2+2n)を示す。一般に、特
に、本願発明の符号化方法の利点は、図1のステップD
において実証された逆複素数変換は、上記式を用いて次
に示すステップEで、非常に能率よく算出され得る。
テップEおよびDに対しての特に利点となる計算手順
は、以下に、図2を参照してさらに詳細に述べられる。
図2において、上述した逆複素数変換は、以下に示す連
続したステップによって計算される。 ステップa:n=0,...,N/4−1に対して(1
2)式を用いて、すなわち、N/2個の現実の加算の有
効な消費で、znを算出する。 ステップb:n=0,...,N/4−1に対して、 (13)… z’n=zn・WNnを算出する。
ため、N/4−1個の複素数乗算を必要とする。 WNn=cos(2π・(n/N))+j・sin(2π・(n/N)) ステップc:次式(14)を得るため、ステップbにお
いて得られた変数z'nに基づく逆高速フーリエ変換を行
なう。 (14)… Z'k=FFT-1・z'n 高速フーリエ変換は、周知の処理技術を用いた、大きさ
N/4のフーリエ変換である。
プdによって続けられる。図2は(15)式を計算する
ことから成る。 (15)… Zk=((-1)k+1・W8-1・W4N4k+1・Z'k) 上記式において、W8は8乗根を示し、かつ、W4N
4k+1は4k+1乗された4N乗根である。上記式は、k
=0,...,N/4−1に対して算出される。 ステップe:次式によって与えられた偶数階数変換係数
を連続して決定する。 (16)… Y2k ε Im[Zk] Zkの虚数部 (17)… Y2k+N/4 ε Re[Zk] Zkの実数部 したがって、偶数階数係数は上述した式のZkの虚数部
および実数部によって得られるであろう。
タルサンプルx(n)の代わりに上記ステップeで得ら
れた変換係数を転送する。この転送は符号化された信号
の形の転送によって行なわれる。偶数階数係数が決定さ
れた後の、奇数階数変換係数を決定するための同一のス
テップは、上記奇数係数が偶数階数係数である (18)… Y2K+1=-Y2k (19)… Y2k+1+N/4=-Y2K+N/4 に基づいて、非対称の関係によって与えられているた
め、必要でないことに気がつくかもしれない。
するために、序文に詳細に記載されている「Implementa
tion of Split-Radix FFT algorithms for complex rea
l and real symmetric data」と名付けられている刊行
物にPierre duhamelによって述べられる演算方法を用い
るのが好ましい。このような状況の下、本願発明の符号
化方法を実行するために要求される演算処理の数は下記
のように評価され得る。乗算の全数は、 N=2qであるので、N/4(q+1)個の実際の乗
算、および加算の全数は、 N/4(3q−1)個の実際の加算となる。
符号化方法を実施するために要求される作業の最小全数
がN・qに厳密に等しいと見ることができる。作業の上
記数を評価するにおいて、従来の方法では、難しい演算
に換算して、複素数乗算は3つの実際の乗算に3つの実
際の加算を加えることに相当する、と推定される。ま
た、それは、さらに複雑な式を生じさせるが、4つの実
際の乗算と2つの実際の加算とを用いることと等価であ
ると評価することが可能である。なお、いかなる場合で
も、作業の全数N・qは変わることがない。
明細書の記載の序文で述べた従来の符号化方法に関し
て、4つの要因によって改善点を提供する。特に、難し
い演算に対して上記定義された項は、最良のFFTコー
ディング方法と両立できる。このタイプのFFTコーデ
ィング処理と、本願発明のコーディング方法のステップ
Eにおいて与えられたZkを定義する(11)式と(1
2)式との間の関係は、両面性があり、かつ、1対1で
ある。ASSPへ移行されたjounal IEEE発行、38
巻、1504−1511頁、1990年9月号、「Algo
rithms meeting the lower bound on the multiplicati
ve complexity oflength2n DFT's and their connectio
n with practical algorithms」と名付けられた論文に
おいて、同一の著者、Pierre Duhamel は、大きさNと
一致するスプリット根型の高速フーリエ変換の活性部
が、大きさN/16の上述した式によって実際に与えら
れること、および、上記活性部が、例えば、W8-1およ
び(−1)k+1のような定数がなければ、おそらく、上
述した主要部を用いて連続して算出され得ることを示し
ている。
ト根型高速フーリエ変換による処理の改善がより良い修
正ディスクリートコサイン変換(MDCT)過程を得る
という事実として理解されなければならない。しかし、
本願発明のコーディング方法に提供されるような、修正
ディスクリートコサイン変換(MDCT)過程の改善
は、現在知られている方法よりも、高速フーリエ変換
(FFT)を用いる処理のより良い方法も提供する。そ
の結果として、本願発明の方法は、このタイプの処理に
関する技術の現在の状態に提供された方法よりも、作業
の数の点から、すなわち実際の乗算および加算の数の点
からより良いものである。
て、ステップdからW8-1による乗算を取ること、およ
びステップbにそれを含めることによって、上述したス
テップbおよびdを変更することが可能である。
インテーブル、すなわち本願発明のコーディング方法を
向上するために要求される値のテーブルにおける対称性
を変更する。ステップbにW8-1による乗算を含めるこ
とは、対応するメモリテーブルの対称性を変更する。そ
れでも、上記テーブルの利用は容易にし得る。加えて、
次式で表されるZkの補助計算のための式を用いて、図
2に示すステップbおよびdは交換され得る。 (12b)… Zk=(-1)k+1・W8-1・WNk・Σ(zn・W4N)(4n+1)・WN/4nk n=0,...,N/4-1 W4N(4n+1)による複素数乗算は、演算の初期または最
終に位置する。最終は選択され、この最終は最小演算時
間の関数として決定される。
した符号化方法を用いられる符号化されたデジタルサン
プルx(n)によって表されるデジタル信号xを復号化
する方法は、図3を参照して説明する。図3において、
ステップA’では、予め転送された実数部YN/4+2
kおよび虚数部Y2kがこの発明の方法によって符号化
された後に授受される。復号化方法は、転送された変換
係数の実数部YN/4+2kおよび虚数部Y2kを置換
するステップB’からなる。本願発明の符号化およびデ
コーディング方法の顕著な利点には、該方法が図3に示
すようにステップC’において復号化している間に存在
する。ステップC’では、上述した置換された係数は、
形式Re[Zk]、Im[Zk]の先のMDCTの複素
数部を計算することによって符号化処理へ服従させられ
る。
より、置換した係数を符号化処理へ服従させることから
成る作業は、図2を参照して説明したステップb,cお
よびdの組み合わせが可逆複素変換を形成することによ
って正当化される。したがって、上記ステップによって
構成された変換の逆変換は、入力信号の実数部および虚
数部が置き換えられた後に、修正ディスクリートコサイ
ン変換または逆MDCTの核を構成する。該MDCTは
対応する前のMDCTによって計算される。上記置き換
えは、形式z*=b+jaの転換された複素信号を発生
するところの複素形式z=a+jbのいかなる入力信号
に対しても効果を有し、これによって、転換された入力
信号が発生する。その後、前のMDCTが、Y*kで書
かれる出力信号を出力すべく、転換された入力信号に関
して実行される。出力信号Y*kの実数部および虚数部
を置き換えることは、逆出力信号Ykを発生することを
可能にするということに容易に気が付くであろう。これ
は、上述したMDCTの複素数部の実数部と虚数部を置
き換えることに相当する。
ーラップがサンプルx(n)を再構成するために算出さ
れるステップG’によって引き続き行なわれる。本願発
明の復号化方法では、先のMDCTの複素数部が処理さ
れる図3のステップC’、D’およびE’によって構成
される処理は、符号化方法と同様な符号化する間に実行
される処理と同一である。また、上記処理は、下記の演
算からなるが、転送後に授受されると、置き換えられた
実数部および虚数部に基づき実行される。ステップC’
では、n=0,...,N/4−1に対してz’n=z
n・WNnを計算し、ステップD’では、大きさN/4の
逆高速フーリエ変換FFT-1を用いて計算し、すなわ
ち、Z’k=FFT-1・z’nを計算する。そして、ス
テップE’では、k=0,...,N/4−1に対して
次式によって与えられる式を計算する。 (15)… Zk=((-1)k+1・W8-1・W4N4k+1)・Z'k
号化するための装置の詳細は図4を参照して説明する。
図4において、N個の連続したサンプルのブロックに一
体に集められた上述したようなサンプルを用いて、装置
は、各々、N/2個の連続したサンプルを記憶すること
ができる第1入力メモリ111および第2入力メモリ1
12を備える。限界のない例では、第2入力メモリ11
2は、N/2だけ遅延された(分割された?)サンプル
を記憶すべく動作する。他の異なる例においては、第1
入力メモリ111は、N/2個の偶数番目の箇所または
アドレッシングされるサンプルを記憶すべく動作する。
一方、第2入力メモリ112はN/2個の奇数番目のサ
ンプルまたは箇所を記憶すべく動作する。
付けられた演算回路は、各々、以下に示すことを処理す
べく動作する。(12)式に従ってznを算出し、(1
3)式を満足すべくz'nを算出し、(14)式を用い
てN/4の逆高速フーリエ変換FFT-1を計算し、か
つ、(15)式を用いてZkを算出する。当然、限界の
ない例では、サンプルの空間または時間重み付けから成
る上述した窓処理操作が続き、演算回路1,2,3およ
び4は、特に、補助の計算式を用いて、(16)式およ
び(17)式を満足する可逆複素数変換Y2k+N/4
+jY2kを連続して確定すべく、修正ディスクリート
コサイン変換MDCTを実行する。
めの上述した演算回路1,2,3および4は、連続して
カスケードに接続された本願発明の符号化装置であると
ころの、znを算出するための複素数演算回路1、n=
0,...,N/4−1に対してz'n=zn・WNnを
計算するための複素数乗算回路2、FFT-1、すなわ
ち、(14)式を満足するZ'kを算出するための逆高
速フーリエ変換演算回路3およびk=0,...,N/
4−1に対して(15)式を満足するZkを算出するた
めの第3回路4から成る。この第3回路4は、上述した
ように、変換の偶数番目および奇数番目の係数を連続し
て決定すべく動作する。図4に示すように、znを算出
する複素数演算回路1は、第1および第2入力メモリ回
路111および112に接続された第1および第2の加
算回路121および122を備える。
られたサンプル値y2nを、直接、第1入力メモリ11
1から授受し、次に、乗算器123を通して、「−1」
が乗算された第2入力メモリ112から供給される重み
付けられたサンプル値yN/2−1−2nを授受する。
第2加算回路122iは、まず、第1入力メモリ111
からの重み付けサンプル値yN/2+2nを授受し、次
に、第2入力メモリからの重み付けサンプル値yN−1
−2nを授受する。当然、第1および第2加算回路12
1,122は、複素数値znの組を築き上げるために、
実際に、上述した重み付けサンプル値を授受する基本的
な加算回路とともに、基本的な加算回路121i,12
2iによって構成されてもよい。当然、各第1加算回路
121iは、「−1」を乗算するための乗算回路123
iを通して第2入力メモリ112へ接続されている。し
たがって、第1加算回路121は、N/2個の基本加算
回路121iによって構成されてもよい。同様に、第2
加算回路122に適用し、基本加算回路122iで構成
してもよい。
121iおよび122iの数は、処理速度と記憶容量と
の間の妥協点に決定される。さらに、第1複素数演算回
路1は、N/4個の値を記憶するための第1および第2
の補助メモリ131,132を備える。これら第1およ
び第2の補助メモリ131,132は、第1および第2
の加算回路121,122の出力に、特に、それらを構
成している基本加算回路へ接続されている。補助メモリ
131および132は、znの実数部および虚数部の各
々を記憶可能にする。同様に、かつ、図4に示すよう
に、z'nを算出する複素数乗算のための第2の回路2
は、定数WNnと名付けられる一定値を記憶するための
メモリ回路20aおよび20bを備える。これらのメモ
リ回路は、高速読み込みアクセスに適したRAM(リー
ド・オンリ・メモリ)タイプのメモリ回路によって構成
され、かつ、それ自身に上記一定値を記憶している読み
込み専用メモリ(ROM)からロードされてもよい。
aは、まず、zn=a+jbの実数部aを授受すべく用
意され、かつ、次に、WNnの実数部cos2π(n/
N)または虚数部sin2π(n/N)のいずれか一方
を授受すべく用意される。これら第1の乗算回路21a
および22aは、a・cos2π(n/N)およびa・s
in2π(n/N)なる積を送出する。同様にして、第
2の乗算回路21bおよび22bは、第1に、znの虚
数部bを授受し、かつ、第2に、WNnの実数部cos
2π(n/N)または虚数部sin2π(n/N)を授
受する。これら第2の乗算回路21bおよび22bは、
積b・cos2π(n/N)およびb・sin2π(n/
N)を送出する。図4に示すように、第1の加算回路2
3は、第1の乗算回路からのa・cos2π(n/N)
なる積を授受するともに、第2の乗算器22bから「−
1」を乗算する乗算器23aを介して積b・sin2π
(n/N)を授受する。したがって、第1の加算回路2
3は、z'n=a'+jb’の実数部であるa’=a・c
os2π(n/N)−b・sin2π(n/N)を送出
する。
に、第1の乗算器22aからの積a・sin2π(n/
N)を授受し、かつ、第2に、乗算器21bからの積b
・cos2π(n/N)を授受する。この第2の加算回
路24は、znの虚数部であるb’=a・sin2π
(n/N)+b・cos2π(n/N)を送出する。乗
算器および加算器内の符号iは、上記条件の下で、考え
得る最大数N/4に達するまでの多数のそれらがあるこ
とを示す。加えて、第1および第2の補助メモリ25,
26は、各々、z'nの実数部a'と虚数部b'の値を記
憶する。上記補助メモリ25,26は、各々、N/4個
の値を記憶することが可能なRAMタイプのメモリであ
る。
の逆高速フーリエ変換FFT-1を実行する。これは、上
述した文献のPierre Duhamel によって述べられた計算
方法を実行する備えられた回路によって実行されてもよ
い。最後に、同じく図4に示すように、式Zkの第3の
複素数乗算回路4は、定数W8-1およびW4N4k+1の実
数部および虚数部を記憶するための回路40a,40b
を備える。一般に、定数に対する永久的なメモリ、例え
ば、符号20a,20bまたは40a,40bで示され
るようなメモリは、図を押込め過ぎることを避けるため
に、対応する演算回路の図から省略されている。また、
図4に見ることができるように、第3の複素数乗算器4
は、第2の複素数乗算器2のそれらに相似する構成部分
を備えている。特に、それは、第1に、Z'kの実数部
を授受し、第2に、メモリ回路40aおよび40bから
の定数W8-1およびW4N4k+1の各々の実数部および虚
数部を授受する第1の乗算器41aおよび42aを備え
る。第1の乗算器41a,42aは、対応する積を送出
する。
は、第1に、Z'kの虚数部を授受し、第2に、同じメ
モリ回路40aおよび40bからの定数W8-1およびW
4N4k+1の各々の実数部および虚数部を授受する。第2
の乗算器41bおよび42bは、対応する積を送出す
る。さらに、第1および第2の加算回路43,44が備
えられている。第1の加算回路43は、第1の乗算器4
1aによって送出される積および「−1」を乗算する乗
算器43aを介して第2の乗算回路42bによって送出
される積を授受する。上記第1の加算回路43は次式で
示すようなZkの実数部を送出する。 Re[Zk]=YN/4+2k
算回路42aによって送出される積、および、第2に、
第2の乗算回路41bによって送出される積を授受し、
次式におけるZkの虚数部を送出する。 Im[Zk]=Y2k サンプルx(n)によって表されるデジタル信号xを復
号化するための装置のさらに詳細な説明は、図5を参照
して与えられる。なお、上記サンプルは上述した符号化
方法を用いて符号化されている。図5に示すように、復
号化装置は、図4の符号化回路によって実証されたよう
に、転送された変換係数を記憶したり、読み出したりす
るための記憶読出し回路5を備える。上記係数は、YN
/4+2k=Re[Zk] および Y2k=Im[Z
k]と書かれる。記憶したり、読み出したりするための
回路5は、読み出しに関し適切なアドレッシングによっ
て記憶されると、転送された変換係数を置き換えること
を可能とする。上記置き換えは、記憶された実数部およ
び虚数部をクロスさせることによって、図5に示めされ
る。
素数部を算出するための複素数部算出回路6を備えてい
る。上記回路6はカスケード接続された、Z'kを算出
するための複素数乗算回路62および対応する逆フーリ
エ変換FFT-1を実行するための回路63を備える。回
路64は次式で示されるZkを算出する。 Zk=(−1)k+1・W8-1・W4N4k+1・Z'k 先のMDCTの複素数部を算出するための回路6は、符
号化装置のモジュール2,3および4と、各々、同一の
モジュール62,63および64を備える。したがっ
て、上記モジュール62,63および64の詳細な説明
は省略する。加えて、図5に示すように、復号化装置
は、上述した入力モジュール5により実行される置換に
相似する、読み出しのためのアドレッシングによってZ
kの実数部および虚数部を交換または置換した後に以下
に述べるオーバーラップを算出するためのモジュール7
を備える。
によって先のMDCTの複素数部を算出するためのモジ
ュール6は、内部で、データを逆の手順で処理し、最終
的に、次式で示される複素数変数Znを送出することが
理解されるであろう。 Zn=(x2n・h2n−xN/2−2n−1・hN/2−2n−1) +j(xN−2n+1・hN−2n+1−xN/2+2n・hN/2+2n) 上記式において、実数部は、虚数部に関するブロックの
半分の、すなわちN/2個のサンプル分遅延される。上
述した数式において、符号化処理の重み付け窓および復
号化処理に用いられる窓は同一であると仮定される。上
述した遅延に対する補償後、上記式は、下記のように書
き換えることができる。 Zn=(x2n+jxN/2−2n+1)(h2n/+jh(N/2+2n)) 上述した条件は、原信号を完全に再構成し、数式:h2
n−jhN/2+2nによるznに対する上述した式の
複素数乗算によって、サンプルx(n)のシーケンス
(順列)を再構成可能とする。
するためのモジュール7は、上述した式内の対応する虚
数部に関係する実数部のN/2サンプル分の遅延に対し
て補償するために、各々、ブロックmおよびm−1に一
致するサンプルを記憶可能にするN/2個のサンプルを
記憶するための第1および第2のメモリ71aおよび7
1bを備える。
4や、図5に示した複素数乗算回路62,64は、例え
ば、図6に示すモジュール7に示すように、乗算器73
ai,73bi,74ai,74biと、インバータ7
5aiと、2つの加算器75i,76iとから構成され
る。ここで、モジュール7は、オーバラップを計算する
回路部である。
示すように、それぞれメモリ70a,70bに記憶され
る。各乗算器および各加算器の指数iは、単に、図番の
記載が混雑するのを避けるために付与したものに過ぎな
い。また、各乗算器および各加算器の数は、N/2に等
しくなる。これは、N/2値が記憶される操作に対応さ
せた数に合致させるためである。加算器75i,76i
は、それぞれ再構成(復元)されたサンプルx2n,x
N/2−2n−1に対応した信号を出力する。
装置の実施例をさらに改良したものでは、「窓」処理や
「オーバラップ」処理を考慮に入れた直接的な計算プロ
セスが可能となっており、好都合である。こうしたこと
は、「窓」が完全に再構成(復元)される計算プロセス
であることを表わす。「窓」を完全に再構成(復元)す
るには、まず、第1に入力されるサンプルx(n)の内
から偶数階数と奇数階数とからなる2つのサンプル成分
を回転させることで達成される。さらに、オーバラップ
計算に関しては次のようになる。つまり、復号化により
サンプルを再構成(復元)する場合には、サンプルブロ
ック間がN/2サンプルづつ適正に重なり合うという性
質が要求される。このような性質を備えさせるべく、N
/2サンプル分の遅延を補償する。すなわち、h2n−
jhN/2+2nで表現される式に基づき複素数乗算さ
れたサンプルx2nを再構成する順序を得るためにN/
2サンプル分の遅延を補償している。
「窓」を用いて復号化プロセスあるいは符号化プロセス
の簡略化を図る場合、(6)式で定義される演算子は4
N次根の4n+1累乗、すなわち、W4N4n+1に他なら
ない。なお、この(6)式は、出力サンプルx2nを再
構成(復元)する操作、あるいは入力サンプルx(n)
の重み付け操作の期間中に任意サンプルに適用される重
み付け窓として定義される。
において、前述した(15)式は回転に相似した処理操
作によって求められる。この処理操作は、逆高速フーリ
エ変換により係数z’kを算出するモジュール3(図4
参照)に適用される。このように、窓処理部1と、c=
cos2π(4n+1/4N),s=sin2π(4n
+1/4N)で表わされる処理技術によりz’nを算出
する複素数乗算部2(図4参照)とは、下記(18)〜
(20)式によって表現される。
すると、 (20)… z'n=(x2n-j・xN/2-1-2n) +(-(x2n・s2+xN/2-2n-1・cs) +j(x2n・cs+xN/2-1-2n・c2) +zo-N/2・((x2n・s2+xN/2-2n-1・cs) +j(x2n・cs+xN/2-2n-1・c2))
るサンプル値に適用されるN/2サンプル遅延演算子
は、zo-N/2により表現される。このように、前述した
プロセスは、サンプルNの各ブロックに含まれる2つの
量子、すなわち、(x2n・s2+xN/2−2n−1
・cs)と(x2n・cs+xN/2−2n−1・
c2)とで定義される量子に基づく特殊な窓処理が考慮
されている。
照)と、znを算出する窓処理部1とにおいてなされる
回転処理の回数を減少させるための処理が、図7に示す
加算値N/4の記憶により可能となる。ここで、上記2
つの回転処理は、一般に、3N回の操作で行われてお
り、これをN乗算とN加算とからなる2N回の操作を行
うことで処理回数を減少させている。
プロセスを図7を参照して説明する。図示される装置
は、窓処理とz’nを算出する処理とにおける処理回数
を低減するモジュール8を含む。なお、窓処理を行うモ
ジュール1とz’nを複素数乗算して算出するモジュー
ル2とは図4に開示されている。このモジュール8の動
作は、N/4のサイズで逆高速フーリエ変換を行うモジ
ュール3に従う。また、該モジュール3の動作は、偶数
の変換係数Y2kと奇数の変換係数YN/2−2k−1
とを算出するモジュール9に従う。
を示す図である。このような構成において、乗算や加算
などの演算子が指数iを伴って数値的な基準として与え
られる場合、これに対応する複数の要素演算子が生成さ
れる。すなわち、この図に示すモジュール8,9におい
ては、指数iの全てがN/4に対応する複数の要素演算
子として表示される。これにより、メモリ容量を最適化
し、リアルタイム処理が可能となる符号化/復号化装置
が実現される訳である。
出処理を行うモジュール8は、サンプルx(n)の内、
奇数サンブルを記憶する第1入力メモリ81aと偶数サ
ンプルを記憶する第2入力メモリ81bとを備える。第
1乗算器82a,82bは、それぞれメモリ81aに記
憶される奇数サンプルxN/2−1−2n、メモリ81
bに記憶される偶数サンプルx2nが供給される。この
第1乗算器82aには、外部メモリ80aから出力され
る定数Cがインバータ821iを介して供給される。一
方、第1乗算器82b側には、外部メモリ80aから定
数Cが直接供給される。ここで、インバータ821i
は、入力する信号に「−1」を乗じて出力する。
が供給される第2乗算器83bには、外部メモリ80b
から出力される定数Sが乗算器89を介して供給され
る。偶数サンプルx2nが供給される第2乗算器83a
には、外部メモリ80bから定数Sが直接供給される。
なお、上記乗算器89は、入力される定数Sに「−1」
を乗算して出力する。
の出力と第2乗算器83aの出力とを加算して出力す
る。一方、第1加算回路84bは、第1乗算器82bの
出力と第2乗算器83bの出力とを加算して出力する。
第2加算回路86aは、第1加算器84aの出力信号
と、この出力信号を遅延回路85aでN/2値分遅延さ
せると共に、インバータ861iを介した信号とを加算
する。同様に、第2加算回路86bは、第1加算器84
bの出力信号と、この出力信号を遅延回路85bでN/
2値分遅延させた信号とを加算する。
aの出力信号と、インバータ871iを介して反転され
た奇数サンプルxN/2−1−2nとを加算して出力す
る。第3加算回路87bは、第2加算回路86bの出力
信号と、偶数サンプルx2nとを加算して出力する。こ
の第3加算回路87a,87bの各出力は、z’n(複
素数)に対応するものである。すなわち、第3加算回路
87aの出力は、複素数z’nの実数部に相当し、一
方、第3加算回路87bの出力は、複素数z’nの虚数
部に相当する。そして、このz’nの実数部と虚数部と
は、それぞれRAM型の一時記憶メモリ88a,88b
に記憶される。ここで、メモリ88a,88bには、
z’nの実数部と虚数部とがそれぞれN/4値として記
憶される。
3は、図4に示したモジュール3と同一の構成である。
モジュール9は、このモジュール3から出力される逆フ
ーリア変換結果に基づいて偶数の変換係数Y2kと奇数
の変換係数YN/2−2k−1とを算出する。図7に示
すように、モジュール9は第1乗算器91a,91b,
91c,91dを備える。第1乗算器91aには、逆フ
ーリエ変換結果の実数部Re{Z’N/8−k}が供給
される。第1乗算器91bには、逆フーリエ変換結果の
虚数部Im{Z’N/8−k}が供給される。さらに、
第1乗算器91c,91dには、それぞれ逆フーリエ変
換結果の実数部Re{Z’N/8+k}、虚数部Im
{Z’N/8+k}が供給される。また、これら第1乗
算器91a,91b,91c,91dには、外部メモリ
90aから読み出される定数Cが供給される。
dには、外部メモリ90bから読み出される定数Sが供
給される。そして、第2乗算器92aは、上記実数部R
e{Z’N/8−k}にこの定数Sを乗じて出力する。
第2乗算器92bは、上記虚数部Im{Z’N/8−
k}に定数Sを乗じて出力する。さらに、第2乗算器9
2cは、上記実数部Re{Z’N/8+k}に定数Sを
乗じて出力する。第2乗算器92dは、上記虚数部Im
{Z’N/8+k}に定数Sを乗じて出力する。
の出力と第2乗算器92bの出力とを加算する。第2加
算回路92bは、第2乗算器92aの出力をインバータ
94bで反転した信号と、第1乗算器91bの出力とを
加算する。第3加算回路93cは、第1乗算器91cの
出力と、第2乗算器92dの出力をインバータ94cで
反転した信号とを加算する。さらに、第4加算回路93
dは、第1乗算器91dの出力と第2乗算器92cの出
力とを加算する。
93bの出力とは、奇数の変換係数YN/2−2k−1
として一時記憶メモリ95aに記憶される。一方、第3
加算回路93cの出力と第4加算回路93dの出力と
は、偶数の変換係数Y2kとして一時記憶メモリ95b
に記憶される。
ールは、特に、計算時間を短縮する点で有利となる。す
なわち、図4に示した実施例では、窓処理に対応する操
作が3N回要求されるのに対し、図7に示す構成によれ
ば、同様の窓処理が2N回の操作で実現される。
5に基づいて説明した逆変換原理を図7に適用すれば、
同様の構成として実現できる。この場合、復号化モジュ
ールは、単に、図7に示す入出力関係を対応する入力と
出力とに置き換えると共に、前述した回転操作の定数
C,Sによって乗算されるデータ列の出力手順を与える
だけであるから、その詳細に関する説明を省略する。
あるいは復号化する処理プログラムの一例を示すソース
コードリストである。この処理プログラムは、図7に示
す装置によって達成される処理に対応するものであり、
英語のコメントが付与されたC言語によって記述されて
いる。
ば、符号化に必要とされる演算数を実際上、減少させ、
演算の負担を軽減できるという利点が得られる。また、
この発明によれば、処理手順を簡易化でき、リアルタイ
ム処理に適合したコストおよび信頼性を得ることができ
るという利点が得られる。
実行されるかを示す一般のフローチャートである。
行であるところの特に有用なフローチャートである。
の実行のフローカHートである。
ための装置の一実施例の構成を示すブロック図である。
デジタル信号を復号化するための装置の構成を示すブロ
ック図である。
である。
演算処理が複素数積の演算に統合され、これによって、
窓処理またはオーバーラッピングの(処理の数に特有
の)コストを削減する、符号化または復号化するための
装置の一実施例の構成を示すブロック図である。
る処理プログラムの一例を示すソースコードリストであ
る。
Claims (16)
- 【請求項1】 h(N−1−n)が決定された重み付け
法則を表すときの以下に示す式によって表されるブロッ
ク内のそれぞれのサンプルの位置nの関数として、連続
するN個のサンプルからなる複数のブロックにひとまと
めにされた複数のディジタルサンプルx(n)を空間的
あるいは、時間的に重み付けする窓付け操作を行い、 yn=x(n)・h(N−1−n) YN−k=−Yk−1としたときのk∈〔0,・・・,
N/2−1〕に対して、位相シフト項が前記変換係数の
階数kの関数として分布された数式1の形で表される偶
数階数の変換係数Y2kを計算するために、サンプルの
複数のブロックについて修正ディスクリートコサイン変
換を実行し、 【数1】 y’n=x2n・h2n、y’’n=(xN−2n−
1)・(hN−2nー1)、W4N=cos(2π/4
N)+jsin(2π/4N)としたときのk=0,・
・・,N/4に対して数式2に示す可逆複素変換の形と
して前記計算された偶数階数の変換係数Y2kを表し、 【数2】 n=0,・・・,N/4−1に対してzn=(y2n−y
N/2−1−2n)+j(yN-1−2n+yN+2+2n)
であるときの数式3に示す補助計算方程式に基づいて前
記可逆複素変換を計算して、 【数3】 前記複数のディジタルサンプルx(n)で表されるディ
ジタル信号(x)を符号化することを特徴とするディジ
タル信号の符号化方法。 - 【請求項2】 前記可逆複素変換は、 n=0,・・・,N/4−1に対して前記znを計算す
る第1のステップと、 n=0,・・・,N/4−1に対してz’n=zn・W
Nnを計算する第2のステップと、 サイズN/4の逆高速フーリエ変換FFT-1によって以
下に示す式で表される対応する係数Z’kを計算する第
3のステップと、 Z’k=FFT-1・z’n k=0,・・・,N/4−1に対して以下に示す方程式
を計算する第4のステップと、 Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k Zkの虚数部Im〔Zk〕=Y2kと、Zkの実数部R
e〔Zk〕=(Y2k+N/4)とからなる偶数階数の
変換係数を連続的に決定する第5のステップと、 前記ディジタル信号(x)を構成する前記複数のディジ
タルサンプルx(n)の代わりに、前記偶数階数の変換
係数Y2kおよび(Y2k+N/4)を符号化された信
号の形で伝送する第6のステップとからなる連続するス
テップにより計算されることを特徴とする請求項1記載
のディジタル信号の符号化方法。 - 【請求項3】 前記補助計算方程式は、数式4の形で表
され、 【数4】 前記第2のステップと、前記第4のステップとが交換さ
れ、 W4N(4n+1)による虚数乗法が前記計算の最初あるい
は、最後に設けられていることを特徴とする請求項2記
載のディジタル信号の符号化方法。 - 【請求項4】 伝送された変換係数の実数部(YN/4
+2k)および虚数部Y2kを交換し、 Re〔Zk〕およびIm〔Zk〕の形式のフォワードM
DCTの複素数部を計算することによって前記交換され
た変換係数を符号化処理し、 前記複素数部の実数部および虚数部を交換し、 前記複数のディジタルサンプルx(n)の再構成をする
ためにオーバーラップ計算を実行して、 請求項1記載のディジタル信号の符号化方法によって符
号化された前記複数のディジタルサンプルx(n)で表
されるディジタル信号(x)を復号化することを特徴と
するディジタル信号の復号化方法。 - 【請求項5】 前記フォワードMDCTの複素数部を計
算することによる符号化処理と同一の前記符号化処理
は、 n=0,・・・,N/4−1に対してz’n=zn・W
Nnを計算する第2のステップと、 サイズN/4の逆高速フーリエ変換FFT-1によって以
下に示す式で表される係数Z’kを計算する第3のステ
ップと、 Z’k=FFT-1・z’n k=0,・・・,N/4−1に対して以下に示す方程式
を計算する第4のステップと Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k からなることを特徴とする請求項4記載のディジタル信
号の復号化方法。 - 【請求項6】 可逆複素変換を形成する一連の前記第2
ないし第4のステップとともに、前記第2ないし第4の
ステップによって構成された逆変換は、逆MDCTのコ
アを構成し、 前記逆MDCTは、 逆入力信号を生成するために、以下に示す式により入力
信号の実数部と虚数部とを交換する第7のステップと、 z=a+jb→z*=b+ja 出力信号Y*kを生成するために、逆入力信号のフォー
ワードMDCTを計算する第8のステップと、 逆出力信号Ykを生成するために、出力信号Y*kの実
数部と虚数部とを交換する第9のステップとからなる対
応するフォーワードMDCTによって計算されることを
特徴とする請求項5記載のディジタル信号の復号化方
法。 - 【請求項7】 連続するN個のサンプルからなる複数の
ブロックにひとまとめにされた複数のディジタルサンプ
ルx(n)の連続するN/2個のサンプルを記憶する第
1の記憶手段と、 前記複数のディジタルサンプルx(n)の連続するN/
2個のサンプルをN/2遅延して記憶する第2の記憶手
段と、 h(N−1−n)が決定された重み付け法則を表すとき
の以下に示す式によって表される複数のサンプルynの
空間的重み付けを実行することからなる前記複数のディ
ジタルサンプルx(n)の窓付け操作を行い、 yn=x(n)・h(N−1−n) YN−k=−Yk−1としたときのk∈〔0,・・・,
N/2−1〕に対して、位相シフト項が前記変換係数の
階数kの関数として分布された数式5の形で表される偶
数階数の変換係数Y2kを計算するために、サンプルの
複数のブロックについて修正ディスクリートコサイン変
換を実行し、 【数5】 n=0,・・・,N/4−1に対してzn=(y2n−y
N/2−1−2n)+j(yN-1−2n+yN+2+2n)
であるときの数式6に示す補助計算方程式を用いた可逆
複素変換(Y2k+N/4+jY2k)から偶数階数の
変換係数Y2kを得る第1の計算手段と 【数6】 を具備し、前記複数のディジタルサンプルx(n)で表
されるディジタル信号(x)を符号化することを特徴と
するディジタル信号の符号化装置。 - 【請求項8】 前記第1の計算手段は、 前記znを計算する第1の複素数計算手段と、 n=0,・・・,N/4−1に対してz’n=zn・W
Nnを計算する第2の複素数計算手段と、 (Z’k=FFT-1・z’n)を計算する第1の逆高速
フーリエ変換計算手段と、 k=0,・・・,N/4−1に対して以下に示す方程式
を計算し、 Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k Zkの虚数部Im〔Zk〕=Y2kと、Zkの実数部R
e〔Zk〕=(Y2k+N/4)とからなる偶数階数の
変換係数を連続的に決定する第3の複素数計算手段とが
縦続接続されていることを特徴とする請求項7記載のデ
ィジタル信号の符号化装置。 - 【請求項9】 前記第1の複素数計算手段は、 前記第1および第2の記憶手段に接続され、 重み付けされたサンプル値y2nが前記第1の記憶手段
から第1の入力端に直接供給され、重み付けされたサン
プル値(yN/2−1−2n)が(−1)を乗算する乗
算器を介して第2の入力端に供給される第1の加算手段
と、 重み付けされたサンプル値(yN/2+2n)が前記第
1の記憶手段から第1の入力端に供給され、重み付けさ
れたサンプル値(yN−1−2n)が前記第2の記憶手
段から第2の入力端に供給される第2の加算手段と、 それぞれ記憶容量N/4を有し、前記第1および第2の
加算手段のそれぞれの出力端にそれぞれ接続され、zn
=(y2n−yN/2−1−2n)+j(yN-1−2n+yN
+2+2n)の実数部および虚数部をそれぞれ記憶する
第1および第2の補助記憶手段とを具備することを特徴
とする請求項8記載のディジタル信号の符号化装置。 - 【請求項10】 前記第2の複素数計算手段は、 定数(WNn=cos(2πn/N)+jsin(2π
n/N))を記憶する第1の定数記憶手段と、 (zn=a+jb)の実数部aが第1の入力端に供給さ
れ、前記定数WNnの実数部(cos(2πn/N))
が第2の入力端に供給され、積(a・cos(2πn/
N)を出力する第1の乗算手段と、 (zn=a+jb)の実数部aが第1の入力端に供給さ
れ、前記定数WNnの虚数部(sin(2πn/N))
が第2の入力端に供給され、積(a・sin(2πn/
N))を出力する第2の乗算手段と、 (zn=a+jb)の虚数部bが第1の入力端に供給さ
れ、前記WNnの実数部(cos(2πn/N))が第
2の入力端に供給され、積(b・cos(2πn/
N))を出力する第3の乗算手段と、 (zn=a+jb)の虚数部bが第1の入力端に供給さ
れ、前記WNnの虚数部(sin(2πn/N))が第
2の入力端に供給され、積(b・sin(2πn/
N))を出力する第4の乗算手段と、 前記積(a・cos(2πn/N))が前記第1の乗算
手段から第1の入力端に供給され、前記積(b・sin
(2πn/N))が前記第4の乗算手段から(−1)を
乗算する乗算器を介して第2の入力端に供給され、a’
=a・cos(2πn/N)−b・sin(2πn/
N)、b’=a・sin(2πn/N)+b・cos
(2πn/N))としたときのa’を出力する第3の加
算手段と、 前記積(a・sin(2πn/N))が前記第2の乗算
手段から第1の入力端に供給され、前記積(b・cos
(2πn/N))が前記第3の乗算手段から第2の入力
端に供給され、前記(z’n=a’+jb’)の虚数部
b’を出力する第4の加算手段と、 前記z’nの実数部a’および虚数部b’をそれぞれ記
憶する第3および第4の補助記憶手段とを具備すること
を特徴とする請求項8記載のディジタル信号の符号化装
置。 - 【請求項11】 前記第3の複素数計算手段は、 定数W8-1およびW4N4k+1の実数部および虚数部を記
憶する第2の定数記憶手段と、 式Z’kの実数部が第1の入力端に供給され、前記定数
WNnおよびW4N4k+1の実数部が前記第2の定数記憶
手段から第2の入力端に供給され、対応する第1の積を
出力する第5の乗算手段と、 式Z’kの実数部が第1の入力端に供給され、前記定数
WNnおよびW4N4k+1の虚数部が前記第2の定数記憶
手段から第2の入力端に供給され、対応する第2の積を
出力する第6の乗算手段と、 式Z’kの虚数部が第1の入力端に供給され、前記定数
WNnおよびW4N4k+1の実数部が前記第2の定数記憶
手段から第2の入力端に供給され、対応する第3の積を
出力する第7の乗算手段と、 式Z’kの虚数部が第1の入力端に供給され、前記定数
WNnおよびW4N4k+1の虚数部が前記第2の定数記憶
手段から第2の入力端に供給され、対応する第4の積を
出力する第8の乗算手段と、 前記第1の積が前記第5の乗算手段から第1の入力端に
供給され、前記第4の積が前記第8の乗算手段から(−
1)を乗算する乗算器を介して第2の入力端に供給さ
れ、式Zkの実数部Re〔Zk〕=YN/4+2kを出
力する第5の加算手段と、 前記第2の積が前記第2の乗算手段から第1の入力端に
供給され、前記第3の積が前記第7の乗算手段から第2
の入力端に供給され、前記式Zkの虚数部Im〔Zk〕
=Y2kを出力する第6の加算手段とを具備することを
特徴とする請求項8記載のディジタル信号の符号化装
置。 - 【請求項12】 伝送された変換係数の実数部Re〔Z
k〕=YN/4+2kおよび虚数部Im〔Zk〕=Y2
kを記憶する変換係数記憶手段と、 該変換係数記憶手段に記憶された前記伝送された変換係
数を交換して読み出す読出手段と、 前記z’nを計算する第4の複素数計算手段と、逆高速
フーリエ変換(FFT-1・z’n)を計算する第2の逆
高速フーリエ変換計算手段と、以下に示す方程式を計算
する方程式計算手段とが縦続接続され、フォワードMD
CTの複素数部を計算する第2の計算手段と Zk=((−1)k+1・W8-1・W4N4k+1)・Z’k を具備し、連続するN個のサンプルからなる複数のブロ
ックにひとまとめにされ、請求項7記載のディジタル信
号の符号化装置によって符号化された複数のディジタル
サンプルx(n)で表されるディジタル信号(x)を復
号化することを特徴とするディジタル信号の復号化装
置。 - 【請求項13】 前記第2の計算手段は、 請求項7ないし請求項11のいずれかに記載の符号化装
置の複数のモジュールにそれぞれ等しい複数のモジュー
ルから構成されていることを特徴とする請求項12記載
のディジタル信号の復号化装置。 - 【請求項14】 c=cos((2π・(4n+1))
/4N)、s=sin((2π・(4n+1)/N)と
したときの偶数のサンプルx2nおよび奇数のサンプル
(xN/2−1−2n)に基づいて、以下に示す方程式
を虚数乗法によって計算する第3の計算手段と、 z’n=((x2n・c−xN/2−1−2n・s) +jzo-N/2・(x2n・s+xN/2−2n−1・c))・(c+js ) 式Z’(N/8−k)=FFT-1・z’nを計算するた
めに、前記式z’nに対して逆高速フーリエ変換FFT
-1を実行する逆高速フーリエ変換計算手段と、変換の奇
数係数Y2kおよび偶数係数(YN/2−2K−1)を
計算するために、虚数乗法を実行する第4の計算手段と
を具備し、連続するN個のサンプルからなる複数のブロ
ックにひとまとめにされた複数のディジタルサンプルx
(n)で表されるディジタル信号(x)を符号化および
復号化することを特徴とするディジタル信号の符号化お
よび復号化装置。 - 【請求項15】 前記第3の計算手段は、 偶数サンプル(xN/2−1/2n)および奇数サンプ
ルx2Nをそれぞれ記憶する第3および第4の記憶手段
と、 前記値cおよび前記値sをそれぞれ記憶する第5および
第6の記憶手段と、 前記偶数サンプルx2Nが第1の入力端に供給され、前
記値cが(−1)を乗算する乗算器を介して第2の入力
端に供給される第9の乗算手段と、 前記奇数サンプル(xN/2−1/2n)が第1の入力
端に供給され、前記値cが直接第2の入力端に供給され
る第10の乗算手段と、 前記偶数サンプルx2Nが第1の入力端に供給され、前
記値sが前記第6の記憶手段から直接第2の入力端に供
給される第11の乗算手段と、 前記奇数サンプル(xN/2−1/2n)が第1の入力
端に供給され、前記値sが前記第6の記憶手段から(−
1)を乗算する乗算器を介して第2の入力端に供給され
る第12の乗算手段と、 前記第9の乗算手段から出力される積と前記第11の乗
算手段から出力される積とが供給される第7の加算手段
と、 前記第10の乗算手段から出力される積と前記第12の
乗算手段から出力される積とが供給される第8の加算手
段と、 前記第7の加算手段から出力される加算結果が直接第1
の入力端に供給され、前記加算結果がディレイラインお
よびインバータを経て第2の入力端に供給される第9の
加算手段と、 前記第8の加算手段から出力される加算結果が直接第1
の入力端に供給され、前記加算結果がディレイラインを
経て第2の入力端に供給される第10の加算手段と、 前記奇数サンプル(xN/2−1/2n)がインバータ
を経て第1の入力端に供給され、前記第9の加算手段か
ら出力される加算結果が直接第2の入力端に供給され、
前記式z’nの実数部を出力する第11の加算手段と、 前記偶数サンプルx2Nがインバータを経て第1の入力
端に供給され、前記第10の加算手段から出力される加
算結果が直接第2の入力端に供給され、前記式z’nの
虚数部を出力する第12の加算手段と、 前記式z’nの実数部および虚数部をそれぞれ記憶する
RAMタイプの第7および第8の記憶手段とを具備する
ことを特徴とする請求項14記載のディジタル信号の符
号化および復号化装置。 - 【請求項16】 前記第3の計算手段は、 前記値cおよび前記値sをそれぞれ記憶する第9および
第10の記憶手段と、 前記式z’nに対する前記逆高速フーリエ変換FFT-1
・z’nの演算結果の実数部Re〔Z’N/8−k〕が
第1の入力端に供給され、前記値cが第2の入力端に供
給される第13の乗算手段と、 前記式z’nに対する前記逆高速フーリエ変換FFT-1
・z’nの演算結果の虚数部Im〔Z’N/8−k〕が
第1の入力端に供給され、前記値cが第2の入力端に供
給される第14の乗算手段と、 前記式z’nに対する前記逆高速フーリエ変換FFT-1
・z’nの演算結果の実数部Re〔Z’N/8+k〕が
第1の入力端に供給され、前記値cが第2の入力端に供
給される第15の乗算手段と、 前記式z’nに対する前記逆高速フーリエ変換FFT-1
・z’nの演算結果の虚数部Im〔Z’N/8+k〕が
第1の入力端に供給され、前記値cが第2の入力端に供
給される第16の乗算手段と、 前記式z’nに対する前記逆高速フーリエ変換FFT-1
・z’nの演算結果の実数部Re〔Z’N/8−k〕が
第1の入力端に供給され、前記値sが第2の入力端に供
給される第17の乗算手段と、 前記式z’nに対する前記逆高速フーリエ変換FFT-1
・z’nの演算結果の虚数部Im〔Z’N/8−k〕が
第1の入力端に供給され、前記値sが第2の入力端に供
給される第18の乗算手段と、 前記式z’nに対する前記逆高速フーリエ変換FFT-1
・z’nの演算結果の実数部Re〔Z’N/8+k〕が
第1の入力端に供給され、前記値sが第2の入力端に供
給される第19の乗算手段と、 前記式z’nに対する前記逆高速フーリエ変換FFT-1
・z’nの演算結果の虚数部Im〔Z’N/8+k〕が
第1の入力端に供給され、前記値sが第2の入力端に供
給される第20の乗算手段と、 前記第13の乗算手段の乗算結果が第1の入力端に供給
され、前記第14の乗算手段の乗算結果が第2の入力端
に供給される第13の加算手段と、 前記第17の乗算手段の乗算結果がインバータを経て第
1の入力端に供給され、前記第14の乗算手段の乗算結
果が第2の入力端に供給される第14の加算手段と、 前記第15の乗算手段の乗算結果が第1の入力端に供給
され、前記第20の乗算手段の乗算結果がインバータを
経て第2の入力端に供給される第15の加算手段と、 前記第19の乗算手段の乗算結果が第1の入力端に供給
され、前記第16の乗算手段の乗算結果が第2の入力端
に供給される第16の加算手段と、 前記第13および第14の加算手段のそれぞれの加算結
果が変換の奇数係数(YN/2−2k−1)として記憶
される第9の記憶手段と、 前記第15および第16の加算手段のそれぞれの加算結
果が変換の偶数係数Y2kとして記憶される第10の記
憶手段とを具備することを特徴とする請求項14記載の
ディジタル信号の符号化および復号化装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR9105064 | 1991-04-24 | ||
| FR9105064A FR2675969B1 (fr) | 1991-04-24 | 1991-04-24 | Procede et dispositif de codage-decodage d'un signal numerique. |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH05145428A true JPH05145428A (ja) | 1993-06-11 |
| JP3147984B2 JP3147984B2 (ja) | 2001-03-19 |
Family
ID=9412212
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP10705392A Expired - Lifetime JP3147984B2 (ja) | 1991-04-24 | 1992-04-24 | デジタル信号の符号化−復号化方法およびその装置 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US5363096A (ja) |
| EP (1) | EP0511095B1 (ja) |
| JP (1) | JP3147984B2 (ja) |
| CA (1) | CA2066471C (ja) |
| DE (1) | DE69209129T2 (ja) |
| FR (1) | FR2675969B1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101495027B1 (ko) * | 2012-05-11 | 2015-02-24 | 누보톤 테크놀로지 코포레이션 | 오디오 시스템 및 그 안에 클래스 d 증폭기를 가지는 집적 회로 칩 |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5486826A (en) * | 1994-05-19 | 1996-01-23 | Ps Venture 1 Llc | Method and apparatus for iterative compression of digital data |
| US5621403A (en) * | 1995-06-20 | 1997-04-15 | Programmed Logic Corporation | Data compression system with expanding window |
| US5913186A (en) * | 1996-03-25 | 1999-06-15 | Prometheus, Inc. | Discrete one dimensional signal processing apparatus and method using energy spreading coding |
| TW301103B (en) * | 1996-09-07 | 1997-03-21 | Nat Science Council | The time domain alias cancellation device and its signal processing method |
| SE9703849L (sv) | 1997-03-14 | 1998-09-15 | Ericsson Telefon Ab L M | Nedskalning av bilder |
| JP3279228B2 (ja) * | 1997-08-09 | 2002-04-30 | 日本電気株式会社 | 符号化音声復号装置 |
| SE9802286L (sv) * | 1998-06-26 | 1999-12-27 | Ericsson Telefon Ab L M | Effektiv nedskalning av DCT-komprimerade bilder |
| US6430529B1 (en) * | 1999-02-26 | 2002-08-06 | Sony Corporation | System and method for efficient time-domain aliasing cancellation |
| US6460061B1 (en) * | 1999-10-29 | 2002-10-01 | Xilinx Inc. | 2-dimensional discrete cosine transform using a polynomial transform |
| US7203717B1 (en) * | 1999-10-30 | 2007-04-10 | Stmicroelectronics Asia Pacific Pte. Ltd. | Fast modified discrete cosine transform method |
| US7987089B2 (en) * | 2006-07-31 | 2011-07-26 | Qualcomm Incorporated | Systems and methods for modifying a zero pad region of a windowed frame of an audio signal |
| US9653088B2 (en) * | 2007-06-13 | 2017-05-16 | Qualcomm Incorporated | Systems, methods, and apparatus for signal encoding using pitch-regularizing and non-pitch-regularizing coding |
| US8548815B2 (en) | 2007-09-19 | 2013-10-01 | Qualcomm Incorporated | Efficient design of MDCT / IMDCT filterbanks for speech and audio coding applications |
| RU2451998C2 (ru) * | 2007-09-19 | 2012-05-27 | Квэлкомм Инкорпорейтед | Эффективный способ проектирования набора фильтров для mdct/imdct в приложениях для кодирования речи и аудиосигналов |
| US20090099844A1 (en) * | 2007-10-16 | 2009-04-16 | Qualcomm Incorporated | Efficient implementation of analysis and synthesis filterbanks for mpeg aac and mpeg aac eld encoders/decoders |
| US10003377B1 (en) * | 2016-12-19 | 2018-06-19 | Cisco Technology, Inc. | Spread spectrum acoustic communication techniques |
| US10404319B1 (en) | 2018-07-09 | 2019-09-03 | Cisco Technology, Inc. | Fast correlation of prometheus orthonormal sets (PONS) for communications |
| US10396846B1 (en) | 2018-10-12 | 2019-08-27 | Cisco Technology, Inc. | Adaptive decoding of spread spectrum signals using multiple correlator peaks |
| US10601459B1 (en) | 2018-11-02 | 2020-03-24 | Cisco Technology, Inc. | Efficient handling of clock offset in spread spectrum decoders |
-
1991
- 1991-04-24 FR FR9105064A patent/FR2675969B1/fr not_active Expired - Lifetime
-
1992
- 1992-04-21 US US07/871,543 patent/US5363096A/en not_active Expired - Lifetime
- 1992-04-22 EP EP92401139A patent/EP0511095B1/fr not_active Expired - Lifetime
- 1992-04-22 DE DE69209129T patent/DE69209129T2/de not_active Expired - Lifetime
- 1992-04-24 CA CA002066471A patent/CA2066471C/en not_active Expired - Lifetime
- 1992-04-24 JP JP10705392A patent/JP3147984B2/ja not_active Expired - Lifetime
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101495027B1 (ko) * | 2012-05-11 | 2015-02-24 | 누보톤 테크놀로지 코포레이션 | 오디오 시스템 및 그 안에 클래스 d 증폭기를 가지는 집적 회로 칩 |
Also Published As
| Publication number | Publication date |
|---|---|
| US5363096A (en) | 1994-11-08 |
| FR2675969A1 (fr) | 1992-10-30 |
| EP0511095B1 (fr) | 1996-03-20 |
| EP0511095A1 (fr) | 1992-10-28 |
| FR2675969B1 (fr) | 1994-02-11 |
| JP3147984B2 (ja) | 2001-03-19 |
| DE69209129T2 (de) | 1996-11-07 |
| CA2066471A1 (en) | 1992-10-25 |
| CA2066471C (en) | 2002-07-16 |
| DE69209129D1 (de) | 1996-04-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3147984B2 (ja) | デジタル信号の符号化−復号化方法およびその装置 | |
| US5999656A (en) | Overlapped reversible transforms for unified lossless/lossy compression | |
| US6963842B2 (en) | Efficient system and method for converting between different transform-domain signal representations | |
| Stark | Wavelets and signal processing: an application-based introduction | |
| CN101930426B (zh) | 信号处理方法、数据处理方法及装置 | |
| JP4398979B2 (ja) | 変換表現に転換するか、または変換表現を逆変換するための装置および方法 | |
| JP3814611B2 (ja) | 時間離散オーディオサンプル値を処理する方法と装置 | |
| TW200923913A (en) | Encoding and/or decoding digital content | |
| EP0128298B1 (en) | Orthogonal transformer and apparatus operational thereby | |
| KR20070001115A (ko) | 복소수 값 데이터를 이용하는 오디오 신호 디코딩 | |
| JP3541680B2 (ja) | 音声音楽信号の符号化装置および復号装置 | |
| EP0590790A2 (en) | Modified DCT signal transforming system | |
| JPS6037514B2 (ja) | 2次元離散フ−リエ変換計算装置 | |
| RU2530327C2 (ru) | Способ обновления кодера посредством интерполяции фильтра | |
| JP2002245027A (ja) | フィルタリング処理方法およびフィルタリング処理装置 | |
| JP3185371B2 (ja) | 改良dctの順変換計算装置、逆変換計算装置及び改良dctの順変換計算方法 | |
| JPH06149861A (ja) | Dct及び逆dct演算装置並びにその演算方法 | |
| US7031404B2 (en) | Inverse DWT method and apparatus | |
| Dai et al. | Fast algorithm for modulated complex lapped transform | |
| US6292817B1 (en) | Discrete cosine transformation circuit | |
| JP4033906B2 (ja) | 時間領域エイリアシング相殺を用いる効率的奇数積重ね単一側波帯フィルタバンクを有する分析・合成濾波システム | |
| JP2003345779A (ja) | フィルタリング処理方法およびフィルタリング処理装置 | |
| WO1994019791A1 (en) | Improved filter for use in audio compression and decompression systems | |
| JP3547971B2 (ja) | 離散コサイン変換回路及びその動作方法 | |
| JP4293005B2 (ja) | 音声音楽信号の符号化装置および復号装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20001205 |
|
| 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 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090112 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090112 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100112 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110112 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110112 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120112 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130112 Year of fee payment: 12 |
|
| EXPY | Cancellation because of completion of term | ||
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130112 Year of fee payment: 12 |