WO2016133007A1

WO2016133007A1 - 音場補正装置、音場補正方法及び音場補正プログラム

Info

Publication number: WO2016133007A1
Application number: PCT/JP2016/054062
Authority: WO
Inventors: 一智福江; 橋本　武志; 哲生渡邉; 隆富熊谷; 藤田　康弘
Original assignee: クラリオン株式会社
Priority date: 2015-02-16
Filing date: 2016-02-12
Publication date: 2016-08-25
Also published as: JP6541364B2; JP2016152422A; US10270412B2; EP3261362A4; US20180034434A1; CN107409256B; CN107409256A; EP3261362A1; EP3261362B1

Abstract

　音場補正装置を、所定の音声信号に基づいてパラメトリックイコライザによる音場補正の目標となる目標振幅特性及び該目標振幅特性と同等若しくはそれよりもなだらかな振幅特性を持つ複数の副目標振幅特性を算出する目標算出手段と、算出された目標振幅特性及び複数の副目標振幅特性に基づいて音場を補正するための補正振幅特性を算出する補正振幅特性算出手段と、算出された補正振幅特性に基づいてパラメトリックイコライザを設定する設定手段とを備える構成とする。

Description

音場補正装置、音場補正方法及び音場補正プログラム

　本発明は、音場補正装置、音場補正方法及び音場補正プログラムに関する。

　スピーカより放出される音声（以下、「音場測定音」と記す。）の音声信号と、聴取位置に配置されたマイクで収録された音声（以下、「マイク収録音」と記す。）の音声信号との違いに基づいて、音場測定音の音声信号の周波数帯域毎の出力レベルを調整する音場補正装置が知られている。この種の音場補正装置の具体的構成は、例えば特開２００７－２９５５２８号公報（以下、「特許文献１」と記す。）や特開２００８－２４５１２３号公報（以下、「特許文献２」と記す。）に記載されている。

　特許文献１に記載の音場補正装置は、補正対象の周波数範囲に含まれるエリア内における目標特性からのゲイン差が最大となる周波数に基づいて中心周波数を決定し、決定された中心周波数におけるゲイン値をその中心周波数での目標特性とのゲイン差に基づいて決定し、更に、所定の候補の中からＱ値を決定する構成となっている。しかし、特許文献１に記載の構成では、エリア内における目標特性とのゲイン差の形状や重心が考慮されていないため、該エリアを精度良く補正することが難しい場合がある。

　特許文献２に記載の音場補正装置は、マイク収録音から周波数帯域毎の信号補正レベルを算出し、算出された信号補正レベルを周波数帯域順に接続した補正レベルの波形の変極点を基準として、周波数帯域を複数のグループに区分し、区分されたグループ単位でレベル補正を行う。特許文献２に記載の構成では、グループ内の重心位置が考慮されている点ではレベル補正を精度良く行うのに好適であるが、グループ分割処理の結果によっては例えばパラメトリックイコライザの特性が先鋭となり、少ないバンド数では十分な補正効果が得られないという問題が指摘される。

　本発明はこのような事情を鑑みてなされたものであり、その目的とするところは、パラメトリックイコライザのバンド数が少ない場合であっても十分な音場補正効果を得るのに好適な音場補正装置、音場補正方法及び音場補正プログラムを提供することである。

　本発明の一実施形態に係る音場補正装置は、所定の音声信号に基づいてパラメトリックイコライザによる音場補正の目標となる目標振幅特性及び該目標振幅特性と同等若しくはそれよりもなだらかな振幅特性を持つ複数の副目標振幅特性を算出する目標算出手段と、算出された目標振幅特性及び複数の副目標振幅特性に基づいて音場を補正するための補正振幅特性を算出する補正振幅特性算出手段と、算出された補正振幅特性に基づいてパラメトリックイコライザを設定する設定手段とを備える。

　また、本発明の一実施形態において、目標算出手段は、音声信号に基づいてその振幅特性を算出し、算出された振幅特性に対して第一の分解能で平均化処理を施し、平均化処理された振幅特性を基に目標振幅特性を算出する構成としてもよい。また、目標算出手段は、音声信号に基づいてその振幅特性を算出し、算出された振幅特性に対して第一の分解能と同等若しくはそれよりも粗い分解能であってそれぞれ異なる複数種類の分解能で平均化処理を施し、異なる分解能で平均化処理された各振幅特性を基に複数の副目標振幅特性を算出する構成としてもよい。

　また、本発明の一実施形態において、補正振幅特性算出手段は、所定の条件に基づいて各副目標振幅特性を周波数領域において複数のグループに分割し、分割された各グループの優先度を算出し、算出された優先度が最も高いグループを該副目標振幅特性毎に選択し、選択された周波数領域のグループに基づいて補正候補のパラメトリックイコライザのパラメータを副目標振幅特性毎に算出し、算出された副目標振幅特性毎のパラメータに基づいて補正対象の該パラメトリックイコライザのパラメータを取得する構成としてもよい。この場合、設定手段は、補正対象のパラメトリックイコライザのパラメータに基づいて該パラメトリックイコライザを設定する。

　また、本発明の一実施形態において、補正振幅特性算出手段は、補正候補のパラメトリックイコライザのパラメータに基づいて補正候補の振幅特性を副目標振幅特性毎に算出し、算出された各補正候補の振幅特性のうち、目標振幅特性との差分が最も小さいもののパラメータを補正対象のパラメトリックイコライザのパラメータとして取得する構成としてもよい。

