(57)【特許請求の範囲】
【請求項１】
輻射ノイズ変調に基づくプロペラ翼数特徴の抽出方法であって、
（１）船舶の輻射ノイズ信号を採集し、フーリエ変換により変調スペクトル図を取得す
るステップと、
（２）変調スペクトル図にローカルピークを探索し、ローカルピーク位置における共振
周波数を取得するステップと、
（３）変調スペクトル図における１番目ローカルピークの共振周波数を確定し、軸周波
数を確定するとともに、ローカルピークの数に基づいて共振周波数の数Ｎを確定するステ
ップと、
（４）軸周波数と翼周波数間の逓倍関係を確定し、他の共振周波数を確定するステップ
と、
（５）各共振周波数位置の平均スペクトルコヒーレント値を確定し、ナイーブベイズの
推断法によりプロペラ翼数を取得するとともに、最後に翼周波数を確定するステップと、
を含み、
ステップ（２）では、ローカルピークを探索するときは、隣り合う両ローカルピークの
共振周波数の差が下記関係を満たし、
0006836041.tif 000083
ただし、
0006836041.tif 000084
はｎ番目ローカルピーク位置の循環周波数の値であり、
ステップ（３）では、軸周波数を確定するときは、１番目ローカルピーク位置の共振周
波数を軸周波数とし、１番目ローカルピーク位置の共振周波数値が０．９Ｈｚより小さけ
れば、それを除去して、２番目共振周波数を軸周波数とする、
ことを特徴とする輻射ノイズ変調に基づくプロペラ翼数特徴の抽出方法。
【請求項２】
ステップ（４）では、他の共振周波数を確定する公式は、
0006836041.tif 000085
であり、ただし、
0006836041.tif 000086
はｎ番目共振周波数、
0006836041.tif 000087
；
0006836041.tif 000088
は軸周波数である、ことを特徴とする請求項１に記載の輻射ノイズ変調に基づくプロペラ
翼数特徴の抽出方法。
【請求項３】
各共振周波数位置の平均スペクトルコヒーレント値
0006836041.tif 000089
は
0006836041.tif 000090
区間内の平均コヒーレント値から平均を求めて得られ、ただし、
0006836041.tif 000091
はサンプリング点
0006836041.tif 000092
前の5番目サンプリング点から
0006836041.tif 000093
後の5番目サンプリング点までの区間を示す、ことを特徴とする請求項１に記載の輻射ノ
イズ変調に基づくプロペラ翼数特徴の抽出方法。
【請求項４】
ステップ（５）では、ナイーブベイズの推断法に基づいて、取得の共振周波数値と対応の
平均コヒーレント値との間の関係により、該変調スペクトル図はどの類型のプロペラ翼数
の集合に従属するかを確定し、最後にはプロペラ翼数を確定し、
ナイーブベイズ公式は、
0006836041.tif 000094
であり、ただし、Ｙ＝｛
0006836041.tif 000095
｝（Ｙはすべての可能のプロペラ翼数の集合、
0006836041.tif 000096
はその一つのプロペラ翼数、
0006836041.tif 000097
はプロペラ翼数が３、
0006836041.tif 000098
はプロペラ翼数が４、
0006836041.tif 000099
はプロペラ翼数が５、
0006836041.tif 000100
はプロペラ翼数が６、
0006836041.tif 000101
はプロペラ翼数が７とする）、Ｘ＝｛
0006836041.tif 000102
｝（Ｘは分類目標の変調スペクトル、
0006836041.tif 000103
は該変調スペクトルにおける各共振周波数位置の振幅大小関係とする）、Ｐ（Ｘ）は分類
自身の確率（定常数）、
0006836041.tif 000104
は各プロペラ翼数の類型
0006836041.tif 000105
の先験的確率、
0006836041.tif 000106
は所定Ｘの
0006836041.tif 000107
プロペラ翼数類型への従属確率、
0006836041.tif 000108
はプロペラ翼数類型
0006836041.tif 000109
のＸ発生確率、
0006836041.tif 000110
はプロペラ翼数類型
0006836041.tif 000111
の変調スペクトルにおける特徴
0006836041.tif 000112
発生の確率であり、各
0006836041.tif 000113
を算出すると、
0006836041.tif 000114
が最大値であれば、Ｘは類型
0006836041.tif 000115
に従属する、ことを特徴とする請求項１に記載の輻射ノイズ変調に基づくプロペラ翼数特
徴の抽出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、信号の周波数成分の抽出の技術領域に関し、特に輻射ノイズ変調に基づくプロ
ペラ翼数特徴の抽出方法に関する。
【背景技術】
【０００２】
船舶の輻射ノイズは、回転翼（プロペラ）の回転際におけるキャビテーション・ノイズに
より引き起こされる。船舶の輻射ノイズの変調スペクトルは船舶のプロペラ翼数の特徴情
報を含む。各共振周波数と振幅の関係を分析することにより、船舶プロペラの軸周波数、
翼周波数およびプロペラ翼数の特徴量を抽出することができる。
【０００３】
民船のプロペラの軸周波数、k翼周波数およびプロペラ翼数の特徴量情報が獲得可能であ
れば、船舶型番の判断に寄与する。海に沿っては密輸・密航などがある。音声センサーに
より船舶のノイズ情報を捕捉し、フーリエ変換により変調スペクトルを取得し、プロペラ
の軸周波数、翼周波数およびプロペラ翼数の特徴量情報を抽出できれば、監督者が船舶の
種類を判断できる。
【０００４】
従来、世界では以下のような識別方法が使用されている。変調スペクトルを取得し、変調
スペクトルから幅値のピークおよびピーク位置の共振周波数を読み込む。プロペラ翼数特
徴識別規則（表１）に基づいて、変調スペクトルにおける共振周波数および振幅からプロ
ペラ翼数特徴値を抽出する。Ｐ（ｎ）は軸周波数のｎ番目高調波線スペクトルの幅値を示
す。
0006836041.tif 000002
【０００５】
ところが、船舶の構造、動作状況や環境などによっては、変調スペクトルの構成が複雑で
あり、識別規則に基づいては典型的な場合しか識別できなく、すべての変調スペクトル構
成は適用不可能である。例えば
0006836041.tif 000003
になった場合、表１のすべてのプロペラ翼数識別規則の条件は満たされるため、プロペラ
翼数を識別できない。
【０００６】
中国特許出願２０１９１０７９０２１７．Ｘにおいて、回転翼特徴の抽出方法に関して、
回転翼キャビテーション後流微細特徴の多次元統計モデリング方法が開示されたが、抽出
特徴は回転翼の幾何学パラメータと動作状況の特徴に主に限り、軸周波数、翼周波数やプ
ロペラ翼数の特徴に及ばない。
【０００７】
戴衛国、邱家興様等は２０１５年に「ベクトル計算機に多クラス分類適用の船舶プロペラ
翼数識別への研究」が公表されている。目標船舶の輻射ノイズの包絡信号識別スペクトル
に基づいて船舶のプロペラ翼数分類を実行する実験に適応する、誤り訂正コード組合出力
のベクトル計算機に多クラス分類適用の改良アルゴリズムが提案されたが、該方法はベク
トル計算機適用の推断法を利用し、大量の既知サンプル、および変調スペクトルにおける
３３次元まで多くの特徴抽出が必要であり、操作が複雑である。また、本文では軸周波数
や翼周波数の具体的識別法は提出されていなった。
【発明の概要】
【０００８】
従来の技術問題に対して、本発明は、輻射ノイズ変調に基づくプロペラ翼数特徴の抽出方
法を提出する。各種類民船の輻射ノイズ変調スペクトルを分析することができ、軸周波数
、翼周波数やプロペラ翼数の特徴量の抽出に便利である。
【０００９】
この輻射ノイズ変調に基づくプロペラ翼数特徴の抽出方法は、下記ステップを含む。
（１）船舶の輻射ノイズ信号を採集し、フーリエ変換により変調スペクトル図を取得する
。
（２）変調スペクトル図にローカルピークを探索し、ローカルピーク位置における共振周
波数を取得する。
（３）変調スペクトル図における１番目ローカルピークの共振周波数を確定し、軸周波数
を確定するとともに、ローカルピークの数に基づいて共振周波数の数Ｎを確定する。
（４）軸周波数と翼周波数間の逓倍関係を確定し、他の共振周波数を確定する。
（５）各共振周波数位置の平均スペクトルコヒーレント値を確定し、ナイーブベイズの推
断法によりプロペラ翼数を取得するとともに、最後に翼周波数を確定する。
ステップ１では、取得する変調スペクトル図は循環周波数を横軸とし、平均スペクトルコ
ヒーレント値を縦軸とする。
【００１０】
ＭＡＴＬＡＢ（登録商標）ソフトウェアのｆｉｎｄｐｅａｋｓ関数により平均コヒーレン
ト値のローカルピークおよび対応の循環周波数を探索する。ＭＡＴＬＡＢ（登録商標）中
のｆｉｎｄｐｅａｋｓツールボックス関数により元の波形における波ピーク位置を探出す
る。まず、１番目ローカルピークおよびその共振周波数
0006836041.tif 000004
を探出する。
【００１１】
ステップ（2）では、ローカルピークを探索するときは、隣り合う両ローカルピークの共
振周波数の差が下記関係を満たす。
0006836041.tif 000005
ただし、
0006836041.tif 000006
はｎ番目ローカルピーク位置の循環周波数の値である。
【００１２】
ステップ（３）では、軸周波数を確定するときは、１番目ローカルピーク位置の共振周波
数を軸周波数とする。もし１番目ローカルピーク位置の共振周波数値が０．９Ｈｚより小
さければ、それを除去して、２番目共振周波数を軸周波数とする。
【００１３】
ステップ（４）では、他の共振周波数を確定する公式：
0006836041.tif 000007
ただし、
0006836041.tif 000008
はｎ番目共振周波数、
0006836041.tif 000009
；
0006836041.tif 000010
は軸周波数である。
ステップ（５）では、共振周波数によって対応の平均コヒーレント値を確定する。共振周
波数
0006836041.tif 000011
位置の平均コヒーレント値は
0006836041.tif 000012
とする。
0006836041.tif 000013
はｎ番目共振周波数位置の幅値を示す。
各共振周波数位置の平均スペクトルコヒーレント値
0006836041.tif 000014
は
0006836041.tif 000015
区間内の平均コヒーレント値から平均を求めて得られる。ただし、
0006836041.tif 000016
はサンプリング点
0006836041.tif 000017
前の5番目サンプリング点から
0006836041.tif 000018
後の5番目サンプリング点までの区間を示す。
【００１４】
ステップ（5）では、ナイーブベイズの推断法に基づいて、取得の共振周波数値と対応の
平均コヒーレント値との間の関係により、該変調スペクトル図はどの類型のプロペラ翼数
の集合に従属するかを確定する。最後には、プロペラ翼数を確定する。
サンプルの数が十分である場合、直接にナイーブベイズの推断法を使用してもよく、その
結果も非常に確実的である。サンプルの数量が少ない場合、表1のプロペラ翼数識別規則
における各プロペラ翼数時の幅値間の関係に基づいてアナログサンプルを生成して、判断
目標の変調スペクトル図を分類してもよい。その分類結果はプロペラ翼数となる。ナイー
ブベイズ公式：
0006836041.tif 000019
ただし、Ｙ＝｛
0006836041.tif 000020
｝（Ｙはすべての可能のプロペラ翼数の集合、
0006836041.tif 000021
はその一つのプロペラ翼数、
0006836041.tif 000022
はプロペラ翼数が３、
0006836041.tif 000023
はプロペラ翼数が４、
0006836041.tif 000024
はプロペラ翼数が５、
0006836041.tif 000025
はプロペラ翼数が６、
0006836041.tif 000026
はプロペラ翼数が７とする）、Ｘ＝｛
0006836041.tif 000027
｝（Ｘは分類目標の変調スペクトル、
0006836041.tif 000028
は該変調スペクトルにおける各共振周波数位置の振幅大小関係とする）、Ｐ（Ｘ）は分類
自身の確率（定常数）、
0006836041.tif 000029
は各プロペラ翼数の類型
0006836041.tif 000030
の先験的確率、
0006836041.tif 000031
は所定Ｘの
0006836041.tif 000032
プロペラ翼数類型への従属確率、
0006836041.tif 000033
はプロペラ翼数類型
0006836041.tif 000034
のＸ発生確率、
0006836041.tif 000035
はプロペラ翼数類型
0006836041.tif 000036
の変調スペクトルにおける特徴
0006836041.tif 000037
発生の確率である。各
0006836041.tif 000038
を算出すると、
0006836041.tif 000039
が最大値であれば、Ｘは類型
0006836041.tif 000040
に従属する。
【００１５】
ナイーブベイズ推断法は、ベイズ定理（Ｂａｙｅｓ´  Ｔｈｅｏｒｅｍ）に基づき、各特
徴条件は相互独立的であると考えられる。事前に提供される訓練集合サンプルによって入
力から出力までの同時確率分布を学習し、学習により得られたモデルに基づいて、分類目
標Ｘを入力して、後験的確率Ｙが最大にできる出力を求めるものである。
軸周波数およびプロペラ翼数を取得すると、翼周波数はプロペラ翼数×軸周波数である。
従来の技術に比べて、本発明は下記の効果を有している。
本発明は、ローカルピークの探索に限定条件を加味することにより、位置が近い両ローカ
ルピークを識別しないようにできる。軸周波数の確定に限定条件を加味することにより、
ノイズの影響で過低軸周波数を識別する誤り判断を避けることができる。共振周波数位置
のピークを確定するときは、一つの周波数区間内でピーク平均値を求める。プロペラ翼数
を判断するときは、ナイーブベイズの推断法を使用する。その優勢は小サンプルに適用可
能であるとともに、従来の技術に述べた従来の識別規則により解決できない問題を解決で
きる。最後に、様々な変調スペクトル構成から軸周波数、翼周波数やプロペラ翼数特徴量
を抽出することができる。
【図面の簡単な説明】
【００１６】
【図１】輻射ノイズ変調に基づくプロペラ翼数特徴の抽出方法の模式的フローチャート図である。
【図２】本発明の実施例に係る特徴抽出目標の変調スペクトルである。
【図３】本発明の実施例に係るＭＡＴＬＡＢ（登録商標）のｆｉｎｄｐｅａｋｓ関数により得られるピークを示す図である。
【図４】本発明の実施例に係る判断により得られる軸周波数の周波数を示す図である。
【図５】本発明の実施例に係る最後に識別した結果図である。
【発明を実施するための形態】
【００１７】
以下、図面および実施例を参考しながら本発明をさらに詳しく説明する。特筆すべきのは
、実施例の作用は、本発明を容易に理解させることにあり、発明の制限としてはならない
。
【００１８】
図１に示されるように、輻射ノイズ変調に基づくプロペラ翼数特徴の抽出方法は下記ステ
ップを含む。
Ｓ０１：本例では、ある商船からの輻射ノイズを用いる。商船のプロペラ翼数は３、商船
のプロペラ回転数は１１１回転／分、翼周波数は５．５５Ｈｚ、軸周波数は１．８５Ｈｚ
とする。商船の輻射ノイズを短時間フーリエ変換することにより変調スペクトルを取得す
る。変調スペクトル図においてローカルピーク、およびローカルピーク位置の共振周波数
を探索する。
【００１９】
特徴抽出目標の変調スペクトルは図２に示される。該ステップにおいては、Ｍａｔｌａｂ
中のｆｉｎｄｐｅａｋｓ関数により平均コヒーレント値のローカルピーク、および対応の
循環周波数を探索する。
隣り合う両ローカルピークの共振周波数の差は下式の関係を満たす。
0006836041.tif 000041
式中、
0006836041.tif 000042
はｎ番目ローカルピーク位置の循環周波数の値である。ピークの探索結果は図３に示され
る。
【００２０】
Ｓ０２：軸周波数およびその共振周波数を確定する。もし１番目ローカルピーク位置の共
振周波数値が０．９Ｈｚより小さいければ、それを除去し、２番目共振周波数を軸周波数
とする。軸周波数は
0006836041.tif 000043
とする。図４に示されるように、本実施例において
0006836041.tif 000044
とし、１番目ローカルピーク位置の共振周波数が０．９Ｈｚ未満でないという条件を満た
し、１番目位置の共振周波数が軸周波数であり、
0006836041.tif 000045
となっている。ローカルピークの数により、共振周波数の数Ｎを確定する。本例では、７
つのピークがあるから、Ｎ＝７となる。
【００２１】
Ｓ０３：軸周波数および翼周波数の逓倍関係により、他の共振周波数の値を確定する。
軸周波数および翼周波数は逓倍関係：
0006836041.tif 000046
（ただし、
0006836041.tif 000047
はｎ番目の共振周波数、
0006836041.tif 000048
とする）ので、
逓倍関係に基づいて、軸周波数が既定であれば、翼周波数位置の可能の共振周波数を確定
することができる。
【００２２】
Ｓ０４：各共振周波数位置の平均スペクトルコヒーレント値を確定する。
誤差が存在するので、各共振周波数位置のコヒーレント値を計算するときは、
0006836041.tif 000049
区間内の平均コヒーレント値の平均を求めて、
0006836041.tif 000050
位置の平均コヒーレント値
0006836041.tif 000051
を得られる。ただし
0006836041.tif 000052
はサンプリング点
0006836041.tif 000053
前の５番目サンプリング点から
0006836041.tif 000054
後の５番目サンプリング点までの区間を示す。
【００２３】
Ｓ０５：ナイーブベイズ推断法に基づいて、プロペラ翼数を取得する。分類問題を解決す
る。即ちナイーブベイズの推断法に基づいて、得られた共振周波数値と対応の平均コヒー
レント値間の関係により、該変調スペクトル図がどの種類のプロペラ翼数の集合に従属す
るかを確定する。最後には、プロペラ翼数を確定する。本例では、得られた変調スペクト
ル図における幅値間の関係は
0006836041.tif 000055
ある。ナイーブベイズの推断法に基づくナイーブベイズの公式：
0006836041.tif 000056
本例では、Ｙ＝｛
0006836041.tif 000057
｝（Ｙはすべての可能のプロペラ翼数の集合、
0006836041.tif 000058
はその一つのプロペラ翼数、
0006836041.tif 000059
はプロペラ翼数が３、
0006836041.tif 000060
はプロペラ翼数が４、
0006836041.tif 000061
はプロペラ翼数が５、
0006836041.tif 000062
はプロペラ翼数が６、
0006836041.tif 000063
はプロペラ翼数が７である）、Ｘ＝｛
0006836041.tif 000064
｝（Ｘは分類目標の変調スペクトル、
0006836041.tif 000065
は該変調スペクトルにおける各共振周波数位置の振幅の大小関係である。またＰ（Ｘ）は
分類自身の確率である（定乗数）。
0006836041.tif 000066
は各プロペラ翼数類型の先験的確率、即ち
0006836041.tif 000067
の確率である。
0006836041.tif 000068
は所定Ｘの
0006836041.tif 000069
プロペラ翼数類型への従属確率である。
0006836041.tif 000070
はプロペラ翼数類型
0006836041.tif 000071
のＸの発生確率である。
0006836041.tif 000072
はプロペラ翼数類型
0006836041.tif 000073
の変調スペクトルにおける特徴
0006836041.tif 000074
発生の確率である。各
0006836041.tif 000075
を計算すると、
0006836041.tif 000076
が最大値であれば、Ｘが類型
0006836041.tif 000077
に従属すると思われる。
【００２４】
プロペラ翼数の変調スペクトル図の大量サンプルが既知であり、或いは表１の規則により
模擬したサンプルによれば
0006836041.tif 000078
が既定量であるから、
0006836041.tif 000079
を算出できる。本例では、得られる変調スペクトル図における幅値間の関係は
0006836041.tif 000080
がある。該変調スペクトル図における振幅間の関係はＸのある特徴
0006836041.tif 000081
により示される。最後に計算により
0006836041.tif 000082
が最大値となる。したがって、プロペラ翼数が３であると確定できる。最後に翼周波数は
プロペラ翼数×軸周波数となっており、即ち翼周波数は５．５６５Ｈｚとなる。最後の識
別結果は図５に示される。
【００２５】
上記の実施例の説明は、本発明の方法およびその精神を理解するためのものである。当業
者にとって、本発明の原理から逸脱することなく、本発明をいくつかの改良および修正を
加えて実施することもでき、これらの改良および修正はすべて本発明の保護範囲に含まれ
ることは言うまでもない。
【要約】      （修正有）【課題】簡単な操作で船舶の種類を判断するために、音声センサによるプロペラの軸周波数、翼周波数およびプロペラ翼数の特徴量情報を抽出する方法を提供する。【解決手段】（１）船舶の輻射ノイズ信号を採集し、フーリエ変換により変調スペクトル図を取得するステップと、（２）変調図にローカルピークを探索し、ローカルピーク位置における共振周波数を取得するステップＳ０１と、（３）変調図における１番目ローカルピークの共振周波数を確定し、軸周波数を確定するとともに、ローカルピークの数に基づいて共振周波数の数Ｎを確定するステップＳ０２と、（４）軸周波数と翼周波数間の逓倍関係を確定し、他の共振周波数を確定するステップＳ０３と、（５）各高調波周波数位置の平均スペクトルコヒーレント値を確定しＳ０４、ナイーブベイズの推断法によりプロペラ翼数を取得するＳ０５、とともに、最後に翼周波数を確定するステップと、を含む。【選択図】図１【図１】
0006836041.tif 000116
【図２】
0006836041.tif 000117
【図３】
0006836041.tif 000118
【図４】
0006836041.tif 000119
【図５】
0006836041.tif 000120