"音響信号のうち残響成分または残響成分以外の初期音成分の比率に応じた比率評価値を算定する解析処理手段と、前記音響信号の前記残響成分および前記初期音成分の少なくとも一方の比率または強度を示す残響評価画像を前記比率評価値に応じて表示装置に表示させる表示制御手段とを具備する残響解析装置。前記解析処理手段は、前記音響信号の単位期間毎に複数の周波数の各々について前記比率評価値を算定し、前記表示制御手段は、時間軸と周波数軸とが設定された座標領域を含む前記残響評価画像を前記表示装置に表示させ、前記座標領域内の複数の座標点の各々の表示態様を、時間軸上で当該座標点に対応する単位期間と周波数軸上で当該座標点に対応する周波数とについて算定された前記比率評価値に応じて制御する請求項1の残響解析装置。前記解析処理手段は、複数の周波数の各々について前記比率評価値を算定し、前記表示制御手段は、前記音響信号のスペクトルと、当該音響信号に前記比率評価値を乗算した前記残響成分または前記初期音成分のスペクトルとを、相異なる表示態様で相互に重ねて配置した前記残響評価画像を前記表示装置に表示させる請求項1の残響解析装置。前記解析処理手段は、複数の周波数の各々について前記比率評価値を算定し、前記表示制御手段は、前記音響信号に前記各周波数の比率評価値を乗算した前記残響成分または前記初期音成分の強度を前記複数の周波数について平均した平均強度を、前記残響評価画像として前記表示装置に表示させる請求項1の残響解析装置。前記解析処理手段は、複数の周波数の各々について前記比率評価値を算定し、前記表示制御手段は、周波数軸と定位方向を示す定位軸とが設定された座標領域を含む前記残響評価画像を前記表示装置に表示させ、前記座標領域のうち前記音響信号の各周波数成分の周波数と当該周波数成分の定位方向とに対応する位置に、当該周波数の比率評価値に応じた表示態様の音像点を配置する請求項1の残響解析装置。ところで、例えば複数の音響を適宜に合成してひとつの楽曲を作成する場合や、複数の楽曲を選択してひとつのアルバムを作成する場合には、サンプル間または楽曲間で音響の残響感残響の度合を統一させることが望ましい。再生音を聴取して残響感を正確に把握できれば、音響信号の残響成分を適宜に調整して残響感を統一させることも可能であるが、再生音の聴取で残響感を定量的に把握することは実際には困難である。以上の事情を考慮して、本発明は、音響信号に含まれる残響成分の多寡を容易に把握できるようにすることを目的とする。以上の課題を解決するために本発明が採用する手段を説明する。なお、本発明の理解を容易にするために、以下の説明では、本発明の各要素と後述の各実施形態の要素との対応を括弧書で付記するが、本発明の範囲を実施形態の例示に限定する趣旨ではない。本発明の残響解析装置は、音響信号のうち残響成分または残響成分以外の初期音成分の比率に応じた比率評価値例えば初期音比率Gdや残響比率Grを算定する解析処理手段例えば解析処理部34と、音響信号の残響成分および初期音成分の少なくとも一方の比率または強度を示す残響評価画像例えば残響評価画像72,74,76,78を比率評価値に応じて表示装置に表示させる表示制御手段例えば表示制御部60とを具備する。以上の構成では、音響信号の残響成分および初期音成分の少なくとも一方の比率または強度を示す残響評価画像が比率評価値に応じて表示されるから、音響信号に含まれる残響成分の多寡を利用者が容易に把握できるという利点がある。本発明の第1態様において、解析処理手段は、音響信号の単位期間毎に複数の周波数の各々について比率評価値を算定し、表示制御手段は、時間軸と周波数軸とが設定された座標領域例えば第2座標領域722を含む残響評価画像例えば残響評価画像72を表示装置に表示させ、座標領域内の複数の座標点の各々の表示態様を、時間軸上で当該座標点に対応する単位期間と周波数軸上で当該座標点に対応する周波数とについて算定された比率評価値に応じて制御する。第1態様では、時間軸と周波数軸とが設定された座標領域内に残響成分または初期音成分の比率が表示されるから、周波数間で残響成分や初期音成分の多寡を対比したり、特定の周波数における初期音成分や残響成分の時間変化を確認したりすることが可能である。なお、第1態様の具体例は、例えば第1実施形態として後述される。なお、本発明における表示態様とは、観察者が視覚的に識別できる画像の性状を意味する。例えば、画像の階調,色彩,模様,形状,サイズまたはこれらの組合せが、表示態様の概念に包含される。本発明の第2態様において、解析処理手段は、複数の周波数の各々について比率評価値を算定し、表示制御手段は、音響信号のスペクトルと、当該音響信号に比率評価値を乗算した残響成分または初期音成分のスペクトル例えば振幅スペクトルXdとを、相異なる表示態様で相互に重ねて配置した残響評価画像例えば残響評価画像74を表示装置に表示させる。第1態様では、音響信号のスペクトルと残響成分または初期音成分のスペクトルとが相異なる表示態様で相互に重ねて表示されるから、初期音成分または残響成分の周波数毎の比率に加えて初期音成分や音響成分の強度振幅も視覚的に把握できるという利点がある。なお、第2態様の具体例は、例えば第2実施形態として後述される。本発明の第3態様において、解析処理手段は、複数の周波数の各々について比率評価値を算定し、表示制御手段は、音響信号に各周波数の比率評価値を乗算した残響成分または初期音成分の強度を複数の周波数について平均した平均強度例えば平均強度Adや平均強度Arを、残響評価画像として表示装置に表示させる。例えば、表示制御手段は、所定方向の長さが音響信号の強度に応じて設定された強度図形例えば強度図形760を含む残響評価画像例えば残響評価画像76を表示装置に表示させ、所定方向の長さが初期音成分の平均強度に応じて設定された第1部分例えば第1部分761と、所定方向の長さが残響成分の平均強度に応じて設定された第2部分例えば第2部分762とに、強度図形を所定方向に沿って区分する。以上の態様では、初期音成分または残響成分の比率に加えて初期音成分や音響成分の強度振幅も視覚的に把握できるという利点がある。なお、第3態様の具体例は、例えば第3実施形態として後述される。本発明の第4態様において、解析処理手段は、複数の周波数の各々について比率評価値を算定し、表示制御手段は、周波数軸と定位方向を示す定位軸とが設定された座標領域例えば座標領域780を含む残響評価画像例えば残響評価画像78を表示装置に表示させ、座標領域のうち音響信号の各周波数成分の周波数と当該周波数成分の定位方向とに対応する位置に、当該周波数の比率評価値に応じた表示態様の音像点例えば音像点782を配置する。以上の態様では、音響信号の各周波数成分の定位方向を示す音像点の表示態様が当該周波数の比率評価値に応じて制御されるから、各周波数成分における初期音成分または残響成分の比率または強度と各周波数成分の定位方向との関係を利用者が視覚的および直観的に把握できるという利点がある。なお、第4態様の具体例は、例えば第4実施形態として後述される。第2態様から第4態様の好適例において、解析処理手段は、音響信号の単位期間毎に残響評価値を算定し、表示制御手段は、残響評価画像を単位期間毎に順次に更新する。以上の態様では、残響評価画像が単位期間毎に更新されるから、音響信号の残響成分または初期音成分の比率や強度の時間変化を利用者が把握できるという利点がある。本発明の好適な態様において、解析処理手段は、音響信号の時間変化に追従する第1指標値例えば第1指標値Q1と、第1指標値と比較して低い追従性で音響信号の時間変化に追従する第2指標値例えば第2指標値Q2とを算定する指標値算定手段例えば指標値算定部50A,50Bと、第1指標値と第2指標値との相違に応じて比率評価値を算定する評価値算定手段例えば評価値算定部60とを含む。以上の態様では、音響信号の時間変化に追従する第1指標値と第2指標値との相違に応じて比率評価値が算定されるから、例えば予測フィルタ係数の確率モデルを利用して残響成分を推定する特許文献1の技術と比較して簡易な処理で音響信号の残響成分および初期音成分を解析できるという利点がある。ただし、本発明における比率評価値の算定残響成分の解析には公知の技術が任意に採用される。本発明の好適な態様において、指標値算定手段は、音響信号の強度の時系列を平滑化して第1指標値を算定する第1平滑手段例えば第1平滑部51と、第1平滑手段による平滑化の時定数を上回る時定数で音響信号の強度の時系列を平滑化することで第2指標値を算定する第2平滑手段例えば第2平滑部52とを含む。他の態様において、指標値算定手段は、第2指標値の時間変化が第1指標値の時間変化を遅延させた関係となるように、音響信号の信号強度の時系列を平滑化した第1指標値および第2指標値を生成する。本発明の好適な態様において、評価値算定手段は、第2指標値に対する第1指標値の比例えば比Rを算定する比算定手段例えば比算定部62と、比が所定値を上回る場合に当該所定値に設定され、比が所定値例えば1を下回る場合に比に設定される初期音比率を比率評価値として算定する第1処理手段例えば第1処理部64と、初期音比率を所定値から減算することで残響比率を比率評価値として算定する第2処理手段例えば第2処理部66とを含む。以上の態様では、第2指標値に対する第1指標値の比の演算と所定値から初期音比率を減算する演算とを含む簡易な演算で残響比率を算定できるという利点がある。以上の各態様に係る音響処理装置は、音響信号の処理に専用されるDSPDigitalSignalProcessorなどのハードウェア電子回路によって実現されるほか、CPUCentralProcessingUnitなどの汎用の演算処理装置とプログラムとの協働によっても実現される。本発明に係るプログラムは、音響信号のうち残響成分または残響成分以外の初期音成分の比率に応じた比率評価値を算定する解析処理と、音響信号の残響成分および初期音成分の少なくとも一方の比率または強度を示す残響評価画像を比率評価値に応じて表示装置に表示させる表示制御処理とをコンピュータに実行させる。以上のプログラムによれば、本発明に係る残響解析装置と同様の作用および効果が実現される。なお、本発明のプログラムは、コンピュータが読取可能な記録媒体に格納された形態で提供されてコンピュータにインストールされるほか、通信網を介した配信の形態で提供されてコンピュータにインストールされる。<第1実施形態>図1は、本発明の第1実施形態に係る音響処理装置100のブロック図である。音響処理装置100には信号供給装置12と放音装置14と表示装置16とが接続される。信号供給装置12は、音響信号xを音響処理装置100に供給する。音響信号xは、音響の波形を示す時間領域の信号であるt:時間。周囲の音響を収音して音響信号xを生成する収音機器や、可搬型または内蔵型の記録媒体から音響信号xを取得して音響処理装置100に供給する再生装置や、通信網から音響信号xを受信して音響処理装置100に供給する通信装置が信号供給装置12として採用され得る。音響信号xが示す音響は、初期音成分ドライ成分と残響成分ウェット成分とを包含する。残響成分は、音源の発音動作の停止後も経時的に減衰しながら継続する響き成分である。具体的には、音響空間の壁面での反射および散乱後に受音点に到来する音響初期反射音,後部残響音や、自然楽器の響板による共鳴音胴鳴り,箱鳴り等が残響成分に該当する。初期音成分は、残響成分以外の音響成分である。具体的には、音源の発音動作に直接的に起因する音響反射や共鳴を殆ど経ていない音響が初期音成分に該当する。収音時に音響空間内で発生した残響成分や収音後に人為的に生成された残響成分が初期音成分に付加される。なお、以下の説明では、初期音成分に関連する要素に添字dd:dryを付加し、残響成分に関連する要素に添字rr:reverberationを付加する場合がある。音響処理装置100は、音響信号xの残響成分を解析するとともに音響信号xから音響信号yを生成する。音響信号yは、音響信号xの残響成分を抑圧または強調した時間領域の信号である。残響成分の抑圧/強調は例えば利用者から指示される。放音装置14は、音響処理装置100が生成した音響信号yに応じた音波を放射する。なお、音響信号yをデジタルからアナログに変換するDA変換器の図示は便宜的に省略されている。図1の表示装置16例えば液晶表示器は、音響処理装置100から指示された画像を表示する。具体的には、音響処理装置100による残響成分の解析結果残響評価画像が表示装置16に表示される。図1に示すように、音響処理装置100は、演算処理装置22と記憶装置24とを具備するコンピュータシステムで実現される。記憶装置24は、演算処理装置22が実行するプログラムPGMや演算処理装置22が使用する各種のデータを記憶する。半導体記録媒体や磁気記録媒体などの公知の記録媒体や複数種の記録媒体の組合せが記憶装置24として任意に採用され得る。音響信号xを記憶装置24に記憶した構成したがって信号供給装置12は省略されるも好適である。演算処理装置22は、記憶装置24に格納されたプログラムPGMを実行することで複数の機能周波数分析部32,解析処理部34,残響調整部36,波形合成部38,表示制御部40を実現する。なお、演算処理装置22の各機能を複数の集積回路に分散した構成や、専用の電子回路DSPが各機能を実現する構成も採用され得る。周波数分析部32は、音響信号xのスペクトル複素スペクトルXを時間軸上の単位期間フレーム毎に順次に生成する。記号kは、周波数軸上の任意の1個の周波数帯域を指定する変数であり、記号mは、時間軸上の任意の1個の単位期間時間軸上の特定の時点を指定する変数である。スペクトルXの生成には、短時間フーリエ変換等の公知の周波数分析が任意に採用され得る。なお、通過帯域が相違する複数の帯域通過フィルタで構成されるフィルタバンクも周波数分析部32として採用され得る。解析処理部34は、音響信号xの残響成分を解析する要素である。具体的には、解析処理部34は、音響信号xのスペクトルXに応じた初期音比率Gdおよび残響比率Grを複数の周波数の各々について単位期間毎に算定する。初期音比率Gdは、音響信号xのうち初期音成分の比率に応じた変数比率評価値である。概略的には、1個の単位期間のスペクトルXにて初期音成分の強度が高い周波数初期音成分が優勢である周波数の初期音比率Gdほど大きい数値に設定されるという傾向がある。すなわち、初期音比率Gdは、音響信号xのスペクトルXにおける初期音成分の優勢度とも換言され得る。他方、残響比率Grは、音響信号xにおける残響成分の比率に応じた変数比率評価値である。概略的には、スペクトルXにて残響成分の強度が高い周波数残響成分が優勢である周波数の残響比率Grほど大きい数値に設定されるという傾向がある。すなわち、残響比率Grは、音響信号xのスペクトルXにおける残響成分の優勢度とも換言され得る。図2は、解析処理部34のブロック図である。図2に示すように、第1実施形態の解析処理部34は、指標値算定部50Aと評価値算定部60とを具備する。指標値算定部50Aは、音響信号xに応じた第1指標値Q1と第2指標値Q2とを順次に算定する。具体的には、指標値算定部50Aは、第1平滑部51と第2平滑部52とを含んで構成される。第1平滑部51は、音響信号xのパワーX2の時系列を平滑化することで各周波数の第1指標値Q1を単位期間毎に順次に算定する。同様に、第2平滑部52は、音響信号xのパワーX2の時系列を平滑化することで各周波数の第2指標値Q2を単位期間毎に順次に算定する。第1指標値Q1は、以下の数式で定義されるように、相前後するM1個M1は2以上の自然数の単位期間で構成される第1期間内のパワーX2の移動平均単純移動平均である。第1期間は、例えば第m番目の単位期間を最後尾とするM1個の単位期間の集合である。他方、第2指標値Q2は、以下の数式で定義されるように、相前後するM2個M2は2以上の自然数の単位期間で構成される第2期間内のパワーX2の移動平均である。第2期間は、例えば第m番目の単位期間を最後尾とするM2個の単位期間の集合である。以上の説明から理解されるように、第1平滑部51および第2平滑部52はFIRfiniteimpulseresponse型のローパスフィルタに相当する。第2指標値Q2の算定に加味される単位期間の個数M2は、第1指標値Q1の算定に加味される単位期間の個数M1を上回るM2>M1。すなわち、第2期間は第1期間よりも長い。例えば、第1期間は100ミリ秒から300ミリ秒程度の時間に設定され、第2期間は300ミリ秒から600ミリ秒程度の時間に設定される。したがって、第2平滑部52による平滑化の時定数τ2は第1平滑部51による平滑化の時定数τ1を上回るτ2>τ1。第1平滑部51および第2平滑部52をローパスフィルタで実現する場合を想定すると、第2平滑部52の遮断周波数が第1平滑部51の遮断周波数を下回ると換言することも可能である。図3の部分は、音響信号xの任意の周波数について算定される第1指標値Q1および第2指標値Q2の時間変化のグラフである。図3の部分のようにパワーX2パワー密度が指数減衰する室内インパルス応答RIRを音響信号xとして音響処理装置100に供給した場合の第1指標値Q1および第2指標値Q2が図3の部分には図示されている。図3の部分から理解されるように、第1指標値Q1および第2指標値Q2は、音響信号xのパワーX2に追従して経時的に変化する。ただし、第2平滑部52による平滑化の時定数τ2は第1平滑部51による平滑化の時定数τ1を上回るから、第2指標値Q2は、第1指標値Q1と比較して低い追従性変化率で音響信号xのパワーX2の時間変化に追従する。具体的には、図3の部分に示すように、室内インパルス応答の開始の時点t0の直後の区間では、第1指標値Q1が第2指標値Q2を上回る変化率で増加する。そして、第1指標値Q1および第2指標値Q2は、時間軸上の相異なる時点でピークに到達し、第1指標値Q1は第2指標値Q2を上回る変化率で減少する。以上のように第1指標値Q1と第2指標値Q2とは相異なる変化率で変化するから、第1指標値Q1と第2指標値Q2との大小は時間軸上の特定の時点txで反転する。すなわち、時点t0から時点txまでの区間SAでは第1指標値Q1が第2指標値Q2を上回り、時点tx以降の区間SBでは第2指標値Q2が第1指標値Q1を上回る。区間SAは、室内インパルス応答の初期音成分直接音が存在する区間に相当し、区間SBは、室内インパルス応答の残響成分後部残響音が存在する区間に相当する。図2の評価値算定部60は、指標値算定部50Aが算定した第1指標値Q1と第2指標値Q2とに応じた初期音比率Gdおよび残響比率Grを各周波数について単位期間毎に順次に算定する。第1実施形態の評価値算定部60は、比算定部62と第1処理部64と第2処理部66とを含んで構成される。比算定部62は、第1指標値Q1と第2指標値Q2との比Rを算定する。具体的には、比算定部62は、以下の数式で表現される通り、第2指標値Q2に対する第1指標値Q1の比Rを単位期間毎に算定する。図2の第1処理部64は、比算定部62が算定した比Rに応じて初期音比率Gdを各周波数について単位期間毎に順次に算定する。具体的には、第1処理部64は、数式で表現される閾値処理を実行する。第1に、比Rが1閾値を上回る場合、R≧1、第1処理部64は、初期音比率Gdを1に設定する。第2に、比Rが所定の下限値Gminを下回る場合R≦Gmin、第1処理部64は、初期音比率Gdを下限値Gminに設定する。下限値Gminは、0以上かつ1未満の数値に設定される。第3に、比Rが1と下限値Gminとの間の数値である場合Gmin<R<1、第1処理部64は、比Rを初期音比率Gdとして設定する。第1指標値Q1および第2指標値Q2が図3の部分のように変化する場合の初期音比率Gdに図示されている。図3の部分から理解されるように、概略的には、第1指標値Q1が第2指標値Q2を上回る場合区間SAの初期音比率Gdが第2指標値Q2を下回る場合区間SBの初期音比率Gdが第2指標値Q2を上回る区間SA内では比Rが1を上回るから、初期音比率Gdが第2指標値Q2を下回る区間SBのうち比Rが下限値Gminを上回る区間SB1では、初期音比率Gdに設定されて経時的に減少する。そして、区間SBのうち比Rが下限値Gminを下回る区間SB2では、初期音比率Gdにおける初期音成分の比率の指標値として初期音比率Gdを利用することが可能である。図2の第2処理部66は、第1処理部64が算定した初期音比率Gdを各周波数について単位期間毎に順次に算定する。初期音比率Gdが減少するように残響比率Grは算定される。具体的には、第2処理部66は、前掲の数式で算定された初期音比率Gdを算定するGr=1-Gd。したがって、残響比率Grは、初期音成分が存在する区間SAではゼロに維持され、残響成分が存在する区間SBでは所定値1-Gminまで経時的に増加する。すなわち、第1指標値Q1が第2指標値Q2を上回る場合区間SAの残響比率Grは、第1指標値Q1が第2指標値Q2を下回る場合区間SBの残響比率Grよりも小さい数値となる。したがって、音響信号xにおける残響成分の比率の指標値として残響比率Grを利用することが可能である。以上に説明したように、第1実施形態では、音響信号xの時間変化に追従する第1指標値Q1および第2指標値Q2の比Rに応じて初期音比率Gdおよび残響比率Grが算定されるから、予測フィルタ係数の確率モデルを利用して残響成分を推定する特許文献1の技術と比較して簡易な処理で音響信号xの残響成分を解析できるという利点がある。図1の残響調整部36は、解析処理部34が算定した初期音比率Gdまたは残響比率Grを音響信号xのスペクトルXに作用させることで音響信号yのスペクトルYを単位期間毎に生成する。残響成分の抑圧が指示された場合、残響調整部36は、初期音比率GdをスペクトルXに乗算することでスペクトルYを生成するY=GdX。初期音成分が優勢な周波数の初期音比率Gdほど大きい数値に設定されるから、音響信号xの残響成分を抑圧したスペクトルYが生成される。他方、残響成分の強調が指示された場合、残響調整部36は、残響比率GrをスペクトルXに乗算することでスペクトルYを生成するY=GrX。残響成分が優勢な周波数の残響比率Grほど大きい数値に設定されるから、音響信号xの残響成分を強調抽出したスペクトルYが生成される。波形合成部38は、残響調整部36が残響調整部36が生成したスペクトルYから時間領域の音響信号yを生成する。具体的には、波形合成部38は、各単位期間のスペクトルYの短時間逆フーリエ変換で生成した時間領域信号を相互に連結することで音響信号yを生成する。図1の表示制御部40は、解析処理部34が音響信号xの残響成分および初期音成分を解析した結果を表示装置16に表示させる。第1実施形態の表示制御部40は、図4の残響評価画像72を表示装置16に表示させる。残響評価画像72は、第1座標領域721と第2座標領域722とを含んで構成される。第1座標領域721と第2座標領域722とは、利用者が両者を視覚的に対比できる位置およびサイズで表示される。第1座標領域721および第2座標領域722の各々は、時間軸Tと周波数軸Fとが設定された座標平面である。第1座標領域721には音響信号xのスペクトログラムスペクトルXの時系列が表示される。他方、第2座標領域722には、各周波数の初期音比率Gdの時系列が表示される。具体的には、時間軸Tのうち第m番目の単位期間の時点と、周波数軸Fのうち第k番目の周波数とに対応する座標点時間軸T上の各時点と周波数軸F上の各周波数とに対応して格子状に配列する点は、第m番目の単位期間と第k番目の周波数とに対応する初期音比率Gdの数値に応じた表示態様色彩や階調に制御される。図4では、初期音比率Gdが大きい初期音成分が残響成分と比較して優勢であるほど座標点を低階調濃い階調で表示する場合が例示されている。すなわち、第2座標領域722のうち低階調の領域に対応する時点および周波数では音響信号xが初期音成分を豊富に包含し、第2座標領域722のうち高階調の領域に対応する時点および周波数では音響信号xが残響成分を豊富に包含することを、利用者は視覚的および直観的に把握することが可能である。また、第1座標領域721と第2座標領域722とを対比することで、第1座標領域721内で音響信号xの強度が高い位置時間および周波数に対応する第2座標領域722内の領域ほど低階調初期音成分が優勢であるになるという傾向を把握することもできる。以上に説明したように、第1実施形態では、初期音比率Gdが残響評価画像72として表示装置16に表示されるから、音響信号x内の初期音成分や残響成分の多寡音響信号xの全体に対する比率を利用者が視覚的および直観的に把握できるという利点がある。第1実施形態では特に、各周波数の初期音比率Gdの時系列が残響評価画像72として表示されるから、音響信号xの複数の周波数成分の間で残響成分や初期音成分の多寡を対比したり、音響信号xの特定の周波数について初期音成分や残響成分の時間変化を確認したりすることが可能である。さらに、音響信号xのスペクトログラムを表現する第1座標領域721と初期音比率Gdを表現する第2座標領域722とが対比的に配置されるから、音響信号xの各周波数の強度の経時変化と初期音成分や残響成分の経時変化との関係を容易に把握できるという利点もある。ただし、第1座標領域721を省略することも可能である。なお、以上の例示では第2座標領域722内の各座標点の表示態様を初期音比率Gdに応じて制御したが、残響比率Grに応じて各座標点の表示態様を制御する構成も採用され得る。また、初期音比率Gdまたは残響比率Grおよび音響信号xの強度の双方に応じて第2座標領域722の各座標点の表示態様を制御する構成も好適である。例えば、各座標点の階調および色彩の一方を初期音比率Gdに応じて設定し、各座標点の階調および色彩の他方を音響信号xの強度に応じて設定する構成が採用される。以上の構成によれば、音響信号xの強度と初期音成分または残響成分の比率との双方を第2座標領域722から把握できるという利点がある。<第2実施形態>本発明の第2実施形態を以下に説明する。なお、以下に例示する各形態において作用や機能が第1実施形態と同等である要素については、第1実施形態の説明で参照した符号を流用して各々の詳細な説明を適宜に省略する。第2実施形態の表示制御部40は、音響信号xの特定の単位期間以下対象単位期間というについて図5の残響評価画像74を表示装置16に表示させる。例えば、音響信号xのうち利用者が入力装置図示略に対する操作で指定した時点の単位期間や、音響信号xを順次に処理して音響信号yを再生する場合の現在の再生点に対応する単位期間が対象単位期間として選択される。対象単位期間の変更毎に例えば音響信号yの再生進行に連動して、残響評価画像74は、変更後の対象単位期間に対応する内容に更新される。図5に示すように、第2実施形態の残響評価画像74は、音響信号xの振幅スペクトルXを含んで構成される。振幅スペクトルXは、周波数軸上に配列された複数の強度線740で構成される。各強度線740は、周波数軸に直交する線分である。各周波数に対応する強度線740の長さLAは、その周波数での音響信号xの振幅強度に応じて設定される。図5に拡大して示すように、各強度線740は、初期音成分に対応する第1部分741と残響成分に対応する第2部分742とに区分される。第1部分741と第2部分742とは相異なる表示態様色彩や階調で表示される。第1部分741の長さLA1は初期音比率Gdに応じて設定され、第2部分742の長さLA2は残響比率Grに応じて設定される。具体的には、第1部分741の長さLA1は強度線740の長さLAと初期音比率Gdとの乗算値に設定されLA1=GdLA、第2部分742の長さLA2は強度線740の長さLAと残響比率Grとの乗算値に設定されるLA2=GrLA。すなわち、初期音成分の強度が各強度線740の第1部分741の長さLA1で表現され、残響成分の強度が各強度線740の第2部分742の長さLA2で表現される。例えば、音響信号xの対象単位期間において残響成分が優勢であるほど強度線740のうち第2部分742の割合LA2/LAが増加し、初期音成分が優勢であるほど強度線740のうち第1部分741の割合LA1/LAが増加する。以上の説明から理解されるように、第1部分741の集合は、音響信号xの対象単位期間における初期音成分の振幅スペクトルXdを意味する。すなわち、第2実施形態の残響評価画像74は、音響信号xの振幅スペクトルXと音響信号xの初期音成分の振幅スペクトルXdとを相異なる表示態様で相互に重ねて配置した画像と換言される。音響信号xの振幅スペクトルXと初期音成分の振幅スペクトルXdとの差分各強度線740の第2部分742の集合が残響成分に相当する。第2実施形態では、初期音比率Gdおよび残響比率Grが残響評価画像74として表示装置16に表示されるから、第1実施形態と同様に、音響信号x内の初期音成分や残響成分の多寡音響信号xの全体に対する比率を利用者が視覚的および直観的に把握できるという利点がある。また、第2実施形態では、音響信号xにおける初期音成分や残響成分の比率に加えて初期音成分や残響成分の強度振幅も視覚的に把握できるという利点がある。なお、以上の例示では、強度線740のうち周波数軸に近い第1部分741を初期音成分に対応させ、第2部分742を残響成分に対応させたが、第1部分741を残響成分に対応させるとともに第2部分742を初期音成分に対応させることも可能である。具体的には、第1部分741の長さLA1は強度線740の長さLAと残響比率Grとの乗算値に設定され、第2部分742の長さLA2は強度線740の長さLAと初期音比率Gdとの乗算値に設定される。したがって、第1部分741の集合は、音響信号xのうち残響成分の振幅スペクトルに相当する。また、振幅スペクトルXおよび振幅スペクトルXdの各々の表示/非表示を利用者からの指示に応じて切替える構成も採用され得る。図5に例示した残響評価画像74では、1個の対象単位期間について振幅スペクトルXおよび振幅スペクトルXdを表示したが、図6の例示のように時間軸Tを設定し、振幅スペクトルXおよび振幅スペクトルXdの時間変化を立体的に表示することも可能である。音響信号xの振幅スペクトルXの時間軌跡となる曲面を透過画像で表示し、初期音成分の振幅スペクトルXdの時間軌跡となる曲面を非透過画像で表示すれば、振幅スペクトルXdおよび振幅スペクトルXの双方の時間変化を視認することが可能である。<第3実施形態>第3実施形態の表示制御部40は、音響信号xの初期音成分の平均強度Adと残響成分の平均強度Arとを単位期間毎に順次に算定する。初期音成分の平均強度Adは、以下の数式で表現されるように、音響信号xの振幅スペクトルXと初期音比率Gdとの乗算値の自乗で算定される強度X2を周波数軸上のK個の周波数について平均した数値平均パワーである。他方、残響成分の平均強度Arは、以下の数式で表現されるように、音響信号xの振幅スペクトルXと残響比率Grとの乗算値の自乗で算定される強度X2を周波数軸上のK個の周波数について平均した数値である。表示制御部40は、音響信号xの対象単位期間について図7の残響評価画像76を表示装置16に表示させる。対象単位期間の変更毎に残響評価画像76が更新される構成は第2実施形態と同様である。残響評価画像76は、方向αに長尺な帯状長方形状の強度図形760を含んで構成される。方向αにおける強度図形760の長さLBは、対象単位期間での音響信号xの強度振幅やパワーに応じて選定される。したがって、対象単位期間の変更毎に強度図形760の長さLBは刻々と変化する。図7に示すように、強度図形760は、初期音成分に対応する第1部分761と残響成分に対応する第2部分762とに区分される。第1部分761と第2部分762とは相異なる表示態様色彩や階調で表示される。方向αにおける第1部分761の長さLB1は平均強度Adに応じて設定され、方向αにおける第2部分762の長さLB2は平均強度Arに応じて設定される。具体的には、第1部分761の長さLB1と第2部分762の長さLB2との比LB1/LB2が、平均強度Adと平均強度Arとの比Ad/Arと合致するように、第1部分761の長さLB1と第2部分762の長さLB2とが選定される。すなわち、初期音成分の平均強度Adが第1部分761の長さLB1で表現され、残響成分の平均強度Arが第2部分762の長さLB2で表現される。例えば、音響信号xの対象単位期間にて残響成分が優勢であるほど強度図形760のうち第2部分762の割合LB2/LBが増加し、初期音成分が優勢であるほど強度図形760のうち第1部分761の割合LB1/LBが増加する。第3実施形態では、初期音成分の平均強度Adと残響成分の平均強度Arとが残響評価画像76として表示装置16に表示されるから、第1実施形態と同様に、音響信号xの初期音成分や残響成分の多寡を利用者が視覚的および直観的に把握できるという利点がある。また、第3実施形態では、第2実施形態と同様に、音響信号xにおける初期音成分や残響成分の比率に加えて初期音成分や残響成分の平均強度も視覚的に把握できるという利点がある。なお、以上の例示では、初期音成分の平均強度Adと残響成分の平均強度Arとを図形的に表示したが、平均強度Adおよび平均強度Arの表示方法は任意である。例えば、平均強度Adおよび平均強度Arの各々の数値を残響評価画像76として表示する構成や、平均強度Adと平均強度Arとの比例えばAd/Arを残響評価画像76として表示する構成も採用され得る。<第4実施形態>第4実施形態の音響処理装置100には、右チャネルの音響信号xRと左チャネルの音響信号xLとが信号供給装置12から供給される。音響信号xRと音響信号xLとは、再生音の受聴者が特定の方向に音像を知覚するように生成されたステレオ信号である。音響処理装置100の周波数分析部32は、音響信号xRのスペクトルXRと音響信号xLのスペクトルXLとを単位期間毎に順次に生成する。解析処理部34は、音響信号xRと音響信号xLとを加算した音響信号xを第1実施形態と同様に解析することで初期音比率Gdと残響比率Grとを算定する。残響調整部36は、音響信号xRのスペクトルXRに初期音比率Gdまたは残響比率Grを作用させることで右チャネルの音響信号yRのスペクトルYRを生成し、音響信号xLのスペクトルXLに初期音比率Gdまたは残響比率Grを作用させることで左チャネルの音響信号yLのスペクトルYLを生成する。波形合成部38は、スペクトルYRに応じた音響信号yRとスペクトルYLに応じた音響信号yLとを生成する。表示制御部40は、音響信号xの複数の周波数成分の各々について、その周波数成分が定位する方向以下定位方向というθを単位期間毎に算定する。定位方向θの算定には公知の技術が任意に採用され得るが、例えば、音響信号xRの振幅スペクトルXRと音響信号xLの振幅スペクトルXLとを適用した以下の数式の演算が好適である。なお、数式については、例えば、MVinyesJBonadaALoscosDemixingCommercialMusicProductionswiaHumanAssistedTimeFrequencyMasking,AudioEngineeringSociety120thConventionFrance2006にも開示されている。第4実施形態の表示制御部40は、音響信号xの対象単位期間について図8の残響評価画像78を表示装置16に表示させる。対象単位期間の変更毎に残響評価画像78が更新される構成は第2実施形態と同様である。残響評価画像78は、周波数軸Fと定位軸Pとが設定された座標領域780を含んで構成される。座標領域780内には、音響信号xの各周波数成分に対応する複数の音像点782が配置される。音像点782は、例えば円形状の画像である。各周波数成分に対応する音像点782は、周波数軸Fのうちその周波数成分の周波数と、定位軸Pのうちその周波数成分について算定された定位方向θとに対応する座標に位置する。したがって、複数の音像点782の分布を視認することで、利用者は、音響信号xの各周波数成分について定位方向θを把握することが可能である。例えば、相異なる音域に対応する複数の楽器音の混合音を音響信号xが示す場合、利用者は楽器毎の定位方向θを把握することができる。表示制御部40は、各周波数成分に対応する音像点782のサイズ直径を、その周波数成分の強度振幅に応じて可変に設定する。具体的には、強度が高い周波数成分の音像点782ほど大径に設定される。また、表示制御部40は、各音像点782の表示態様階調や色彩を初期音比率Gdに応じて可変に設定する。例えば、初期音比率Gdが大きい初期音成分が残響成分と比較して優勢であるほど音像点782は低階調濃い階調で表示される。以上に説明したように、第4実施形態では、初期音比率Gdに応じて音像点782の表示態様が制御された残響評価画像78が表示装置16に表示されるから、第1実施形態と同様に、音響信号x内の初期音成分や残響成分の多寡を利用者が視覚的および直観的に把握できるという利点がある。また、第4実施形態では、音響信号xの各周波数成分に対応する音像点782の表示態様が初期音比率Gdに応じて可変に制御されるから、利用者は、各定位方向θに定位する周波数成分において初期音成分および残響成分の何れが優勢であるか、あるいは、初期音成分または残響成分が豊富な定位方向θを、視覚的および直観的に把握することが可能である。なお、以上の例示では、座標領域780内の各音像点782の表示態様を初期音比率Gdに応じて制御したが、残響比率Grに応じて各音像点782の表示態様を制御する構成も採用され得る。また、以上の例示では、1個の対象単位期間のみについて音像点782を表示したが、複数の単位期間にわたる音像点782を相互に重ねて表示することも可能である。具体的には、現在の単位期間の音像点782と過去または未来の単位期間の音像点782とを相異なる表示態様例えば階調や色彩で表示することが可能である。<第5実施形態>図9は、第5実施形態における解析処理部34のブロック図である。第5実施形態の解析処理部34は、図2に例示した第1実施形態の指標値算定部50Aを指標値算定部50Bに置換した構成である。指標値算定部50Bは、第1指標値Q1および第2指標値Q2を単位期間毎に順次に算定する要素であり、第1平滑部51と第2平滑部52と遅延部54とを含んで構成される。なお、評価値算定部60の構成および動作は第1実施形態と同様である。第1平滑部51は、第1実施形態と同様に、音響信号xのパワーX2の時系列を平滑化することで第1指標値Q1を単位期間毎に順次に算定する。遅延部54は、音響信号xのスペクトルXを単位期間のd個分dは自然数に相当する時間だけ遅延させる記憶回路である。第2平滑部52は、遅延部54による遅延後のスペクトルXのパワーX2の時系列を平滑化することで第2指標値Q2を単位期間毎に順次に算定する。ただし、第2実施形態では、第2平滑部52による平滑化の時定数τ2は第1平滑部51による平滑化の時定数τ1と同等であるτ2=τ1。したがって、第2指標値Q2の時間変化は、第1指標値Q1の時間変化を単位期間のd個分だけ遅延させた関係にあるQ2=Q1。なお、第1平滑部51による平滑化の時定数τ1と第2平滑部52による平滑化の時定数τ2とを相違させることも可能である。また、第1平滑部51が算定した第1指標値Q1を遅延させることで第2指標値Q2を算定する構成第2平滑部52を省略した構成も採用され得る。図10の部分は、図3の部分と同様の室内インパルス応答図10の部分を音響信号xとして第5実施形態の音響処理装置100に供給した場合の第1指標値Q1および第2指標値Q2の時間変化のグラフである。図10の部分から理解されるように、第1指標値Q1と第2指標値Q2とで時間変化の態様波形は共通するが、第2指標値Q2の時間変化は第1指標値Q1の時間変化に対して単位期間のd個分だけ遅延する。すなわち、第2指標値Q2は、第1指標値Q1と比較して低い追従性で音響信号xのパワーX2に追従する。したがって、第1実施形態と同様に、第1指標値Q1と第2指標値Q2との大小は時間軸上の特定の時点txで反転する。すなわち、時点txまでの区間SAでは第1指標値Q1が第2指標値Q2を上回り、時点tx以降の区間SBでは第2指標値Q2が第1指標値Q1を上回る。比算定部62による比Rの算定数式や第1処理部64による初期音比率Gdの算定や第2処理部66による残響比率Grの算定は第1実施形態と同様である。したがって、図10の部分に示すように、初期音比率Gdは、初期音成分が存在する区間SAにて1に設定され、残響成分が存在する区間SBでは下限値Gminまで経時的に減少する。したがって、第5実施形態においても第1実施形態と同様の効果が実現される。なお、第2実施形態から第4実施形態に第5実施形態を適用することも可能である。<変形例>前述の各形態は多様に変形され得る。具体的な変形の態様を以下に例示する。以下の例示から任意に選択された2以上の態様は適宜に併合され得る。1前述の各形態では、音響信号xのパワーX2の単純移動平均を第1指標値Q1および第2指標値Q2として算定したが、第1指標値Q1および第2指標値Q2の算定方法は以上の例示に限定されない。例えば、以下の数式および数式で表現されるように、音響信号xのパワーX2の指数平均指数移動平均を第1指標値Q1および第2指標値Q2として算定することも可能である。すなわち、第1平滑部51および第2平滑部52は、IIRinfiniteimpulseresponse型のローパスフィルタに相当する。数式の記号α1および数式の記号α2は平滑化係数忘却係数である。具体的には、平滑化係数α1は、過去の第1指標値Q1に対する現在のパワーX2の重みを意味し、平滑化係数α2は、過去の第2指標値Q2に対する現在のパワーX2の重みを意味する。平滑化係数α2は、平滑化係数α1を下回る数値に設定されるα2<α1。したがって、第1実施形態と同様に、第2平滑部52による平滑化の時定数τ2は第1平滑部51による平滑化の時定数τ1を上回るτ2>τ1。すなわち、第2指標値Q2は、第1指標値Q1と比較して低い追従性で音響信号xのパワーX2に追従する。また、以下の数式および数式で表現されるように、音響信号xのパワーX2の加重移動平均を第1指標値Q1および第2指標値Q2として算定することも可能である。数式の記号w1および数式の記号w2は、第m番目の単位期間からみて前方の第i番目に位置する単位期間に対する加重値を意味する。第2期間が第1期間よりも長いという条件N2>N1は前掲の例示と同様である。2初期音比率Gdや残響比率Grの算定方法は任意である。例えば、前述の各形態の説明から理解されるように、残響比率Grを以下の数式の演算で算定することも可能である。残響比率Grを数式で算定する場合、初期音比率Gdを前掲の数式で算定する構成すなわち、初期音比率Gdと残響比率Grとを並列に算定する構成や、数式で算定された残響比率Grを1から減算することで初期音比率GdGd=1-Grを算定する構成すなわち、第1処理部64と第2処理部66とを入替えた構成も採用され得る。また、例えば、第1指標値Q1および第2指標値Q2を変数とする所定の演算や比Rを変数とする所定の演算により初期音比率Gdや残響比率Grを算定する構成も採用され得る。また、前述の各形態では初期音比率Gdから残響比率Grを算定したが、例えば第1指標値Q1に対する第2指標値Q2の比Rを数式の演算に適用することで残響比率Grを直接的にすなわち初期音比率Gdの算定を経ずに算定することも可能である。以上の説明から理解されるように、初期音比率Gdや残響比率Grは、音響信号のうち残響成分または初期音成分の比率に応じた変数比率評価値として包括され、算定方法の如何は不問である。3前述の各形態では、初期音比率Gdおよび残響比率Grの双方を算定したが、初期音比率Gdおよび残響比率Grの一方のみを算定する構成も採用され得る。例えば第1実施形態や第4実施形態では、残響評価画像に表示される初期音比率Gdのみを算定する残響比率Grの算定を省略することも可能である。4前述の各形態では、第1平滑部51による平滑化の時定数τ1と第2平滑部52による平滑化の時定数τ2との各々を複数の周波数にわたり共通させたが、時定数τ1と時定数τ2とを周波数毎帯域毎に個別に設定することも可能である。また、時定数τ1および時定数τ2の一方または双方を経時的に変化させることも可能である。5前述の各形態では、音響信号xのパワーX2の時系列を平滑化することで第1指標値Q1および第2指標値Q2を算定したが、第1平滑部51や第2平滑部52による平滑化の対象はパワーX2に限定されない。例えば、音響信号xの振幅Xや振幅の4乗X4を平滑化することで第1指標値Q1や第2指標値Q2を算定する構成も採用され得る。すなわち、前述の各形態における第1平滑部51や第2平滑部52は、音響信号xの強度の時系列を平滑化する要素として包括され、音響信号xの強度は、パワーX2のほかに振幅Xや振幅の4乗X4を包含する。6前述の各形態における残響調整部36および波形合成部38を省略することも可能である。すなわち、本発明は、音響信号xの残響を解析して解析結果を表示する残響解析装置としても実施され得る。"