　また、本発明の一実施形態において、補正振幅特性算出手段は、補正対象のパラメトリックイコライザのパラメータが取得されると、副目標振幅特性毎に算出された補正候補の振幅特性の中で目標振幅特性との差分が最も小さい補正対象の振幅特性と、該目標振幅特性との差分を新たな目標振幅特性として算出し、且つ該補正対象の振幅特性と各副目標振幅特性との差分を新たな複数の副目標振幅特性として算出し、取得された補正対象のパラメトリックイコライザの数が所定数に達していなければ、新たな目標振幅特性及び新たな複数の副目標振幅特性に基づいて、該補正対象のパラメトリックイコライザのパラメータを更に取得する構成としてもよい。

　また、本発明の一実施形態において、パラメトリックイコライザのパラメータは、例えば、パラメトリックイコライザの中心周波数、ゲイン及び周波数帯域幅を含む。

　また、本発明の一実施形態に係る音場補正方法は、所定の音声信号に基づいてパラメトリックイコライザによる音場補正の目標となる目標振幅特性及び該目標振幅特性と同等若しくはそれよりもなだらかな振幅特性を持つ複数の副目標振幅特性を算出する目標算出ステップと、算出された目標振幅特性及び複数の副目標振幅特性に基づいて音場を補正するための補正振幅特性を算出する補正振幅特性算出ステップと、算出された補正振幅特性に基づいてパラメトリックイコライザを設定する設定ステップとを含む方法である。

　また、本発明の一実施形態に係る音場補正プログラムは、上記の音場補正方法をコンピュータに実行させるためのプログラムである。

　本発明の一実施形態によれば、パラメトリックイコライザのバンド数が少ない場合であっても十分な音場補正効果を得るのに好適な音場補正装置、音場補正方法及び音場補正プログラムが提供される。

本発明の一実施形態に係る音響システムの構成を示すブロック図である。本発明の一実施形態に係る音響システムにおいて実行される音場測定処理のフローチャートを示す図である。本発明の一実施形態に係る音響システムにおいて実行される音場補正処理のフローチャートを示す図である。本発明の一実施形態に係る音響システムにおいて実行される音場補正処理のフローチャートを示す図である。図３Ａの処理ステップＳ２２（インパルス応答の振幅特性の算出）にて算出されるインパルス応答の振幅特性の一例を示す図である。図３Ａの処理ステップＳ２３（目標振幅特性の生成）にて算出されるインパルス応答の正規化振幅特性の一例を示す図である。図３Ａの処理ステップＳ２３（目標振幅特性の生成）にて生成される目標振幅特性の一例を示す図である。図３Ａの処理ステップＳ２４（複数の副目標振幅特性の生成）にて生成される副目標振幅特性の一例を示す図である。図３Ａの処理ステップＳ２５（各副目標振幅特性に対するグループ分割処理）のフローチャートを示す図である。図３Ａの処理ステップＳ２５（各副目標振幅特性に対するグループ分割処理）のフローチャートを示す図である。図８Ａの処理ステップＳ２５ｂにて選択される正の暫定グループの一例を模式的に示す図である。図３Ａの処理ステップＳ２５（各副目標振幅特性に対するグループ分割処理）の実行結果（各副目標振幅特性に対するグループ分割処理の結果）を例示する図である。図３Ａの処理ステップＳ２６（各グループの優先度の算出）にて優先度が算出される本グループの振幅特性を模式的に示す図である。図３Ａの処理ステップＳ２７（優先度に基づくグループの選択）にて選択された補正対象候補グループの一例を示す図である。各副目標振幅特性と各補正候補ＰＥＱの振幅特性とを示す図である。各目標振幅特性と各補正候補ＰＥＱの振幅特性とを示す図である。図３Ｂの処理ステップＳ３２（新たな副目標振幅特性の算出）にて算出される新たな副目標振幅特性を示す図である。図３Ｂの処理ステップＳ３３（新たな目標振幅特性の算出）にて算出される新たな目標振幅特性を示す図である。本発明の一実施形態に係る制御部の内部メモリに記録される各補正候補ＰＥＱのパラメータを示す図（図１７（ａ））及び目標振幅特性及びこれに近似する補正振幅特性を示す図（図１７（ｂ））である。従来（特許文献２）、本実施形態のそれぞれにおいて、ＰＥＱバンド数と目標振幅特性との誤差を比較した結果を示す図である。

　以下、本発明の実施形態について図面を参照しながら説明する。なお、以下においては、本発明の一実施形態として、車室内に配置された音響システムを例に取り説明する。

［音響システム１の概要説明］
　図１は、本発明の一実施形態に係る音響システム１の構成を示すブロック図である。本実施形態に係る音響システム１は、車室内という聴取環境における音場を補正、より詳細には、車室内に配置された複数のスピーカに対して出力される音声信号の周波数帯域毎の信号レベルを調整する、音場補正機能（音場補正装置）を搭載したシステムである。

　なお、音響システム１における各種処理は、音響システム１に備えられるソフトウェアとハードウェアとが協働することにより実行される。音響システム１に備えられるソフトウェアのうち少なくともＯＳ（Operating System）部分は、組み込み系システムとして提供されるが、それ以外の部分、例えば、音場補正を実行するためのソフトウェアモジュールについては、ネットワーク上で配布可能な又はメモリカード等の記録媒体にて保持可能なアプリケーションとして提供されてもよい。すなわち、本実施形態に係る音場補正機能は、音響システム１に予め組み込まれた機能であっても、ネットワーク経由や記録媒体経由で追加可能な機能であってもよい。

