JPS59632A - 信号分析装置 - Google Patents

信号分析装置

Info

Publication number
JPS59632A
JPS59632A JP9942383A JP9942383A JPS59632A JP S59632 A JPS59632 A JP S59632A JP 9942383 A JP9942383 A JP 9942383A JP 9942383 A JP9942383 A JP 9942383A JP S59632 A JPS59632 A JP S59632A
Authority
JP
Japan
Prior art keywords
autocorrelation function
pitch
linear prediction
signal
resonance frequency
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
JP9942383A
Other languages
English (en)
Other versions
JPH0139056B2 (ja
Inventor
Hisashi Nishiyama
久司 西山
Kageyoshi Katakura
景義 片倉
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.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP9942383A priority Critical patent/JPS59632A/ja
Publication of JPS59632A publication Critical patent/JPS59632A/ja
Publication of JPH0139056B2 publication Critical patent/JPH0139056B2/ja
Granted legal-status Critical Current

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔発明の利用分野〕 本発明は・、音響等の連続信号から帯域幅、共振周波数
等を抽出して出力する信号分析装置に関し。
−1:、とじて船舶の航走音などの1機械音から対象に
特有の帯域幅、共振周波数、インパルス応答等を特徴と
して捕えることを目的とする信号分析装置に関するもの
である。
〔発明の背景〕
この種の信号分析装置として、線形予測係数を用いるも
のが提案(特開昭51−112377参照)されている
が、その場合9時間差零付近の自己相関関数から線形予
測係数を求めていたので、雑音除去が不充分で、抽出精
度が悪かった。
〔発明の目的〕
本発明は9時間差零付近の自己相関関数を用し・る代り
に2時間差が信号ピークの1ピツチだけ離れた時点を中
心とするある領域内の自己相関関数を用いて線形予測係
数を求めることで、雑音除去を可能にし、抽出精度の向
上化を図ろうとするも〔発明の概要〕、 本発明の特徴は、上記目的を達成するために。
入力信号をサンプリングして得られ名デジタル信号から
一定時間間隔ごとに自己相関関数を求める自己相関計と
、この自己相関関数をある所定数だけ累加する累加装置
と、この累加自己相関関数から入力信号のピークのピッ
チ概略値を抽出するピッチ抽出装置と9時間差が時間零
時点から上記ピッチ概略値だけ離れた点を中心とする前
後の、ある一定時間幅領域内の自己相関関数を補間で求
める補間装置と、補間された自己相関関数から時間差が
正確に1ピツチ離れた点の自己相関関数を再生する再生
装置と、この再生自己相関関数からAIC(情報量規準
)による最適次数の線形予測係数を求める線形予測係数
抽出装置と、この線形予測係数を用いてスペクトル」ユ
で意味のある帯域幅共振周波数あるいはインパルス応答
を求めてこれをスペクトル表示装置、インパルス応答聴
音器にそれぞれ出力する回路手段とを備えた構成とする
にある。
〔発明の実施例〕
まず9本発明の原理について説明する。
十、述のように2本発明においては9時間差が1ピツチ
だけ離れた時点付近の自己相関関数から線形予測係数を
求める。そのために、シャノンのサンプリング定理によ
り、ピッチ非同期の自己相関関数から、補間により、ピ
ッチに同期した自己相関関数を再生する。即ち、ピッチ
非同期の相関関数V5が1時間中+b−lzでの信号波
形をfl。
1=o、 ]、 2.・・・Nとして ■、−Σfifi−Ls+  8=0.1 、2.・・
・−二〇 で与えられるとき、再生自己相関関数R1はn= o 
+ 1+ ・・・+ N+ (Nll+ 1) 1とな
る。ただし、Δは信号のサンプリング時間間・・1隔〔
第1図fbl参照L ”+は1ピッチ当りのデータ数で
あり、予め求めた概略ピッチをP+ [第1図tal参
照〕とすれば N1=Pi/1 の関係があり、また再生自己相関関数におけるサンプリ
ング時間T(ピッチ分解能とも呼ばれる)は T−Δ/(No+1) となる。
このとき、再生自己相関関数列R0のピークの存在する
時刻をnmaXとすれば、推定ピッチP1 はP(””
 2 N1”Δ+(n、、、、、x−1)Tであり、補
間再生された自己相関関数V5′はVs””R((nm
ax 1)T+sΔ)+  S=0+ 1+”’+ P
で与えられる。ただし、 R(nT)=Rnである。
従って、正規方程式を解いて線形予測係数を求める際に
は、第1図talに示すように1ピッチ付近の雑音に乱
されない、信号だけの自己相関関数V:を用いることに
より、雑音の影響のない線形予測係数か得られる3゜ 予測次数の決定に際しては、AIC(情報量規準)を用
いる。即ち AIC=Ntogeσ、+2(p+1)で定義され、p
は予測次数、Nは全データ数、σ2は残差自乗平均であ
る。各次数について正規方程式を解き、線形予測係数 (A1)       ・ p=1 (A+、A2)     1  p=2p (A+ l A21・・・+ Ap) +  p””p
を求め、そのときのAICの値が最小となる予測次数を
最適パラメータ数とする。
最適パラメータ数をp。とすれば、96次の線形予測係
数(A1.A2.・・・+ Apg )を係数とする代
数方程式を解くことにより、98個の(共振周波数fi
 +帯域幅B、)の組(f+、 BIL (f2. B
4C”””l (rpo。
B、。)が得られる。このとき、対象となる物理系の極
の数が既知であれば、帯域幅の小さい方から。
その次数分だけの共振周波数、帯域幅を残し、残りは捨
て去る。あるいは、帯域幅の極端に大きいものは、雑音
や計算誤差を近似している次数であるとみなせるので、
取り込む必要はない。
なお、第1図falは時間差零付近に雑音の相関が加わ
って本来の自己相関関数波形が損なわれるが。
1ピツチ付近では雑音の影響のない、信号だけの自己相
関関数が現われることを示す図である。図中の“′Δ′
′印はピッチ非同期のサンプリング点を示し、“0”′
印はそれから補間再生されたピッチ同期のサンプリング
点を示している。第1図tblはサンプリング時間間隔
ΔにNl(個の補間点をとることを示す図である。
以下9本発明の一実施例を、第2図に示すブロック構成
図に従って説明する。入力された連続信号をA/D変換
器1でデジタル符号列に変換し。
自己相関計2で、制御装置12の制御により、一定・時
間間隔ごとに自己相関関数を求める。この自己相関関数
は、累加装置3で累加される。累加装置3において累加
される自己相関関数の個数は制御装置12により制御さ
れる。累加された自己相関関数はピッチ抽出装置4に送
られ、ここで、ピッチ1 。    3 概略値が抽出され、〜ヒツチからlピッチまでの時間領
域1りにおける自己相関関数が、補間装置5により補間
で求められる。このときの補間の個数は制御装置12で
制御される。次の自己相関関数再生装置6において、」
ユ記の補間された自己相関関数から、1ピツチロのピー
クを正確に求め、それを推定ピッチとし、この推定ピッ
チを中心にして。
自己相関関数を元のサンプリング時間でサンプリングす
ることにより、ピッチに同期した自己相関関数が再生さ
れる。この再生自己相関関数をAIOによる最適次数の
線形予測係数抽出装置7に加え。
得られた最適次数の線形予測係数から帯域幅、共振周波
数を求めるが、帯域幅、共振周波数選択装置8において
、スペクトル上で意味のある帯域幅。
共振周波数をより分ける。選択された帯域幅、共・振周
波数はスペクトル波形としてスペクトル装置9で表示さ
れる。また、最適次数の線形予測係数はインパルス応答
装置10にも印加され、ここでインパルス応答が求めら
れ、これをインパルス応答聴音器11で聴くことも可能
である。これらの全体の動作は制御装置12により制御
される。制御装置12は、操作パネル13を介して、使
用者からの指示により、動作モードを変更する。
〔発明の効果〕
以下9本発明の効果を、シミュレーション実験結果によ
り説明する。従来の信号分析装置においては、SN比が
悪いとき、帯域幅が真値よりかなり広くなる傾向がある
。即ち、単一共振モデルの信号として5 exp (−
t ) sin 2πtの減衰正弦波をピッチ周期5秒
で繰り返した第3図波形を、サンプリノブ時間0.2秒
でサンプリングし、それに雑音として正規乱数を加法的
に加えることで第4図波形が得られる。この場合、ピッ
チに同期したサンプリング時間なので、補間数を零とし
た。そのとき。
でSN比を定義すれば、SN比= −2,75dBの減
衰正弦波信号である。この第4図の波形信号を、従来の
信号分析装置を用いて、予測次数は単共振であるので2
とし、自己相関関数の累加数を500回。
積分区間は75データ数として2分析した。その結果、
帯域幅は1.29 Hz 、共振周波数は1.14Hz
と求められた。この結果は、自己相関関数の累加数及び
積分区間をさらに増大させても、はとんど変化しなかっ
た。
」1記の減衰正弦波から得られる帯域幅、共振周波数の
真値は、それぞれ0.318 Hz 、 1.00 H
zである。
従って、従来方式での分析結果は、帯域幅では真値の約
4倍となる過大推定であり、共振周波数で10%以」二
の誤差となる。これに対し1本考案の方式によるときは
、帯域幅及び共振周波数がそれぞれ0.310 Hz 
、 1.01 Hzと求めらし、抽出精度カ著シく改善
されることが確認された。
なお、第2図実施例回路において自己相関関数を求める
積分時間が雑音に比較して長いときは。
累加装置3の設置は省略しても差支えない。またインパ
ルス応答聴音器11とスペクトル表示装置9とは、必要
に応じてその一方を省略することも。
あるいはその機能の一部のみを使用することも可能であ
る。
【図面の簡単な説明】
第1図falは自己相関関数と時間差との関係を説明す
る図、第1図fblはサンプリング時間間隔Δ内・にN
14個の補間点をとることを示す図、第2図は本発明の
一実施例のブロック構成図、第3図は効果の比較実験用
に想定した物理系から得られる信号波形図、第4図は第
3図の信号に雑音を加法的に加えて得られる信号波形図
である。 符号の説明 1・・・A/D変換器    2・・・自己相関計3・
・・累加装置     4・・・ピッチ抽出装置5・・
・補間装置 6・・・自己相関関数再生装置 7・・・線形予測係数抽出装置 8・・・帯域幅、共振周波数選択装置 9−・・スペクトル表示装置 10・・・インパルス応答装置 11・・・インパルス応答聴音器 12・・・制御装置     13・・・操作、(ネル
代理人弁理士 中村純之助 11図 (引

Claims (1)

  1. 【特許請求の範囲】 1 音響等の連続信号の信号分析装置において人力信号
    から自己相関関数に対応した量を計算する計算手段と、
    該自己相関関数に対応した量から入力信号のピークのピ
    ッチ概略値を算出するピッチ算出手段と2時間差が時間
    零時点から上記ピッチ概略値だけ離れた概略時点を基準
    とする所定時間領域内の自己相関関数を補間演算により
    求める補間演算手段と、補間演算により求められた自己
    相関関数にもとづき上記時間零時点から1ピツチだけ離
    れた時点の自己相関関数を再生する再生手段と、再生さ
    れた自己相関関数から線形予測係数を算出する線形予測
    係数算出手段とを備えたことを特徴とする信号分析装置
    。 2、 上記自己相関関数に対応した量はフレーム周期ご
    とに計算された自己相関関数を累加した量であることを
    特徴とする特許請求の範囲第1項の信号分析装置。 3、上記所定領域は上記概略時点を中心とする前後を含
    む所定時間幅の領域であることを特徴とする特許請求の
    範囲第1項または第2項の信号分析装置。 4、上記所定領域は」ユ記概略時点の前または後におけ
    る所定時間幅の領域であることを特徴とする特許請求の
    範囲第1項または第2項の信号分析装置。 5、上記線形予測係数算出手段は再生された自己相関関
    数から情報量規準により最適次数の線形予測係数を求め
    ることを特徴とする特許請求の範囲第1項乃至第4項い
    ずれか1項の信号分析装置。 6、上記線形予測係数算出手段は線形予測係数から求め
    られた帯域幅と共振周波数とスペクトル包絡またはイン
    パルス応答とを表示するスペクトル表示手段と、インパ
    ルス応答聴音手段とを含むことを特徴とする特許請求の
    範囲第1項乃至第5項いずれか1項の信号分析装置。
JP9942383A 1983-06-06 1983-06-06 信号分析装置 Granted JPS59632A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP9942383A JPS59632A (ja) 1983-06-06 1983-06-06 信号分析装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP9942383A JPS59632A (ja) 1983-06-06 1983-06-06 信号分析装置

Publications (2)

Publication Number Publication Date
JPS59632A true JPS59632A (ja) 1984-01-05
JPH0139056B2 JPH0139056B2 (ja) 1989-08-17

Family

ID=14247050

Family Applications (1)

Application Number Title Priority Date Filing Date
JP9942383A Granted JPS59632A (ja) 1983-06-06 1983-06-06 信号分析装置

Country Status (1)

Country Link
JP (1) JPS59632A (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006259789A (ja) * 2005-03-15 2006-09-28 National Institute Of Advanced Industrial & Technology 周期判定装置、周期判定方法および周期判定プログラム
JP2014081352A (ja) * 2012-09-27 2014-05-08 Daihen Corp 周波数解析装置、当該周波数解析装置を用いた信号処理装置、および、当該信号処理装置を用いた高周波測定装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006259789A (ja) * 2005-03-15 2006-09-28 National Institute Of Advanced Industrial & Technology 周期判定装置、周期判定方法および周期判定プログラム
JP2014081352A (ja) * 2012-09-27 2014-05-08 Daihen Corp 周波数解析装置、当該周波数解析装置を用いた信号処理装置、および、当該信号処理装置を用いた高周波測定装置

Also Published As

Publication number Publication date
JPH0139056B2 (ja) 1989-08-17

Similar Documents

Publication Publication Date Title
US8233636B2 (en) Method, apparatus, and computer program for suppressing noise
US20090047003A1 (en) Playback apparatus and method
KR970001166B1 (ko) 언어 처리 방법 및 장치
US5576978A (en) High resolution frequency analyzer and vector spectrum analyzer
JP2004505304A (ja) デジタルオーディオ信号の連続可変時間スケール変更
JP2001298368A (ja) 圧縮方法及び装置、伸長方法及び装置、圧縮伸長システム、記録媒体
JP3402748B2 (ja) 音声信号のピッチ周期抽出装置
EP1289152A1 (en) Compression method and apparatus, expansion method and apparatus, compression and expansion system, recorded medium, program
JP2625622B2 (ja) 信号プロセッサ及びアナログ信号を変換するための方法
JPS59632A (ja) 信号分析装置
JPS62229200A (ja) ピツチ検出器
JP5035815B2 (ja) 周波数測定装置
JP3139803B2 (ja) インパルス応答測定装置
JPS642960B2 (ja)
JP3919359B2 (ja) 楽音信号のアタック位置検出装置
JPH0229239B2 (ja)
JPH0758734A (ja) ディジタルデータ検出器
CN114387989B (zh) 语音信号处理方法、装置、系统及存储介质
Katsumata et al. An audio sampling frequency conversion using digital signal processors
JPH10503299A (ja) 改良されたピッチ検出を備えた音声符号化用伝送システム及び方法
KR100745133B1 (ko) 계산량을 줄일 수 있는 샘플링 레이트 컨버젼 방법 및샘플링 레이트 컨버젼 시스템
JPS5876891A (ja) 音声ピツチ抽出方法
KR0147185B1 (ko) 영상가요재생시스템의 점수산정제어장치
CN114465628A (zh) 一种精确亚音频ctcss频率搜索系统及方法
JPH07234672A (ja) ピッチ検出装置