JP3379624B2 - 波形同期方法 - Google Patents

波形同期方法

Info

Publication number
JP3379624B2
JP3379624B2 JP03181497A JP3181497A JP3379624B2 JP 3379624 B2 JP3379624 B2 JP 3379624B2 JP 03181497 A JP03181497 A JP 03181497A JP 3181497 A JP3181497 A JP 3181497A JP 3379624 B2 JP3379624 B2 JP 3379624B2
Authority
JP
Japan
Prior art keywords
waveform
time
time series
phase
periodicity
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
JP03181497A
Other languages
English (en)
Other versions
JPH10229336A (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.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone 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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP03181497A priority Critical patent/JP3379624B2/ja
Publication of JPH10229336A publication Critical patent/JPH10229336A/ja
Application granted granted Critical
Publication of JP3379624B2 publication Critical patent/JP3379624B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】この発明は、音響信号処理、
通信、およびその他の分野における一次元信号処理にお
いて、処理対象の波形と、基準波形とが与えられたとき
に、処理対象の波形の位相を基準波形の位相に同期させ
る方法に関する。
【0002】
【従来の技術】従来、位相同期方法に関しては、位相ロ
ックループ、すなわちローパスフィルタと位相比較器と
による閉ループ回路を用いて対象波形の位相を制御する
方法が知られている。しかし、この方法は、正弦波およ
び方形波など単純な波形の位相を制御することを前提と
して、波形の基本周期の取得やFM検波等の目的で設計
されたものであって、例えばあるピアノと他のピアノと
の同一高さの音の波形のような複雑な成分をもつ波形に
対して局所的な伸縮を行うことで位相を制御することは
できないという欠点があった。
【0003】また、線形予測分析、自己相関分析など信
号処理の手法を用いて基準波形の基本周期を抽出し、こ
れに同期させて対象波形の位相を適応的に制御する方法
が知られている。しかしながら、これらの方法では、基
準波形が複数の基本周期からなる波形の和であった場合
に、適切な位相同期処理ができないという欠点があっ
た。
【0004】従って、上記の各方法は、基準波形が複数
の相異なる基本周期をもつ波形の和で構成されている場
合、すなわち例えばピアノの和音のような複数の周波数
成分からなる複雑な波形である場合に、その中の任意の
周波数成分の位相を基準として対象波形の位相を同期さ
せたい場合にあっては、十分な位相同期処理が期待し難
い。
【0005】
【発明が解決しようとする課題】この発明は、基準波形
が複数の相異なる基本周期をもつ波形の和で構成されて
いる場合、すなわち複数の周波数成分からなる場合で
も、その中の任意の周波数成分の位相を基準として対象
波形に対して局所的な伸縮を行うことで位相を同期させ
ることができ、公知の方法と比較して自由度の高い波形
位相同期方法を提供することを目的としている。
【0006】
【課題を解決するための手段】この発明によれば、基準
波形時系列および対象波形時系列に含まれるほぼ同一の
周波数成分の位相差の時系列を取得し、その各位相差を
時間方向の移動量に換算して時間移動量の時系列を得、
その時間移動量の時系列の各時間移動量ごとに、対象波
形の時系列の対応するものを遅らせ又は進ませる。
【0007】
【発明の実施の形態】次に、この発明の実施形態につい
て図面を用いて説明する。図1は、この発明による波形
位相同期方法を適用した波形位相同期装置の機能構成例
である。この波形位相同期装置10は、その基準入力端
子11からの基準波形と少なくとも1つの対象入力端子
12からの少なくとも一つの対象波形とを入力とし、対
象波形に対して位相同期処理を行った波形を出力端子1
3に出力する。
【0008】基準波形及び対象波形は例えば48kHz
とか96kHzとかでサンプリングされ、デジタル値の
時系列としてそれぞれ入力され、周期性取得手段14,
15でそれぞれ基準波形および対象波形に含まれている
周期性が周波数成分として抽出される。位相差取得手段
16でこれら抽出された基準波形の周期性すなわち周波
数成分の周波数と、対象波形の周期性すなわち周波数成
分の周波数とを比較して、同期させるべき周期性(周波
数)を決定乃至取得、その周波数を通過域とする狭帯域
通過フィルタなどの位相抽出手段によって、同期させる
べき周波数成分の位相差の時系列が抽出される。この位
相差時系列の抽出処理については後で詳細に説明する。
【0009】この抽出された角度を単位とする位相差の
時系列は時間移動量算出手段17で時間を単位とする時
間差、つまり時間移動量の時系列に換算される。この時
間差(移動量)の時系列を参照しながら、波形遅延手段
18で対象波形について時間シフトが行われる。負の遅
延量を考えれば、遅れだけではなく波形の進みも考慮さ
れる。
【0010】次に、上述した各手段14〜18における
各処理を具体的に説明する。周期性取得手段14,15
では、図2に示すように入力される波形、すなわち基準
波形と対象波形のそれぞれを読み込んで(ステップ30
1)、狭帯域通過フィルタバンクなどの方法で周波数解
析を行う(ステップ302)。例えば1オクターブ当た
り40個とか90個程度の狭帯域通過フィルタを用いて
何れのフィルタから出力が得られたかにより解析する。
あるいはいわゆるFFT(高速フーリエ変換)処理など
により周波数解析を行う。
【0011】次に、この周波数解析によって時間周波数
平面上でのパワー表現が得られるので、周波数方向での
パワーの極大点(ローカルピーク)を見出す(ステップ
303)。続いて、時間的に連続するローカルピークを
接続して、一続きのローカルピークとする(ステップ3
04)。ローカルピークを時間的に接続したものは、周
波数成分と呼ばれ、もとの波形に存在する色々な周期性
を表現したものである。この周波数成分を周期性情報と
して出力する(ステップ305)。
【0012】位相差取得手段16では、図3に示すよう
に、まず、基準波形と対象波形に存在する周期性の中
で、同期をとるものを選択して、その周期性に合わせて
狭帯域通過フィルタの中心周波数を設定する(ステップ
401)。同期をとる周期性すなわち周波数成分の組の
選択、つまり基準波形中の同期基準の周波数成分と対象
波形中の同期をとりたい周波数成分との選択は、外部か
らの選択によるか、あるいはステップ401において、
対象波形の基本周波数を選択したり(基準波形にその選
択した基本周波数に近い成分がある)、基準波形と対象
波形とで最も近接した周波数成分を選択したりするなど
任意の設計が可能である。何れの場合でも選択した基本
波形の周波数と対象波形の周波数とは共に前記狭帯域通
過フィルタを通過するようにする。このように、波形に
含まれる複数種類の周期性を任意に選択して同期をとる
ことができるのがこの発明の利点である。狭帯域通過フ
ィルタの中心周波数を設定すると、基準波形に狭帯域通
過フィルタを適用して出力波形を得て(ステップ40
2)、その出力波形の位相の時系列を求めて記憶する
(ステップ403)。狭帯域通過フィルタの出力は正弦
波に近いので、各時刻での位相の計算は容易である。
【0013】前述したように入力波形はデジタル値の時
系列であり、フィルタ処理もデジタル処理、つまり演算
により行われ、フィルタ出力としての正弦波もデジタル
値の時系列として得られ、その時系列におけるデジタル
値の符号(サインビット)が正から負、又は負から正に
それぞれ変化した隣接処理(演算)時刻kとk+1の間
にその正弦波のゼロクロス点があり、時刻k、k+1の
各サンプル値で、時刻k、k+1間のゼロクロス点を、
例えば線形補間して、正確なゼロクロス点の時刻を求め
る。この二つの隣接したゼロクロス点の時刻間の差がそ
の正弦波の半周期、つまりπラジアンであるから、この
周期と、ゼロクロス点時刻とから正弦波のデジタル値時
系列の各デジタル値における位相値(角度)を求めて位
相の時系列を得る。
【0014】続いて対象波形にもステップ402で用い
たものと同じ狭帯域通過フィルタを適用して出力波形を
得て(ステップ404)、その出力波形の位相の時系列
を求めて記憶する(ステップ405)。次に、ステップ
403とステップ405とでそれぞれ記憶した位相の時
系列の対応する時刻のものを減算することによって、位
相差の時系列を得て出力する(ステップ406)。
【0015】時間移動量算出手段17では、図4Aに示
すように位相差取得手段16から出力される位相差の時
系列と、適用された狭帯域通過フィルタの中心周波数を
読み込み(ステップ501)、時間差の時系列に換算す
る(ステップ502)。この換算は、次式(1)によっ
て行う。 δt(k)=(1/(2πf))・δp(k) (1) ただし、kは時刻、δt(k)は時間差時系列の時刻k
での時間差、fはステップ501で読み込まれた狭帯域
通過フィルタの中心周波数、δp(k)はステップ50
1で読み込まれた位相差時系列の時刻kでの位相差であ
る。式(1)によって計算された時間差の時系列は、逐
次出力する(ステップ503)。
【0016】波形遅延手段18では図4Bに示すよう
に、まず対象波形を読み込み(ステップ601)、次に
時間移動量算出手段17から出力された時間差の時系列
を読み込む(ステップ602)。つぎに、対象波形の各
時点において、ステップ602で得られた対応時点の時
間差だけ時間をずらす(遅らせるまたは進ませる)処理
を行う(ステップ603)。ステップ603の出力が、
基本波形のステップ401で設定した中心周波数成分に
位相同期した対象波形である。
【0017】図5にこの発明を適用した波形例を示す。
基準波形Aの基本波成分に、対象波形Bを位相同期させ
ると波形Cとなる。この発明では各時刻ごとに位相差に
応じて対象波形のサンプルを移動させ、対象波形が伸縮
され、基準波形Aのゆらぎに対象波形Bのゆらぎが一致
したものとなる。この例では基準波Aの第1、第2、第
8、第9山又は対象波形Bの第1、第2、第6、第7山
とほぼ同一時刻に生じ、対象波形Bの第1、第2、第
6、第7山はそれぞれほぼそのまま出力されるが、対象
波形Bの第3山は基準波の第3、第4と対応して位相が
進まされたものと遅らされたものの2回現われる。
【0018】基準波形、対象波形の周波数成分が予めわ
かっている場合は、図1中の周期性取得手段14,15
を省略できる。
【0019】
【発明の効果】以上、説明したように、この発明によれ
ば、基準波形が複数の相異なる基本周期をもつ波形の和
で構成されている場合、すなわち複数の周波数成分から
なる場合に、その中の任意の周波数成分の位相を基準と
して対象波形の位相を同期させることができるという利
点がある。
【図面の簡単な説明】
【図1】この発明方法を適用した波形位相同期装置の機
能構成例を示すブロック図。
【図2】周期性取得手段14,15の処理手順を示す流
れ図。
【図3】位相差取得手段16の処理手順を示す流れ図。
【図4】Aは、時間移動量算出手段17の処理手順を示
す流れ図、Bは波形遅延手段18の処理手順を示す流れ
図。
【図5】基準波形と対象波形と、この発明により同期化
した波形との例を示す図。

Claims (2)

    (57)【特許請求の範囲】
  1. 【請求項1】 基準波形中の同期基準の周波数成分と、
    対象波形中の同期をとりたい周波数成分との位相差の時
    系列を取得する位相差取得過程と、 上記取得した位相差時系列の各位相差を極性を含む時間
    差に変換した時間移動量の時系列を算出する時間移動量
    算出過程と、 上記対象波形に対し、上記時間移動量時系列の各時間移
    動量だけ時間的に移動させる波形遅延過程とを有し、 上記位相差取得過程は、 上記同期基準の周波数成分の周波数を中心周波数とする
    狭帯域通過フィルタに基準波形の時系列を通す過程と、
    そのフィルタを通過した基準波形時系列から、その各時
    刻におけるその位相値を計算して位相値時系列を得る過
    程と、 上記狭帯域通過フィルタに上記対象波形の時系列を通す
    過程と、そのフィルタを通過した対象波形時系列から、
    その各時刻における位相値を計算して位相値時系列を得
    る過程と、 上記両位相値時系列の対応する時刻における位相値の差
    をとって上記位相差時系列を得る過程とよりなることを
    特徴とする波形同期方法。
  2. 【請求項2】 上記基準波形から周期性を取得する基準
    周期性取得過程と、上記対象波形から周期性を取得する
    対象周期性取得過程と、上記取得した基準波形の周期性
    と、上記取得した対象波形の周期性とから上記同期基準
    の周波数成分と上記同期をとりたい周波数成分を決定す
    る過程とを有することを特徴とする請求項1記載の波形
    同期方法。
JP03181497A 1997-02-17 1997-02-17 波形同期方法 Expired - Fee Related JP3379624B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP03181497A JP3379624B2 (ja) 1997-02-17 1997-02-17 波形同期方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP03181497A JP3379624B2 (ja) 1997-02-17 1997-02-17 波形同期方法

Publications (2)

Publication Number Publication Date
JPH10229336A JPH10229336A (ja) 1998-08-25
JP3379624B2 true JP3379624B2 (ja) 2003-02-24

Family

ID=12341570

Family Applications (1)

Application Number Title Priority Date Filing Date
JP03181497A Expired - Fee Related JP3379624B2 (ja) 1997-02-17 1997-02-17 波形同期方法

Country Status (1)

Country Link
JP (1) JP3379624B2 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MXPA03010751A (es) * 2001-05-25 2005-03-07 Dolby Lab Licensing Corp Segmentacion de senales de audio en eventos auditivos.
WO2006046632A1 (ja) 2004-10-27 2006-05-04 Nec Corporation ディジタル無線受信装置
CN106130543A (zh) * 2016-06-21 2016-11-16 中国电子科技集团公司第五十八研究所 一种高分辨率时钟相移架构与算法的实现方法

Also Published As

Publication number Publication date
JPH10229336A (ja) 1998-08-25

Similar Documents

Publication Publication Date Title
Dolson The phase vocoder: A tutorial
JP4940588B2 (ja) ビート抽出装置および方法、音楽同期画像表示装置および方法、テンポ値検出装置および方法、リズムトラッキング装置および方法、音楽同期表示装置および方法
US6657567B2 (en) Compressing method and device, decompression method and device, compression/decompression system, and recorded medium
EP2401740B1 (en) Apparatus and method for determining a plurality of local center of gravity frequencies of a spectrum of an audio signal
EP2881944B1 (en) Audio signal processing apparatus
US4937868A (en) Speech analysis-synthesis system using sinusoidal waves
WO2002086866A1 (fr) Procede et dispositif de compression, procede et dispositif de decompression, systeme de compression/decompression, procede de detection de crete, programme et support d'enregistrement
JP2006293230A (ja) 音響信号処理装置、音響信号処理プログラム及び音響信号処理方法
JP3379624B2 (ja) 波形同期方法
JPH05281996A (ja) ピッチ抽出装置
EP1288912A1 (en) Speech recognition method and device, speech synthesis method and device, recording medium
CA1319994C (en) Speech analysis method
JP2001136073A (ja) 圧縮方法及び装置、圧縮伸長システム、記録媒体
JP2000180484A (ja) 高調波測定装置
Suzuki et al. Instantaneous frequencies of signals obtained by the analytic signal method
JP2007225476A (ja) 波形分析装置
JP3019603B2 (ja) 音声の基本周波数の抽出装置
JPH07218565A (ja) ジッタ周波数成分の検出方法
JP3438376B2 (ja) 周期信号処理方法および装置
JP2009186762A (ja) 拍タイミング情報生成装置およびプログラム
JP3223564B2 (ja) ピッチ抽出方法
JPS5858625B2 (ja) 信号処理装置
JP2543271B2 (ja) 時間軸補正装置
JPS6091227A (ja) 音響解析装置の合成装置
Goh Noise robust AM_FM demodulation using least-squares truncated power series approximation

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20071213

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20081213

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20091213

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20101213

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20101213

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20111213

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20111213

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20121213

Year of fee payment: 10

LAPS Cancellation because of no payment of annual fees