　図１に示されるように、音響システム１は、音場装置１０、マイクロフォン１１６、スピーカＦＣ、ＦＲ、ＦＬ、ＲＲ、ＲＬ、ＳＷを備えている。スピーカＦＣは、車室内の前部中央位置に配置されたセンタスピーカであり、スピーカＦＲは、車室内の前部右側方位置に配置された右フロントスピーカであり、スピーカＦLは、車室内の前部左側方位置に配置された左フロントスピーカであり、スピーカＲＲは、車室内の後部右側方位置に配置された右リアスピーカであり、スピーカＲＬは、車室内の後部左側方位置に配置された左リアスピーカであり、スピーカＳＷは、車室内の後部中央位置に配置されたサブウーファである。

　音場装置１０は、制御部１００、表示部１０２、操作部１０４、測定信号発生部１０６、記録媒体再生部１０８、ＰＥＱ（Parametric Equalizer）部１１０、Ｄ／Ａコンバータ１１２、パワーアンプ１１４、マイクアンプ１１８、Ａ／Ｄコンバータ１２０、信号収録部１２２及び計算部１２４を備えている。

［音場測定処理］
　図２は、本発明の一実施形態に係る音響システム１において実行される音場測定処理のフローチャートを示す図である。本フローチャートに示される音場測定処理をはじめとする、音響システム１内での各種処理は、制御部１００の制御下で実行される。制御部１００は、表示部１０２に対する所定のタッチ操作又は操作部１０４に対する所定の操作を受けると、音場測定条件の入力画面を表示部１０２に表示する。これにより、本フローチャートに示される音場測定処理が開始される。

［図２のＳ１１（音場測定条件の設定）］
　本処理ステップＳ１１では、表示部１０２に表示された入力画面に対してユーザによる音場測定条件の入力が行われると、入力された音場測定条件が設定される。ここで入力される音場測定条件は、例えばチャンネル数（又は測定対象のスピーカの指定）である。本実施形態に係る音響システム１には６つのスピーカが搭載されているため、入力可能なチャンネル数は最大「６」である。本実施形態では、６つのチャンネル数が入力されたものとして説明を行う。

［図２のＳ１２（音場測定音の再生）］
　本処理ステップＳ１２では、測定信号発生部１０６が所定の測定信号を発生させる。発生された測定信号は、例示的には、Ｍ系列符号（Maximal length sequence）やＴＳＰ信号（Time Stretched Pulse）であり、Ｄ／Ａコンバータ１１２及びパワーアンプ１１４を介して、処理ステップＳ１１（音場測定条件の設定）にて設定された各スピーカＦＣ、ＦＲ、ＦＬ、ＲＲ、ＲＬ、ＳＷに所定の時間間隔を空けて順次入力される。これにより、所定の音場測定音が所定の時間間隔を空けて各スピーカＦＣ、ＦＲ、ＦＬ、ＲＲ、ＲＬ、ＳＷから順次再生される。

［図２のＳ１３（インパルス応答の記録）］
　本実施形態において、車室内には、４つの座席（運転席、助手席、左右一対の後部座席）が設置されている。マイクロフォン１１６は、これら４つの座席の各乗員に適切な音場補正を行うため、各座席から等距離の位置（４つの座席の中央の位置）に配置されている。マイクロフォン１１６は、音場補正の効果を与えたい乗員（換言すると、車室内の位置）に応じて設置される位置が変わる。例えば運転手に最適な音場補正を行いたい場合、マイクロフォン１１６は運転席に設置される。

　本処理ステップＳ１３では、処理ステップＳ１２（音場測定音の再生）にて再生された音場測定音がマイクロフォン１１６にて収音されて、マイクアンプ１１８及びＡ／Ｄコンバータ１２０を介して信号収録部１２２に入力される。信号収録部１２２では、インパルス応答が演算される。インパルス応答は、例えば、入力された音場測定音とリファレンスの測定信号（ＴＳＰ信号等）を時間軸上で逆転させた逆特性のリファレンス信号をフーリエ変換して周波数上で乗算し、乗算された値を逆フーリエ変換することで求まる。演算されたインパルス応答は、制御部１００の内部メモリ１００Ｍに記録される。

［図２のＳ１４（音場測定処理の完了判定）］
　処理ステップＳ１２（音場測定音の再生）及びＳ１３（インパルス応答の記録）は、処理ステップＳ１１（音場測定条件の設定）にて設定されたスピーカ毎に実行される。処理ステップＳ１１（音場測定条件の設定）にて設定された全てのスピーカについて、処理ステップＳ１２（音場測定音の再生）及びＳ１３（インパルス応答の記録）が実行されると（Ｓ１４：ＹＥＳ）、本フローチャートに示される音場測定処理が完了する。

［音場補正処理］
　図３（図３Ａ及び図３Ｂ）は、本発明の一実施形態に係る音響システム１において実行される音場補正処理のフローチャートを示す図である。図３に示される音場補正処理は、図２に示される音場測定処理の完了次第開始され、且つ各スピーカに対して実行される。

［図３ＡのＳ２１（補正条件の設定）］
　本処理ステップＳ２１では、補正条件を設定するための設定画面が表示部１０２に表示される。表示部１０２に表示された設定画面に対してユーザによる補正条件の入力が行われると、入力された補正条件が設定される。ここで設定される補正条件は、ＰＥＱバンド数及び補正周波数範囲である。ＰＥＱバンド数は、一台のスピーカに割り当てられるパラメトリックイコライザの数を示しており、本実施形態では「７」である。補正周波数範囲は、補正される周波数の範囲を示しており、スピーカの再生可能周波数等に基づいてスピーカ毎に設定される。

