(57)【要約】【課題】受信電波に雑音を含んでいる場合でも、不法無線局のどの無線機から発信されているかまで特定できるようにする。【解決手段】この発明の無線機同定装置は、受信した電波を用いてその発信元の無線機を同定する装置であり、受信電波からトリガ位置を決定するトリガ位置決定手段１００Ｂと、受信電波の読み取りをトリガ位置で開始し、受信電波の立上がり波形を捕捉する立上がり波形捕捉手段１１，１５と、立上がり波形からスペクトログラム・パターンを演算するスペクトログラム・パターン演算手段２４と、無線機のメーカおよびモデル毎に予め求めたスペクトログラム・パターンを基準スペクトログラム・パターンとして記憶する基準情報記憶手段２２と、今回の受信電波の立上がり波形から演算したスペクトログラム・パターンと、基準スペクトログラム・パターンとを対比することで、今回の受信電波の発信元である無線機のモデルを同定する同定手段６，２４，２２と、を具備することを特徴としている。【選択図】  図１
2004214817.tif 000002
【特許請求の範囲】
【請求項１】
受信した電波を用いてその発信元の無線機を同定する無線機同定装置において、
上記受信電波からトリガ位置を決定するトリガ位置決定手段と、
上記受信電波の読み取りをトリガ位置で開始し、受信電波の立上がり波形を捕捉する立上がり波形捕捉手段と、
上記立上がり波形からスペクトログラム・パターンを演算するスペクトログラム・パターン演算手段と、
上記無線機のメーカおよびモデル毎に予め求めたスペクトログラム・パターンを基準スペクトログラム・パターンとして記憶する基準情報記憶手段と、
今回の受信電波の立上がり波形から演算したスペクトログラム・パターンと、基準スペクトログラム・パターンとを対比することで、今回の受信電波の発信元である無線機のモデルを同定する同定手段と、
を具備することを特徴とする無線機同定装置。
【請求項２】
上記トリガ位置決定手段は、受信電波から変換して得られたスペクトルが予め設定したレベル以上となる最初のポイントをトリガ位置とする、請求項１に記載の無線機同定装置。
【請求項３】
上記立上がり波形捕捉手段は、受信電波を中間周波数帯信号に変換し、その後、所定の帯域幅に制限することで、立上がり波形の捕捉処理を行う、請求項１に記載の無線機同定装置。
【請求項４】
上記スペクトログラム・パターンは、立上がり波形から演算されたスペクトログラムの信号レベルをさらに等高線で結んで得られる無線機に固有のパターンである、請求項１に記載の無線機同定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、不法無線局が送信する電波を捕らえて、発信元の無線機のモデル（機種）および個体の同定を行う無線機同定装置に関する。
【０００２】
【従来の技術】
不法無線局は、重要な無線通信や一般業務用無線へ混信・妨害の原因となったり違法に出力を強くして道路沿いのテレビやラジオなどに妨害を与えている。近年、不法無線局は増加をたどり根絶へ向け適切な対応が強く求められている。そして、不法無線局を取り締まる方法としては、従来から方向探知器により不法電波の発射源を求める方法が一般的に知られている。
【０００３】
ところが、この方向探知器により得られるのは電波の発射場所であり無線機自体でない。例えば、そこに同じ周波数帯の無線機が複数台あった場合、どの無線機から電波を発射したのかを特定することはできなかった。そこで、受信した電波からその発信元である無線機を特定する手法が求められており、例えば無線機個々から発射される電波をデータベース化し、不法電波を受信したときに、指紋を照合するように、予め保有する基準データと照合し、無線機を特定できるようにした無線機同定装置が開発されつつある（例えば、特許文献１参照）。
【０００４】
【特許文献１】
特開２００２−２６８２６号公報
【０００５】
【発明が解決しようとする課題】
しかし、上記特許文献１で開示されている無線機同定装置の場合、受信電波に雑音が入っていると、受信電波の立上がり時点を的確に捕捉できない場合があり、そのようなときには基準データとの照合を精度良く行うことができず、不法電波を発信する無線機を特定できないこともあった。
【０００６】
この発明は上記に鑑み提案されたもので、受信電波に雑音を含んでいる場合でも、不法無線局のどの無線機から発信されているかまで特定することができる無線機同定装置を提供することを目的とする。
【０００７】
【課題を解決するための手段】
上記目的を達成するために、請求項１に記載の発明は、受信した電波を用いてその発信元の無線機を同定する無線機同定装置において、上記受信電波からトリガ位置を決定するトリガ位置決定手段と、上記受信電波の読み取りをトリガ位置で開始し、受信電波の立上がり波形を捕捉する立上がり波形捕捉手段と、上記立上がり波形からスペクトログラム・パターンを演算するスペクトログラム・パターン演算手段と、上記無線機のメーカおよびモデル毎に予め求めたスペクトログラム・パターンを基準スペクトログラム・パターンとして記憶する基準情報記憶手段と、今回の受信電波の立上がり波形から演算したスペクトログラム・パターンと、基準スペクトログラム・パターンとを対比することで、今回の受信電波の発信元である無線機のモデルを同定する同定手段と、を具備することを特徴としている。
【０００８】
また、請求項２に記載の発明は、上記した請求項１に記載の発明の構成に加えて、上記トリガ位置決定手段は、受信電波から変換して得られたスペクトルが予め設定したレベル以上となる最初のポイントをトリガ位置とする、ことを特徴としている。
【０００９】
また、請求項３に記載の発明は、上記した請求項１に記載の発明の構成に加えて、上記立上がり波形捕捉手段は、受信電波を中間周波数帯信号に変換し、その後、所定の帯域幅に制限することで、立上がり波形の捕捉処理を行う、ことを特徴としている。
【００１０】
さらに、請求項４に記載の発明は、上記した請求項１に記載の発明の構成に加えて、上記スペクトログラム・パターンは、立上がり波形から演算されたスペクトログラムの信号レベルをさらに等高線で結んで得られる無線機に固有のパターンである、ことを特徴としている。
【００１１】
【発明の実施の形態】
以下にこの発明の実施の形態を図面に基づいて詳細に説明する。
【００１２】
図１は本発明の無線機同定装置の構成図である。図において、本発明の無線機同定装置は、装置メイン部１００Ａとトリガ位置決定部１００Ｂとから構成されている。装置メイン部１００Ａのアンテナ２は無線機１からの電波を受信し、その立上がり状態を捕捉する。なお、以下の説明では、無線機１からは不法電波が送信されているものとする。このアンテナ２は駆動部３により任意の方向に回転可能である。また、方向探知器４は駆動部３を制御することで不法電波の方向にアンテナ２を自動的に追尾させる機能を有している。受信部５はアンテナ２で捕捉した不法電波を受信する。電波解析部７は受信部５で受信した電波から、その電波が不法電波に類するものかを分析する。設定部８は外部から操作員が設定可能なキーボード及びマウスなどである。表示部６は、捕捉した不法電波のスペクトルや波形等の表示、後述の不法電波を発生している無線機のモデルや個体の番号等を表示するために使用する。
【００１３】
なお、受信部５には、空中を伝搬する電波が受信されるが、無線機自体の試験のためには空中を伝搬する電波を捕捉するのではなく、直接無線機の送信端子から有線で受信信号を取得するように構成してもかまわない。
【００１４】
受信部５は、ＦＭ波（例えば３０〜３０００ＭＨｚ）としての受信信号Ｓ１を出力し、ダウンコンバータ１０に入力する。この受信信号Ｓ１は、ダウンコンバータ１０において例えば２０．４８ＭＨｚでサンプリングされて中間周波数帯信号（ＩＦ信号）に変換され、その後、Ａ／Ｄ変換器１１に出力される。
【００１５】
このＩＦ信号は、その一部が分岐されてトリガ位置決定部１００Ｂに入力される。トリガ位置決定部１００Ｂでは、ＩＦ信号に対して、所定の演算処理を施すことで受信電波に対するトリガ位置を決定し、そのトリガ位置ＴｒをＡ／Ｄ変換器１１に出力する。その詳細は後述する。
【００１６】
Ａ／Ｄ変換器１１では、今回の受信電波に対し、そのトリガ位置決定部１００Ｂから受け取ったトリガ位置Ｔｒから読み取りを開始し、その読み取った信号をデジタル化しデジタルＩＦ信号Ｓ２として出力する。
【００１７】
デジタルＩＦ信号Ｓ２はデジタル直交検波部１５に入力され、その前段部において、処理前Ｉ信号Ｓ３と処理前Ｑ信号Ｓ４に変換される。処理前Ｉ信号Ｓ３と処理前Ｑ信号Ｓ４は間引き部１６で間引きされ、また必要に応じて帯域幅を制限された後、一度メモリ１７にバッファリングされる。すなわち、デジタルＩＦ信号Ｓ２はデジタル直交検波部１５で検波され、最終的に立上がり波形の振幅情報としてのＩ信号Ｓ５と、位相情報としてのＱ信号Ｓ６となってメモり１７から出力される。
【００１８】
図２はアンテナ経由での受信電波波形を示す図、図３は同一無線機からの受信信号を有線で試験的に受信部で直接取得した場合の受信電波波形を示す図である。この受信電波波形は、図１では図示されていない振幅−位相表示部で表示され、各図において、上段は振幅（Ｉ信号Ｓ５）を、下段は位相（Ｑ信号Ｓ６）をそれぞれ示している。
【００１９】
この図２と図３との比較から分かるように、アンテナ経由での測定デ−タ（図２）は全体的に雑音が重畳して立上り波形の特徴も識別できないほど信号が歪んでいる様子が分かる（図２の上段参照）。そして、このような雑音が重畳する受信電波に対しては、その立上がり波形を捕捉するのが困難である。一方、無線機の特定（同定）は、従来から立上がり波形に基づいて行うようになっており、したがって、雑音を含む受信電波に対してはその機能を十分に発揮できず、無線機特定が極めて困難であった。
【００２０】
これに対し、本発明では、（１）受信電波の立上がり時点を高精度で把握し、その立上がり時点にトリガを掛けることで受信電波を捕捉したこと、および（２）受信電波の立上がり波形に対し、本発明に特有の演算処理を施すことで、無線機の各々に固有のスペクトログラム・パターンを高精度で抽出できるようにしたことの２点により、雑音のひどい状況であっても無線機を特定できるようにしている。
【００２１】
以下の説明では、先ず（２）のスペクトログラム・パターン抽出について図４〜図１０を用いて説明し、続いて（１）トリガ位置決定の手法について図１１を用いて説明する。
【００２２】
上記のように、本発明では、取得デ−タからスペクトログラムに変換し、更に、無線機の特徴を適切に表すスペクトログラム・パタ−ンに変換している。パタ−ンの照合を効率よく行うためには、適切なスペクトログラム・パターン表示手法が重要である。そこで、迅速にスペクトログラム・パターンの全体像を把握できるツールとして、所定のソフトウェア処理を試みた。このソフトウェア処理は、図１の時間−周波数スペクトル解析部２４で実行され、デジタル直交検波部１５で得られたＩ信号Ｓ５とＱ信号Ｓ６とに対して行われる。その結果得られたスペクトログラム・パターンは表示部で表示される。一方、無線機同定データ部２２には、無線機のメーカおよびモデル毎に予め求めたスペクトログラム・パターンが基準スペクトログラム・パターンとして記憶されており、今回の受信電波のＩ信号Ｓ５等の立上がり波形から演算したスペクトログラム・パターンと、基準スペクトログラム・パターンとを対比することで、今回の受信電波の発信元である無線機のモデルを同定するようになっている。
【００２３】
図４はＡ社の無線機ａからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図である。上段に振幅の絶対値が表示されている。下段のスペクトログラム・パターンは、データ長Ｎｆｆｔ＝１０２４の高速フーリエ変換（ＦＦＴ）を行い、スペクトログラムを計算し、その信号レベルをさらに等高線で結んでいる。すなわち、データ取得時間全体にわたり最大信号レベルＰｍａｘの９５％から１０％までを０．０５スッテプの等高線でパターン化して表示している。デシベル表示選択の場合は、Ｐｍａｘの１ｄＢから２ｄＢステップで−２２ｄＢまでの表示となる。データ長Ｎｆｆｔ＝１０２４を採用した理由は、周波数分解能と時間分解能およびスペクトラム密度の広がり具合から決めた。縦軸は中心周波数からの周波数偏差をこの場合は±１５ｋＨｚ範囲で示し、横軸はトリガ開始ポイントからの経過時間を１０２ｍｓに亘ってパターン化して表示してある。このツールは対話形式で実行され、表示はリニヤとデシベルの選択が可能で、周波数のレンジ等は適切な大きさのパターンが得られるように工夫されている。
【００２４】
このように、スペクトログラム・パターンの信号レベルを等高線で結ぶことで、無線機に固有のパターンを把握することができ、図４下段のスペクトログラム・パターンでは、無線機ａに固有の等高線パターンが、丁度指紋のように完璧に捕らえられている。
【００２５】
図５は無線機ａからの電波が強い干渉を受けた場合の振幅波形の絶対値とスペクトログラム・パターンを示す図である。下段のスペクトログラム・パターンでの縦軸（周波数偏差）はスケールダウンして表示されている。図４と図５の振幅波形を比較してみると前者（図４）では振幅波形は若干、雑音の影響と隣接局からの干渉を受けている様子が伺える。しかし、後者（図５）の振幅波形になると干渉波と雑音の影響が支配的になり信号波形の原型が確認できないほどに劣化している。
【００２６】
また、スペクトログラム・パターンを比較してみると、周波数に段差が発生するポイントに着目すると、図５ではパターン全体が約４μｓ見かけ上早く立ち上がっている。これはトリガ設定の変更とＳ／Ｎ比変化のためであると思われるが、図４と図５のスペクトログラム・パターンを重ね合わせて一致具合を調べた結果、スペクトログラム・パターンの変形等の影響はほとんどなく両者のパターンはよい一致を示した。すなわち、スペクトログラム・パターンは、雑音の影響をほとんど受けず、無線機に固有のパターンを保持できることが分かった。
【００２７】
図６および図７はＢ社の無線機ｂからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図である。この無線機ｂは０〜３０ｍｓにかけて立上り時の周波数変化が激しく振動しその後は設定周波数に収斂する特性を示している。測定毎の立上がり時の変動を明示すために、測定を１０回連続して行い、その中で周波数変化の最大と最小の測定例を選んで示してある。図６が最小の場合、図７が最大の場合である。
【００２８】
このように、測定毎の立上がり時の変動が大きい場合、確かに立上がりの周波数のバラツキは大きいものの、全体としは、両者のパターンはよい一致を示しており、立上り時の周波数変化が激しくなる無線機に対しても、その無線機同定に本発明を十分に適用できることが分かった。
【００２９】
図８はＣ社の無線機ｃからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図である。この無線機ｃの場合、最初大きく周波数が振動しその後は減衰振動する傾向示す。他の無線機で行った結果では、減衰振動のパターンは複雑ではあるが、示すパターンは機別に異なり無線機の識別が可能であった。
【００３０】
上記図４〜図８に示した測定例は全てＩＦ帯域幅に関しては２５０ｋＨｚに設定し測定を行い、振幅波形の絶対値表示とスペクトログラム・パターンの処理段階で４０ｋＨｚのローパスフィルタを通して得られたものであるが、このスペクトログラム・パターンに対してＩＦ帯域幅と処理過程のローパスフィルタとが及ぼす影響について説明する。なお、この帯域幅制限は、図１の間引き部１６で施されるようになっている。
【００３１】
図９はＢ社の無線機ｂからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図である。この図９では、ＩＦ帯域幅５０ｋＨｚで測定し、測定したＩＱデータに遮断周波数２０ｋＨｚのローパスフィルタを通して処理した点が、上記の図６および図７とは相違している。この図９と、ＩＦ帯域幅２５０ｋＨｚでの測定例である図６，７と較べてみる。先ず、各図の上段側の振幅波形の比較では、図６，７の両図とも最初オーバシュートしてその後、一定のレベルに落ち着いている。図９では、オーバシュートの区間の振幅波形に歪みが生じている。この歪みの原因は、雑音除去に使用したローパスフィルタの影響によるものである。この区間のデータからローパスフィルタの特性を越える周波数成分が取り除かれたことによるものである。しかしながら、パターンの有効性を損なうような影響は観測されていない。
【００３２】
図１０はＢ社の無線機ｂからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図であって、ＩＦ帯域幅１２ｋＨｚで測定し、測定したＩＱデータに遮断周波数１０ｋＨｚのローパスフィルタを通して処理した結果を示している。この図１０での振幅波形の歪みは一層激しく、図６，７でオーバシュートを示した区間では、帯域制限を越えた周波数成分がフィルタの特性により削り取られ、振幅波形が大きく陥没している。この現象を下のパターンと関係つけてみれば、この部分の周波数成分が除去されていることが明らかになっている。しかし、下段側のスペクトログラム・パターン全体におけるパターンの傾向はそのまま保存され無線機の識別が可能なことが分かる。
すなわち、この図９および図１０から次の結果が得られた。すなわち、スペクトログラム・パターンと帯域幅との関係は、帯域幅が狭くなればなるほど雑音を抑圧してＳ／Ｎ比を改善する効果に期待できる。しかし、その反面、必要とする信号の情報が失われるおそれがある。時間領域の振幅波形はこの影響が大きく、振幅波形の特性が低下する。これに反して、周波数領域のスペクトログラム・パターンでは影響を受けた部分がはっきり確認できる。雑音の影響が大きい状況下において、無線機の識別を行う場合、雑音除去のために狭帯域化の信号処理して得たパターンであっても有効に活用できることが分かった。
【００３３】
次に、上記したトリガ位置決定の手法について説明する。アンテナ経由での受信電波は、図２の上段に示したように、全体的に雑音が重畳して立上り波形の特徴も識別できないほど信号が歪んでいる。このため従来、電波を受信しても、その電波の立上がり時点を捕らえるのが極めて困難であった。これに対し、本発明では、受信電波をダウンコンバータ１０でＩＦ信号に変換後、分岐させ、トリガ位置決定部１００Ｂに入力し、所定のデータ処理を施すことにより、トリガ位置を決定している。
【００３４】
このトリガ位置決定部１００Ｂは、上記の装置メイン部１００Ａと同様に、Ａ／Ｄ変換器１１０とデジタル直交検波部１５０とを有している。ダウンコンバータ１０から出力されたＩＦ信号は、サンプルレート１０．２４ＭＨｚでＡ／Ｄ変換器１１０に印加され、Ａ／Ｄ変換器１１０においてデジタルＩＦ信号Ｓ２０に変換される。続いてデジタル直交検波部１５で処理前Ｉ信号Ｓ３０と処理前Ｑ信号Ｓ４０に変換され、さらに間引き部１６０で間引きされ、また帯域幅を、例えば３１．２５ｋＨｚ、１５．６３ｋＨｚ、７．８１ｋＨｚの３種類の内の一つに制限された後、一旦メモリ１７０にバッファリングされる。その後振幅情報としてのＩ信号Ｓ５０と、位相情報としてのＱ信号Ｓ６０となって出力される。そして、直交検波部１５０の後段に設けられたスペクトル演算部１８０は、このＩ信号Ｓ５０とＱ信号Ｓ６０とをスペクトルに変換し、そのスペクトルはトリガ検出部１９０に入力される。
【００３５】
図１１はスペクトル演算部で得られたフリーランでのスペクトルを示す図である。この図に示すように、トリガ検出部１９０では、計算されたスペクトルの最大値と、設定部８で設定されたトリガレベルＴ０とを比較し、スペクトルが予め設定したレベル以上となる最初のポイントをトリガ位置Ｔｒとし、そのトリガ位置情報を取り込むように、装置メイン部１００ＡのＡ／Ｄ変換器１１に指示を送る。
【００３６】
このように、本発明では、トリガ位置を決定する際に、帯域幅を制限し、スペクトルに変換後、そのスペクトルが予め設定したレベル以上となる最初のポイントをトリガ位置としたので、受信電波がホワイトノイズ等でその立上がり時点を把握できないような場合でも、受信電波に対する読み取り開始位置を的確に検出することができる。
【００３７】
また、上記のトリガ位置から読みとった受信電波に対し、発明に特有の演算処理を施すことで、無線機の各々に固有のスペクトログラム・パターンを高精度で抽出できるようにしたので、雑音のひどい状況であっても無線機を特定できるようになった。
【００３８】
【発明の効果】
以上述べたように、この発明では、トリガ位置を決定する際に、帯域幅を制限し、スペクトルに変換後、そのスペクトルが予め設定したレベル以上となる最初のポイントをトリガ位置としたので、受信電波がホワイトノイズ等でその立上がり時点を把握できないような場合でも、受信電波に対する読み取り開始位置を的確に検出することができる。
【００３９】
また、上記のトリガ位置から読みとった受信電波に対し、発明に特有の演算処理を施すことで、無線機の各々に固有のスペクトログラム・パターンを高精度で抽出でき、そのパターンを基準スペクトログラム・パターンと比較するようにしたので、雑音のひどい状況であっても無線機を特定することができる。したがって、受信電波に雑音を含んでいる場合でも、不法無線局のどの無線機から発信されているかまで特定することができ、不法電波局の取り締まり等に大きく貢献することができる。
【図面の簡単な説明】
【図１】本発明の無線機同定装置の構成図である。
【図２】アンテナ経由での受信電波波形を示す図である。
【図３】同一無線機からの受信信号を有線で試験的に受信部で直接取得した場合の受信電波波形を示す図である。
【図４】図４はＡ社の無線機ａからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図である。
【図５】無線機ａからの電波が強い干渉を受けた場合の振幅波形の絶対値とスペクトログラム・パターンを示す図である。
【図６】Ｂ社の無線機ｂからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図である。
【図７】Ｂ社の無線機ｂからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図である。
【図８】Ｃ社の無線機ｃからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示す図である。
【図９】Ｂ社の無線機ｂからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示し、ＩＦ帯域幅５０ｋＨｚで測定し、測定したＩＱデータに遮断周波数２０ｋＨｚのローパスフィルタを通して処理した結果を示す図である。
【図１０】Ｂ社の無線機ｂからの受信電波にこのソフトウェアを適応して得られた振幅波形とスペクトログラム・パターンを示し、ＩＦ帯域幅１２ｋＨｚで測定し、測定したＩＱデータに遮断周波数１０ｋＨｚのローパスフィルタを通して処理した結果を示す図である。
【図１１】スペクトル演算部で得られたフリーランでのスペクトルを示す図である。
【符号の説明】
１      無線機
２      アンテナ
３      駆動部
４      方向探知器
５      受信部
６      表示部
７      電波解析部
８      設定部
１０      ダウンコンバータ
１１      Ａ／Ｄ変換器
１５      デジタル直交検波部
１６      間引き部
１７      メモリ
２２      無線機同定データ部
２４      時間−周波数スペクトル解析部
１００Ａ  装置メイン部
１００Ｂ  トリガ位置決定部
１１０    Ａ／Ｄ変換器
１５０    デジタル直交検波部
１６０    間引き部
１７０    メモリ
１８０    スペクトル演算部
１９０    トリガ検出部
Ｓ１      受信信号
Ｓ２      デジタルＩＦ信号
Ｓ３      処理前Ｉ信号
Ｓ４      処理前Ｑ信号
Ｓ５      Ｉ信号
Ｓ６      Ｑ信号
Ｓ２０    デジタルＩＦ信号
Ｓ３０    処理前Ｉ信号
Ｓ４０    処理前Ｑ信号
Ｓ５０    Ｉ信号
Ｓ６０    Ｑ信号
【図１】
2004214817.tif 000003
【図２】
2004214817.tif 000004
【図３】
2004214817.tif 000005
【図４】
2004214817.tif 000006
【図５】
2004214817.tif 000007
【図６】
2004214817.tif 000008
【図７】
2004214817.tif 000009
【図８】
2004214817.tif 000010
【図９】
2004214817.tif 000011
【図１０】
2004214817.tif 000012
【図１１】
2004214817.tif 000013