JP5145733B2 - 音声信号処理装置および音声信号処理方法ならびにプログラム - Google Patents

Description

本発明は、音声信号処理装置および音声信号処理方法ならびにプログラムに係り、特に、音声信号に対して倍音を付加する技術に係る。

従来、信号音に対し信号音の倍音を付加して音質を向上する技術が知られている。例えば、特許文献１には、音源からの信号音に含まれる周波数成分を利用し、信号音を乗算器あるいは非直線特性素子により高調波成分を増加させ、帯域濾波器に接続し、源信号音の人工倍音を作り出し、この人工倍音を源信号音に混合するように構成した音質改善装置が開示されている。このような装置によれば、電気信号経路で失われがちな倍音成分を増強することで、通過する信号音の明瞭性を上げることができる。

なお、関連する技術として、基本波の倍音を音響信号に付加して、再生音における低音感を向上させる音響信号処理装置が、特許文献２および３に記載されている。

特開昭５７−５４９９３号公報 特開２００６−２２２８６７号公報 特開２００６−３２４７８６号公報

携帯電話端末装置等の受話信号にあって、信号帯域が狭い場合などでは、音声の明瞭性が損なわれることがある。このような場合、信号音に対し倍音を付加することによって、音声品質を向上させる技術が一般的に用いられている。しかしながら、従来の倍音付加技術では、（１）常に生成される倍音が付加されるため、音質の変化が大きい、（２）音声信号を認識しやすくするために音声信号の一部を倍音で強調すると、それ以外の部分も強調されてしまうため、倍音が強すぎて耳障りな音となる虞がある、といった弊害があった。

したがって、本発明の目的は、音声信号をより聞き取りやすい音質に加工する音声信号処理装置および音声信号処理方法ならびにプログラムを提供することにある。

本発明者は、携帯電話端末装置等において、聴取する音声信号の明瞭性が損なわれる大きな要因として、音声信号の子音部に着目した。すなわち、子音部は、音声信号のアタック部（立ち上り部の短期間）に存在し、子音部を明瞭に聞き取れるか否かで聴取する音声信号の明瞭性が概ね決定されてしまうことに注目した。そして、入力信号の音量変化に応じて倍音を付与する量を調整すると共に、必要な部分（アタック部分）にのみ選択的に倍音を付加することで、より聞き取りやすい音質となることを見出し、本発明を創案するに至った。

本発明の１つのアスペクトに係る音声信号処理装置は、入力信号の絶対値に対する時間方向のエンベロープにおける立ち上り部分を検出するアタック検出部と、入力信号から倍音成分を生成する倍音生成部と、検出した立ち上り部分の信号の大きさに応じて倍音成分を増幅する増幅部と、増幅した倍音成分を入力信号に加算して出力する加算部と、を備える。

本発明の音声信号処理装置において、アタック検出部は、入力信号の絶対値に対する時間方向のエンベロープにおいて立ち上りに第１の時定数で追従して立ち下りに該第１の時定数より大きな第２の時定数で下降する第１の信号と、該立ち上りに該第１の時定数より大きく該第２の時定数より小さい第３の時定数で追従して該立ち下りに該第２の時定数で下降する第２の信号との差分信号を求めるアタック検出器と、差分信号における高域成分を除去して立ち上り部分の信号成分とするエンベロープ生成器と、を備えるようにしてもよい。

本発明の音声信号処理装置において、倍音生成部は、入力信号に対して高調波成分を生成し、所定の周波数の帯域制限を行って出力するようにしてもよい。

本発明の他のアスペクトに係る音声信号処理方法は、入力信号の絶対値に対する時間方向のエンベロープにおける立ち上り部分を検出するステップと、入力信号から倍音成分を生成するステップと、検出した立ち上り部分の信号の大きさに応じて倍音成分を増幅するステップと、増幅した倍音成分を入力信号に加算して出力するステップと、を含む。

本発明のさらに他のアスペクトに係るプログラムは、音声信号処理装置を構成するコンピュータに、入力信号の絶対値に対する時間方向のエンベロープにおける立ち上り部分を検出する処理と、入力信号から倍音成分を生成する処理と、検出した立ち上り部分の信号の大きさに応じて倍音成分を増幅する処理と、増幅した倍音成分を入力信号に加算して出力する処理と、を実行させる。

本発明によれば、入力信号のパワー増加区間に含まれる期間の信号成分の大きさに応じて倍音成分を付加するようにする。この期間の信号成分は、音声信号の明瞭性において重要な役割を果たす子音部を含む場合があるため、この期間の信号成分のみを強調することでより聞き取りやすい音声信号を作り出すことができる。

本発明の実施形態に係る音声信号処理装置は、入力信号のパワー増加区間に含まれる期間（アタック部分）の信号成分（アタック成分）を検出するアタック検出部（図２の２１、２２）と、入力信号から倍音成分を生成する倍音生成部（図２の２６）と、検出したアタック成分の大きさに応じて倍音成分を増幅する増幅部（図２の２７）と、増幅した倍音成分を入力信号に加算して出力する加算部（図２の２９）と、を備える。アタック検出部は、入力信号の絶対値に対する時間方向のエンベロープにおいて立ち上りに第１の時定数で追従して立ち下りに該第１の時定数より大きな第２の時定数で下降する第１の信号と、該立ち上りに該第１の時定数より大きく該第２の時定数より小さい第３の時定数で追従して該立ち下りに該第２の時定数で下降する第２の信号との差分信号を求めるアタック検出器（図２の２１）と、差分信号における高域成分を除去してアタック成分の信号とするエンベロープ生成器（図２の２２）と、を備え、入力信号の絶対値に対する時間方向のエンベロープにおける立ち上り部分をアタック成分として検出する。倍音生成部は、入力信号に対して高調波成分を生成し、所定の周波数の帯域制限を行って出力する。

このような構成の音声信号処理装置によれば、音声信号の明瞭性において重要な役割を果たす子音部を含む必要な部分（アタック部分）にのみ選択的に倍音を付加することで、より聞き取りやすい音質の音声信号が得られる。以下、実施例に即し、図面を参照して詳細に説明する。

図１は、本発明の実施例に係る携帯電話端末装置の構成を示すブロック図である。図１において、携帯電話端末装置１０は、マイク１１、スピーカ１２、音声入出力回路１３、符号器（エンコーダ）１４、送信無線回路１５、アンテナ１６、受信無線回路１７、復号器（デコーダ）１８、信号処理回路２０を備える。マイク１１から入力される送話信号は、音声入出力回路１３によって増幅され、符号器（エンコーダ）１４で符号化され、送信無線回路１５によって無線信号に変換され、アンテナ１６から送信される。

一方、アンテナ１６によって受信された無線信号は、受信無線回路１７によって受信信号に変換され、さらに復号器（デコーダ）１８によって復号され、受信音声信号として、信号処理回路２０に入力される。信号処理回路２０は、受信音声信号に対して入力信号のパワー増加区間に含まれる期間の信号成分（アタック成分）を検出し、アタック部分に倍音を付加する。倍音を付加された受信音声信号は、音声入出力回路１３を介して増幅されスピーカ１２から再生され、携帯電話端末装置１０の使用者によって聴取される。

図２は、本発明の実施例に係る信号処理回路の構成を示すブロック図である。図２において、信号処理回路２０は、大きくアタック検出セクションと倍音付加セクションとによって構成される。

まず、アタック検出セクションについて説明する。アタック検出セクションには、アタック検出器２１、エンベロープ生成器２２、遅延回路２３が含まれる。

アタック検出器２１は、信号処理回路２０の入力信号のアタック成分を検出する。アタック検出器２１の詳細については後述する。エンベロープ生成器２２は、アタック検出器２１で得られたアタック成分の大きさを計算する。エンベロープ生成器２２の詳細については後述する。アタック検出器２１とエンベロープ生成器２２とによって処理された信号は、アタック成分の大きさを示す信号として後段の倍音付加セクションに送られる。

遅延回路２３は、信号処理回路２０の入力信号を一定時間遅らせる回路であって、一般的な信号遅延回路で構成される。アタック成分の大きさを示す信号（アタック検出器２１とエンベロープ生成器２２で生成）は、信号処理回路２０の入力信号に対してアタック検出器２１とエンベロープ生成器２２との処理による遅延を生じる。そこで、この遅延と入力信号の位相を合わせるために、遅延回路２３は、入力信号を遅延させて倍音付加セクションに送る。