［図３ＡのＳ２２（インパルス応答の振幅特性の算出）］
　本処理ステップＳ２２では、処理ステップＳ１３（インパルス応答の記録）にて記録されたインパルス応答が読み出されて、計算部１２４によりフーリエ変換されて、周波数領域におけるインパルス応答の振幅特性が算出される。図４に、ここで算出されるインパルス応答の振幅特性の一例を示す。なお、図４中、縦軸、横軸はそれぞれ、信号レベル（Power（単位：ｄＢ））、周波数（Frequency（単位：Ｈｚ））を示す。Powerは、振幅を自乗したものである。また、人間の聴覚特性は、周波数に対して対数的である。横軸の周波数は、人間の聴覚特性に合わせて対数表示となっている。

［図３ＡのＳ２３（目標振幅特性の生成）］
　本処理ステップＳ２３では、計算部１２４により、処理ステップＳ２２（インパルス応答の振幅特性の算出）にて算出された振幅特性に対し、１サンプルずつシフトしながら所定のサンプル数内の平均値が算出されて、スムージング（平均化処理）が行われる。ここで、平均化処理は、聴覚の周波数分解能として知られている１／３オクターブバンド幅と同程度の分解能で行われる。

　次いで、基準バンド幅（本実施形態では５００Ｈｚ～３０００Ｈｚ）内における信号レベルを基に、インパルス応答の正規化振幅特性が算出される。図５に、ここで算出されるインパルス応答の正規化振幅特性の一例を示す。

　算出されたインパルス応答の正規化振幅特性の信号レベルの符号が反転されて、所定の重み付け（例えば音響システム１にて生成すべき音場の振幅特性の重み付け）が行われる。重み付けされた振幅特性のうち、処理ステップＳ２１（補正条件の設定）にて設定された補正周波数範囲における振幅特性が、音場補正の目標となる目標振幅特性として取得される。図６に、ここで生成され取得される目標振幅特性の一例を示す。

［図３ＡのＳ２４（複数の副目標振幅特性の生成）］
　本処理ステップＳ２４では、計算部１２４により、処理ステップＳ２３（目標振幅特性の生成）にて生成された目標振幅特性と同様の方法で複数の副目標振幅特性が生成される。本実施形態では、３種類の副目標振幅特性が生成される。１つは、１/１オクターブバンド幅と同程度の分解能で平均化処理を施すことによって生成されたものであり、説明の便宜上、「第一の副目標振幅特性」と記す。他の１つは、１/２オクターブバンド幅と同程度の分解能で平均化処理を施すことによって生成されたものであり、説明の便宜上、「第二の副目標振幅特性」と記す。残りの１つは、目標振幅特性と同じく１/３オクターブバンド幅と同程度の分解能で平均化処理を施すことによって生成されたものであり、説明の便宜上、「第三の副目標振幅特性」と記す。広いオクターブバンド幅の分解能で平均化処理を施すことによって生成された副目標振幅特性ほど、その振幅特性がなだらかになる。

　図７（ａ）に、第一の副目標振幅特性の一例を示し、図７（ｂ）に、第二の副目標振幅特性の一例を示し、図７（ｃ）に、第三の副目標振幅特性の一例を示す。なお、図７（ａ）～図７（ｃ）の各図は、図面内に３つのグラフを記載する便宜上、縦軸の縮尺が他の図（図４～図６等）と大きく異なっている。第三の副目標振幅特性は、図７（ｃ）に示されるように、目標振幅特性と同等の振幅特性を有している。これに対し、第二の副目標振幅特性は、図７（ｂ）に示されるように、目標振幅特性（又は第三の副目標振幅特性）よりもなだらかな振幅特性を有している。また、第一の副目標振幅特性は、図７（ａ）に示されるように、第二の副目標振幅特性よりもなだらかな振幅特性を有している。すなわち、本処理ステップＳ２４では、目標振幅特性と同等若しくはそれよりもなだらかな振幅特性を持つ複数の副目標振幅特性が生成される。

［図３ＡのＳ２５（各副目標振幅特性に対するグループ分割処理）］
　本処理ステップＳ２５では、計算部１２４により、処理ステップＳ２４（複数の副目標振幅特性の生成）にて生成された各副目標振幅特性に対してグループ分割処理が行われる。図８（図８Ａ及び図８Ｂ）に、本処理ステップＳ２５にて行われるグループ分割処理のフローチャートを示す。

・図８ＡのＳ２５ａ
　本処理ステップＳ２５ａでは、処理ステップＳ２４（複数の副目標振幅特性の生成）にて生成された副目標振幅特性毎に、信号レベルの符号による暫定的なグループ分割が行われる。より詳細には、副目標振幅特性は、正又は負の信号レベルが連続する周波数領域（すなわち、信号レベルの反転の無い周波数領域）毎にグループ分割される。

・図８ＡのＳ２５ｂ
　本処理ステップＳ２５ｂでは、処理ステップＳ２５ａにて暫定的に分割された暫定グループ群の中から１つの暫定グループが選択される。

・図８ＡのＳ２５ｃ
　本処理ステップＳ２５ｃでは、処理ステップＳ２５ｂにて選択された暫定グループ内の所定の特異点について検出処理が行われる。本処理ステップＳ２５ｃにて特異点が検出されなかった場合（Ｓ２５ｃ：ＮＯ）は、処理ステップＳ２５ｂにて選択された暫定グループが本グループとして決定されて、処理ステップＳ２５ｊに進む。本処理ステップＳ２５ｃにて特異点が検出された場合（Ｓ２５ｃ：ＹＥＳ）は、処理ステップＳ２５ｄに進む。この場合、正の暫定グループ（正の信号レベルが連続する周波数領域）では、少なくとも１つの極小値が検出され、負の暫定グループ（負の信号レベルが連続する周波数領域）では、少なくとも１つの極大値が検出される。

・図８ＡのＳ２５ｄ
　本処理ステップＳ２５ｄでは、処理ステップＳ２５ｃにて検出された特異点のうち絶対値が最も小さいものが選択される。

・図８ＡのＳ２５ｅ
　本処理ステップＳ２５ｅでは、処理ステップＳ２５ｄにて選択された特異点を境に、暫定グループが仮分割される。

・図８ＡのＳ２５ｆ
　本処理ステップＳ２５ｆでは、処理ステップＳ２５ｅにて仮分割された各仮分割グループ内で最大となる絶対値が検出される。

・図８ＢのＳ２５ｇ
　本処理ステップＳ２５ｇでは、処理ステップＳ２５ｄにて選択された特異点の絶対値（暫定グループ内で最小となる絶対値）と、処理ステップＳ２５ｆにて検出された各仮分割グループ内で最大となる絶対値との差分が所定の閾値以上か否かが判定される。

・図８ＢのＳ２５ｈ
　本処理ステップＳ２５ｈは、処理ステップＳ２５ｇにて上記差分が所定の閾値以上と判定された場合（Ｓ２５ｇ：ＹＥＳ）に実行される。本処理ステップＳ２５ｈでは、暫定グループが、処理ステップＳ２５ｄにて選択された特異点を境に分割されて、２つの本グループとして決定される。

・図８ＢのＳ２５ｉ
　本処理ステップＳ２５ｉは、処理ステップＳ２５ｇにて上記差分が所定の閾値未満と判定された場合（Ｓ２５ｇ：ＮＯ）に実行される。本処理ステップＳ２５ｉでは、処理ステップＳ２５ｄにて選択された特異点が処理ステップＳ２５ｃにおける検出対象から除外される。次いで、本フローチャートに示されるグループ分割処理は、処理ステップＳ２５ｃに戻る。

・図８ＢのＳ２５ｊ
　本処理ステップＳ２５ｊでは、処理ステップＳ２５ａにて分割された全ての暫定グループに対して、処理ステップＳ２５ｃ以降の処理（本グループとして決定されるための処理）が実行されたか否かが判定される。本フローチャートに示されるグループ分割処理は、未処理の暫定グループが残っている場合（Ｓ２５ｊ：ＮＯ）、処理ステップＳ２５ｂに戻り、全ての暫定グループに対する処理が完了している場合（Ｓ２５ｊ：ＹＥＳ）、終了する。

　図９は、処理ステップＳ２５ｂにて選択される正の暫定グループの一例を模式的に示す図である。図９の例では、処理ステップＳ２５ｃにて極小値ｍｉｎ１及びｍｉｎ２が検出される。処理ステップＳ２５ｄでは、絶対値の小さい極小値ｍｉｎ２が選択される。処理ステップＳ２５ｅでは、極小値ｍｉｎ２を境に暫定グループが仮分割される。処理ステップＳ２５ｆでは、仮分割された各仮分割グループＧ１内、Ｇ２内のそれぞれで、最大となる極大値ｍａｘ１、ｍａｘ２が検出される。処理ステップＳ２５ｇでは、極小値ｍｉｎ２と極大値ｍａｘ１との差分値、極小値ｍｉｎ２と極大値ｍａｘ２との差分値の何れもが所定の閾値以上か否かが判定される。図９の例では、何れの差分値も所定の閾値以上であるため、仮分割グループＧ１、Ｇ２が本グループとして決定される。

　図１０に、処理ステップＳ２５ｇにて判定に用いられる閾値を１とした場合における、各副目標振幅特性に対するグループ分割処理の結果を例示する。図１０（ａ）～（ｃ）の各図中、実線で示される振幅特性は、信号レベルが正の本グループであり、破線で示される振幅特性は、信号レベルが負の本グループである。図１０（ａ）、図１０（ｂ）、図１０（ｃ）はそれぞれ、第一、第二、第三の副目標振幅特性に対するグループ分割処理の結果を示す。図１０（ａ）、図１０（ｂ）、図１０（ｃ）の各図に示されるように、第一、第二、第三の副目標振幅特性はそれぞれ、９個、１４個、１５個の本グループに分割されている。

［図３ＡのＳ２６（各グループの優先度の算出）］
　本処理ステップＳ２６では、計算部１２４により、本グループの信号レベルに基づいて本グループの優先度が副目標振幅特性毎に算出される。図１１に、優先度が算出される本グループの振幅特性を模式的に示す。図１１では、各信号レベルｇ_ｎ－１、ｇ_０、・・・ｇ_ｋ、ｇ_ｋ＋１に対応する周波数がｆ_ｎ－１、ｆ_０、・・・ｆ_ｋ、ｆ_ｋ＋１と定義される。この場合、各副目標振幅特性の各本グループの優先度は、次式１により算出される。
（式１）

