(57)【要約】【課題】周期的な波形要素からなる音の音質を評価する場合、その波形要素の性質を考慮して、実際の変動感に近づけられるようにする。【解決手段】周期的な波形要素を含む音波の変動強度を推定する場合、受け付けた音波の振幅包絡を周波数解析して最もスペクトルの大きい基本周波数を算出する。そして、その基本周波数の逆数から前記波形要素の周期を求め、当該波形要素の周期を用いて波形要素の振幅が最大となるオンセット時刻、振幅が最小となるオフセット時刻、スペクトル重心、音響パワー、正規化された音響パワーのうち少なくとも一つを含む三次変動の特徴量を求める。また、各波形要素の振幅逸脱量や時間逸脱量を考慮した二次変動や、従来の変動成分である一次変動を考慮して最終的な変動強度を推定する。【選択図】図７
2012159443.tif 000002
【特許請求の範囲】
【請求項１】
周期的な波形要素を含む音波の変動強度を推定する音質評価システムにおいて、
前記音波の振幅包絡を周波数解析して最もスペクトルの大きい基本周波数を算出する周波数解析手段と、
当該算出された基本周波数の逆数を算出して前記波形要素の周期を求め、当該波形要素の周期を用いて波形要素の振幅が最大となる時刻であるオンセット時刻、振幅が最小となる時刻であるオフセット時刻、スペクトル重心、音響パワー、正規化された音響パワーのうち少なくとも一つを含む要素特徴量を求める要素特徴量算出手段と、
当該算出された要素特徴量を用いて変動強度を推定する変動強度推定手段と、
を備えたことを特徴とする音質評価システム。
【請求項２】
請求項１における音質評価システムにおいて、
前記算出された波形要素の１／４周期ごとに振幅包絡の傾きを求め、当該傾きの符号が変化する時刻間の長さである時間領域を求めて当該時間領域内で振幅が最大となる時刻を探索し、当該時刻における振幅の逸脱量である振幅逸脱量を算出する振幅逸脱量算出手段を備え、
前記変動強度推定手段が、当該振幅逸脱量および前記要素特徴量を用いて変動強度を推定するものである音質評価システム。
【請求項３】
請求項２における音質評価システムにおいて、
前記探索された時刻ごとの間隔を抽出し、各時間間隔の逸脱量である時間逸脱量を算出する時間逸脱量算出手段を備え、
前記変動強度推定手段が、当該時間逸脱量および前記要素特徴量を用いて変動強度を推定するものである音質評価システム。
【請求項４】
周期的な波形要素を含む音波の変動強度を推定する音質評価方法において、
前記音波の振幅包絡を周波数解析して最もスペクトルの大きい基本周波数を算出するステップと、
当該算出された基本周波数の逆数から前記波形要素の周期を求め、当該波形要素の周期を用いて波形要素の振幅が最大となる時刻であるオンセット時刻、振幅が最小となる時刻であるオフセット時刻、スペクトル重心、音響パワー、正規化された音響パワーのうち少なくとも一つを含む要素特徴量を求めるステップと、
当該算出された要素特徴量を用いて変動強度を推定するステップと、
を備えたことを特徴とする音質評価方法。
【請求項５】
請求項４における音質評価方法において、
前記算出された波形要素の１／４周期ごとに振幅包絡の傾きを求め、当該傾きの符号が変化する時刻間の長である時間領域を求めて当該時間領域内で振幅が最大となる時刻を探索し、当該時刻における振幅の逸脱量である振幅逸脱量を算出するステップを備え、
当該振幅逸脱量および前記要素特徴量を用いて変動強度を推定するようにした音質評価方法。
【請求項６】
請求項５における音質評価方法において、
前記探索された時刻ごとの間隔を抽出し、各時間間隔の逸脱量である時間逸脱量を算出するステップを備え、
当該時間逸脱量および前記要素特徴量を用いて変動強度を推定するようにした音質評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば、エンジンの排気音などのように断続的に発せられる音の性質を評価する音質評価システムおよびその評価方法に関するものである。
【背景技術】
【０００２】
一般に、人間が音を聴く場合、その音の大きさ（ラウドネス）や音の鋭さ（シャープネス）、音の変動強度（Fluctuation Strength）、音の粗さ（ラフネス）などに基づいて刺激を受けるとされている（非特許文献１）。このため、これらの音を心理音響評価量として数値化することによって（非特許文献２）、例えば、家電製品の騒音を最小限度に抑えたり、自動車の騒音を最小限度に抑えるようにすることなどが考えられている。
【０００３】
ところで、近年では騒音の低減とは逆に、静粛な電気自動車の接近音を歩行者に知らせるように警告音を報知することも考えられている。このような電気自動車に使われる警告音としては、例えば、50dBの排気音を使うことなどが提案されているが（非特許文献３）、現実的には、どのような音を発すれば最も効果的であるのかは判っていない。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】「小野測器技術レポート  音質評価とは」（http://www.onosokki.co.jp/HP-WK/c_support/newreport/soundquality/index.htm）H23年1月19日確認
【非特許文献２】「人に優しい家電製品や自動車などの音響デザインに関する研究」（章忠  豊橋技術科学大学）（www.sangetsu.co.jp/hibizaidan/pdf/hibi_17/syou.pdf）H23年1月19日確認
【非特許文献３】「電気自動車の接近通知と歩行者の安全確保に関する考察と提案」（川崎医療福祉学会誌）Vol.19 No.2 2010 269-275
【発明の概要】
【発明が解決しようとする課題】
【０００５】
このような排気音を発する場合、一般的には、その音の大きさ（ラウドネス）や音の鋭さ（シャープネス）などを大きくすればよいが、このように音の大きさや鋭さを大きくすると、逆に周囲に対する騒音が大きくなってしまう。
【０００６】
一方、エンジンの排気音などのような振幅変動音を発するとともに、その変調周波数を4Hz〜8Hz程度にすると、最も歩行者に刺激を与えることができる。しかしながら、従来の変動強度の推定方法を用いて最適な振幅変動音を評価する場合、発せられる衝撃音の波形形状が異なっている場合であっても、これを評価することができなかった。すなわち、従来における変動強度の推定方法では、その音のラウドネスの時間履歴にどれだけ変動成分が含まれているかを調べるだけであるため、断続的に発せられるそれぞれの音のエンベロープの立ち上がりが速い場合と遅い場合を比較した場合、実際には、立ち上がりの速い音の方が変動感が強く感じられるにもかかわらず、これらの違いをうまく評価することができなかった。
【０００７】
また、４サイクルエンジン音の場合、「吸気、圧縮、爆発、排気」の４つのサイクルを繰り返すため、排気音もこれらの４つのサイクルに合わせた波形要素を含むことになるが、実際には、爆発音の最大振幅がずれたり、あるいは、発音時刻がずれたりする場合がある。これらの違いは、人間に変動感の違いとして感じられるが、従来の方法では、これを大きく評価に組み込むことができなかった。
【０００８】
そこで、本発明は、上記課題に着目して、周期的な波形要素からなる音の音質を評価する場合、その波形要素の性質を考慮して、実際の変動感に近づけられるようにした音質評価システムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
すなわち、本発明は上記課題を解決するために、周期的な波形要素を含む音波の変動強度を推定する音質評価システムにおいて、前記音波の振幅包絡を周波数解析して最もスペクトルの大きい基本周波数を算出する周波数解析手段と、当該算出された基本周波数の逆数を算出して前記波形要素の周期を求め、当該波形要素の周期を用いて波形要素の振幅が最大となる時刻であるオンセット時刻、振幅が最小となる時刻であるオフセット時刻、スペクトル重心、音響パワー、正規化された音響パワーのうち少なくとも一つを含む要素特徴量を求める要素特徴量算出手段と、当該算出された要素特徴量を用いて変動強度を推定する変動強度推定手段とを備えるようにしたものである。
【００１０】
このようにすれば、周期的な波形要素の特徴を考慮することによって、振幅変動の値に波形要素の値を組み込むことができ、より人間が聴いたときの変動感に近づけることができるようになる。
【００１１】
また、このような発明において、波形要素の１／４周期ごとに振幅包絡の傾きを求め、当該傾きの符号が変化した時刻間の長さである時間領域を求めて当該時間領域内で振幅が最大となる時刻を探索し、当該時刻における振幅の逸脱量である振幅逸脱量を算出する振幅逸脱量算出手段を備え、当該振幅逸脱量および前記要素特徴量を用いて変動強度を推定するようにする。
【００１２】
このようにすれば、エンジンの不規則な爆発音の大きさのずれなども考慮して変動強度を求めることができるようになる。
【００１３】
さらには、前記探索された時刻ごとの間隔を抽出し、各時間間隔の逸脱量である時間逸脱量を算出する時間逸脱量算出手段を備え、当該時間逸脱量および前記要素特徴量を用いて変動強度を推定するようにする。
【００１４】
このようにすれば、エンジンの不規則な爆発音の時間のずれなども考慮して変動強度を求めることができるようになる。
【発明の効果】
【００１５】
本発明によれば、周期的な波形要素を含む音波の変動強度を推定する音質評価システムにおいて、前記音波の振幅包絡を周波数解析して最もスペクトルの大きい基本周波数を算出する周波数解析手段と、当該算出された基本周波数の逆数から前記波形要素の周期を求め、当該算出された要素周期を用いて波形要素の振幅が最大となるオンセット時刻、振幅が最小となるオフセット時刻、スペクトル重心、音響パワー、正規化された音響パワーのうち少なくとも一つを含む要素特徴量を求める要素特徴量算出手段と、当該算出された要素特徴量を用いて変動強度を推定する変動強度推定手段とを備えるようにしたので、振幅変動の値に波形要素の値を組み込むことができ、より人間が聴いたとき変動感に近づけることができるようになる。
【図面の簡単な説明】
【００１６】
【図１】本発明の一実施の形態における音質評価システムの機能ブロック図
【図２】同形態における入力された音響波形から振幅包絡を求めて周波数解析をする処理を示す図
【図３】同形態における波形要素を示す図
【図４】同形態における一次変動の特徴量を算出する場合のフローチャート
【図５】同形態における二次変動の特徴量を算出する場合のフローチャート
【図６】同形態における三次変動の特徴量を算出する場合のフローチャート
【図７】同形態における変動強度を推定する場合のフローチャート
【図８】同形態における一次変動と二次変動の特徴量をまとめた図
【図９】同形態における三次変動の特徴量をまとめた図
【発明を実施するための最良の形態】
【００１７】
以下、本発明の一実施の形態について図面を参照しながら説明する。この実施の形態における音質評価システム１００は、パーソナルコンピューターなどによって構成されるものであって、マイクを介して入力された音の振幅変動を算出して出力するようにしたものである。より具体的には、音波の入力を受け付ける入力受付手段１と、その受け付けた音の振幅包絡を算出する振幅包絡算出手段２と、その算出された振幅包絡を周波数解析して基本周波数を算出する周波数解析手段３とを備えている。そして、特徴的には、その基本周波数の逆数から波形要素の周期を求め、その周期ごとに波形要素を推定して、その波形要素のオンセット時刻、オフセット時刻、スペクトル重心、音響パワー、正規化された音響パワーなどを求めて変動強度を求めるようにしたものである。以下、本実施の形態における音質評価システム１００について詳細に説明する。
【００１８】
まず、入力受付手段１は、評価の対象となる音源からの音波の入力を受け付ける。ここで入力される音波としては、バイクの排気音を想定する。この排気音は、エンジン音を構成する音の一つであり、シリンダー内で燃料と空気を混合させて爆発させる衝撃波の一種である。通常、エンジンは１秒間に何回も爆発や排気などを繰り返すため、その排気音も断続的なものになる。一般に、インジェクターを搭載したエンジンでは、燃料と空気をコンピューターで電子的に制御して混合させるのに対して、キャブレターを搭載したエンジンは、吸気口内の気圧が下がることによって燃料と空気を混合させるようになっており、このようなエンジンでは、その排気音の発音時刻の間隔や爆発音の大きさが不均一になって、図２に示すような複数の波形要素を含む波形形状となる。
【００１９】
図２（ａ）のうち、縦方向の太い破線で示された符号１１は波形要素の発音時刻、横線で示される符号１２は発音時刻における振幅、符号１３は発音時刻の間隔Inter-Onset-Interval（以下、IOIと称する）である。図２（ａ）により、単一の爆発音である波形要素の最大振幅および発音時刻は不均一であり、等間隔なIOIからの逸脱や爆発音の最大振幅の逸脱が含まれていることがわかる。それらの逸脱量は、知覚印象に大きく影響を与えるため、変動感を音響波形から推定する場合は、それらを考慮する必要がある。また、排気音は、波形要素である爆発音の繰り返しで構成されているため、その繰り返し回数が変動感に影響を及ぼすと考えられる。また、排気音は車種によっても大きく変わるため、図３に示すような排気音の波形要素から、スペクトル重心や音響パワー、正規化された音響パワー、立ち上がり時間および立ち下がり時間の長さ、立ち上がり時間および立ち下がり時間における傾きやその高さなどを求め、変動強度を推定するものとする。
【００２０】
振幅包絡算出手段２は、この入力された波形を平滑化して振幅包絡を算出する。この振幅包絡を算出する場合は、まず、入力された音響波形を二乗してパワー関数を算出し、その後、所定の窓幅の移動平均を数回繰り返して算出する。すると、図２（ｂ）に示すような平滑化された振幅包絡が得られる。
【００２１】
次に、周波数解析手段３は、このようにして得られた振幅包絡を周波数解析する（図２（ｃ））。この周波数解析を行う場合は、全波整流によって得られた振幅包絡に対してFast Fourier Transform(以下、ＦＦＴと称する)を行い、得られた周波数から最もスペクトルの大きな周波数である基本周波数を抽出する。この得られた基本周波数は、入力された音響波形に含まれる波形要素の周波数とほぼ等しい周波数となる。
【００２２】
一次変動算出手段４は、このようにして得られた値を用いて変動成分を抽出する。この変動成分を抽出する場合は、周波数解析手段３によって解析された周波数に対して、1Hz〜20Hzを通す重み付き帯域通過フィルタを通し、そのフィルタに通した波形に対し逆フーリエ変換（Inverse Fast Fourier Transform）を行う。そして、逆フーリエ変換した後の波形の実効値（RMS値）を求めて、そのRMS値を一次変動パラメタの値とする。
【００２３】
二次変動算出手段５は、各波形要素の最大振幅やIOIがどれだけ逸脱しているかを算出する。この二次変動の値を算出する場合、各波形要素のオンセット時刻およびオフセット時刻を推定する。このオンセット時刻やオフセット時刻の推定においては、同様に、音響波形を二乗して得られたパワー関数に所定の窓幅の移動平均を数回行い、得られた波形の基本周波数を求めた後にそれの逆数を算出して波形要素の周期を求める。次に、その波形要素の１／４周期ごとにその傾きを最小二乗法を用いて推定し、その傾きの符号が正から負に変化した場合の長さ、すなわち１／２周期の長さを持つ時間領域を求め、その領域内で振幅が最大となる時刻を「オンセット時刻」とし、また、傾きの符号が負から正に変化した場合の長さ、すなわち１／２周期の長さを持つ時間領域を求め、その領域内で振幅が最小となる時刻を「オフセット時刻」とする。次に、得られたオンセット時刻に対して、逸脱がない所定の状態のオンセット時刻からのずれを抽出し、時間逸脱量ODk(k=1〜N)とした。ここでkは波形要素のID、Nは波形要素の総数を示すものとする。なお、k=1の波形要素におけるオンセット時刻を基準とし、各波形要素におけるODkを算出している。また、「オンセット時刻における音圧レベル」を音圧レベルの最大値からのずれとして振幅逸脱量ADkとした。そして、得られた各逸脱量xkに対して移動平均を行い、傾向曲線を求め、次の＜数１＞〜＜数５＞を用いて表１における特徴量pi,j(i=1-2, j=0-4)を求める。
【００２４】
【数１】
2012159443.tif 000003
【００２５】
【数２】
2012159443.tif 000004
【００２６】
【数３】
2012159443.tif 000005
【００２７】
【数４】
2012159443.tif 000006
【００２８】
【数５】
2012159443.tif 000007
【００２９】
このように算出される一次変動と二次変動の特徴量をまとめると図８のようになり、排気音の変動成分や、時間逸脱量や振幅逸脱量における標準偏差や傾向からの差分量の和、傾向に対する全体の変動幅、傾向の階差の輪、逸脱がない状態における値からの誤差の和などを求めることができる。
【００３０】
三次変動算出手段６は、音響波形から波形要素を切り出した後に、各波形要素のスペクトル重心、音響パワー、正規化音響パワー、立ち上がり時間および立ち下がり時間の長さ、立ち上がり時間および立ち下がり時間における傾きやその高さを算出する。
【００３１】
この三次変動の値を算出する場合、まず、二次変動の算出の場合と同様に、オンセット時刻およびオフセット時刻を抽出し、これらの時刻を用いて波形要素を切り出す。そして、その切り出された波形要素から、スペクトル重心、音響パワー、正規化音響パワー、立ち上がり時間および立ち下がり時間の長さ、立ち上がり時間および立ち下がり時間における傾きやその高さの９種類の特徴量を算出する。
【００３２】
スペクトル重心を算出する場合は＜数６＞、音響パワーの算出をする場合は＜数７＞、正規化音響パワーを算出する場合は＜数８＞、立ち上がり時間における高さの算出を行う場合は＜数９＞、立ち下がり時間における高さを算出する場合は＜数１０＞を用いる。また、各傾きは最小二乗法を用いて得られた一次回帰式の傾きであり、その傾きを用いて各高さを算出する。そして、１つの波形要素に対して算出した各特徴におけるそれらの平均、変動幅、標準偏差および絶対値の総和を計算し、図９に示すように、９種類の特徴×４＝３６個の特徴量を抽出する。
【００３３】
【数６】
2012159443.tif 000008
【００３４】
【数７】
2012159443.tif 000009
【００３５】
【数８】
2012159443.tif 000010
【００３６】
【数９】
2012159443.tif 000011
【００３７】
【数１０】
2012159443.tif 000012
【００３８】
変動強度推定手段７は、このように算出された一次変動の特徴量から三次変動の特徴量までを用いて最終的な変動強度を推定する。この変動強度を推定する場合、抽出した一次変動の特徴量から三次変動までの特徴量について主成分分析を行い、それらの特徴量を無相関となる特徴量に修正する。次に、実験によって得られた主観評価値と修正後の特徴量を用いて重回帰分析を行うことで変動感を推定する重回帰式を求めるようにする。
【００３９】
次に、このように構成された音質評価システム１００を用いて一次変動の特徴量から三次変動までの特徴量を抽出し、入力された音の変動強度を求める方法について説明する。
【００４０】
まず、一次変動の特徴量を抽出する場合は、図４に示すように、入力された音響波形に対して（ステップＳ１）全波整流を行って振幅包絡を算出し（ステップＳ２）、周波数解析であるＦＦＴを行う（ステップＳ３）。次に、1Hz〜20Hzを通す重み付き帯域通過フィルタを通し（ステップＳ４）、そのフィルタに通した波形に対して逆フーリエ変換(IFFT)を行い（ステップＳ５）、IFFT後の波形の実効値（RMS値）である一次変動の特徴量を求める（ステップＳ６）。
【００４１】
次に、二次変動の特徴量を求める場合は、図５に示すように、入力された音響波形（ステップＴ１）を二乗して得られたパワー関数に（ステップＴ２）、例えば、窓幅20msの移動平均を３回行って振幅包絡を算出し（ステップＴ３）、その得られた振幅包絡の波形を周波数解析して基本周波数を算出する（ステップＴ４）。そして、その基本周波数の逆数から波形要素の周期を求め（ステップＴ５）。次に、排気音の１／４周期ごとにその傾きを最小二乗法を用いて推定し、その傾きの符号が正から負に変化した場合の長さ、すなわち１／２周期の長さを持つ時間領域を求め、その領域内で振幅が最大となる時刻を「オンセット時刻」として抽出し、また、傾きの符号が負から正に変化した場合の長さ、すなわち１／２周期の長さを持つ時間領域を求め、その領域内で振幅が最小となる時刻を「オフセット時刻」として抽出する（ステップＴ６）。
【００４２】
次に、得られたオンセット時刻に対して逸脱がない状態のオンセット時刻からのずれを抽出し、時間逸脱量を算出する（ステップＴ７）。また、オンセット時刻における音圧レベルを、音響波形における音圧レベルの最大値からのずれとして音響波形から抽出し、振幅逸脱量を抽出する（ステップＴ８）。そして、このように得られた逸脱量に対して移動平均を行って傾向曲線を求め、上記＜数１＞〜＜数５＞を用いて二次変動の特徴量を求める（ステップＴ９）。
【００４３】
三次変動の特徴量を求める場合は、図６に示すように、二次変動の特徴量の算出の場合と同様に、各波形要素のオンセット時刻およびオフセット時刻を抽出し（ステップＵ１）、その各時刻を用いて音響波形から波形要素を切り出す（ステップＵ２）。そして、その切り出された波形要素におけるスペクトル重心、音響パワー、正規化された音響パワー、立ち上がり時間の長さや立ち下がり時間の長さ、立ち上がり時間および立ち下がり時間における傾きやその高さの９種類の特徴量を＜数６＞〜＜数１０＞を用いて算出する（ステップＵ３）。そして、このように算出された各特徴量における平均、変動幅、標準偏差および絶対値の総和を算出し、三次変動の特徴量を抽出する（ステップＵ４）。
【００４４】
次に、このように一次変動から三次変動の特徴量を算出した後、図７に示すように、各特徴量を主成分分析（ステップＳＴＵ１）を行うことによって無相関となる特徴量に修正し、重回帰式による推定を行った後（ステップＳＴＵ２）、変動強度の推定を行う（ステップＳＴＵ３）。
【００４５】
このように上記実施の形態によれば、周期的な波形要素を含む音波の変動強度を推定する音質評価システム１００において、当該音波の振幅包絡を周波数解析して最もスペクトルの大きい基本周波数を算出する周波数解析手段３と、当該算出された基本周波数の逆数から前記波形要素の周期を求め、当該要素周期を用いて波形要素のオンセット時刻、オフセット時刻、スペクトル重心、音響パワー、正規化された音響パワーのうち少なくとも一つを含む特徴量を求める三次変動算出手段６と、当該算出された要素特徴量を用いて変動強度を推定する変動強度推定手段７とを備えるようにしたので、振幅変動の値に波形要素の値を組み込むことができ、より人間が聴いた際の変動感に近づけることができるようになる。
【００４６】
なお、上記実施の形態では、振幅変動を有する音波の変動強度を推定する場合について説明したが、例えば、救急自動車のサイレンなどのように周波数変動を有する音波の変動強度も推定する場合にも適用することができる。このような場合としては、周波数変動音を周波数解析して基本周波数を算出し、その基本周波数を有する時間関数を上記実施の形態における振幅包絡とみなして同様の処理をすることができる。この場合、スペクトル重心は求めずに、音響パワー、正規化音響パワー、立ち上がり時間および立ち下がり時間の長さ、立ち上がり時間および立ち下がり時間における傾きやその高さの合計８個のパラメータを計算する。なお、音響パワーは基本周波数の時間関数から１周期の波形の面積を用いて求めるものとする。
【実施例１】
【００４７】
一次変動と二次変動を含む排気音に対する変動感の主観評価実験を行うためには、様々な二次変動を持つ排気音が必要となるが、収録した現実の排気音が持つ二次変動は網羅的でないという問題がある。そこで、収録した排気音から時間および振幅の逸脱を抽出した後に、得られた各逸脱を強調または抑制することにより、様々な二次変動を持つ排気音を合成する。排気音は、爆発音の繰り返しによって構成されるので、記録した排気音から爆発音を抽出し、それを時間軸上でつなぎ合わせることで排気音を合成する。しかし、そのような方法で合成した排気音はオートバイの音色を持っているため、変動感のみだけでなく、別の印象や感覚を感じてしまう可能性がある。そこで、オートバイの音色をもたない排気音を合成するために、記録した排気音から抽出した爆発音をそのまま用いるのではなく、その音の振幅包絡をもつピンクノイズを用いるようにした。
【００４８】
合成した音を用いて、防音室にて正常な聴力を持つ健聴者５名に対して主観評価実験を行った。用いた刺激音は上述の排気音とする。なお、用いた二次変動の傾向は、キャブレターを搭載したエンジン駆動車３台、インジェクターを搭載したエンジン駆動車３台の計６台に含まれるものであり、逸脱量を強調または抑制する割合は、0%、150%、300%の３通りとし、合成した音は、時間および振幅逸脱量において強調する割合を組み合わせた９通りの方法によって強調された音である。各オートバイの排気音は、量子化ビット数16bit、サンプリング周波数48kHzで録音された排気音である。主観評価実験では、シェッフェの一対比較法の浦変法を用いた。聴取者には、２つの音を提示し、最初に聞いた音に比べて後に聞いた音がどの程度変動しているかを７段階評価で回答させた。
【００４９】
抽出した特徴量と、実験によって得られた主観評価値に対して重回帰分析を行い、ＦＳ推定で用いる重回帰式を求める。重回帰分析では、目的変数を主観評価値、説明変数を抽出したパラメタとする。重回帰分析において、複数の説明変数間に強い相関があると、偏回帰係数の推定量が不安定になるため、求められた予測精度の信頼性が低くなる。そこで、抽出した特徴に対して主成分分析を行い、それらの特徴を無相関な変数に修正し、それらを用いて重回帰分析を行う。なお、新たに用いる説明変数は累積寄与率がほぼ100%までの第５主成分までとする。回帰分析によって得られた回帰方程式の予測精度は決定係数（Ｒ2）で表され、決定係数が0.8以上でかなり良い精度、決定係数が0.5以上且つ0.8未満で良い精度であると言える。重回帰分析では、説明変数を増やせば増やすほど決定係数が高くなると言われており、そのような決定係数の欠点を補うために定義されている自由度調整済み決定係数（adjusted R2）を算出する。
【００５０】
１次変動の特徴のみで変動感を推定した場合のR2と、抽出した全ての特徴を考慮に入れた場合のadjustR2を比較した。その結果、抽出したすべての特徴を考慮に入れた場合の方が、一次変動の特徴のみで変動感を推定した場合よりも予測精度が高いことが確認された。
【符号の説明】
【００５１】
１００・・・音質評価システム
１・・・入力受付手段
２・・・振幅包絡算出手段
３・・・周波数解析手段
４・・・一次変動算出手段
５・・・二次変動算出手段
６・・・三次変動算出手段
７・・・変動強度推定手段
【図１】
2012159443.tif 000013
【図２】
2012159443.tif 000014
【図３】
2012159443.tif 000015
【図４】
2012159443.tif 000016
【図５】
2012159443.tif 000017
【図６】
2012159443.tif 000018
【図７】
2012159443.tif 000019
【図８】
2012159443.tif 000020
【図９】
2012159443.tif 000021