[go: up one dir, main page]

JP4019824B2 - 波形生成装置及び方法並びに復号装置 - Google Patents

波形生成装置及び方法並びに復号装置 Download PDF

Info

Publication number
JP4019824B2
JP4019824B2 JP2002199070A JP2002199070A JP4019824B2 JP 4019824 B2 JP4019824 B2 JP 4019824B2 JP 2002199070 A JP2002199070 A JP 2002199070A JP 2002199070 A JP2002199070 A JP 2002199070A JP 4019824 B2 JP4019824 B2 JP 4019824B2
Authority
JP
Japan
Prior art keywords
sine wave
waveform
signal
output
value
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.)
Expired - Fee Related
Application number
JP2002199070A
Other languages
English (en)
Other versions
JP2004048101A (ja
Inventor
千智 劔持
晃 井上
正之 西口
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.)
Sony Corp
Original Assignee
Sony Corp
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 Sony Corp filed Critical Sony Corp
Priority to JP2002199070A priority Critical patent/JP4019824B2/ja
Priority to US10/516,819 priority patent/US7405499B2/en
Priority to EP03733581A priority patent/EP1538740A4/en
Priority to CNB038161567A priority patent/CN100477486C/zh
Priority to PCT/JP2003/008228 priority patent/WO2004006425A1/ja
Publication of JP2004048101A publication Critical patent/JP2004048101A/ja
Application granted granted Critical
Publication of JP4019824B2 publication Critical patent/JP4019824B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H7/00Instruments in which the tones are synthesised from a data store, e.g. computer organs
    • G10H7/08Instruments in which the tones are synthesised from a data store, e.g. computer organs by calculating functions or polynomial approximations to evaluate amplitudes at successive sample points of a tone waveform
    • G10H7/12Instruments in which the tones are synthesised from a data store, e.g. computer organs by calculating functions or polynomial approximations to evaluate amplitudes at successive sample points of a tone waveform by means of a recursive algorithm using one or more sets of parameters stored in a memory and the calculated amplitudes of one or more preceding sample points
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H7/00Instruments in which the tones are synthesised from a data store, e.g. computer organs
    • G10H7/08Instruments in which the tones are synthesised from a data store, e.g. computer organs by calculating functions or polynomial approximations to evaluate amplitudes at successive sample points of a tone waveform
    • G10H7/10Instruments in which the tones are synthesised from a data store, e.g. computer organs by calculating functions or polynomial approximations to evaluate amplitudes at successive sample points of a tone waveform using coefficients or parameters stored in a memory, e.g. Fourier coefficients
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03BGENERATION OF OSCILLATIONS, DIRECTLY OR BY FREQUENCY-CHANGING, BY CIRCUITS EMPLOYING ACTIVE ELEMENTS WHICH OPERATE IN A NON-SWITCHING MANNER; GENERATION OF NOISE BY SUCH CIRCUITS
    • H03B28/00Generation of oscillations by methods not covered by groups H03B5/00 - H03B27/00, including modification of the waveform to produce sinusoidal oscillations
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H2250/00Aspects of algorithms or signal processing methods without intrinsic musical character, yet specifically adapted for or used in electrophonic musical processing
    • G10H2250/055Filters for musical processing or musical effects; Filter responses, filter architecture, filter coefficients or control parameters therefor
    • G10H2250/111Impulse response, i.e. filters defined or specified by their temporal impulse response features, e.g. for echo or reverberation applications
    • G10H2250/121IIR impulse
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10HELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H2250/00Aspects of algorithms or signal processing methods without intrinsic musical character, yet specifically adapted for or used in electrophonic musical processing
    • G10H2250/131Mathematical functions for musical analysis, processing, synthesis or composition
    • G10H2250/215Transforms, i.e. mathematical transforms into domains appropriate for musical signal processing, coding or compression
    • G10H2250/221Cosine transform; DCT [discrete cosine transform], e.g. for use in lossy audio compression such as MP3
    • G10H2250/225MDCT [Modified discrete cosine transform], i.e. based on a DCT of overlapping data

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Algebra (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Power Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Electrophonic Musical Instruments (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、正弦波信号等の波形を生成するための波形生成装置及び方法並びに復号装置に関するものであり、例えば、正弦波テーブルを読み出すためのモジュロアドレッシングを行わずに正弦波等の周期的波形を生成するような波形生成装置及び方法並びに復号装置に関するものである。
【0002】
【従来の技術】
ディジタル的に正弦波等の波形を生成する技術は、信号符号化や信号合成等の種々の分野において用いられている。例えば、正弦波を決定する特徴量が周期的に与えられて一定長Nサンプルの正弦波を生成する技術として、予め正弦波の波高値が記憶されたテーブルを作成しておき、DSP(Digital Signal Processor:デジタル信号処理装置)を用いて、上記特徴量に基づきモジュロアドレッシングを行いながらテーブルを読み込み、振幅補正を行う方法が知られている。
【0003】
以下、周期的に変化する正弦波生成装置の従来の構成例について、図面を参照しながら説明する。図13は従来の正弦波生成装置のブロック図、図14は該正弦波生成装置の動作を説明するためのフローチャートの一例である。
【0004】
図13において、周期的に変化する特徴量s101が入力手段101により取り込まれ、特徴量検出手段102へ渡される。特徴量検出手段102では、図14のステップS11に示すように、渡された特徴量をもとに正弦波テーブル読み込み方法特徴量s102と正弦波振幅補正量s103が算出される。例えば、読み込み方法特徴量s102は、Ω、Ωの関係とΦにより、読み込み開始点、及び何点置きにテーブルを読み込むかが算出される。同様に振幅補正量s103もA、Aの割合より算出される。特徴量検出手段102で算出される読み込み方法特徴量とは、生成する正弦波と正弦波テーブルの間の「位相のずれ」、「周波数の関係」などである。
【0005】
特徴量検出手段102で算出された読み込み方法特徴量s2が読み込み位置算出手段103に送られて正弦波テーブル読み込み開始位置が指示され、読み込み位置算出手段103により正弦波テーブル読み込み位置信号s104として正弦波テーブル読み込み手段104に送られ、ステップS12に示すように、正弦波テーブルより値を読み込む。
【0006】
正弦波テーブルより読み込んだ値(読み込み信号s105)を正弦波振幅補正手段105に渡し、ステップS13において、正弦波振幅補正手段105は、正弦波テーブル読み込み手段4より渡された読み込み信号s105と特徴量検出手段102より指示される正弦波振幅補正量s103により出力する正弦波の振幅を補正する。また、ステップS14で、得られた正弦波信号を蓄積する。
【0007】
ステップS15では、処理すべき全データの読み込みが完了したか否かを判別し、NOのときはステップS16に進んで読み込み位置(readIndex)を更新した後、上記ステップS12に戻る。
【0008】
ここで、上記読み込み位置(readIndex)の更新は、一時点前の読み込み位置信号s104(readIndex)に、特徴量検出手段102より指示される読み込み位置更新量(readInc)をモジュロ加算し、新たな読み込み位置信号s104(readIndex)を算出することにより行われ、これは、上記正弦波テーブルの読み込み位置最大値をindexMaxとするとき、
readIndex= (readIndex+readInc)%indexMax(1)
の式により表される。なお、a%bは、aをbで割った余り(剰余)を示す。この更新された新たな読み込み位置(readIndex)は、読み込み位置算出手段103へ渡される。すなわち、読み込み位置信号s104は、一時点前の読み込み位置信号s104と読み込み更新量と正弦波テーブルのサンプル数Nより、算出される。
【0009】
図14に戻って、上記ステップS12からS14までの操作を、出力正弦波のサンプル数Nが揃うまで、すなわちステップS15でYES(全データ読み込み完了)と判別されるまで繰り返し行う。出力サンプル数が揃った時点で出力手段6へ正弦波信号をわたす。
【0010】
【発明が解決しようとする課題】
ところで、上述したような方法により正弦波を生成する場合、テーブルの読み込みを行うためにモジュロアドレッシングを行う必要、すなわちアドレス計算のためにモジュロ演算を行う必要がある。しかしながら、DSP(デジタル信号処理装置)には種々の仕様のものが提供されており、モジュロアドレッシングを行うのに制限があるデジタル信号処理装置も存在している。このようなモジュロアドレッシングを行うのに制限があるデジタル信号処理装置であっても、他の機能が優れているために採用されることも多い。
【0011】
例えば、上述のような正弦波生成装置は、オーディオ信号の復号装置等に用いられることが多く、上記DSPは、正弦波生成のみに用いられるよりも、復号処理全体に用いられることが多い。このような復号装置において、例えば高い積和演算の処理能力が求められる場合、積和演算能力が優れたDSPであれば、他の例えばモジュロアドレッシングの能力に制限があっても採用されることは十分考えられる。
【0012】
このような場合には、モジュロアドレッシングを行わずに正弦波生成を行うことが必要とされる。
【0013】
本発明は、上述のような実情に鑑みて提案されたものであって、デジタルオーディオ信号処理装置等において例えば周期的に正弦波等の周期的波形を決定する特徴量が与えられる場合に、正弦波テーブルを参照する際のモジュロアドレッシングを行わずに正弦波等の周期的波形の生成を行うこと可能とするような波形生成装置及び方法、並びにこのような波形生成が適用される復号装置を提供することを目的とする。
【0014】
【課題を解決するための手段】
本発明に係る波形生成装置及び方法は、上述の課題を解決するために、入力された特徴量に基づいて周期的波形を生成する際に、上記入力された特徴量を検出し、検出された上記特徴量に基づき少なくとも2点のサンプル点を初期値として漸化式を演算することにより上記周期的波形を発生し、発生された上記周期的波形を出力するものであり、上記特徴量として、複数の波形に対応する複数の特徴量が与えられ、各特徴量ごとに周期的波形の信号を生成してそれらを加算したものを最終出力波形として出力し、上記複数の特徴量のうち、上記最終出力波形への寄与分が少ないと判断された波形の生成、加算を行わないことを特徴としている。
【0015】
ここで、上記周期的波形は正弦波であり、任意のn時点の正弦波信号の値をy[n]とするとき、上記特徴量として出力正弦波の位相Φ及び周波数Ωが与えられるとき、初期値y[0]、y[1]として、
[0]=Asin(Φ)
[1]=Asin(Ω+Φ)
を用い、n+2の時点の正弦波信号の値y[n+2]を、n+1の時点の正弦波信号の値 y[n+1]と、n時点の正弦波信号の値 y[n]とにより、
[n+2]=2×Acos(Ω)×y[n+1]−y[n]
と表される式を上記漸化式として用いることが挙げられる。
【0016】
また、時系列信号を符号化フレーム毎に正弦波分析することで得られた特徴量を含む符号化データが供給される復号装置の正弦波合成に適用される波形生成において、上記特徴量は上記符号化フレーム毎に周期的に与えられ、上記生成された正弦波信号を蓄積手段に上記1フレーム分蓄積し、蓄積された正弦波信号を上記出力手段を介して出力することが挙げられる。
【0017】
また、本発明に係る復号装置は、時系列信号を符号化フレーム毎に正弦波分析することで得られた特徴量を含む符号化データが供給される復号装置において、上記符号化データ中の特徴量を分離する分離手段と、上記分離された特徴量を検出する検出手段と、上記検出手段により検出された上記特徴量に基づき少なくとも2点のサンプル点を初期値として漸化式を演算することにより上記周期的波形を発生する発振手段と、上記発振手段から発生された上記周期的波形を出力する出力手段とを備える波形合成部を有し、上記特徴量として、複数の波形に対応する複数の特徴量が与えられ、各特徴量ごとに周期的波形の信号を生成してそれらを加算したものを最終出力波形として出力し、上記複数の特徴量のうち、上記最終出力波形への寄与分が少ないと判断された波形の生成、加算を行わないことを特徴とする。
【0018】
ここで、本発明によれば、周期的に正弦波を決定する特徴量が変化する正弦波を生成する際に、特徴量を取り込み、取り込まれた特徴量を検出し、検出された正弦波の特徴量と、上記特徴量を初期設定として発振器(正弦波生成手段)を用いて正弦波信号を生成し、生成された正弦波信号を蓄積し、蓄積された正弦波を出力する。
【0019】
【発明の実施の形態】
以下、本発明に係る好ましい実施の形態について、図面を参照しながら説明する。
【0020】
図1は、本発明に係る第1の実施の形態となる正弦波生成装置の概略構成を示すブロック図である。また、図2は、図1に示す正弦波生成装置の動作を説明するためのフローチャートである。
【0021】
先ず、図1に示す入力手段1により、周期的に変化する特徴量s1が取り込まれ、特徴量検出手段2へ送られて、図2のステップS21に示すように、特徴量検出が行われる。この特徴量検出手段2に渡された特徴量s1をもとに、生成する正弦波の振幅、位相及び周波数が算出され、これらの振幅、位相及び周波数を示す信号s2が初期設定手段3に送られる。初期設定手段3では、特徴量検出手段2より指示された振幅、位相及び周波数により正弦波の最初の2点を算出し、上記指示された値と併せて、初期値(発振器初期化信号s4)として発振器(正弦波生成手段)4に送られることにより、ステップS22に示す初期値設定がなされる。この初期値設定は、例えば、正弦波テーブルを参照することにより行われる。またこのとき、出力波形のサンプル点の位置を示す変数nが初期化(n=0)される。
【0022】
発振器(正弦波生成手段)4は、ステップS23にて、初期設定手段3より指示される初期値に従い正弦波信号を生成する。この正弦波信号の生成は、後述するように、例えばIIRフィルタ演算を行うことにより、モジュロアドレッシングを行わずに実現している。
【0023】
発振器(正弦波生成手段)4により生成された出力信号は、正弦波信号s5として正弦波蓄積手段5に送られる。正弦波蓄積手段5は、ステップS24において、発振器(正弦波生成手段)4からの正弦波信号s5を蓄積する。
【0024】
次のステップS25では、上記変数nが出力正弦波信号のサンプル点数N以上となったか否かを判別しており、NOと判別されたときは、ステップS26で上記変数nをインクリメント(n=n+1)した後にステップS23に戻り、上記ステップS23、S24の処理をnがNになるまで(生成波形データが一定長Nサンプル集まるまで)繰り返す。ステップS25でYES(nがNになった)と判別されると、正弦波蓄積手段5により蓄積された正弦波信号s6を出力手段6より出力する。
【0025】
次に、発振器(正弦波生成手段)4によるモジュロアドレッシングを行わない正弦波生成の一例について説明する。以下の説明では、出力正弦波信号のn時点(nサンプル目)の値をy[n]と表す。
【0026】
図3は、発振器(正弦波生成手段)4の構成例としてのIIRによる正弦波生成回路を示す。すなわち、例えば、2次のIIRフィルタ処理を用いるとき、出力正弦波信号の初期値y[0]、y[1]として、
[0]=Asin(Φ) (2)
[1]=Asin(Ω+Φ) (3)
が指示される。これらの式(2)、(3)において、Φは出力正弦波の位相、Ωは出力正弦波の周波数であり、上記特徴量として例えば符号化フレーム周期毎に与えられる。また、n+2の時点の正弦波信号の値y[n+2]は、n+1の時点の正弦波信号の値y[n+1]と、n時点の正弦波信号の値y[n]とを用いて、
[n+2]=2×Acos(Ω)×y[n+1]−y[n] (4)
の漸化式により表され、この式(4)により正弦波信号の各サンプル点の値が逐次算出される。算出された信号は正弦波信号s5として正弦波蓄積手段5に渡される。
【0027】
上記式(4)は以下のように導かれる。先ず、sin関数のz変換は、次の式(5)のように与えられる。
【0028】
【数1】
Figure 0004019824
【0029】
この(5)式を展開して、逆z変換を行うと、
[n+2]−2cos(ωT)×y[n+1]+y[n]
=2sin(ωT)×x[n] (6)
ここで、図3のIIRフィルタには入力がないことから、x[n]の項は常に0となり、これを整理すると、上記式(4)が導かれる。
【0030】
次に、上述したような正弦波合成装置が用いられる信号復号装置の例について説明する。この復号装置の例は、オーディオ信号の少なくとも一部を正弦波合成により復号するものであり、符号化側では、例えばGHA(Generalized Harmonic Analysis:一般調和解析)により符号化を行うものである。以下、符号化装置、復号装置の順に説明する。
【0031】
図4は、GHA(Generalized Harmonic Analysis:一般調和解析)を用いたオーディオ信号(音響時系列信号)の符号化装置の構成例を示すブロック図である。この図4において、ステレオオーディオ信号の左右チャンネル(Lch、Rch)の信号が入力端子121L、121Rにそれぞれ入力されている。これらの左右チャンネルの構成は同様であるため、左チャンネルについて説明し、右チャンネルの説明を省略する。入力端子121L(あるいは121R)に供給された左(あるいは右)チャンネルのオーディオ信号は、所定の符号化フレーム周期(例えば2048サンプル)毎に符号化処理が行われる。
【0032】
入力端子121L(あるいは121R)からの入力信号は、例えば16-band帯域分割部122に送られて、16のバンド(周波数帯域)に周波数分割され、それぞれの帯域の信号(上記符号化フレーム周期毎に例えば1バンド当たり128サンプル)がGHA分析部123に送られる。各バンドのGHA分析部123では、音声の母音部やドラム等のように周期性の強い信号を抽出・分離して、別途符号化を行い、特徴量(例えば周波数、振幅、位相等)の情報を適応ビットアロケーション/ステレオコーディング部126に送っている。また、各バンド(周波数帯域)のGHA分析部123で抽出された周期性の強い信号を、16-band帯域分割部122からの対応するバンドの信号から減算しており、その残差信号を、ゲイン制御部124を介してMDCT(Modified Discrete Cosine Transform:変形離散コサイン変換)部125に送ってMDCT処理している。各バンドのMDCT部125からの出力(MDCT係数データ)は、適応ビットアロケーション/ステレオコーディング部126に送られ、この適応ビットアロケーション/ステレオコーディング部126においては、上記各バンドのGHA分析部123からの特徴量情報と共に、上記左右チャンネルのそれぞれ16バンドについて、適応ビット割り当てを伴うステレオコーディング処理が施される。適応ビットアロケーション/ステレオコーディング部126からの出力は、量子化/符号化部127に送られて、量子化及び符号化が施され、出力端子128を介してエンコード出力信号(符号列)が取り出される。
【0033】
図5は、上記図4に示したような符号化装置からのエンコード出力信号を復号(デコード)する復号装置の構成例を示すブロック図である。この図5の入力端子131には、上記図4の出力端子128からのエンコード出力信号(符号列)が入力されているものとする。
【0034】
図5の入力端子131からの信号は、復号化/逆量子化部132に送られて、上記図4の量子化/符号化部127における処理の逆処理となる復号処理及び逆量子化処理が施される。復号化/逆量子化部132からは、上述したステレオ左右チャンネルの各バンド毎に、上記特徴量及びMDCT係数が取り出され、上記特徴量は各バンドのGHA合成部133に、上記MDCT係数は各バンドのIMDCT(逆MDCT)部134にそれぞれ送られる。GHA合成部133では、上記特徴量(例えば周波数、振幅、位相等)に応じた正弦波を生成し、当該バンド内で複数組の特徴量が与えられて複数の正弦波が生成される場合にはこれらを合成し、出力する。IMDCT部134では、上記MDCT係数を逆変換して上記図4のMDCT部125の入力に相当する信号を出力し、ゲイン補償部135に送って上記図4のゲイン制御部124の逆の処理を施して出力する。各バンドのゲイン制御部124からの出力(上記残差信号)と当該バンドの上記GHA合成部133からの正弦波出力とを加算し、各バンド毎に得られた加算出力の16バンド分を16-band帯域合成部136に送って合成することにより、左チャンネル(あるいは右チャンネル)のオーディオ信号を出力端子137L(あるいは137R)から取り出している。
【0035】
本発明に係る波形生成装置及び方法は、例えば図5の復号装置のGHA部133に適用することができる。なお、この他、正弦波等の周期的波形の合成を伴う種々の復号装置の波形合成部等に本発明に係る波形生成装置及び方法を適用することができることは勿論である。
【0036】
次に、本発明の第2の実施の形態について、図6及び図7を参照しながら説明する。図6は、本発明の第2の実施の形態となる正弦波生成装置の概略構成を示すブロック図であり、図7は図6の装置の動作を説明するためのフローチャートである。この第2の実施の形態では、複数(例えばk個)の発振器(正弦波生成手段)を用いて並列処理により正弦波生成を行っている。
【0037】
図6の入力手段1から特徴量s1が入力され、特徴量検出手段2へ渡されて、図7のステップS31に示す特徴量検出が行われる。特徴量検出手段2により、生成する正弦波の振幅、位相及び周波数が算出され、信号s2として初期設定手段3へと渡される。初期設定手段3は、特徴量検出手段2より指示された振幅、位相及び周波数により、正弦波の最初の(2×k)点を算出し、上記指示された値と併せて発振器に初期値s4として渡す。ここでkは並列化する本数である(ステップS32)。
【0038】
例えば、2並列(k=2)で処理を行う場合は、最初の(2×2)=4点を正弦波テーブルより参照する。
【0039】
次に、初期設定手段3より指示される上記初期値に従い、複数(例えばk個)の発振器(正弦波生成手段)4を発振させて正弦波信号を生成する(ステップS33)。ここで、複数の各発振器4からは、1つの正弦波信号の時系列順の各サンプル点が順次発生される。
【0040】
例えば、2次のIIRフィルタ処理を2並列で用いるとき初期値として
[0]=Asin(Φ) (7)
[1]=Asin(Ω+Φ) (8)
[2]=Asin(2×Ω+Φ) (9)
[3]=Asin(3×Ω+Φ) (10)
が指示され
[n+4]=2×Acos(2×Ω)×y[n+2]−y[n] (11)
[n+5]=2×Acos(2×Ω)×y[n+3]−y[n+1] (12)
と2並列で生成される。ここでは生成される正弦波を入れ子で作成したが特にこのようにする必要はない。理論的には上記式(5)において、z=ejωTであることを考慮して、z→z、ω→2×ω の変換を行うことにより、sin関数のz変換式は、次の式(13)のように与えられる。
【0041】
【数2】
Figure 0004019824
【0042】
この(13)式を展開して、逆z変換を行うと、
[n+4]−2cos(2ωT)×y[n+2]+y[n]
=2sin(2ωT)×x[n+2] (14)
ここで、入力がないことから、x[n+2]の項は常に0となり、これを整理すると上記式(11)が導かれる。また、nを1進めると上記式(12)が導かれる。算出された信号は正弦波信号s5として正弦波蓄積手段5に渡される。
【0043】
次に、図7のステップS34では、発振器(正弦波生成手段)4により指示された正弦波信号s5を正弦波蓄積手段により蓄積する。
【0044】
次のステップS35では、上記変数nが出力正弦波信号のサンプル点数N以上となったか否かを判別し、NOのときはステップS36でnをkだけインクリメント(n=n+k)した後にステップS33に戻り、上記ステップS33、S34の処理をnがNになるまで(生成波形データが一定長Nサンプル集まるまで)繰り返す。ステップS35でYES(nがNになった)と判別されると、ステップS37に進んで、正弦波蓄積手段5により蓄積された正弦波信号s7を正弦波合成手段7に送り、一つ或いは複数の正弦波に合成する。合成した正弦波信号s6は出力手段6より出力する。
【0045】
この第2の実施の形態に示すような正弦波生成装置は、例えばDSP(Digital Signal Processor:デジタル信号処理装置)が複数の積和演算部を有している場合に、各積和演算部で各発振器の動作を行わせることによって容易に実現でき、処理量や処理速度を向上できるのみならず、漸化式で次のサンプル点を順次算出することによる蓄積誤差を小さくできるという利点もある。
【0046】
次に、本発明の第3の実施の形態について、図8及び図9を参照しながら説明する。図8は、本発明の第3の実施の形態となる正弦波生成装置の概略構成を示すブロック図であり、図9は図8の装置の動作を説明するためのフローチャートである。この第3の実施の形態では、発振器(正弦波生成手段)を時分割的に用いて複数種類(複数本)の正弦波を生成し、合成している。
【0047】
図8の構成において、正弦波生成部20は、上述した図1に示すような正弦波生成装置と同様な動作を行うものであり、対応する部分に同じ指示符号を付している。
【0048】
図8の入力手段1から特徴量s1が入力され、特徴量検出手段2へ渡されて、図9のステップS41に示す特徴量検出が行われる。特徴量検出手段2により、生成する正弦波の振幅、位相及び周波数が算出され、信号s2として正弦波生成部20の初期設定手段3へと渡される。正弦波生成部20において、発振器(正弦波生成手段)4で生成された正弦波信号s5は、正弦波蓄積手段5に蓄積され、正弦波信号s6として出力される。正弦波生成部20としては、上記図1(あるいは図6)に示した構成の装置を用いることができるため、詳細な説明を省略する。また、ここまでの動作は、図9のステップS41からS46までに相当し、これは、上記図2のステップS21からS26までと同様であるため、詳細な説明を省略する。
【0049】
ここで、この第3の実施の形態の場合には、生成する正弦波の本数(種類数)が複数(N)の場合を想定しており、図9の特徴量検出のステップS41に先行して、生成する正弦波の本数を示す変数n’を初期化(n’=0)するステップS40を設けている。また、図8の正弦波蓄積手段5に蓄積された正弦波信号s6は、波形合成手段8に送られ、波形合成手段8からの出力が波形蓄積手段9に送られると共に、波形蓄積手段9からの出力が波形合成手段8に送られるようになっている。すなわち、波形合成手段8は、正弦波生成部20の正弦波蓄積手段5からの正弦波信号s6と、波形蓄積手段9より渡される蓄積合成波信号s9とを、適当な合成法により合成し、得られた合成波信号s8を波形蓄積手段9に送っている。波形合成手段8において、例えば、重み付け無しに合成する場合は、合成波形出力z[k]を、
[k]=z[k]+y[k] (0≦k<N
の式により求めればよい。図9のステップS47は、これらの波形合成手段8及び波形蓄積手段9による波形合成・蓄積の動作を示している。
【0050】
次のステップS25では、上記変数n’が上記生成する正弦波の本数N以上か(n'≧N)否かを判別し、NOと判別されたときは、ステップS49で上記変数n'をインクリメント(n'=n'+1)した後にステップS41に戻り、上記ステップS41〜S47の処理を、上記変数n'がNになるまで(生成された正弦波が上記N本揃うまで)繰り返す。ステップS48でYES(n'がNに達した)と判別されると、正弦波蓄積手段9により蓄積された合成正弦波を出力手段6より出力する。
【0051】
次に、本発明の第4の実施の形態について、図10及び図11を参照しながら説明する。図10は、本発明の第4の実施の形態となる正弦波生成装置の概略構成を示すブロック図であり、図11は図10の装置の動作を説明するためのフローチャートである。この第4の実施の形態は、上記図8及び図9と共に説明した第3の実施の形態の一部を変更したものであり、複数種類(複数本)の正弦波を生成し、合成する場合に、合成波形の振幅に比べて充分小さい振幅の正弦波については、当該正弦波の生成・合成を省略している。
【0052】
すなわち、図10において、上記図8の構成に比べて、特徴量検出手段2と正弦波生成部20の初期設定手段3との間に特徴量判別手段10が設けられている点が異なっており、この特徴量判別手段10による判別結果により正弦波生成部20での動作を制御している。また、図11のフローチャートにおいては、上記ず9のフローチャートと比べて、ステップS41の特徴検出とステップS42の初期値設定との間に、現時点で生成しようとする正弦波の寄与が充分小さいか否かを判別するステップS50が挿入され、NOのときはステップS42に進み、YESのときはステップS48に飛ぶ点が異なっている。図10、図11の他の部分は、上記図8、図9と同様であり、対応する部分には同じ指示符号を付して説明を省略する。
【0053】
このような構成の第4の実施の形態において、複数の正弦波を生成して合成する場合に、図10の特徴量判別手段10は、特徴量に基づく特徴量検出手段2からの1つの正弦波についての振幅、位相、周波数情報の信号s2を分析し、図11のステップS50に示すように当該正弦波を生成・合成するかしないかを判別して、信号s2'を正弦波生成部20の初期設定手段3に送っている。これは、当該正弦波が、合成して得られる最終的な合成波への寄与が充分小さい場合、例えば合成波の振幅に比べて当該正弦波の振幅が充分小さい場合には、当該正弦波の生成・合成を行わない。このとき、特徴量判別手段10としては、特徴量検出手段2からの上記信号s2を正弦波生成部20の初期設定手段3に送らないようにしたり、あるいは正弦波生成部20での正弦波生成動作を停止させるように制御すること等が挙げられる。図11においては、ステップS50でYES(当該正弦波を生成・合成しない)と判別されたとき、ステップS42からS47までの正弦波生成動作を飛ばして、ステップS48に進んでいる。
【0054】
このような第4の実施の形態によれば、不要な処理を省略できるため、処理能力に余裕が生じ、速度の向上や他の部分への処理パワーの振り分け等が可能となる。
【0055】
ところで、上述した実施の形態において、発振器(正弦波生成手段)4は一定長Nを出力する間に充分な精度が保証されるものであれば任意のものを用いることが可能である。例えば、長さ128点、振幅1、初期位相0で1kHzの正弦波を上記図3に示したような2次のIIRフィルタ処理を用いて求めると、初期値はそれぞれ
[0]=sin(0)
[1]=sin(2π/128)
となり、求める正弦波は、
[n+2]=2cos(2π/128)×y[n+1]−y[n]
と求められる。この方法で求めた正弦波を図12に示す。なお、sin関数を用いてテーブル参照等により算出された正弦波も、図12の正弦波曲線に略々重なって示され、実用上充分な上記2次のIIRフィルタ処理を用いて求めた正弦波は、実用上充分な精度が得られている。
【0056】
なお、本発明は上述した実施の形態に限定されるものではなく、例えば、正弦波以外の各種周期的波形の生成にも容易に本発明を適用可能である。
【0057】
【発明の効果】
本発明によれば、入力された特徴量に基づいて周期的波形を生成する際に、上記入力された特徴量を検出し、検出された上記特徴量に基づき少なくとも2点のサンプル点を初期値として漸化式を演算することにより上記周期的波形を発生し、発生された上記周期的波形を出力することにより、モジュロアドレッシングをすることなく正弦波等の周期的波形を円滑に生成することができる。
【0058】
また、本発明によれば、時系列信号を符号化フレーム毎に正弦波分析することで得られた特徴量を含む符号化データが供給される復号装置において、上記符号化データ中の特徴量を分離する分離手段と、上記分離された特徴量を検出する検出手段と、上記検出手段により検出された上記特徴量に基づき少なくとも2点のサンプル点を初期値として漸化式を演算することにより上記周期的波形を発生する発振手段と、上記発振手段から発生された上記周期的波形を出力する出力手段とを備える波形合成部を有することにより、復号装置を実現するためのDSP等の能力に制限があっても、特にモジュロアドレッシング能力が低くても、正弦波等の周期的波形の生成が円滑に行えるため、モジュロアドレッシングを不得意とする装置においてモジュロアドレッシングをすることなく正弦波の生成を実現でき、復号装置全体の機能実現に好適なDSPを採用することができ、設計の自由度が高まる。
【図面の簡単な説明】
【図1】本発明の第1の実施の形態となる正弦波生成装置の概略構成を示すブロック図である。
【図2】本発明の第1の実施の形態の動作を説明するためのフローチャートである。
【図3】発振器(正弦波生成手段)の構成例を示すブロック図である。
【図4】オーディオ信号(音響時系列信号)の符号化装置の一例の概略構成を示すブロック図である。
【図5】本発明の実施の形態が適用される復号装置の一例の概略構成を示すブロック図である。
【図6】本発明の第2の実施の形態となる正弦波生成装置の概略構成を示すブロック図である。
【図7】本発明の第2の実施の形態の動作を説明するためのフローチャートである。
【図8】本発明の第3の実施の形態となる正弦波生成装置の概略構成を示すブロック図である。
【図9】本発明の第3の実施の形態の動作を説明するためのフローチャートである。
【図10】本発明の第4の実施の形態となる正弦波生成装置の概略構成を示すブロック図である。
【図11】本発明の第4の実施の形態の動作を説明するためのフローチャートである。
【図12】発振器(正弦波生成手段)により生成される正弦波を示す図である。
【図13】従来の正弦波生成装置の構成例を示すブロック図である。
【図14】従来の正弦波生成装置の動作を説明するためのフローチャートである。
【符号の説明】
1 入力手段、 2 特徴量検出手段、 3 初期設定手段、 4 発振器(正弦波生成手段)、 5 正弦波蓄積手段、 6 出力手段、 7 正弦波合成手段、 8 波形合成手段、 9 波形蓄積手段、 20 正弦波生成部、 133 GHA合成部

Claims (9)

  1. 入力された特徴量に基づいて周期的波形を生成する波形生成装置において、
    上記入力された特徴量を検出する検出手段と、
    上記検出手段により検出された上記特徴量に基づき少なくとも2点のサンプル点を初期値として漸化式を演算することにより上記周期的波形を発生する発振手段と、
    上記発振手段から発生された上記周期的波形を出力する出力手段とを有し、
    上記特徴量として、複数の波形に対応する複数の特徴量が与えられ、各特徴量ごとに周期的波形の信号を生成してそれらを加算したものを最終出力波形として出力し、上記複数の特徴量のうち、上記最終出力波形への寄与分が少ないと判断された波形の生成、加算を行わないことを特徴とする波形生成装置。
  2. 上記周期的波形は正弦波であり、任意のn時点の正弦波信号の値をy[n]とするとき、上記特徴量として出力正弦波の位相Φ及び周波数Ωが与えられるとき、初期値y[0]、y[1]として、
    [0] = Asin(Φ)
    [1] = Asin(Ω+Φ)
    を用い、n+2の時点の正弦波信号の値y[n+2]を、n+1の時点の正弦波信号の値 y[n+1]と、n時点の正弦波信号の値 y[n]とにより、
    [n+2] = 2×Acos(Ω)×y[n+1]−y[n]
    と表される式を上記漸化式として用いることを特徴とする請求項1記載の波形生成装置。
  3. 時系列信号を符号化フレーム毎に正弦波分析することで得られた特徴量を含む符号化データが供給される復号装置の正弦波合成に用いられる波形生成装置であって、
    上記特徴量は上記符号化フレーム毎に周期的に与えられ、上記発振手段からの正弦波信号を蓄積手段に上記1フレーム分蓄積し、蓄積された正弦波信号を上記出力手段を介して出力することを特徴とする請求項1記載の波形生成装置。
  4. 上記発振手段を複数用い、上記周期的波形の時系列順の各サンプル点を順次上記複数の発振手段により発生させることを特徴とする請求項1記載の波形生成装置。
  5. 入力された特徴量に基づいて周期的波形を生成する波形生成方法において、
    上記入力された特徴量を検出する検出工程と、
    上記検出手段により検出された上記特徴量に基づき少なくとも2点のサンプル点を初期値として漸化式を演算することにより上記周期的波形を発生する波形生成工程と、
    上記波形生成工程により発生された上記周期的波形を出力する出力工程とを有し、
    上記特徴量として、複数の波形に対応する複数の特徴量が与えられ、各特徴量ごとに周期的波形の信号を生成してそれらを加算したものを最終出力波形として出力し、上記複数の特徴量のうち、上記最終出力波形への寄与分が少ないと判断された波形の生成、加算を行わないことを特徴とする波形生成方法。
  6. 上記周期的波形は正弦波であり、任意のn時点の正弦波信号の値をy[n]とするとき、上記特徴量として出力正弦波の位相Φ及び周波数Ωが与えられるとき、初期値y[0]、y[1]として、
    [0] = Asin(Φ)
    [1] = Asin(Ω+Φ)
    を用い、n+2の時点の正弦波信号の値y[n+2]を、n+1の時点の正弦波信号の値 y[n+1]と、n時点の正弦波信号の値 y[n]とにより、
    [n+2] = 2×Acos(Ω)×y[n+1]−y[n]
    と表される式を上記漸化式として用いることを特徴とする請求項記載の波形生成方法。
  7. 時系列信号を符号化フレーム毎に正弦波分析することで得られた特徴量を含む符号化データが供給される復号装置の正弦波合成に用いられる波形生成方法であって、
    上記特徴量は上記符号化フレーム毎に周期的に与えられ、上記波形生成工程からの正弦波信号を上記1フレーム分蓄積し、蓄積された正弦波信号を出力することを特徴とする請求項記載の波形生成方法。
  8. 時系列信号を符号化フレーム毎に正弦波分析することで得られた特徴量を含む符号化データが供給される復号装置において、
    上記符号化データ中の特徴量を分離する分離手段と、
    上記分離された特徴量を検出する検出手段と、
    上記検出手段により検出された上記特徴量に基づき少なくとも2点のサンプル点を初期値として漸化式を演算することにより上記周期的波形を発生する発振手段と、
    上記発振手段から発生された上記周期的波形を出力する出力手段と
    を備える波形合成部を有し、
    上記特徴量として、複数の波形に対応する複数の特徴量が与えられ、各特徴量ごとに周期的波形の信号を生成してそれらを加算したものを最終出力波形として出力し、上記複数の特徴量のうち、上記最終出力波形への寄与分が少ないと判断された波形の生成、加算を行わないことを特徴とする復号装置。
  9. 上記周期的波形は正弦波であり、任意のn時点の正弦波信号の値をy[n]とするとき、上記特徴量として出力正弦波の位相Φ及び周波数Ωが与えられるとき、初期値y[0]、y[1]として、
    [0] = Asin(Φ)
    [1] = Asin(Ω+Φ)
    を用い、n+2の時点の正弦波信号の値y[n+2]を、n+1の時点の正弦波信号の値 y[n+1]と、n時点の正弦波信号の値 y[n]とにより、
    [n+2] = 2×Acos(Ω)×y[n+1]−y[n]
    と表される式を上記漸化式として用いることを特徴とする請求項記載の復号装置。
JP2002199070A 2002-07-08 2002-07-08 波形生成装置及び方法並びに復号装置 Expired - Fee Related JP4019824B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2002199070A JP4019824B2 (ja) 2002-07-08 2002-07-08 波形生成装置及び方法並びに復号装置
US10/516,819 US7405499B2 (en) 2002-07-08 2003-06-27 Waveform generating apparatus, waveform generating method, and decoder
EP03733581A EP1538740A4 (en) 2002-07-08 2003-06-27 DEVICE AND METHOD FOR GENERATING WAVEFORMS AND DECODER
CNB038161567A CN100477486C (zh) 2002-07-08 2003-06-27 波形生成装置及方法和解码装置
PCT/JP2003/008228 WO2004006425A1 (ja) 2002-07-08 2003-06-27 波形生成装置及び方法並びに復号装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002199070A JP4019824B2 (ja) 2002-07-08 2002-07-08 波形生成装置及び方法並びに復号装置

Publications (2)

Publication Number Publication Date
JP2004048101A JP2004048101A (ja) 2004-02-12
JP4019824B2 true JP4019824B2 (ja) 2007-12-12

Family

ID=30112446

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002199070A Expired - Fee Related JP4019824B2 (ja) 2002-07-08 2002-07-08 波形生成装置及び方法並びに復号装置

Country Status (5)

Country Link
US (1) US7405499B2 (ja)
EP (1) EP1538740A4 (ja)
JP (1) JP4019824B2 (ja)
CN (1) CN100477486C (ja)
WO (1) WO2004006425A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4640321B2 (ja) * 2006-11-20 2011-03-02 ヤマハ株式会社 波形発生回路
ES2401817T3 (es) * 2008-01-31 2013-04-24 Agency For Science, Technology And Research Procedimiento y dispositivo de distribución/truncado de la velocidad de transmisión de bits para codificación de audio escalable
JP5242478B2 (ja) * 2009-03-26 2013-07-24 株式会社エヌ・ティ・ティ・ドコモ 無線通信装置、無線通信システム及び無線通信方法
JP6284298B2 (ja) * 2012-11-30 2018-02-28 Kddi株式会社 音声合成装置、音声合成方法および音声合成プログラム
CN103795377B (zh) * 2014-01-14 2016-08-24 深圳市文鼎创数据科技有限公司 正弦波生成方法和装置以及电子设备
CN105974997B (zh) * 2016-04-29 2018-08-14 洛阳理工学院 一种正弦波信号的数字实现方法
CN108010537A (zh) * 2017-12-30 2018-05-08 长沙迪普美医疗科技有限公司 一种声音报警处理方法、装置、系统
EP4155767A1 (en) * 2021-09-24 2023-03-29 SuperSonic Imagine Method and system for processing a set of signals received by a transducer element

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5850802A (ja) 1981-09-19 1983-03-25 Anritsu Corp 多周波デイジタル正弦波発生装置
AT399236B (de) * 1986-02-03 1995-04-25 Siemens Ag Oesterreich Digitaler sinusgenerator
JPH01278104A (ja) 1988-04-30 1989-11-08 Oki Electric Ind Co Ltd ディジタル発振器
JPH10171462A (ja) 1996-12-16 1998-06-26 Yamaha Corp 正弦波分析合成方法及び装置
JP2000347678A (ja) 1999-06-08 2000-12-15 Nec Corp トーン発生回路
US6591230B1 (en) 1999-11-12 2003-07-08 Texas Instruments Incorporated Coprocessor for synthesizing signals based upon quadratic polynomial sinusoids
JP2002111503A (ja) 2000-09-27 2002-04-12 Victor Co Of Japan Ltd 復号化装置

Also Published As

Publication number Publication date
US7405499B2 (en) 2008-07-29
EP1538740A1 (en) 2005-06-08
US20050251709A1 (en) 2005-11-10
WO2004006425A1 (ja) 2004-01-15
CN1666407A (zh) 2005-09-07
CN100477486C (zh) 2009-04-08
EP1538740A4 (en) 2006-04-12
JP2004048101A (ja) 2004-02-12

Similar Documents

Publication Publication Date Title
JP6849847B2 (ja) 信号生成システム及び信号生成方法
US7649135B2 (en) Sound synthesis
KR20190104457A (ko) 고조파 전위에 기초하여 개선된 서브밴드 블록
JP4019824B2 (ja) 波形生成装置及び方法並びに復号装置
KR20020031103A (ko) 압축 방법 및 장치, 압축 해제 방법 및 장치, 압축·압축해제 시스템, 기록 매체
EP1851752B1 (en) Sound synthesis
EP1840871B1 (en) Audio waveform processing device, method, and program
RU2008105555A (ru) Синтез аудиосигнала
US5185491A (en) Method for processing a waveform
JP2759163B2 (ja) ペダル効果付加装置を内蔵した電子楽器
JPS642960B2 (ja)
JPH0422275B2 (ja)
JP2011090189A (ja) 音響信号の符号化方法および装置
JP3419005B2 (ja) 楽音信号の合成に用いられる原楽音信号の発生装置
JP3112462B2 (ja) 音声符号化装置
JPH1020886A (ja) 波形データに存在する調和波形成分の検出方式
KR0181587B1 (ko) 엠펙-1 오디오 복호기의 합성 필터링 장치 및 방법
JPS60140299A (ja) 素片編集型音声分析装置
JPH0644199B2 (ja) 可変長フレ−ム音声分析合成方式
JPH0795010A (ja) ディジタルフィルタ
JP2000357969A (ja) オーディオ信号の符号化装置
JP2002073067A (ja) オーディオ信号の復号方法、及びオーディオ信号の復号装置
KR960015132A (ko) 디지탈 정현파 발생방법 및 그 회로
JPH0519767A (ja) 楽音合成装置
JPH02239296A (ja) 音声合成方法および装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040419

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070327

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070510

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: 20070904

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070917

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101005

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101005

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees