(57)【要約】  低ビットレートの音声符号化においても少ない計算量で符号化性能を落とさずに固定音源符号帳の探索を行う符号化装置。この符号化装置では、閾値算出部（２１１）は、各チャネルでソーティングされた候補位置の相関値の大きさを、複数チャネルに跨って評価して各チャネルの候補位置を特定し、特定した各チャネルの候補位置の相関値を加算して閾値を求める。候補位置ソーティング部（２１２）は、各チャネルのパルス位置に配置されたパルスの合成音と量子化ターゲットとの相関に基づいてパルス候補位置のソーティングを行う。探索制御部（２２１）は、ソーティングされたパルス候補位置で探索を行うように制御を行い、探索対象のチャネルで既に探索された相関値の和が閾値未満の場合には、制御信号を探索打ち切り部（２２３）に出力する。探索打ち切り部（２２３）は、探索制御部（２２１）から入力された制御信号により、そのパルス位置候補における探索を打ち切る。
2008108077.tif 000002
【特許請求の範囲】
【請求項１】
少数のパルスにより固定音源が表現され、チャネル毎にパルスの候補位置が設定される符号帳を用い、チャネル数分の深さのループにより各チャネルの候補位置から最適なパルスの位置を探索する符号化装置であって、
各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて、探索ループの中で探索される各チャネルの候補位置を並び替える候補位置ソーティング手段と、
各チャネルで並び替えられた候補位置の相関値の大きさを複数チャネルに跨って評価し、評価結果に基づいて各チャネルの基準候補位置を特定し、前記各チャネルの基準候補位置の相関値を用いて閾値を求める閾値算出手段と、
並び替えられた候補位置の順でパルスの位置の探索を行い、探索対象のチャネルで既に探索された相関値を用いて設定された代表値と前記閾値との比較結果に基づいて探索を打ち切る探索制御手段と、
を有する符号化装置。
【請求項２】
前記候補位置ソーティング手段は、探索ループの中で探索される各チャネルの候補位置を、前記相関値が大きい順番に並び替える、請求項１に記載の符号化装置。
【請求項３】
前記閾値算出手段は、各チャネルに跨って相関値が大きい順番に候補位置を探索し、探索した候補位置の数の和が予め指定した候補数になった時点で基準候補位置を特定する、請求項２に記載の符号化装置。
【請求項４】
前記閾値算出手段は、各チャネルに跨って相関値が大きい順番に候補位置を探索し、探索した候補位置の数の積が予め指定した積になった時点で基準候補位置を特定する、請求項２に記載の符号化装置。
【請求項５】
前記閾値算出手段は、前記各チャネルの基準候補位置の相関値を加算して前記閾値を求め、
前記探索制御手段は、探索対象のチャネルで既に探索された相関値の和を算出して前記代表値とし、前記代表値が前記閾値未満の場合には探索を打ち切る、請求項１に記載の符号化装置。
【請求項６】
少数のパルスにより固定音源が表現され、チャネル毎にパルスの候補位置が設定される符号帳を用い、チャネル数分の深さのループにより各チャネルの候補位置から最適なパルスの位置を探索する符号化方法であって、
各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて、探索ループの中で探索される各チャネルの候補位置を並び替える候補位置ソーティングステップと、
各チャネルで並び替えられた候補位置の相関値の大きさを複数チャネルに跨って評価し、評価結果に基づいて各チャネルの基準候補位置を特定し、前記各チャネルの基準候補位置の相関値を用いて閾値を求める閾値算出ステップと、
並び替えられた候補位置の順でパルスの位置の探索を行い、探索対象のチャネルで既に探索された相関値を用いて設定された代表値と前記閾値との比較結果に基づいて探索を打ち切る探索制御ステップと、
を有する符号化方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音声信号やオーディオ信号を符号化する符号化装置および符号化方法に関する。
【背景技術】
【０００２】
移動体通信においては、電波などの伝送路容量や記憶媒体の有効利用を図るため、音声や画像のディジタル情報に対して圧縮符号化を行うことが必須であり、これまでに多くの符号化／復号方式が開発されてきた。【０００３】
その中で、音声符号化技術は、音声の発声機構をモデル化してベクトル量子化を巧みに応用した基本方式「ＣＥＬＰ」（Code Excited Linear Prediction）によって性能が大きく向上した。また、オーディオ符号化等の楽音符号化技術は、変換符号化技術（ＭＰＥＧ標準ＡＣＣやＭＰ３等）により性能が大きく向上した。【０００４】
また、音声符号化技術は、代数的符号帳（非特許文献１に記載）の様な少数パルスにより固定音源を表現する発明により、一段とその性能を向上させ、計算量が比較的少なくても良好な性能を得ることができる様になった。【０００５】
適応音源符号帳あるいは固定音源符号帳の音源を探索するアルゴリズムとして、ＣＥＬＰ方式では、計算量を削減するために一般にオープンループ探索が用いられる。【０００６】
以下、固定音源符号帳の従来の探索方法について代数的符号帳を例に説明する。まず、音源の符号の導出は下記式（１）の符号化歪を最小化する音源を探索することにより行われる。【０００７】
【数１】
2008108077.tif 000003
【０００８】
適応音源はオープンループで探索されるので、代数的符号帳の符号の導出は、下記式（２）の符号化歪を最小化する固定音源を探索することにより行われる。【０００９】
【数２】
2008108077.tif 000004
【００１０】
ここで、ゲインｐ、ｑは、音源の符号を探索した後で決定するので、ここでは最適ゲインで探索を進めることとする。すると、上記式（２）は、下記式（３）とすることができる。【００１１】
【数３】
2008108077.tif 000005
【００１２】
そして、この符号化歪の式を最小化することは、下記式（４）の関数Ｃを最大化することと同値であることがわかる。【００１３】
【数４】
2008108077.tif 000006
【００１４】
そこで、代数符号帳の音源のような少数パルスからなる音源の探索の場合は、ｙＨとＨＨを予め計算しておけば、少ない計算量で上記関数Ｃを算出することができる。【００１５】
ここで、代数的符号帳のパルス探索アルゴリズムについて説明する。代数的符号帳の音源は、少数の振幅１で極性（＋、−）のあるパルスで構成されている。以後、サブフレーム３２でパルス本数４（チャネル数４）の場合を例として説明する。【００１６】
代数的符号帳のパルス位置は、下記式（５）のように、お互いが重ならないような配置をとる。【００１７】
【数５】
2008108077.tif 000007
【００１８】
したがって、極性も含めると、この代数的符号帳は（３＋１）×４で１６ビットの符号により符号化されることがわかる。【００１９】
ここで、探索の手順を以下に示す。まず、前処理として、ベクトルｙＨとマトリクスＨＨを算出する。ｙＨはベクトルｙを逆順にしてマトリクスＨを畳み込み、更にその結果を逆順にすることにより求める。ＨＨはマトリクス同士の掛け算により求める。【００２０】
ベクトルｙＨの要素の極性（＋−）から、事前にパルスの極性を決める。具体的には、各位置に立つパルスの極性をｙＨのその位置の値に合わせることとし、ｙＨの値の極性を別の配列に格納しておく。各位置の極性を別の配列に格納した後、ｙＨの値は全て絶対値をとり正の値に変換しておく。また、その極性にあわせてＨＨの値も極性を乗ずることによって変換しておく。【００２１】
４重ループ（後述の実施の形態ではパルスが４本（チャネルが４つ）であるからである）を用いて、ｙＨとＨＨの値を加算することにより関数Ｃを求め、この値が最も大きくなる位置を探索する。各パルスの位置の符号と極性の符号を合わせて固定音源の符号とする。【００２２】
ここで、上記の４重ループを用いた探索について図１及び図２を用いて説明する。図１及び図２は、従来の代数的符号帳探索アルゴリズムのフロー図である。【００２３】
まず、ステップ（以下、ＳＴと省略する）１１で初期化を行い、まず、代数的符号帳から１つ目のパルスｉ０を一つずつ出力してｙＨ及びＨＨから値を取り出し、それぞれ相関値ｓｙ０、パワーｓｈ０とする（ＳＴ１３）。この計算をｉ０が８（パルス位置候補数）になるまで行う（ＳＴ１２〜ＳＴ１４）。【００２４】
一つのｉ０における計算において、代数的符号帳から２つ目のパルスｉ１を一つずつ出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ０、パワーｓｈ０にそれぞれ加算し、相関値ｓｙ１、パワーｓｈ１とする（ＳＴ１６）。この計算をｉ１が８（パルス位置候補数）になるまで行う（ＳＴ１５〜ＳＴ１７）。【００２５】
一つのｉ１における計算において、代数的符号帳から３つ目のパルスｉ２を一つずつ出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ１、ｓｈ１にそれぞれ加算し、相関値ｓｙ２、パワーｓｈ２とする（ＳＴ１９）。この計算をｉ２が８（パルス位置候補数）になるまで行う（ＳＴ１８〜ＳＴ２０）。【００２６】
一つのｉ２における計算において、代数的符号帳から４つ目のパルスｉ３を一つずつ出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ２、ｓｈ２にそれぞれ加算し、相関値ｓｙ３、パワーｓｈ３とする（ＳＴ２２）。この計算をｉ３が８（パルス位置候補数）になるまで行う（ＳＴ２１〜ＳＴ２３）。【００２７】
このようにして得られた相関値ｓｙ３、パワーｓｈ３を用いて、それまでの最大値と比較し（ＳＴ２４）、最も大きな値になる位置を探索する（ＳＴ２５）。このアルゴリズムによって代数的符号帳の探索がなされる。
【非特許文献１】Algebraic Codebook：Salami,Laflamme,Adoul,”8kbit/s ACELP Coding of Speech with 10ms Speech-Frame:a Candidate for CCITT Standardization”,IEEE Proc. ICASSP94,pp.II-97n
【発明の開示】
【発明が解決しようとする課題】
【００２８】
上記のように、代数的符号帳を用いることにより比較的少ない計算量で固定音源符号帳の探索を行うことができるものの、パルスの本数分の深さのループ探索（上記の例では、パルスが４本なので４重ループ）が必要で、パルスの本数によって計算量は指数的に増大する（上記の例では、８の４乗）。特に低ビットレートの音声符号化ではより多くのパルスが必要になるためより多くの計算量が必要になるという問題を有している。【００２９】
本発明の目的は、低ビットレートの音声符号化においても少ない計算量で符号化性能を落とさずに固定音源符号帳の探索を行うことができる符号化装置および符号化方法を提供することである。
【課題を解決するための手段】
【００３０】
本発明の符号化装置は、少数のパルスにより固定音源が表現され、チャネル毎にパルスの候補位置が設定される符号帳を用い、チャネル数分の深さのループにより各チャネルの候補位置から最適なパルスの位置を探索する符号化装置であって、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて、探索ループの中で探索される各チャネルの候補位置を並び替える候補位置ソーティング手段と、各チャネルで並び替えられた候補位置の相関値の大きさを複数チャネルに跨って評価し、評価結果に基づいて各チャネルの基準候補位置を特定し、前記各チャネルの基準候補位置の相関値を用いて閾値を求める閾値算出手段と、並び替えられた候補位置の順でパルスの位置の探索を行い、探索対象のチャネルで既に探索された相関値を用いて設定された代表値と前記閾値との比較結果に基づいて探索を打ち切る探索制御手段と、を有する構成を採る。【００３１】
本発明の符号化方法は、少数のパルスにより固定音源が表現され、チャネル毎にパルスの候補位置が設定される符号帳を用い、チャネル数分の深さのループにより各チャネルの候補位置から最適なパルスの位置を探索する符号化方法であって、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて、探索ループの中で探索される各チャネルの候補位置を並び替える候補位置ソーティングステップと、各チャネルで並び替えられた候補位置の相関値の大きさを複数チャネルに跨って評価し、評価結果に基づいて各チャネルの基準候補位置を特定し、前記各チャネルの基準候補位置の相関値を用いて閾値を求める閾値算出ステップと、並び替えられた候補位置の順でパルスの位置の探索を行い、探索対象のチャネルで既に探索された相関値を用いて設定された代表値と前記閾値との比較結果に基づいて探索を打ち切る探索制御ステップと、を有する方法を採る。
【発明の効果】
【００３２】
本発明によれば、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値が小さい場合の探索処理にかかる計算量を節約することができ、符号化性能を落とさずに平均的計算量を削減することができる。
【図面の簡単な説明】
【００３３】
【図１】従来の代数的符号帳探索アルゴリズムのフロー図
【図２】従来の代数的符号帳探索アルゴリズムのフロー図
【図３】本発明の一実施の形態に係る音声符号化装置の構成を示すブロック図
【図４】本発明の一実施の形態に係る音声符号化装置の音源探索回路の構成を示すブロック図
【図５】本発明の一実施の形態に係るソーティングのアルゴリズムを示すフロー図
【図６】本発明の一実施の形態に係る代数的符号帳探索アルゴリズムのフロー図
【図７】本発明の一実施の形態に係る代数的符号帳探索アルゴリズムのフロー図
【図８】本発明の一実施の形態に係る閾値を求めるアルゴリズムを示すフロー図
【図９】本発明の一実施の形態に係る閾値を求めるアルゴリズムを示すフロー図
【図１０】本発明の一実施の形態に係るソーティングされた相関値の例を示す図
【図１１】本発明の一実施の形態に係るソーティングされた相関値の例を示す図
【図１２】本発明の一実施の形態に係る閾値を求めるアルゴリズムを示すフロー図
【図１３】本発明の一実施の形態に係るソーティングされた相関値の例を示す図
【発明を実施するための最良の形態】
【００３４】
本実施の形態においては、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて探索中断の判断のための閾値を設け、ループ探索の途中で相関値の和が閾値を下回る場合にそれ以上深いループ探索を行わない。さらに、本実施の形態では、各チャネルでソーティングされた候補位置の相関値の大きさを、複数チャネルに跨って評価して各チャネルの候補位置を特定し、特定した各チャネルの候補位置の相関値を加算して閾値を求める。【００３５】
以下、本発明の一実施の形態について、図面を用いて説明する。【００３６】
図３は、本実施の形態に係る音声符号化装置の構成を示すブロック図である。【００３７】
前処理部１０１は、入力音声信号に対し、ＤＣ成分を取り除くハイパスフィルタ処理や後続する符号化処理の性能改善につながるような波形整形処理やプリエンファシス処理を行い、これらの処理後の信号（Xin）をＬＰＣ分析部１０２および加算部１０５に出力する。【００３８】
ＬＰＣ分析部１０２は、Xinを用いて線形予測分析を行い、分析結果（線形予測係数）をＬＰＣ量子化部１０３に出力する。ＬＰＣ量子化部１０３は、ＬＰＣ分析部１０２から出力された線形予測係数（ＬＰＣ）の量子化処理を行い、量子化ＬＰＣを合成フィルタ１０４に出力するとともに量子化ＬＰＣを表す符号（Ｌ）を多重化部１１４に出力する。【００３９】
合成フィルタ１０４は、量子化ＬＰＣに基づくフィルタ係数により、後述する加算部１１１から出力される駆動音源に対してフィルタ合成を行うことにより合成信号を生成し、合成信号を加算部１０５に出力する。【００４０】
加算部１０５は、合成信号の極性を反転させてXinに加算することにより誤差信号を算出し、誤差信号を聴覚重み付け部１１２に出力する。【００４１】
適応音源符号帳１０６は、過去に加算部１１１によって出力された駆動音源をバッファに記憶し、パラメータ決定部１１３から出力された信号により特定される過去の駆動音源から１フレーム分のサンプルを適応音源ベクトルとして切り出して乗算部１０９に出力する。【００４２】
ゲイン符号帳１０７は、パラメータ決定部１１３から出力された信号によって特定される適応音源ベクトルのゲインと固定音源ベクトルのゲインとをそれぞれ乗算部１０９と乗算部１１０とに出力する。【００４３】
固定音源符号帳１０８は、所定形状のパルス音源ベクトルをバッファに複数記憶し、パラメータ決定部１１３から出力された信号によって特定される形状を有するパルス音源ベクトルを固定音源ベクトルとして乗算部１１０に出力する。なお、パルス音源ベクトルに拡散ベクトルを乗算して得られたものを固定音源ベクトルとして乗算部１１０に出力しても良い。【００４４】
乗算部１０９は、ゲイン符号帳１０７から出力されたゲインを、適応音源符号帳１０６から出力された適応音源ベクトルに乗じて、加算部１１１に出力する。乗算部１１０は、ゲイン符号帳１０７から出力されたゲインを、固定音源符号帳１０８から出力された固定音源ベクトルに乗じて、加算部１１１に出力する。【００４５】
加算部１１１は、利得乗算後の適応音源ベクトルと利得乗算後の固定音源ベクトルとをそれぞれ乗算部１０９と乗算部１１０とから入力し、これらをベクトル加算し、加算結果である駆動音源を合成フィルタ１０４および適応音源符号帳１０６に出力する。なお、適応音源符号帳１０６に入力された駆動音源は、バッファに記憶される。【００４６】
聴覚重み付け部１１２は、加算部１０５から出力された誤差信号に対して聴覚的な重み付けをおこない符号化歪みとしてパラメータ決定部１１３に出力する。【００４７】
パラメータ決定部１１３は、聴覚重み付け部１１２から出力された符号化歪みを最小とする適応音源ベクトル、固定音源ベクトル及び量子化利得の符号を探索し、探索された適応音源ベクトルを表す符号（Ａ）、固定音源ベクトルを表す符号（Ｆ）及び量子化利得を表す符号（Ｇ）を多重化部１１４に出力する。【００４８】
多重化部１１４は、ＬＰＣ量子化部１０３から量子化ＬＰＣを表す符号（Ｌ）を入力し、パラメータ決定部１１３から適応音源ベクトルを表す符号（Ａ）、固定音源ベクトルを表す符号（Ｆ）および量子化利得を表す符号（Ｇ）を入力し、これらの情報を多重化して符号化情報として出力する。【００４９】
図４は、本実施の形態に係る音声符号化装置の音源探索回路の構成を示すブロック図である。この音源探索回路は、図３に示した音声符号化装置のパラメータ決定部１１３に設けられる。【００５０】
図４に示す音源探索回路１５０は、適応音源符号帳１０６の適応音源を探索する適応音源探索部１５１と、固定音源符号帳１０８の固定音源を探索する固定音源探索部１５２とから構成されている。【００５１】
固定音源探索部１５２は、前処理部２０１と、探索部２０２とから構成されている。前処理部２０１は、閾値算出部２１１と、候補位置ソーティング部２１２と、カウンタクリア部２１３と、制限回数設定部２１４と、を有する。また、探索部２０２は、探索制御部２２１と、カウンタ２２２と、探索打ち切り部２２３と、を有する。【００５２】
この固定音源探索部１５２において、前処理部２０１の閾値算出部２１１は、探索を打ち切るための閾値、すなわちこれ以上深い探索を行うか否かを判定するための相関値の閾値を算出する。なお、本実施の形態に係る閾値の算出方法の詳細については後述する。【００５３】
前処理部２０１の候補位置ソーティング部２１２は、チャネル毎にすべての要素の値が正の値に変換されたベクトルｙＨの要素の大きさに基づいて、要素の大きい順に、パルス候補位置を各チャネルから出力する順番を並び替える。この並び替えの順序は、閾値算出部２１１および探索部２０２の探索制御部２２１に出力される。【００５４】
前処理部２０１のカウンタクリア部２１３は、探索部２０２のカウンタ２２２をリセットする。【００５５】
前処理部２０１の制限回数設定部２１４は、カウンタ２２２の値を参照して探索を中断するための制限回数Ｒを定める。この制限回数Ｒは、設定された後にカウンタ２２２に出力され、カウンタ２２２は、カウンタ数が制限回数設定部２１４で設定された制限回数Ｒを超えた場合には、探索打ち切り部２２３に制御信号を送る。【００５６】
探索部２０２の探索制御部２２１は、候補位置ソーティング部２１２にてソーティングされたパルス候補位置で探索を行うように制御を行う。また、探索制御部２２１は、探索対象のチャネルで既に探索された相関値の和に対して閾値判定を行う。そして、探索制御部２２１は、閾値以上の場合には、制御信号をカウンタ２２２に出力し、閾値未満の場合には、制御信号を探索打ち切り部２２３に出力する。【００５７】
探索部２０２のカウンタ２２２は、閾値以上と判定された数をカウントし、このカウント数と制限回数Ｒとの比較を行い、Ｒを超えたときには、全探索を終了させる制御信号を探索打ち切り部２２３に出力する。【００５８】
探索部２０２の探索打ち切り部２２３は、探索制御部２２１あるいはカウンタ２２２から入力された制御信号により、そのパルス位置候補における探索を打ち切り、全探索を終了させる。【００５９】
上記構成を有する固定音源探索部１５２における固定音源符号帳の探索方法について説明する。まず、本発明の代数的符号帳のパルス探索アルゴリズムについて説明する。【００６０】
代数的符号帳の音源は、少数の振幅１で極性（＋、−）のある複数のパルスで構成される。以後、サブフレーム長３２でパルス本数４の場合を例として説明する。代数的符号帳のパルス位置は式（５）のようにお互いが重ならないような配置をとる。したがって、極性も含めると、この代数的符号帳は（３＋１）×４で１６ビットの符号により符号化されることがわかる。【００６１】
探索の手順を以下に示す。まず、前処理として、ベクトルｙＨとマトリクスＨＨを算出する。ｙＨはベクトルｙを逆順にしてマトリクスＨを畳み込み、更にその結果を逆順にすることにより求める。ＨＨはマトリクス同士の掛け算により求める。【００６２】
次いで、ベクトルｙＨの要素の極性（＋−）から、事前にパルスの極性を決める。具体的には、各位置に立つパルスの極性をｙＨのその位置の値に合わせることとし、ｙＨの極性を別の配列に格納する（後にこの配列の極性を参照して極性の符号を求める）。各位置の極性を別の配列に格納したら、ｙＨの値は全て絶対値をとり正の値に変換しておく。また、その極性にあわせてＨＨの値に極性を乗ずることによって変換しておく。【００６３】
次いで、候補位置ソーティング部２１２において、各パルスの候補位置のｙＨの値を参照して、値の大きい順にソーティングを行い、候補位置を別配列に格納する。この並べ替えは、８つの並べ替えなので少ない計算量で実現できる。【００６４】
ベクトルｙＨのある位置の値が大きいということは、ターゲットとその位置のパルスの合成音との相関が高いということを示しており、最適なパルスの組み合わせの一つとして選ばれる可能性が高いということが言える。したがって、ベクトルｙＨの値の大きい順に候補位置を並べ替えることにより、相関の高いパルス候補位置から探索することになり、探索を途中で打ち切ったとしても、最適な組み合わせについての探索が漏れてしまう確率を小さくすることが可能となる。その結果、少ない演算量で高い確率で最適なパルス位置の組み合わせを探索することができる。【００６５】
このソーティングアルゴリズムについて図５を用いて説明する。なお、例としてチャネル０について示す。【００６６】
まず、ｉを初期化し（ＳＴ３０１）、８つのパルス候補位置（０，４，８，１２，１６，２０，２４，２８）のｙＨの値を求め（ＳＴ３０３）、その値をＺ[ｉ]に書き込む（ＳＴ３０２〜ＳＴ３０４）。【００６７】
次いで、再度ｉを初期化し（ＳＴ３０５）、Ｚｍａｘを選択されない十分に小さい値とする（ＳＴ３０７）。そして、ｊを増やしながらＺｍａｘとＺ[ｊ]とを比較し（ＳＴ３１０）、ＺｍａｘよりもＺ[ｊ]が大きければ、ＺｍａｘをＺ[ｊ]に書きかえる（ＳＴ３１２）。Ｚｍａｘの初期値は小さく設定しているので、最初の比較では必ず書きかえが起こることとなる。そして、この書きかえられたＺ[ｊ]のｊをｊｊとして保存する。これをｊが８になるまで行う（ＳＴ３０９〜ＳＴ３１２）。【００６８】
これにより、最大のＺｍａｘとその時のインデクスｊｊが決定される。そして、Ｚ[ｊｊ]を選択されない小さい値とする（ＳＴ３１３）。これにより、一度決定したｊｊは２度と選択されないこととなり、残りの７つについて上記と同じような処理を行って、順に大きい値と、そのときのインデクスが決定されていくことになる（ＳＴ３０６，ＳＴ３０８）。同様にして、チャネル１〜３についても並べ替えを行う。【００６９】
次に、上記ソーティングを行ったパルス候補位置で探索を行う場合について図６及び図７を用いて説明する。図６及び図７は、本実施の形態に係る代数的符号帳探索アルゴリズムのフロー図である。【００７０】
まず、ＳＴ４０１で初期化を行い、固定音源符号帳１０８から１つ目のパルスの位置を出力してｙＨ及びＨＨから値を取り出して、それぞれ相関値ｓｙ０、パワーｓｈ０とする（ＳＴ４０３）。この計算をｉ０が８（パルス位置候補数）になるまで行う（ＳＴ４０２〜ＳＴ４０４）。【００７１】
一つのｉ０における計算において、固定音源符号帳１０８から２つ目のパルスの位置を出力してｙＨ及びＨＨから値を取り出して相関値ｓｙ０、パワーｓｈ０にそれぞれ加算し、相関値ｓｙ１、パワーｓｈ１とする（ＳＴ４０６）。この計算をｉ１が８（パルス位置候補数）になるまで行う（ＳＴ４０５〜ＳＴ４０７）。【００７２】
一つのｉ１における計算において、固定音源符号帳１０８から３つ目のパルスの位置を出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ１、ｓｈ１にそれぞれ加算し、相関値ｓｙ２、パワーｓｈ２とする（ＳＴ４０９）。ここで、探索を打ち切るかどうかを判断する（ＳＴ４１０）。このように本実施の形態では、第３ループが終わった段階で判定する場合について説明を行うこととする。なお、他のループ終了時に判定しても同様の効果が得られる。【００７３】
ＳＴ４１０において、相関値ｓｙ２が閾値Ｔｈ以上であれば、カウンタ２２２のカウント数を１インクリメントする（ＳＴ４１２）。相関値ｓｙ２が閾値Ｔｈ以上であれば、相関値が高いということであり、更に深い探索ループに進むことになる。すなわち、この計算をｉ２が８（パルス位置候補数）になるまで行う（ＳＴ４０８〜ＳＴ４１１，ＳＴ４１４）。【００７４】
次いで、一つのｉ２における計算において、固定音源符号帳１０８から４つ目のパルスの位置を出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ２、ｓｈ２にそれぞれ加算し、相関値ｓｙ３、パワーｓｈ３とする（ＳＴ４１６）。この計算をｉ３が８（パルス位置候補数）になるまで行う（ＳＴ４１５〜ＳＴ４１８）。【００７５】
このようにして得られた相関値ｓｙ３とパワーｓｈ３を用いて、それまでの最大値と比較し（ＳＴ４１７）、最も大きな値になる位置を探索する（ＳＴ４１９）。このアルゴリズムによって代数的符号帳の探索がなされる。【００７６】
すなわち、上記処理においては、予め指定された深さのループまで探索した後で、それまでに加算した相関値の和が上記閾値を上回る場合は更に深いループの探索を行い、下回る場合はそれより深いループの探索を行わない。【００７７】
また、カウンタ２２２では、カウンタ数と制限回数Ｒとを比較し（ＳＴ４１３）、カウンタ数Ｃが制限回数Ｒを超えた場合には、探索打ち切り部２２３に制御信号が送られ、探索打ち切り部２２３で探索が打ち切られる（探索中断）。一方、カウンタ数Ｃが制限回数Ｒ未満である場合には、上述のように更に深い探索ループに進む。【００７８】
すなわち、ここでの処理においては、予め指定された所定の深さのループを探索する数を数えるカウンタを用い、カウンタ２２２により、ある深さより深いループを探索すると判定された場合の数をカウントし、カウンタの値が予め指定された制限回数を超えた場合に、探索を中断する。【００７９】
この４重ループ（パルスが４本であるからである）を用いて、ｙＨとＨＨの値から相関値とパワーを算出し、相関が最も大きくなるパルス位置を探索する。そして各パルスの位置の符号と極性の符号（極性は前述したように、別配列に格納されている）を合わせて固定音源の符号とする。【００８０】
次に、本実施の形態に係る閾値の算出方法の詳細について説明する。本実施の形態では、各チャネルでソーティングされた候補位置の相関値の大きさを複数チャネルに跨って評価し、各チャネルの閾値算出に用いる候補位置（以下、「基準候補位置」という）を特定し、各チャネルの基準候補位置の相関値を加算して閾値を求める。【００８１】
本実施の形態では、各チャネルに跨って相関値が大きい順番に候補位置を探索し、探索した候補位置の数の和が予め指定した候補数になった時点で基準候補位置を特定する。また、この時、各チャネルのエントリ数を越えないようにソーティングすることが必要になる。【００８２】
図８は、本実施の形態に係る閾値を求めるアルゴリズムを示すフロー図であり、３パルスまでの結果で４パルス目の探索を行うかどうかを判定する場合を示す。なお、図８のフロー図において、ｉ０，ｉ１，ｉ２は閾値を求める各チャネルの順位、ｃはカウンタ、Ｔｈは閾値、Ｍは閾値の強度を決める候補数、を示す。Ｍは定数で１２〜１６程度に設定される。【００８３】
また、図９は、本実施の形態に係る閾値を求めるアルゴリズムを示すフロー図であり、２パルスまでの結果で４パルス目の探索を行うかどうかを判定する場合を示す。なお、図９のフロー図において、ｉ０，ｉ１は閾値を求める各チャネルの順位、ｃはカウンタ、Ｔｈは閾値、Ｍは閾値の強度を決める候補数、を示す。Ｍは定数で８〜１２程度に設定される。【００８４】
以下、本実施の形態の閾値の算出方法を、具体例を用いて説明する。例えば、ソーティングされた相関値が以下の図１０に示すものであり、候補数Ｍが「１０」であるとする。この場合、図８に示したアルゴリズムを実行することにより、候補位置の探索は、図１１の枠で囲った部分にまで進む。【００８５】
よって、本例の場合、各チャネルの基準候補位置は、ｉ０＝５、ｉ１＝２、ｉ２＝３となり、閾値Ｔｈは、Ｔｈ＝７＋７＋６＝２０となる。【００８６】
なお、本実施の形態では、各チャネルに跨って相関値が大きい順番に候補位置を探索し、探索した候補位置の数の積が予め指定した積になった時点で基準候補位置を特定しても良い。【００８７】
図１２は、この方法により閾値を求めるアルゴリズムを示すフロー図である。なお、図１２のフロー図において、ｉ０，ｉ１，ｉ２は閾値を求める各チャネルの順位、ｃはカウンタ、Ｔｈは閾値、Ｎは閾値の強度を決める積、を示す。Ｎは定数で３０〜６０程度に設定される。【００８８】
この場合において、例えば、ソーティングされた相関値が以下の上記図１０に示すものであり、積Ｎが「３５」であるとする。この場合、図１２に示したアルゴリズムを実行することにより、「６×２×３＝３６＞３５」で、候補位置の探索は、図１３の枠で囲った部分にまで進む。【００８９】
よって、本例の場合、各チャネルの基準候補位置は、ｉ０＝６、ｉ１＝２、ｉ２＝３となり、閾値Ｔｈは、Ｔｈ＝５＋７＋６＝１８となる。【００９０】
このように、本実施の形態では、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて探索中断の判断のための閾値を設け、ループ探索の途中で相関値の和が閾値を下回る場合にそれ以上深いループ探索を行わない。しかも、本実施の形態では、各チャネルでソーティングされた候補位置の相関値の大きさを、複数チャネルに跨って評価して各チャネルの候補位置を特定し、特定した各チャネルの候補位置の相関値を加算して閾値を求める。【００９１】
これにより、相関値が小さい場合の探索処理にかかる計算量を節約することができ、符号化性能を落とさずに平均的計算量を削減することができる。【００９２】
なお、本発明は上記実施の形態に限定されず、種々変更して実施することが可能である。例えば、上記実施の形態では、第３ループの深さで探索の打ち切りをするか否かを判断した場合について説明しているが、本発明では、他のループの深さで探索の打ち切りをするか否かを判断するようにしても良い。【００９３】
また、上記実施の形態に係るパルス候補位置のソーティングや探索は、音声符号化装置において行うように説明しているが、このソーティングや探索をソフトウェアとして構成しても良い。例えば、上記ソーティングや探索のプログラムをＲＯＭに格納し、そのプログラムにしたがってＣＰＵの指示により動作させるように構成しても良い。また、ソーティングや探索のプログラムをコンピュータで読み取り可能な記憶媒体に格納し、この記憶媒体のソーティングや探索のプログラムをコンピュータのＲＡＭに記録して、このプログラムにしたがって動作させるようにしても良い。このような場合においても、上記実施の形態と同様の作用、効果を呈する。【００９４】
また、本実施の形態では、固定音源符号帳として代数的符号帳を用いた場合の例を示したが、本発明はこれに限られず、マルチパルス符号帳や、固定波形がＲＯＭに書かれた固定符号帳でも有効である。この場合、各チャネルのそれぞれの番号が固定波形のベクトルに対応する。【００９５】
また、本実施の形態では、ＣＥＬＰに対して用いる場合について説明したが、本発明はこれに限られず、音源の符号帳が存在する他の符号化にも有効である。これは、本発明の所在は固定音源符号帳ベクトル内のものであり、適応音源符号帳の有無や、スペクトル包絡の分析方法に依存しないからである。【００９６】
また、本実施の形態では、固定音源符号帳に用いた例を示したが、本発明はこれに限られず、スペクトル符号化の多段ベクトル量子化、画像符号化におけるベクトルの多段ベクトル量子化にも有効である。これは、本発明が、多重ループの探索における候補数の削減に貢献するものであり、符号化対象には限定されないからである。【００９７】
また、本発明に係る信号は、音声信号だけでなく、オーディオ信号でも良い。また、入力信号の代わりに、ＬＰＣ予測残差信号に対して本発明を適用する構成であっても良い。【００９８】
また、本発明に係る符号化装置は、移動体通信システムにおける通信端末装置および基地局装置に搭載することが可能であり、これにより上記と同様の作用効果を有する通信端末装置、基地局装置、および移動体通信システムを提供することができる。【００９９】
また、ここでは、本発明をハードウェアで構成する場合を例にとって説明したが、本発明をソフトウェアで実現することも可能である。例えば、本発明に係るアルゴリズムをプログラミング言語によって記述し、このプログラムをメモリに記憶しておいて情報処理手段によって実行させることにより、本発明に係る符号化装置と同様の機能を実現することができる。【０１００】
また、上記実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されても良いし、一部または全てを含むように１チップ化されても良い。【０１０１】
また、ここではＬＳＩとしたが、集積度の違いによって、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩ等と呼称されることもある。【０１０２】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現しても良い。ＬＳＩ製造後に、プログラム化することが可能なＦＰＧＡ（Field Programmable Gate Array）や、ＬＳＩ内部の回路セルの接続もしくは設定を再構成可能なリコンフィギュラブル・プロセッサを利用しても良い。【０１０３】
さらに、半導体技術の進歩または派生する別技術により、ＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行っても良い。バイオ技術の適用等が可能性としてあり得る。【０１０４】
２００７年３月２日出願の特願２００７−０５３５０１の日本出願に含まれる明細書、図面および要約書の開示内容は、すべて本願に援用される。
【産業上の利用可能性】
【０１０５】
本発明は、音声信号やオーディオ信号を符号化する符号化装置等に用いるに好適である。
【図１】
2008108077.tif 000008
【図２】
2008108077.tif 000009
【図３】
2008108077.tif 000010
【図４】
2008108077.tif 000011
【図５】
2008108077.tif 000012
【図６】
2008108077.tif 000013
【図７】
2008108077.tif 000014
【図８】
2008108077.tif 000015
【図９】
2008108077.tif 000016
【図１０】
2008108077.tif 000017
【図１１】
2008108077.tif 000018
【図１２】
2008108077.tif 000019
【図１３】
2008108077.tif 000020
【手続補正書】
【提出日】平成21年8月12日(2009.8.12)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
少数のパルスにより固定音源が表現され、チャネル毎にパルスの候補位置が設定される符号帳を用い、チャネル数分の深さのループにより各チャネルの候補位置から最適なパルスの位置を探索する符号化装置であって、
各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて、探索ループの中で探索される各チャネルの候補位置を並び替える候補位置ソーティング手段と、
各チャネルで並び替えられた候補位置の相関値の大きさを複数チャネルに跨って評価し、評価結果に基づいて各チャネルの基準候補位置を特定し、前記各チャネルの基準候補位置の相関値を用いて閾値を求める閾値算出手段と、
並び替えられた候補位置の順でパルスの位置の探索を行い、探索対象のチャネルで既に探索された相関値を用いて設定された代表値と前記閾値との比較結果に基づいて探索を打ち切る探索制御手段と、
を有する符号化装置。
【請求項２】
前記候補位置ソーティング手段は、探索ループの中で探索される各チャネルの候補位置を、前記相関値が大きい順番に並び替える、請求項１に記載の符号化装置。
【請求項３】
前記閾値算出手段は、各チャネルに跨って相関値が大きい順番に候補位置を探索し、探索した候補位置の数の和が予め指定した候補数になった時点で基準候補位置を特定する、請求項２に記載の符号化装置。
【請求項４】
前記閾値算出手段は、各チャネルに跨って相関値が大きい順番に候補位置を探索し、探索した候補位置の数の積が予め指定した積になった時点で基準候補位置を特定する、請求項２に記載の符号化装置。
【請求項５】
前記閾値算出手段は、前記各チャネルの基準候補位置の相関値を加算して前記閾値を求め、
前記探索制御手段は、探索対象のチャネルで既に探索された相関値の和を算出して前記代表値とし、前記代表値が前記閾値未満の場合には探索を打ち切る、請求項１に記載の符号化装置。
【請求項６】
少数のパルスにより固定音源が表現され、チャネル毎にパルスの候補位置が設定される符号帳を用い、チャネル数分の深さのループにより各チャネルの候補位置から最適なパルスの位置を探索する符号化方法であって、
各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて、探索ループの中で探索される各チャネルの候補位置を並び替える候補位置ソーティングステップと、
各チャネルで並び替えられた候補位置の相関値の大きさを複数チャネルに跨って評価し、評価結果に基づいて各チャネルの基準候補位置を特定し、前記各チャネルの基準候補位置の相関値を用いて閾値を求める閾値算出ステップと、
並び替えられた候補位置の順でパルスの位置の探索を行い、探索対象のチャネルで既に探索された相関値を用いて設定された代表値と前記閾値との比較結果に基づいて探索を打ち切る探索制御ステップと、
を有する符号化方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音声信号やオーディオ信号を符号化する符号化装置および符号化方法に関する。
【背景技術】
【０００２】
移動体通信においては、電波などの伝送路容量や記憶媒体の有効利用を図るため、音声や画像のディジタル情報に対して圧縮符号化を行うことが必須であり、これまでに多くの符号化／復号方式が開発されてきた。
【０００３】
その中で、音声符号化技術は、音声の発声機構をモデル化してベクトル量子化を巧みに応用した基本方式「ＣＥＬＰ」（Code Excited Linear Prediction）によって性能が大きく向上した。また、オーディオ符号化等の楽音符号化技術は、変換符号化技術（ＭＰＥＧ標準ＡＣＣやＭＰ３等）により性能が大きく向上した。
【０００４】
また、音声符号化技術は、代数的符号帳（非特許文献１に記載）の様な少数パルスにより固定音源を表現する発明により、一段とその性能を向上させ、計算量が比較的少なくても良好な性能を得ることができる様になった。
【０００５】
適応音源符号帳あるいは固定音源符号帳の音源を探索するアルゴリズムとして、ＣＥＬＰ方式では、計算量を削減するために一般にオープンループ探索が用いられる。
【０００６】
以下、固定音源符号帳の従来の探索方法について代数的符号帳を例に説明する。まず、音源の符号の導出は下記式（１）の符号化歪を最小化する音源を探索することにより行われる。
【０００７】
【数１】
2008108077.tif 000021
【０００８】
適応音源はオープンループで探索されるので、代数的符号帳の符号の導出は、下記式（２）の符号化歪を最小化する固定音源を探索することにより行われる。
【０００９】
【数２】
2008108077.tif 000022
【００１０】
ここで、ゲインｐ、ｑは、音源の符号を探索した後で決定するので、ここでは最適ゲインで探索を進めることとする。すると、上記式（２）は、下記式（３）とすることができる。
【００１１】
【数３】
2008108077.tif 000023
【００１２】
そして、この符号化歪の式を最小化することは、下記式（４）の関数Ｃを最大化することと同値であることがわかる。
【００１３】
【数４】
2008108077.tif 000024
【００１４】
そこで、代数符号帳の音源のような少数パルスからなる音源の探索の場合は、ｙＨとＨＨを予め計算しておけば、少ない計算量で上記関数Ｃを算出することができる。
【００１５】
ここで、代数的符号帳のパルス探索アルゴリズムについて説明する。代数的符号帳の音源は、少数の振幅１で極性（＋、−）のあるパルスで構成されている。以後、サブフレーム３２でパルス本数４（チャネル数４）の場合を例として説明する。
【００１６】
代数的符号帳のパルス位置は、下記式（５）のように、お互いが重ならないような配置をとる。
【００１７】
【数５】
2008108077.tif 000025
【００１８】
したがって、極性も含めると、この代数的符号帳は（３＋１）×４で１６ビットの符号により符号化されることがわかる。
【００１９】
ここで、探索の手順を以下に示す。まず、前処理として、ベクトルｙＨとマトリクスＨＨを算出する。ｙＨはベクトルｙを逆順にしてマトリクスＨを畳み込み、更にその結果を逆順にすることにより求める。ＨＨはマトリクス同士の掛け算により求める。
【００２０】
ベクトルｙＨの要素の極性（＋−）から、事前にパルスの極性を決める。具体的には、各位置に立つパルスの極性をｙＨのその位置の値に合わせることとし、ｙＨの値の極性を別の配列に格納しておく。各位置の極性を別の配列に格納した後、ｙＨの値は全て絶対値をとり正の値に変換しておく。また、その極性にあわせてＨＨの値も極性を乗ずることによって変換しておく。
【００２１】
４重ループ（後述の実施の形態ではパルスが４本（チャネルが４つ）であるからである）を用いて、ｙＨとＨＨの値を加算することにより関数Ｃを求め、この値が最も大きくなる位置を探索する。各パルスの位置の符号と極性の符号を合わせて固定音源の符号とする。
【００２２】
ここで、上記の４重ループを用いた探索について図１及び図２を用いて説明する。図１及び図２は、従来の代数的符号帳探索アルゴリズムのフロー図である。
【００２３】
まず、ステップ（以下、ＳＴと省略する）１１で初期化を行い、まず、代数的符号帳から１つ目のパルスｉ０を一つずつ出力してｙＨ及びＨＨから値を取り出し、それぞれ相関値ｓｙ０、パワーｓｈ０とする（ＳＴ１３）。この計算をｉ０が８（パルス位置候補数）になるまで行う（ＳＴ１２〜ＳＴ１４）。
【００２４】
一つのｉ０における計算において、代数的符号帳から２つ目のパルスｉ１を一つずつ出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ０、パワーｓｈ０にそれぞれ加算し、相関値ｓｙ１、パワーｓｈ１とする（ＳＴ１６）。この計算をｉ１が８（パルス位置候補数）になるまで行う（ＳＴ１５〜ＳＴ１７）。
【００２５】
一つのｉ１における計算において、代数的符号帳から３つ目のパルスｉ２を一つずつ出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ１、ｓｈ１にそれぞれ加算し、相関値ｓｙ２、パワーｓｈ２とする（ＳＴ１９）。この計算をｉ２が８（パルス位置候補数）になるまで行う（ＳＴ１８〜ＳＴ２０）。
【００２６】
一つのｉ２における計算において、代数的符号帳から４つ目のパルスｉ３を一つずつ出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ２、ｓｈ２にそれぞれ加算し、相関値ｓｙ３、パワーｓｈ３とする（ＳＴ２２）。この計算をｉ３が８（パルス位置候補数）になるまで行う（ＳＴ２１〜ＳＴ２３）。
【００２７】
このようにして得られた相関値ｓｙ３、パワーｓｈ３を用いて、それまでの最大値と比較し（ＳＴ２４）、最も大きな値になる位置を探索する（ＳＴ２５）。このアルゴリズムによって代数的符号帳の探索がなされる。
【非特許文献１】Algebraic Codebook：Salami,Laflamme,Adoul,”8kbit/s ACELP Coding of Speech with 10ms Speech-Frame:a Candidate for CCITT Standardization”,IEEE Proc. ICASSP94,pp.II-97n
【発明の開示】
【発明が解決しようとする課題】
【００２８】
上記のように、代数的符号帳を用いることにより比較的少ない計算量で固定音源符号帳の探索を行うことができるものの、パルスの本数分の深さのループ探索（上記の例では、パルスが４本なので４重ループ）が必要で、パルスの本数によって計算量は指数的に増大する（上記の例では、８の４乗）。特に低ビットレートの音声符号化ではより多くのパルスが必要になるためより多くの計算量が必要になるという問題を有している。
【００２９】
本発明の目的は、低ビットレートの音声符号化においても少ない計算量で符号化性能を落とさずに固定音源符号帳の探索を行うことができる符号化装置および符号化方法を提供することである。
【課題を解決するための手段】
【００３０】
本発明の符号化装置は、少数のパルスにより固定音源が表現され、チャネル毎にパルスの候補位置が設定される符号帳を用い、チャネル数分の深さのループにより各チャネルの候補位置から最適なパルスの位置を探索する符号化装置であって、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて、探索ループの中で探索される各チャネルの候補位置を並び替える候補位置ソーティング手段と、各チャネルで並び替えられた候補位置の相関値の大きさを複数チャネルに跨って評価し、評価結果に基づいて各チャネルの基準候補位置を特定し、前記各チャネルの基準候補位置の相関値を用いて閾値を求める閾値算出手段と、並び替えられた候補位置の順でパルスの位置の探索を行い、探索対象のチャネルで既に探索された相関値を用いて設定された代表値と前記閾値との比較結果に基づいて探索を打ち切る探索制御手段と、を有する構成を採る。
【００３１】
本発明の符号化方法は、少数のパルスにより固定音源が表現され、チャネル毎にパルスの候補位置が設定される符号帳を用い、チャネル数分の深さのループにより各チャネルの候補位置から最適なパルスの位置を探索する符号化方法であって、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて、探索ループの中で探索される各チャネルの候補位置を並び替える候補位置ソーティングステップと、各チャネルで並び替えられた候補位置の相関値の大きさを複数チャネルに跨って評価し、評価結果に基づいて各チャネルの基準候補位置を特定し、前記各チャネルの基準候補位置の相関値を用いて閾値を求める閾値算出ステップと、並び替えられた候補位置の順でパルスの位置の探索を行い、探索対象のチャネルで既に探索された相関値を用いて設定された代表値と前記閾値との比較結果に基づいて探索を打ち切る探索制御ステップと、を有する方法を採る。
【発明の効果】
【００３２】
本発明によれば、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値が小さい場合の探索処理にかかる計算量を節約することができ、符号化性能を落とさずに平均的計算量を削減することができる。
【発明を実施するための最良の形態】
【００３３】
本実施の形態においては、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて探索中断の判断のための閾値を設け、ループ探索の途中で相関値の和が閾値を下回る場合にそれ以上深いループ探索を行わない。さらに、本実施の形態では、各チャネルでソーティングされた候補位置の相関値の大きさを、複数チャネルに跨って評価して各チャネルの候補位置を特定し、特定した各チャネルの候補位置の相関値を加算して閾値を求める。
【００３４】
以下、本発明の一実施の形態について、図面を用いて説明する。
【００３５】
図３は、本実施の形態に係る音声符号化装置の構成を示すブロック図である。
【００３６】
前処理部１０１は、入力音声信号に対し、ＤＣ成分を取り除くハイパスフィルタ処理や後続する符号化処理の性能改善につながるような波形整形処理やプリエンファシス処理を行い、これらの処理後の信号（Xin）をＬＰＣ分析部１０２および加算部１０５に出力する。
【００３７】
ＬＰＣ分析部１０２は、Xinを用いて線形予測分析を行い、分析結果（線形予測係数）をＬＰＣ量子化部１０３に出力する。ＬＰＣ量子化部１０３は、ＬＰＣ分析部１０２から出力された線形予測係数（ＬＰＣ）の量子化処理を行い、量子化ＬＰＣを合成フィルタ１０４に出力するとともに量子化ＬＰＣを表す符号（Ｌ）を多重化部１１４に出力する。
【００３８】
合成フィルタ１０４は、量子化ＬＰＣに基づくフィルタ係数により、後述する加算部１１１から出力される駆動音源に対してフィルタ合成を行うことにより合成信号を生成し、合成信号を加算部１０５に出力する。
【００３９】
加算部１０５は、合成信号の極性を反転させてXinに加算することにより誤差信号を算出し、誤差信号を聴覚重み付け部１１２に出力する。
【００４０】
適応音源符号帳１０６は、過去に加算部１１１によって出力された駆動音源をバッファに記憶し、パラメータ決定部１１３から出力された信号により特定される過去の駆動音源から１フレーム分のサンプルを適応音源ベクトルとして切り出して乗算部１０９に出力する。
【００４１】
ゲイン符号帳１０７は、パラメータ決定部１１３から出力された信号によって特定され
る適応音源ベクトルのゲインと固定音源ベクトルのゲインとをそれぞれ乗算部１０９と乗算部１１０とに出力する。
【００４２】
固定音源符号帳１０８は、所定形状のパルス音源ベクトルをバッファに複数記憶し、パラメータ決定部１１３から出力された信号によって特定される形状を有するパルス音源ベクトルを固定音源ベクトルとして乗算部１１０に出力する。なお、パルス音源ベクトルに拡散ベクトルを乗算して得られたものを固定音源ベクトルとして乗算部１１０に出力しても良い。
【００４３】
乗算部１０９は、ゲイン符号帳１０７から出力されたゲインを、適応音源符号帳１０６から出力された適応音源ベクトルに乗じて、加算部１１１に出力する。乗算部１１０は、ゲイン符号帳１０７から出力されたゲインを、固定音源符号帳１０８から出力された固定音源ベクトルに乗じて、加算部１１１に出力する。
【００４４】
加算部１１１は、利得乗算後の適応音源ベクトルと利得乗算後の固定音源ベクトルとをそれぞれ乗算部１０９と乗算部１１０とから入力し、これらをベクトル加算し、加算結果である駆動音源を合成フィルタ１０４および適応音源符号帳１０６に出力する。なお、適応音源符号帳１０６に入力された駆動音源は、バッファに記憶される。
【００４５】
聴覚重み付け部１１２は、加算部１０５から出力された誤差信号に対して聴覚的な重み付けをおこない符号化歪みとしてパラメータ決定部１１３に出力する。
【００４６】
パラメータ決定部１１３は、聴覚重み付け部１１２から出力された符号化歪みを最小とする適応音源ベクトル、固定音源ベクトル及び量子化利得の符号を探索し、探索された適応音源ベクトルを表す符号（Ａ）、固定音源ベクトルを表す符号（Ｆ）及び量子化利得を表す符号（Ｇ）を多重化部１１４に出力する。
【００４７】
多重化部１１４は、ＬＰＣ量子化部１０３から量子化ＬＰＣを表す符号（Ｌ）を入力し、パラメータ決定部１１３から適応音源ベクトルを表す符号（Ａ）、固定音源ベクトルを表す符号（Ｆ）および量子化利得を表す符号（Ｇ）を入力し、これらの情報を多重化して符号化情報として出力する。
【００４８】
図４は、本実施の形態に係る音声符号化装置の音源探索回路の構成を示すブロック図である。この音源探索回路は、図３に示した音声符号化装置のパラメータ決定部１１３に設けられる。
【００４９】
図４に示す音源探索回路１５０は、適応音源符号帳１０６の適応音源を探索する適応音源探索部１５１と、固定音源符号帳１０８の固定音源を探索する固定音源探索部１５２とから構成されている。
【００５０】
固定音源探索部１５２は、前処理部２０１と、探索部２０２とから構成されている。前処理部２０１は、閾値算出部２１１と、候補位置ソーティング部２１２と、カウンタクリア部２１３と、制限回数設定部２１４と、を有する。また、探索部２０２は、探索制御部２２１と、カウンタ２２２と、探索打ち切り部２２３と、を有する。
【００５１】
この固定音源探索部１５２において、前処理部２０１の閾値算出部２１１は、探索を打ち切るための閾値、すなわちこれ以上深い探索を行うか否かを判定するための相関値の閾値を算出する。なお、本実施の形態に係る閾値の算出方法の詳細については後述する。
【００５２】
前処理部２０１の候補位置ソーティング部２１２は、チャネル毎にすべての要素の値が
正の値に変換されたベクトルｙＨの要素の大きさに基づいて、要素の大きい順に、パルス候補位置を各チャネルから出力する順番を並び替える。この並び替えの順序は、閾値算出部２１１および探索部２０２の探索制御部２２１に出力される。
【００５３】
前処理部２０１のカウンタクリア部２１３は、探索部２０２のカウンタ２２２をリセットする。
【００５４】
前処理部２０１の制限回数設定部２１４は、カウンタ２２２の値を参照して探索を中断するための制限回数Ｒを定める。この制限回数Ｒは、設定された後にカウンタ２２２に出力され、カウンタ２２２は、カウンタ数が制限回数設定部２１４で設定された制限回数Ｒを超えた場合には、探索打ち切り部２２３に制御信号を送る。
【００５５】
探索部２０２の探索制御部２２１は、候補位置ソーティング部２１２にてソーティングされたパルス候補位置で探索を行うように制御を行う。また、探索制御部２２１は、探索対象のチャネルで既に探索された相関値の和に対して閾値判定を行う。そして、探索制御部２２１は、閾値以上の場合には、制御信号をカウンタ２２２に出力し、閾値未満の場合には、制御信号を探索打ち切り部２２３に出力する。
【００５６】
探索部２０２のカウンタ２２２は、閾値以上と判定された数をカウントし、このカウント数と制限回数Ｒとの比較を行い、Ｒを超えたときには、全探索を終了させる制御信号を探索打ち切り部２２３に出力する。
【００５７】
探索部２０２の探索打ち切り部２２３は、探索制御部２２１あるいはカウンタ２２２から入力された制御信号により、そのパルス位置候補における探索を打ち切り、全探索を終了させる。
【００５８】
上記構成を有する固定音源探索部１５２における固定音源符号帳の探索方法について説明する。まず、本発明の代数的符号帳のパルス探索アルゴリズムについて説明する。
【００５９】
代数的符号帳の音源は、少数の振幅１で極性（＋、−）のある複数のパルスで構成される。以後、サブフレーム長３２でパルス本数４の場合を例として説明する。代数的符号帳のパルス位置は式（５）のようにお互いが重ならないような配置をとる。したがって、極性も含めると、この代数的符号帳は（３＋１）×４で１６ビットの符号により符号化されることがわかる。
【００６０】
探索の手順を以下に示す。まず、前処理として、ベクトルｙＨとマトリクスＨＨを算出する。ｙＨはベクトルｙを逆順にしてマトリクスＨを畳み込み、更にその結果を逆順にすることにより求める。ＨＨはマトリクス同士の掛け算により求める。
【００６１】
次いで、ベクトルｙＨの要素の極性（＋−）から、事前にパルスの極性を決める。具体的には、各位置に立つパルスの極性をｙＨのその位置の値に合わせることとし、ｙＨの極性を別の配列に格納する（後にこの配列の極性を参照して極性の符号を求める）。各位置の極性を別の配列に格納したら、ｙＨの値は全て絶対値をとり正の値に変換しておく。また、その極性にあわせてＨＨの値に極性を乗ずることによって変換しておく。
【００６２】
次いで、候補位置ソーティング部２１２において、各パルスの候補位置のｙＨの値を参照して、値の大きい順にソーティングを行い、候補位置を別配列に格納する。この並べ替えは、８つの並べ替えなので少ない計算量で実現できる。
【００６３】
ベクトルｙＨのある位置の値が大きいということは、ターゲットとその位置のパルスの
合成音との相関が高いということを示しており、最適なパルスの組み合わせの一つとして選ばれる可能性が高いということが言える。したがって、ベクトルｙＨの値の大きい順に候補位置を並べ替えることにより、相関の高いパルス候補位置から探索することになり、探索を途中で打ち切ったとしても、最適な組み合わせについての探索が漏れてしまう確率を小さくすることが可能となる。その結果、少ない演算量で高い確率で最適なパルス位置の組み合わせを探索することができる。
【００６４】
このソーティングアルゴリズムについて図５を用いて説明する。なお、例としてチャネル０について示す。
【００６５】
まず、ｉを初期化し（ＳＴ３０１）、８つのパルス候補位置（０，４，８，１２，１６，２０，２４，２８）のｙＨの値を求め（ＳＴ３０３）、その値をＺ[ｉ]に書き込む（ＳＴ３０２〜ＳＴ３０４）。
【００６６】
次いで、再度ｉを初期化し（ＳＴ３０５）、Ｚｍａｘを選択されない十分に小さい値とする（ＳＴ３０７）。そして、ｊを増やしながらＺｍａｘとＺ[ｊ]とを比較し（ＳＴ３１０）、ＺｍａｘよりもＺ[ｊ]が大きければ、ＺｍａｘをＺ[ｊ]に書きかえる（ＳＴ３１２）。Ｚｍａｘの初期値は小さく設定しているので、最初の比較では必ず書きかえが起こることとなる。そして、この書きかえられたＺ[ｊ]のｊをｊｊとして保存する。これをｊが８になるまで行う（ＳＴ３０９〜ＳＴ３１２）。
【００６７】
これにより、最大のＺｍａｘとその時のインデクスｊｊが決定される。そして、Ｚ[ｊｊ]を選択されない小さい値とする（ＳＴ３１３）。これにより、一度決定したｊｊは２度と選択されないこととなり、残りの７つについて上記と同じような処理を行って、順に大きい値と、そのときのインデクスが決定されていくことになる（ＳＴ３０６，ＳＴ３０８）。同様にして、チャネル１〜３についても並べ替えを行う。
【００６８】
次に、上記ソーティングを行ったパルス候補位置で探索を行う場合について図６及び図７を用いて説明する。図６及び図７は、本実施の形態に係る代数的符号帳探索アルゴリズムのフロー図である。
【００６９】
まず、ＳＴ４０１で初期化を行い、固定音源符号帳１０８から１つ目のパルスの位置を出力してｙＨ及びＨＨから値を取り出して、それぞれ相関値ｓｙ０、パワーｓｈ０とする（ＳＴ４０３）。この計算をｉ０が８（パルス位置候補数）になるまで行う（ＳＴ４０２〜ＳＴ４０４）。
【００７０】
一つのｉ０における計算において、固定音源符号帳１０８から２つ目のパルスの位置を出力してｙＨ及びＨＨから値を取り出して相関値ｓｙ０、パワーｓｈ０にそれぞれ加算し、相関値ｓｙ１、パワーｓｈ１とする（ＳＴ４０６）。この計算をｉ１が８（パルス位置候補数）になるまで行う（ＳＴ４０５〜ＳＴ４０７）。
【００７１】
一つのｉ１における計算において、固定音源符号帳１０８から３つ目のパルスの位置を出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ１、ｓｈ１にそれぞれ加算し、相関値ｓｙ２、パワーｓｈ２とする（ＳＴ４０９）。ここで、探索を打ち切るかどうかを判断する（ＳＴ４１０）。このように本実施の形態では、第３ループが終わった段階で判定する場合について説明を行うこととする。なお、他のループ終了時に判定しても同様の効果が得られる。
【００７２】
ＳＴ４１０において、相関値ｓｙ２が閾値Ｔｈ以上であれば、カウンタ２２２のカウント数を１インクリメントする（ＳＴ４１２）。相関値ｓｙ２が閾値Ｔｈ以上であれば、相
関値が高いということであり、更に深い探索ループに進むことになる。すなわち、この計算をｉ２が８（パルス位置候補数）になるまで行う（ＳＴ４０８〜ＳＴ４１１，ＳＴ４１４）。
【００７３】
次いで、一つのｉ２における計算において、固定音源符号帳１０８から４つ目のパルスの位置を出力してｙＨ及びＨＨから値を取り出して、相関値ｓｙ２、ｓｈ２にそれぞれ加算し、相関値ｓｙ３、パワーｓｈ３とする（ＳＴ４１６）。この計算をｉ３が８（パルス位置候補数）になるまで行う（ＳＴ４１５〜ＳＴ４１８）。
【００７４】
このようにして得られた相関値ｓｙ３とパワーｓｈ３を用いて、それまでの最大値と比較し（ＳＴ４１７）、最も大きな値になる位置を探索する（ＳＴ４１９）。このアルゴリズムによって代数的符号帳の探索がなされる。
【００７５】
すなわち、上記処理においては、予め指定された深さのループまで探索した後で、それまでに加算した相関値の和が上記閾値を上回る場合は更に深いループの探索を行い、下回る場合はそれより深いループの探索を行わない。
【００７６】
また、カウンタ２２２では、カウンタ数と制限回数Ｒとを比較し（ＳＴ４１３）、カウンタ数Ｃが制限回数Ｒを超えた場合には、探索打ち切り部２２３に制御信号が送られ、探索打ち切り部２２３で探索が打ち切られる（探索中断）。一方、カウンタ数Ｃが制限回数Ｒ未満である場合には、上述のように更に深い探索ループに進む。
【００７７】
すなわち、ここでの処理においては、予め指定された所定の深さのループを探索する数を数えるカウンタを用い、カウンタ２２２により、ある深さより深いループを探索すると判定された場合の数をカウントし、カウンタの値が予め指定された制限回数を超えた場合に、探索を中断する。
【００７８】
この４重ループ（パルスが４本であるからである）を用いて、ｙＨとＨＨの値から相関値とパワーを算出し、相関が最も大きくなるパルス位置を探索する。そして各パルスの位置の符号と極性の符号（極性は前述したように、別配列に格納されている）を合わせて固定音源の符号とする。
【００７９】
次に、本実施の形態に係る閾値の算出方法の詳細について説明する。本実施の形態では、各チャネルでソーティングされた候補位置の相関値の大きさを複数チャネルに跨って評価し、各チャネルの閾値算出に用いる候補位置（以下、「基準候補位置」という）を特定し、各チャネルの基準候補位置の相関値を加算して閾値を求める。
【００８０】
本実施の形態では、各チャネルに跨って相関値が大きい順番に候補位置を探索し、探索した候補位置の数の和が予め指定した候補数になった時点で基準候補位置を特定する。また、この時、各チャネルのエントリ数を越えないようにソーティングすることが必要になる。
【００８１】
図８は、本実施の形態に係る閾値を求めるアルゴリズムを示すフロー図であり、３パルスまでの結果で４パルス目の探索を行うかどうかを判定する場合を示す。なお、図８のフロー図において、ｉ０，ｉ１，ｉ２は閾値を求める各チャネルの順位、ｃはカウンタ、Ｔｈは閾値、Ｍは閾値の強度を決める候補数、を示す。Ｍは定数で１２〜１６程度に設定される。
【００８２】
また、図９は、本実施の形態に係る閾値を求めるアルゴリズムを示すフロー図であり、２パルスまでの結果で４パルス目の探索を行うかどうかを判定する場合を示す。なお、図
９のフロー図において、ｉ０，ｉ１は閾値を求める各チャネルの順位、ｃはカウンタ、Ｔｈは閾値、Ｍは閾値の強度を決める候補数、を示す。Ｍは定数で８〜１２程度に設定される。
【００８３】
以下、本実施の形態の閾値の算出方法を、具体例を用いて説明する。例えば、ソーティングされた相関値が以下の図１０に示すものであり、候補数Ｍが「１０」であるとする。この場合、図８に示したアルゴリズムを実行することにより、候補位置の探索は、図１１の枠で囲った部分にまで進む。
【００８４】
よって、本例の場合、各チャネルの基準候補位置は、ｉ０＝５、ｉ１＝２、ｉ２＝３となり、閾値Ｔｈは、Ｔｈ＝７＋７＋６＝２０となる。
【００８５】
なお、本実施の形態では、各チャネルに跨って相関値が大きい順番に候補位置を探索し、探索した候補位置の数の積が予め指定した積になった時点で基準候補位置を特定しても良い。
【００８６】
図１２は、この方法により閾値を求めるアルゴリズムを示すフロー図である。なお、図１２のフロー図において、ｉ０，ｉ１，ｉ２は閾値を求める各チャネルの順位、ｃはカウンタ、Ｔｈは閾値、Ｎは閾値の強度を決める積、を示す。Ｎは定数で３０〜６０程度に設定される。
【００８７】
この場合において、例えば、ソーティングされた相関値が以下の上記図１０に示すものであり、積Ｎが「３５」であるとする。この場合、図１２に示したアルゴリズムを実行することにより、「６×２×３＝３６＞３５」で、候補位置の探索は、図１３の枠で囲った部分にまで進む。
【００８８】
よって、本例の場合、各チャネルの基準候補位置は、ｉ０＝６、ｉ１＝２、ｉ２＝３となり、閾値Ｔｈは、Ｔｈ＝５＋７＋６＝１８となる。
【００８９】
このように、本実施の形態では、各チャネルの候補位置に配置されたパルスの合成音と量子化ターゲットとの相関値に基づいて探索中断の判断のための閾値を設け、ループ探索の途中で相関値の和が閾値を下回る場合にそれ以上深いループ探索を行わない。しかも、本実施の形態では、各チャネルでソーティングされた候補位置の相関値の大きさを、複数チャネルに跨って評価して各チャネルの候補位置を特定し、特定した各チャネルの候補位置の相関値を加算して閾値を求める。
【００９０】
これにより、相関値が小さい場合の探索処理にかかる計算量を節約することができ、符号化性能を落とさずに平均的計算量を削減することができる。
【００９１】
なお、本発明は上記実施の形態に限定されず、種々変更して実施することが可能である。例えば、上記実施の形態では、第３ループの深さで探索の打ち切りをするか否かを判断した場合について説明しているが、本発明では、他のループの深さで探索の打ち切りをするか否かを判断するようにしても良い。
【００９２】
また、上記実施の形態に係るパルス候補位置のソーティングや探索は、音声符号化装置において行うように説明しているが、このソーティングや探索をソフトウェアとして構成しても良い。例えば、上記ソーティングや探索のプログラムをＲＯＭに格納し、そのプログラムにしたがってＣＰＵの指示により動作させるように構成しても良い。また、ソーティングや探索のプログラムをコンピュータで読み取り可能な記憶媒体に格納し、この記憶媒体のソーティングや探索のプログラムをコンピュータのＲＡＭに記録して、このプログ
ラムにしたがって動作させるようにしても良い。このような場合においても、上記実施の形態と同様の作用、効果を呈する。
【００９３】
また、本実施の形態では、固定音源符号帳として代数的符号帳を用いた場合の例を示したが、本発明はこれに限られず、マルチパルス符号帳や、固定波形がＲＯＭに書かれた固定符号帳でも有効である。この場合、各チャネルのそれぞれの番号が固定波形のベクトルに対応する。
【００９４】
また、本実施の形態では、ＣＥＬＰに対して用いる場合について説明したが、本発明はこれに限られず、音源の符号帳が存在する他の符号化にも有効である。これは、本発明の所在は固定音源符号帳ベクトル内のものであり、適応音源符号帳の有無や、スペクトル包絡の分析方法に依存しないからである。
【００９５】
また、本実施の形態では、固定音源符号帳に用いた例を示したが、本発明はこれに限られず、スペクトル符号化の多段ベクトル量子化、画像符号化におけるベクトルの多段ベクトル量子化にも有効である。これは、本発明が、多重ループの探索における候補数の削減に貢献するものであり、符号化対象には限定されないからである。
【００９６】
また、本発明に係る信号は、音声信号だけでなく、オーディオ信号でも良い。また、入力信号の代わりに、ＬＰＣ予測残差信号に対して本発明を適用する構成であっても良い。
【００９７】
また、本発明に係る符号化装置は、移動体通信システムにおける通信端末装置および基地局装置に搭載することが可能であり、これにより上記と同様の作用効果を有する通信端末装置、基地局装置、および移動体通信システムを提供することができる。
【００９８】
また、ここでは、本発明をハードウェアで構成する場合を例にとって説明したが、本発明をソフトウェアで実現することも可能である。例えば、本発明に係るアルゴリズムをプログラミング言語によって記述し、このプログラムをメモリに記憶しておいて情報処理手段によって実行させることにより、本発明に係る符号化装置と同様の機能を実現することができる。
【００９９】
また、上記実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されても良いし、一部または全てを含むように１チップ化されても良い。
【０１００】
また、ここではＬＳＩとしたが、集積度の違いによって、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩ等と呼称されることもある。
【０１０１】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現しても良い。ＬＳＩ製造後に、プログラム化することが可能なＦＰＧＡ（Field Programmable Gate Array）や、ＬＳＩ内部の回路セルの接続もしくは設定を再構成可能なリコンフィギュラブル・プロセッサを利用しても良い。
【０１０２】
さらに、半導体技術の進歩または派生する別技術により、ＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行っても良い。バイオ技術の適用等が可能性としてあり得る。
【０１０３】
２００７年３月２日出願の特願２００７−０５３５０１の日本出願に含まれる明細書、図面および要約書の開示内容は、すべて本願に援用される。
【産業上の利用可能性】
【０１０４】
本発明は、音声信号やオーディオ信号を符号化する符号化装置等に用いるに好適である。
【図面の簡単な説明】
【０１０５】
【図１】従来の代数的符号帳探索アルゴリズムのフロー図
【図２】従来の代数的符号帳探索アルゴリズムのフロー図
【図３】本発明の一実施の形態に係る音声符号化装置の構成を示すブロック図
【図４】本発明の一実施の形態に係る音声符号化装置の音源探索回路の構成を示すブロック図
【図５】本発明の一実施の形態に係るソーティングのアルゴリズムを示すフロー図
【図６】本発明の一実施の形態に係る代数的符号帳探索アルゴリズムのフロー図
【図７】本発明の一実施の形態に係る代数的符号帳探索アルゴリズムのフロー図
【図８】本発明の一実施の形態に係る閾値を求めるアルゴリズムを示すフロー図
【図９】本発明の一実施の形態に係る閾値を求めるアルゴリズムを示すフロー図
【図１０】本発明の一実施の形態に係るソーティングされた相関値の例を示す図
【図１１】本発明の一実施の形態に係るソーティングされた相関値の例を示す図
【図１２】本発明の一実施の形態に係る閾値を求めるアルゴリズムを示すフロー図
【図１３】本発明の一実施の形態に係るソーティングされた相関値の例を示す図
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】全図
【補正方法】変更
【補正の内容】
【図１】
2008108077.tif 000026
【図２】
2008108077.tif 000027
【図３】
2008108077.tif 000028
【図４】
2008108077.tif 000029
【図５】
2008108077.tif 000030
【図６】
2008108077.tif 000031
【図７】
2008108077.tif 000032
【図８】
2008108077.tif 000033
【図９】
2008108077.tif 000034
【図１０】
2008108077.tif 000035
【図１１】
2008108077.tif 000036
【図１２】
2008108077.tif 000037
【図１３】
2008108077.tif 000038
【国際調査報告】
2008108077.tif 000039
2008108077.tif 000040