次に、倍音付加セクションについて説明する。倍音付加セクションには、倍音生成器２６、ＶＣＡ（Ｖｏｌｔａｇｅ Ｃｏｎｔｒｏｌｌｅｄ Ａｍｐｌｉｆｉｅｒ）２７、遅延回路２４、２８、加算器２９が含まれる。

倍音生成器２６は、アタック検出セクションから出力された（遅延回路２３で遅延された）入力信号から、倍音成分を合成する。詳細については後述する。ＶＣＡ（Ｖｏｌｔａｇｅ Ｃｏｎｔｒｏｌｌｅｄ Ａｍｐｌｉｆｉｅｒ）２７は、倍音生成器２６の出力信号を制御信号の大きさ（一般的には電圧）に応じて増幅する回路である。ここで制御信号としては、エンベロープ生成器２２から送出されたアタック成分の大きさを示す信号を遅延回路２４で遅延した信号が使用される。

遅延回路２４、２８の構成は、遅延回路２３と同様、一般的な信号遅延回路である。倍音生成器２６は、遅延回路２３の出力に対して処理遅延を生じる。そこで遅延回路２４は、倍音生成器２６の出力信号とエンベロープ生成器２２の出力信号との位相を合わせるためにエンベロープ生成器２２の出力信号を遅延させてＶＣＡ２７に与える。

加算器２９は、ＶＣＡ２７の出力信号と、遅延回路２３の出力に対して遅延回路２８によって適切に遅延が加えられた信号と、を加算して信号処理回路２０の出力信号として出力する。

次に、アタック検出器２１の詳細について説明する。図３は、アタック検出器２１の構成例を示すブロック図である。アタック検出器２１は、第１のセクション（ｓｈａｐｅ０）と第２のセクション（ｓｈａｐｅ１）と第３のセクション（０−１）とで構成される。以下、第１のセクションにおいて出力される信号をｓｈａｐｅ０、第２のセクションにおいて出力される信号をｓｈａｐｅ１と称する。ｓｈａｐｅ０は、信号処理回路２０の入力信号のピークと、それに続くゆるやかな減衰成分からなる信号であり、ｓｈａｐｅ１は、ｓｈａｐｅ０のピークをなまらせた信号である。

まず、第１のセクションにおけるｓｈａｐｅ０の生成について説明する。第１のセクションには、絶対値回路３１、追従回路３２、大信号選択回路３３、１サンプル遅延回路３４が含まれる。

絶対値回路３１は、アタック検出器２１における入力信号の振幅の絶対値を演算して追従回路３２および大信号選択回路３３に出力する。

追従回路３２は、絶対値回路３１の出力に対し、直前の大信号選択回路３３の出力信号を一定の比率だけ近づけた値を計算する。この計算は、１サンプル遅延回路３４による１サンプル分の遅延（直前の出力信号を生成する）と追従計算により行われる。ここで追従計算は、以下のような式で表現される。これはデジタル信号処理を行う場合の式となる。
Ｙ［ｉ］＝ｙ［ｉ−１］＋（ｘ［ｉ］−ｙ［ｉ−１］）＊ａ
ただし、ｘ［ｉ］は時刻ｉの入力、ｙ［ｉ］は時刻ｉの出力、Ｙ［ｉ］は追従信号、ａは時定数である。なお、ａは、以下の式に基づいて決定される。
ａ＝１−（１／２）＾（１／（Ｔ＊Ｆｓ））
ただし、「＾」はべき乗記号、Ｆｓは信号のサンプリングレートである。また、Ｔは、入力のピーク成分が入力されてからｓｈａｐｅ０がその１／２になるまでの時間を表し、通常５０〜５００ｍｓ前後の値を用いる。

大信号選択回路３３は、絶対値回路３１の出力信号と追従回路３２の出力である追従信号との大きさを比較し、大きい方の信号をｓｈａｐｅ０として採用し出力する。これにより、ピークがあった場合は、大信号選択回路３３によって入力信号（の絶対値）が選択され、ピーク以外の箇所では追従計算により得られたゆるやかに減衰する信号が得られる。

以上をまとめると、ｓｈａｐｅ０は、以下の式で生成される。
ｙ［ｉ］＝ｍａｘ｛ａｂｓ（ｘ［ｉ］），ｙ［ｉ−１］＋（ｘ［ｉ］−ｙ［ｉ−１］）＊ａ｝
ただし、ｍａｘは最大値選択を行う関数、ａｂｓは絶対値演算を行う関数である。

次に、第２のセクションにおけるｓｈａｐｅ１の生成について説明する。第２のセクションには、小信号選択回路３５、追従回路３６、１サンプル遅延回路３７が含まれる。追従回路３６、１サンプル遅延回路３７は、それぞれ追従回路３２、１サンプル遅延回路３４と同等の回路である。ただし、追従回路３６における時定数は、追従回路３２とは異なり、ピークが出現してから、ｓｈａｐｅ１がその１／２に到達するまでの時間をＴとして計算する。計算式は、ｓｈａｐｅ０の場合と同様である。Ｔとしてはｓｈａｐｅ０の算出におけるＴより短い時間が採用され、通常１ｍｓ〜１００ｍｓ程度である。

小信号選択回路３５は、大信号選択回路３３の出力信号（ｓｈａｐｅ０）と追従回路３６の出力信号との大きさを比較し、小さい方の信号を出力する。

第２のセクションでは、ｓｈａｐｅ１として小さい方の信号を採用するので、瞬間的なピークが発生している区間では追従信号が採用され、ｓｈａｐｅ１は、なだらかな立ち上りとなる。ｓｈａｐｅ０がゆるやかに減衰している区間では、ｓｈａｐｅ１も追いつくが、減衰している間は、入力信号が採用されるためｓｈａｐｅ０と同様にゆるやかに減衰する信号となる。

次に、第３のセクション（０−１）について説明する。第３のセクションにおいて、減算器３８は、ｓｈａｐｅ０からｓｈａｐｅ１を減じた信号（０−１信号）を算出し、アタック検出器２１の出力信号とする。ｓｈａｐｅ１の説明で述べたように、瞬間的なピークが発生している区間では、ｓｈａｐｅ１がｓｈａｐｅ０に追いつかないため大きな差分が発生する。逆にピーク以外の部分では、ｓｈａｐｅ１は、ｓｈａｐｅ０と同一の信号となるため差を生じない。これにより、ピーク部分のみで値を持つような信号がアタック検出器２１から出力される。

次に、エンベロープ生成器２２の詳細について説明する。図４は、エンベロープ生成器２２の構成例を示すブロック図である。図４において、エンベロープ生成器２２は、絶対値回路４１、追従回路４２、大信号選択回路４３、１サンプル遅延回路４４、小信号選択回路４５、追従回路４６、１サンプル遅延回路４７を、備え、それぞれ、図３における絶対値回路３１、追従回路３２、大信号選択回路３３、１サンプル遅延回路３４、小信号選択回路３５、追従回路３６、１サンプル遅延回路３７に相当する。すなわち、エンベロープ生成器２２の構成は、アタック検出器２１における第１および第２のセクションと同じ構成であり、ｓｈａｐｅ１をそのまま出力する構成となっている。ただし、時定数はアタック検出器２１と異なっており、エンベロープ生成器２２におけるｓｈａｐｅ０／１の算出に用いるＴは、共に５ｍｓ前後の短い値を用いる。

次に、倍音生成器２６の詳細について説明する。図５は、倍音生成器２６の構成例を示すブロック図である。図５において、倍音生成器２６は、ＨＰＦ（Ｈｉｇｈ−ｐａｓｓ ｆｉｌｔｅｒ、高域通過フィルタ）回路５１、歪み生成回路５２、ＨＰＦ回路５３を備える。倍音生成器２６の入力信号は、まずＨＰＦ回路５１によって高域成分のみが抽出される。これは対象とする信号を選択するためのものであるので、帯域通過型や低域通過型であってもよい。

歪み生成回路５２は、ＨＰＦ回路５１の出力に対して歪み生成処理を行う。歪み生成処理では、過剰な増幅による歪み生成や位相変調等が使用され、入力信号に含まれる周波数の整数倍の信号成分を多く含んだ信号を生成する。さらに、ＨＰＦ回路５３は、歪み生成回路５２の出力信号の不要な歪み成分を除去し、元信号の高域成分のみを抽出した信号を得る。

なお、図５に示す構成は、一例であって、一般的にエンハンサ、エキサイタなどの名称で呼ばれる各種の倍音生成回路を用いることができる。また、場合によっては、倍音生成器２６自体は、何の処理をしなくてもよい（この場合、アタック成分は、単に音量が増強されるだけになる）。