［図３ＡのＳ２７（優先度に基づくグループの選択）］
　本処理ステップＳ２７では、計算部１２４により、処理ステップＳ２６（各グループの優先度の算出）での算出結果に基づいて優先度の最も高い本グループが各副目標振幅特性で選択される。以下、説明の便宜上、第一の副目標振幅特性で優先度の最も高い本グループに符号Ｇ_ｒ１／１を付し、第二の副目標振幅特性で優先度の最も高い本グループに符号Ｇ_ｒ１／２を付し、第三の副目標振幅特性で優先度の最も高い本グループに符号Ｇ_ｒ１／３を付す。

　本処理ステップＳ２７にて選択された本グループＧ_ｒ１／１、Ｇ_ｒ１／２、Ｇ_ｒ１／３が補正対象候補のグループとなる。図１２（ａ）、図１２（ｂ）、図１２（ｃ）にそれぞれ、本処理ステップＳ２７にて選択された補正対象候補グループＧ_ｒ１／１、Ｇ_ｒ１／２、Ｇ_ｒ１／３の一例を示す。

［図３ＢのＳ２８（中心周波数の算出）］
　本処理ステップＳ２８では、計算部１２４により、処理ステップＳ２７（優先度に基づくグループの選択）にて選択された補正対象候補グループに基づいて補正候補ＰＥＱの中心周波数（重心位置）が副目標振幅特性毎に算出される。補正対象候補グループ内の信号レベルと周波数とが図１１のように定義される場合、補正候補ＰＥＱの中心周波数は、次式２により算出される。
（式２）

［図３ＢのＳ２９（ゲインの算出）］
　本処理ステップＳ２９では、計算部１２４により、処理ステップＳ２７（優先度に基づくグループの選択）にて選択された補正対象候補グループに基づいて補正候補ＰＥＱのゲインが副目標振幅特性毎に算出される。補正対象候補グループ内の信号レベルと周波数とが図１１のように定義される場合、補正候補ＰＥＱのゲインは、次式３により算出される。
（式３）

［図３ＢのＳ３０（周波数帯域幅の算出）］
　本処理ステップＳ３０では、計算部１２４により、処理ステップＳ２７（優先度に基づくグループの選択）にて選択された補正対象候補グループに基づいて補正候補ＰＥＱの周波数帯域幅が副目標振幅特性毎に算出される。補正対象候補グループ内の信号レベルと周波数とが図１１のように定義される場合、補正候補ＰＥＱの周波数帯域幅は、次式４により算出される。
（式４）

　補正候補ＰＥＱの中心周波数、ゲイン及び周波数帯域幅が算出されることにより、補正候補ＰＥＱの振幅特性が得られる。図１３（ａ）に、第一の副目標振幅特性（太実線）と、第一の副目標振幅特性から生成された第一の補正候補ＰＥＱの振幅特性（細実線）とを示し、図１３（ｂ）に、第二の副目標振幅特性（太実線）と、第二の副目標振幅特性から生成された第二の補正候補ＰＥＱの振幅特性（細実線）とを示し、図１３（ｃ）に、第三の副目標振幅特性（太実線）と、第三の副目標振幅特性から生成された第三の補正候補ＰＥＱの振幅特性（細実線）とを示す。

［図３ＢのＳ３１（補正対象ＰＥＱの選択）］
　本処理ステップＳ３１では、計算部１２４により、目標振幅特性と第一から第三の各補正候補ＰＥＱの振幅特性とが比較されて、目標振幅特性との差分が最も小さい補正候補ＰＥＱが補正対象のＰＥＱの１つとして選択される。補正対象のＰＥＱの１つとして選択された補正対象ＰＥＱのパラメータ（中心周波数、ゲイン及び周波数帯域幅）は、制御部１００の内部メモリ１００Ｍに記録される。

　図１４（ａ）、図１４（ｂ）、図１４（ｃ）にそれぞれ、目標振幅特性（太実線）と、第一、第二、第三の補正候補ＰＥＱの振幅特性（細実線）とを示す。図１４の例において、目標振幅特性との差分が最も小さい補正候補ＰＥＱは、第二の補正候補ＰＥＱである。従って、ここでは、第二の補正候補ＰＥＱが補正対象ＰＥＱの１つとして選択される。

［図３ＢのＳ３２（新たな副目標振幅特性の算出）］
　本処理ステップＳ３２では、計算部１２４により、処理ステップＳ３１（補正対象ＰＥＱの選択）にて選択された補正対象ＰＥＱの振幅特性と各副目標振幅特性との差分が新たな副目標振幅特性として算出される。図１５（ａ）に、補正対象ＰＥＱの振幅特性と第一の副目標振幅特性との差分（新たな第一の副目標振幅特性）を示し、図１５（ｂ）に、補正対象ＰＥＱの振幅特性と第二の副目標振幅特性との差分（新たな第二の副目標振幅特性）を示し、図１５（ｃ）に、補正対象ＰＥＱの振幅特性と第三の副目標振幅特性との差分（新たな第三の副目標振幅特性）を示す。

［図３ＢのＳ３３（新たな目標振幅特性の算出）］
　本処理ステップＳ３３では、計算部１２４により、処理ステップＳ３１（補正対象ＰＥＱの選択）にて選択された補正対象ＰＥＱの振幅特性と目標振幅特性との差分が新たな目標振幅特性として算出される。図１６に、補正対象ＰＥＱの振幅特性と目標振幅特性との差分（新たな目標振幅特性）を示す。

［図３ＢのＳ３４（終了判定）］
　本処理ステップＳ３４では、計算部１２４により、処理ステップＳ３１（補正対象ＰＥＱの選択）の実行によりパラメータが記録された補正対象ＰＥＱの数が処理ステップＳ２１（補正条件の設定）にて設定されたＰＥＱバンド数に達したか否かが判定される。本フローチャートに示される音場補正処理は、制御部１００の内部メモリ１００Ｍに記録された補正対象ＰＥＱの数がＰＥＱバンド数に達したと判定された場合（Ｓ３４：ＹＥＳ）終了し、該補正対象ＰＥＱの数がＰＥＱバンド数に達していないと判定された場合（Ｓ３４：ＮＯ）、処理ステップＳ２５（各副目標振幅特性に対するグループ分割処理）に戻り、処理ステップＳ３２（新たな副目標振幅特性の算出）にて算出された新たな副目標振幅特性及び処理ステップＳ３２（新たな目標振幅特性の算出）にて算出された新たな目標振幅特性を用いて、処理ステップＳ２５以降の処理が繰り返される。

　処理ステップＳ２５～Ｓ３４がループする毎に、目標振幅特性及び各副目標振幅特性が更新されつつ、補正対象ＰＥＱのパラメータ（中心周波数、ゲイン及び周波数帯域幅）が制御部１００の内部メモリ１００Ｍに順次記録される。

　ＰＥＱ部１１０は、ＩＩＲ（Infinite Impulse Response）フィルタであり、パラメータ（中心周波数、ゲイン及び周波数帯域幅）を調整することが可能なイコライザを複数備えている。制御部１００の内部メモリ１００Ｍに記録された各補正対象ＰＥＱのパラメータがＰＥＱ部１１０に設定されることにより、目標振幅特性に近似する補正振幅特性が設定される。ＰＥＱ部１１０は、記録媒体再生部１０８にて読み出されるＣＤ（Compact Disc）やＤＶＤ（Digital Versatile Disc）等の音声信号に対し、補正振幅特性に応じた周波数帯域毎の信号レベルの調整を行って音場を補正する。

　図１７（ａ）に、制御部１００の内部メモリ１００Ｍに記録される各補正対象ＰＥＱのパラメータを示し、図１７（ｂ）に、目標振幅特性（太実線）及びこれに近似する補正振幅特性（細実線）を示す。図１７（ｂ）に示されるように、補正振幅特性は、目標振幅特性に近似しており、精度良く補正されたものとなっている。

　このように、本発明の一実施形態によれば、大雑把な（マクロ的な）補正から細かい（ミクロ的な）補正が順に行われることにより、ＰＥＱのバンド数が少ない場合であっても十分な音場補正効果が得られる。また、音場補正が簡易な構成で行われるため、音場補正の処理時間の短縮化に好適である。

　図１８（ａ）に、目標振幅特性の一例を示す。また、図１８（ｂ）に、従来（特許文献２）、本実施形態のそれぞれにおいて、図１８（ａ）に示されるものを目標振幅特性とした場合における、ＰＥＱバンド数と目標振幅特性との誤差を比較した結果を示す。図１８（ｂ）に示されるように、本実施形態では、従来と比べて音場補正の誤差が１５～２５％程度減少していることが判る。すなわち、本実施形態では、ＰＥＱのバンド数が少ない場合であっても十分な音場補正効果が得られることが判る。

　以上が本発明の例示的な実施形態の説明である。本発明の実施形態は、上記に説明したものに限定されず、本発明の技術的思想の範囲において様々な変形が可能である。例えば明細書中に例示的に明示される実施例等又は自明な実施例等を適宜組み合わせた内容も本願の実施形態に含まれる。

　上記の実施形態では、音響システムは車室内の音場を補正するが、本発明はこれに限らない。別の実施形態では、音響システムは、住宅内など、他の特定空間の音場を補正するものであってもよい。また、音響システムは、上記の実施形態では複数のスピーカを備えるが、別の実施形態では単一のスピーカを備えるものであってもよい。

　また、上記の実施形態では、単一の装置（音場装置１０）が音場測定機能及び音場補正機能を有しているが、本発明はこれに限らない。別の実施形態では、音響システムが複数の装置で構成されており、音場測定機能と音場補正機能とが該システムを構成する別々の装置に備えられていてもよい。一例として、スマートフォン等の情報処理端末が音場測定を行い、その測定結果に基づいて車載器等の装置が音場補正を行う構成が考えられる。

Claims