以上説明した信号処理回路２０は、内蔵するコンピュータにプログラムを実行させて、各機能を実現するようにしてもよい。

次に、信号処理回路２０の動作について、実際の信号の例を用いて説明する。図６は、実際の音声信号の波形例を示す図である。信号振幅を、「−１〜＋１」の範囲に補正して表示している。また、図７は、図６に示す音声信号に対しアタック検出器２１で生成されたｓｈａｐｅ０とｓｈａｐｅ１の波形を示す図である。さらに、図８は、図７における最初の立ち上り部分の波形を拡大して示す図である。

図７、８を参照すると、信号処理回路２０の構成で説明したように、ｓｈａｐｅ０が入力信号のピークを捕捉し、その後なだらかに減衰するような信号となっている。減衰中の信号より大きな信号が再び入力されると、その信号レベルまで回復する。ｓｈａｐｅ１は、ｓｈａｐｅ０の立ち上り（ピーク検出部分）に追従するように変化する。ｓｈａｐｅ０に追いついた後は、ｓｈａｐｅ０の方が小さくなるため、ｓｈａｐｅ１は、ｓｈａｐｅ０と同一の出力となる。

図９は、アタック検出器２１によって生成された各信号（ｓｈａｐｅ０、ｓｈａｐｅ１、０−１信号）と、エンベロープ生成器２２の出力信号（ＬＰＦ＋ｓｈｉｆｔ）との波形を示す図である。エンベロープ生成器２２の出力については、遅延回路２３で遅延させる分、信号をずらして表示している。また、「０−１信号」および「ＬＰＦ＋ｓｈｉｆｔ」については見やすくなるように、それぞれ−０．２、−０．４を加えて表示している。図９に示すように、アタック検出器２１によってピーク成分の大きさのみが検出され、さらにエンベロープ生成器２２によって、なだらかに変化する信号として整形される。

図１０は、元の音声信号に対し、エンベロープ生成器２２の出力信号を重ねて表示した図である。エンベロープ生成器２２の出力が、元の音声信号の立ち上り部分（アタック成分）を検出していることがわかる。

以上のようにして生成されたアタック成分の大きさを示す信号の量に基づいて、倍音付加セクションで生成された倍音の加算量が調整される。これにより、信号中のアタック成分が存在する期間においてのみ倍音の量が増える。アタック成分の期間には音声信号の子音部分がふくまれる場合が多いので、子音部分が強調されることで音声信号の明瞭性が増すこととなる。

以上説明したような音声信号処理装置を携帯電話端末装置における受信側の信号処理回路に用いることで、受信音声信号の明瞭性が向上する。したがって、携帯電話端末装置の使用者にとって受信音声信号がより認識しやすい音として聞こえるようになる。

以上本発明を上記実施例に即して説明したが、本発明は、上記実施例にのみ限定されるものではなく、本願特許請求の範囲の各請求項の発明の範囲内で当業者であればなし得るであろう各種変形、修正を含むことは勿論である。

本発明の実施例に係る携帯電話端末装置の構成を示すブロック図である。 本発明の実施例に係る信号処理回路の構成を示すブロック図である。 アタック検出器の構成例を示すブロック図である。 エンベロープ生成器の構成例を示すブロック図である。 倍音生成器の構成例を示すブロック図である。 実際の音声信号の波形例を示す図である。 図６に示す音声信号に対しアタック検出器で生成されたｓｈａｐｅ０とｓｈａｐｅ１の信号の波形を示す図である。 図７における最初の立ち上り部の波形を拡大して示す図である。 アタック検出器によって生成された各信号と、エンベロープ生成器の出力信号との波形を示す図である。 元の音声信号に対し、エンベロープ生成器の出力信号を重ねて表示した図である。

符号の説明

１０ 携帯電話端末装置
１１ マイク
１２ スピーカ
１３ 音声入出力回路
１４ 符号器（エンコーダ）
１５ 送信無線回路
１６ アンテナ
１７ 受信無線回路
１８ 復号器（デコーダ）
２０ 信号処理回路
２１ アタック検出器
２２ エンベロープ生成器
２３、２４、２８ 遅延回路
２６ 倍音生成器
２７ ＶＣＡ（Ｖｏｌｔａｇｅ Ｃｏｎｔｒｏｌｌｅｄ Ａｍｐｌｉｆｉｅｒ）
２９ 加算器
３１、４１ 絶対値回路
３２、３６、４２、４６ 追従回路
３３、４３ 大信号選択回路
３４、３７、４４、４７ １サンプル遅延回路
３５、４５ 小信号選択回路
３８ 減算器
５１、５３ ＨＰＦ回路
５２ 歪み生成回路

Claims (10)

1. 入力信号の絶対値に対する時間方向のエンベロープにおける立ち上り部分を検出するアタック検出部と、
   前記入力信号から倍音成分を生成する倍音生成部と、
   検出した前記立ち上り部分の信号の大きさに応じて前記倍音成分を増幅する増幅部と、
   増幅した前記倍音成分を前記入力信号に加算して出力する加算部と、
   を備えることを特徴とする音声信号処理装置。
2. 前記アタック検出部は、
   前記入力信号の絶対値に対する時間方向のエンベロープにおいて立ち上りに第１の時定数で追従して立ち下りに該第１の時定数より大きな第２の時定数で下降する第１の信号と、該立ち上りに該第１の時定数より大きく該第２の時定数より小さい第３の時定数で追従して該立ち下りに該第２の時定数で下降する第２の信号との差分信号を求めるアタック検出器と、
   前記差分信号における高域成分を除去して前記立ち上り部分の信号成分とするエンベロープ生成器と、
   を備えることを特徴とする請求項１記載の音声信号処理装置。
3. 前記倍音生成部は、前記入力信号に対して高調波成分を生成し、所定の周波数の帯域制限を行って出力することを特徴とする請求項１記載の音声信号処理装置。
4. 受話音声信号を前記入力信号として処理する請求項１乃至３のいずれか一に記載の音声信号処理装置を備える携帯電話端末装置。
5. 入力信号の絶対値に対する時間方向のエンベロープにおける立ち上り部分を検出するステップと、
   前記入力信号から倍音成分を生成するステップと、
   検出した前記立ち上り部分の信号の大きさに応じて前記倍音成分を増幅するステップと、
   増幅した前記倍音成分を前記入力信号に加算して出力するステップと、
   を含むことを特徴とする音声信号処理方法。
6. 前記立ち上り部分を検出するステップにおいて、
   前記入力信号の絶対値に対する時間方向のエンベロープにおいて立ち上りに第１の時定数で追従して立ち下りに該第１の時定数より大きな第２の時定数で下降する第１の信号と、該立ち上りに該第１の時定数より大きく該第２の時定数より小さい第３の時定数で追従して該立ち下りに該第２の時定数で下降する第２の信号との差分信号を求めるステップと、
   前記差分信号における高域成分を除去して前記立ち上り部分の信号成分とするステップと、
   を含むことを特徴とする請求項５記載の音声信号処理方法。
7. 前記倍音成分を生成するステップにおいて、前記入力信号に対して高調波成分を生成し、所定の周波数の帯域制限を行って出力することを特徴とする請求項５記載の音声信号処理方法。
8. 音声信号処理装置を構成するコンピュータに、
   入力信号の絶対値に対する時間方向のエンベロープにおける立ち上り部分を検出する処理と、
   前記入力信号から倍音成分を生成する処理と、
   検出した前記立ち上り部分の信号の大きさに応じて前記倍音成分を増幅する処理と、
   増幅した前記倍音成分を前記入力信号に加算して出力する処理と、
   を実行させるプログラム。
9. 前記立ち上り部分を検出する処理において、
   前記入力信号の絶対値に対する時間方向のエンベロープにおいて立ち上りに第１の時定数で追従して立ち下りに該第１の時定数より大きな第２の時定数で下降する第１の信号と、該立ち上りに該第１の時定数より大きく該第２の時定数より小さい第３の時定数で追従して該立ち下りに該第２の時定数で下降する第２の信号との差分信号を求める処理と、
   前記差分信号における高域成分を除去して前記立ち上り部分の信号成分とする処理と、
   を実行させる請求項８記載のプログラム。
10. 前記倍音成分を生成する処理において、前記入力信号に対して高調波成分を生成し、所定の周波数の帯域制限を行って出力する処理を実行させる請求項８記載のプログラム。

Images