　所定の音声信号に基づいてパラメトリックイコライザによる音場補正の目標となる目標振幅特性及び該目標振幅特性と同等若しくはそれよりもなだらかな振幅特性を持つ複数の副目標振幅特性を算出する目標算出手段と、
　算出された目標振幅特性及び複数の副目標振幅特性に基づいて音場を補正するための補正振幅特性を算出する補正振幅特性算出手段と、
　算出された補正振幅特性に基づいて前記パラメトリックイコライザを設定する設定手段と、
を備える、
音場補正装置。
　前記目標算出手段は、
　　前記音声信号に基づいてその振幅特性を算出し、算出された振幅特性に対して第一の分解能で平均化処理を施し、平均化処理された振幅特性を基に前記目標振幅特性を算出し、
　　前記音声信号に基づいてその振幅特性を算出し、算出された振幅特性に対して前記第一の分解能と同等若しくはそれよりも粗い分解能であってそれぞれ異なる複数種類の分解能で平均化処理を施し、異なる分解能で平均化処理された各振幅特性を基に前記複数の副目標振幅特性を算出する、
請求項１に記載の音場補正装置。
　前記補正振幅特性算出手段は、
　　所定の条件に基づいて各前記副目標振幅特性を周波数領域において複数のグループに分割し、分割された各グループの優先度を算出し、算出された優先度が最も高いグループを該副目標振幅特性毎に選択し、
　選択された周波数領域のグループに基づいて補正候補の前記パラメトリックイコライザのパラメータを前記副目標振幅特性毎に算出し、算出された副目標振幅特性毎のパラメータに基づいて補正対象の該パラメトリックイコライザのパラメータを取得し、
　前記設定手段は、
　　前記補正対象のパラメトリックイコライザのパラメータに基づいて該パラメトリックイコライザを設定する、
請求項１又は請求項２に記載の音場補正装置。
　前記補正振幅特性算出手段は、
　　前記補正候補のパラメトリックイコライザのパラメータに基づいて補正候補の振幅特性を前記副目標振幅特性毎に算出し、算出された各補正候補の振幅特性のうち、前記目標振幅特性との差分が最も小さいもののパラメータを前記補正対象のパラメトリックイコライザのパラメータとして取得する、
請求項３に記載の音場補正装置。
　前記補正振幅特性算出手段は、
　　前記補正対象のパラメトリックイコライザのパラメータが取得されると、前記副目標振幅特性毎に算出された補正候補の振幅特性の中で前記目標振幅特性との差分が最も小さい補正対象の振幅特性と、該目標振幅特性との差分を新たな目標振幅特性として算出し、且つ該補正対象の振幅特性と前記各副目標振幅特性との差分を新たな複数の副目標振幅特性として算出し、
　　取得された補正対象のパラメトリックイコライザの数が所定数に達していなければ、前記新たな目標振幅特性及び前記新たな複数の副目標振幅特性に基づいて、該補正対象のパラメトリックイコライザのパラメータを更に取得する、
請求項３又は請求項４に記載の音場補正装置。
　前記パラメータは、
　　中心周波数、ゲイン及び周波数帯域幅を含む、
請求項３から請求項５の何れか一項に記載の音場補正装置。
　所定の音声信号に基づいてパラメトリックイコライザによる音場補正の目標となる目標振幅特性及び該目標振幅特性と同等若しくはそれよりもなだらかな振幅特性を持つ複数の副目標振幅特性を算出する目標算出ステップと、
　算出された目標振幅特性及び複数の副目標振幅特性に基づいて音場を補正するための補正振幅特性を算出する補正振幅特性算出ステップと、
　算出された補正振幅特性に基づいて前記パラメトリックイコライザを設定する設定ステップと、
を含む、
音場補正方法。
　前記目標算出ステップにて、
　　前記音声信号に基づいてその振幅特性を算出し、算出された振幅特性に対して第一の分解能で平均化処理を施し、平均化処理された振幅特性を基に前記目標振幅特性を算出し、
　　前記音声信号に基づいてその振幅特性を算出し、算出された振幅特性に対して前記第一の分解能と同等若しくはそれよりも粗い分解能であってそれぞれ異なる複数種類の分解能で平均化処理を施し、異なる分解能で平均化処理された各振幅特性を基に前記複数の副目標振幅特性を算出する、
請求項７に記載の音場補正方法。
　前記補正振幅特性算出ステップにて、
　　所定の条件に基づいて各前記副目標振幅特性を周波数領域において複数のグループに分割し、分割された各グループの優先度を算出し、算出された優先度が最も高いグループを該副目標振幅特性毎に選択し、
　選択された周波数領域のグループに基づいて補正候補の前記パラメトリックイコライザのパラメータを前記副目標振幅特性毎に算出し、算出された副目標振幅特性毎のパラメータに基づいて補正対象の該パラメトリックイコライザのパラメータを取得し、
　前記設定ステップにて、
　　前記補正対象のパラメトリックイコライザのパラメータに基づいて該パラメトリックイコライザを設定する、
請求項７又は請求項８に記載の音場補正方法。
　前記補正振幅特性算出ステップにて、
　　前記補正候補のパラメトリックイコライザのパラメータに基づいて補正候補の振幅特性を前記副目標振幅特性毎に算出し、算出された各補正候補の振幅特性のうち、前記目標振幅特性との差分が最も小さいもののパラメータを前記補正対象のパラメトリックイコライザのパラメータとして取得する、
請求項９に記載の音場補正方法。
　前記補正振幅特性算出ステップにて、
　　前記補正対象のパラメトリックイコライザのパラメータが取得されると、前記副目標振幅特性毎に算出された補正候補の振幅特性の中で前記目標振幅特性との差分が最も小さい補正対象の振幅特性と、該目標振幅特性との差分を新たな目標振幅特性として算出し、且つ該補正対象の振幅特性と前記各副目標振幅特性との差分を新たな複数の副目標振幅特性として算出し、
　　取得された補正対象のパラメトリックイコライザの数が所定数に達していなければ、前記新たな目標振幅特性及び前記新たな複数の副目標振幅特性に基づいて、該補正対象のパラメトリックイコライザのパラメータを更に取得する、
請求項９又は請求項１０に記載の音場補正方法。
　前記パラメータは、
　　中心周波数、ゲイン及び周波数帯域幅を含む、
請求項９から請求項１１の何れか一項に記載の音場補正方法。
　請求項７から請求項１２の何れか一項に記載の音場補正方法をコンピュータに実行させるための音場補正プログラム。