(57)【特許請求の範囲】
【請求項１】
セルラ通信システムの少なくとも第1及び第2の基地局から到来する信号の受信タイミングを使って端末の位置を推定する無線端末位置測定方法において、
上記端末の記憶回路に上記第1及び第2の基地局からの受信信号を蓄積する第１のステップと、
該第2の基地局からの受信信号よりも受信強度の強い上記第1の基地局からの受信信号の受信タイミングを測定する第２のステップと、
上記蓄積された受信信号から上記第１の基地局からの受信信号のレプリカを作成する第３のステップと、
上記蓄積された受信信号から上記レプリカを差し引いた信号を上記蓄積された受信信号に上書きする第４のステップとを有し、
上記第４のステップの後に、該上書きされた信号に基づいて上記第２の基地局からの受信信号について上記第２のステップを行うことを特徴とする無線端末位置測定方法。
【請求項２】
請求項１において、
符号分割多元接続（ＣＤＭＡ）方式を用いて、上記第1の基地局からの受信信号成分を逆拡散−復調−再拡散の手順で処理することにより上記レプリカを作成する第５のステップを有することを特徴とする無線端末位置測定方法。
【請求項３】
請求項２において、上記第５のステップは、逆拡散−復調−再拡散の後に受信信号を増幅する第６のステップを有することを特徴とする無線端末位置測定方法。
【請求項４】
請求項３において、上記第６のステップは、上記第1の基地局からの信号伝搬路による振幅変動及び位相回転の少なくとも一方を補正する第７のステップを含むことを特徴とする無線端末位置測定方法。
【請求項５】
セルラ通信システムの少なくとも第1及び第2の基地局から到来する信号の受信タイミングを使って端末の位置を推定する無線端末位置測定システムにおける無線端末装置であって、
第１の基地局からの受信信号と第２の基地局からの受信信号とを含む信号を蓄積する記憶回路と、
上記記憶回路に蓄積された信号から上記第２の基地局からの受信信号よりも受信強度の強い上記第１の基地局からの受信信号の受信タイミングを測定する受信タイミング測定部と、
上記記憶回路に蓄積された信号から上記第1の基地局からの受信信号をキャンセルする信号処理回路とを有し、
上記信号処理回路は、該記憶回路に蓄積された信号から上記第1の基地局の受信信号のレプリカを作成するレプリカ信号発生回路と、上記蓄積した信号から該レプリカを差し引く引き算回路とを有し、
上記記憶回路に蓄積された信号は、上記レプリカを差し引いた後の信号で上書きされ、
上記受信タイミング測定部は、上記上書きされた後の記憶回路に蓄積された信号に基づいて上記第２の基地局からの受信信号について受信タイミングを測定することを特徴とする無線端末装置。
【請求項６】
請求項５において、
上記レプリカ信号発生回路は、符号分割多元接続（ＣＤＭＡ）方式を用いて、上記第1の基地局からの受信信号成分を逆拡散する逆拡散回路と、該逆拡散回路の出力信号から復調する復調回路と、該復調回路の出力信号を再度拡散する再拡散回路とを有することを特徴とする無線端末装置。
【請求項７】
請求項６において、上記再拡散回路の出力信号を増幅する増幅回路を有することを特徴とする無線端末装置。
【請求項８】
請求項７において、上記第1の基地局からの信号伝搬路による振幅変動及び位相回転の少なくとも一方を補正する補正回路を含むことを特徴とする無線端末装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、セルラ通信を使って端末位置を測定するシステムに関する。
【０００２】
【従来の技術】
特開平７−１８１２４２は、セルラ方式で各基地局がＧＰＳ同期したＣＤＭＡ方式を用いているシステムにおいて、各基地局から届くパイロット信号の受信タイミングを観測することで端末位置を測定することを開示している。
【０００３】
図２、３は特開平７−１８１２４２に記載された技術を示す。図２にＧＰＳ衛星２０に同期した基地局２１、２２、２３から送信されている基準信号２５、２６、２７を受信する端末２４が示されている。図３に相関演算結果の例が示されている。
【０００４】
端末２４は、受信信号に対し特定符号パタンの相関演算を行うことで、各基地局から送信される基準信号の受信タイミングを検出することができる。ＣＤＭＡ方式では、各基地局が出す共通パイロット信号が特定パタンの信号であり、端末はこの信号に相関演算を行うことで受信タイミングを検出することが可能である。各基地局は、ＧＰＳ同期により基準時間を合わせ、合わせたシステム時間の特定の送信タイミングにおいてパイロット信号を送信している。また上記の特定の送信タイミングは、オフセット時間と呼ばれ、シンクチャネルを通じて情報が送信されており、端末は自由に入手することができる。端末は、測定した受信タイミングから既知である送信タイミングの差を計算することで、電波伝搬にかかった遅延時間を知ることができる。
【０００５】
図３は相関演算結果の例を示している。図は遅延プロファイルと呼ばれるもので、どのような遅延パスが観測されるかを示している。横軸は遅延時間、すなわち送信タイミングによる補正済みの受信タイミングである。単位は拡散符号のチップに対応している。縦軸は相関演算出力で、相関値が大きい部分は信号の受信があったこと、すなわちパスの存在を示している。この結果を使えば、電波が基地局から端末に到来するまでにかかった相対遅延時間を求めることができる。ここで遅延時間が相対値であるのは、端末では絶対時間がわからないからである。得られた相対遅延時間に光速をかけることで、相対伝搬距離差を求めることができる。少なくとも３局以上の基地局について相対伝搬距離差が求まれば、３点測位により端末位置を推定することが可能である。
【０００６】
図４、５を用いてセルラ通信を測位に使う場合のもう１つの課題である遠近問題について説明する。セルラ通信では基地局−端末位置に依存して受信電力が変化し、必要とされるダイナミックレンジは１００ｄＢを越える。このため通常端末には自動利得調整（ＡＧＣ）機能が持たされていて、受信信号の強さに応じて自動的にＡＤ変換前の信号強度が調整される。図４は、端末において近傍の３つ基地局から到来する信号を受信している場合の受信信号を周波数領域で示したものである。帯域内の総受信電力６４は雑音電力６０に３つの基地局からの信号６１、６２、６３を加えたものである。最も近くにある基地局からの信号６１は、伝搬距離が短いため電力が強く、全受信信号の支配項となっている。図５は、２つの基地局Ａ，Ｂ間の受信状態を示している。雑音電力７１は熱雑音が原因として発生するもので、端末位置に依存せず一定である。ＡＧＣは総受信信号レベル７０に適応するように設定されている。したがって例えば端末が基地局Ａ付近にあるときには、基地局Ａからの信号が支配項となってＡＧＣが働く。量子化雑音電力７４はＡＧＣに応じて変動するため、基地局Ａ付近では値が高くなっている。しかもこの時、基地局Ｂからの信号電力は伝搬距離が長くなるため低下しており、基地局Ｂからの信号は等価的にＳ／Ｉ（信号対干渉電力比）＝（Ｂからの信号）／（Ａからの信号＋雑音）が非常に小さな値になってしまい、信号受信が困難になる。しかし、無線測位を行うには、端末が基地局Ａ付近であっても基地局Ｂからの信号を判別する必要がある。
【０００７】
また他の課題について説明する。移動体通信の端末では、低価格化のために周波数精度があまり高くない局部発信機を用いられていて、最寄りの基地局に同期することでキャリア周波数偏差を削減している（ＡＦＣ機能）。しかしＡＦＣ機能の限界から端末−基地局間には数Ｈｚの周波数差が残り、例えフェージングがない場合でも完全に同期することができない。このため、受信した信号の位相は数Ｈｚのゆっくりした回転運動をしている。このため、たとえ端末のユーザが静止あるいは歩行速度程度のゆっくりとした動きしかしていない場合にも、受信信号はどんどん回転していてしまい、長期間の同相加算を行うことが困難である。遠方の基地局からの信号は等価的にＳ／Ｉが低下しているため、同相加算回数を増加させたいが、上記現象のため加算回数はある程度以上は長くすることができない。例えばキャリヤ周波数８００ＭＨｚ、ＡＦＣ後の周波数安定度が０．０１ｐｐｍの端末があった場合、位相回転の周波数は８Ｈｚになる。同相加算に必要な位相回転許容値を３６度以内とすると、１／８×３６／３６０＝０．０１２５［秒］以内であれば同相加算可能なことがわかる。同相加算回数をこれよりも長い時間行った場合、信号ベクトルが回転しまい、かえってＳ／Ｉが劣化する現象が発生する。
【０００８】
【課題を解決するための手段】
上記課題は、セルラ通信システムの少なくとも第１及び第２の基地局から到来する信号を使って端末の位置を推定する無線端末位置測定方法において、該第２の基地局からの受信信号よりも受信強度の強い上記第１の基地局からの受信信号をキャンセルする第１のステップと、該第１のステップの後に上記第２の基地局からの受信信号を上記端末にて処理する第２のステップとを有することを特徴とする無線端末位置測定方法により解決される。
【０００９】
少なくとも２つの基地局から到来する電波が、同一の端末に受信される場合、伝搬距離差やフェージングの状況により、両者の信号強度には、大きな差が発生している。したがって信号強度の弱い基地局からの信号を受信する際に、信号強度の強い基地局の信号は干渉波として働いている。この干渉波をキャンセルして信号受信することにより、受信信号強度の大きな基地局からの干渉の影響を殆どなくすことができる。
【００１０】
【発明の実施の形態】
図１は、移動端末の構成例を示す。
【００１１】
移動通信端末は通常の通信動作を可能とするため、受信機と送信機が備わっている。本実施例はＣＤＭＡセルラを採用している。アンテナ１が受信した信号は、ＲＦユニット２にてベースバンドの信号に変換される。変換された信号はＡ／Ｄ変換器４においてディジタル信号に変換される。変換された信号から、逆拡散ユニット６，７によってパイロット信号と、通話チャネル信号が取り出される。パイロット信号を逆拡散するブロック７の出力により、伝搬路による位相回転と基地局−端末間のキャリヤ周波数差による位相回転が検出される。これをもとにＡＦＣ９の機能が動作する。パイロット信号の逆拡散出力である伝搬路推定結果に対して共役をとり、これに通話チャネルを逆拡散するブロック６の出力をかけることで復調８が行われる。復調回路では伝搬路符号化も復調される。復調された信号は音声コーデック１３において音声信号に変換され、スピーカー１４から出力される。他方、マイク１５から取り込まれた音声信号は音声コーデック１３においてディジタル信号に変換される。変換された信号は変調器１０で伝送路符号化と変調が行われ、拡散部１１において拡散後、Ｄ／Ａ変換器１２によりアナログ信号に変換される。アナログ信号はＲＦ部２でキャリヤ周波数にアップコンバートされてから、アンテナ１より送信される。
【００１２】
第１の実施例を、図６を用いて説明する。図６は干渉除去を行うための干渉除去装置の構成を示している。
【００１３】
端末アンテナ１で受信された信号には、希望信号（再生したい信号）と干渉信号（他の信号）の両方が含まれる。ＣＤＭＡ方式では、両信号は同一タイミングで、同一周波数帯に含まれる。Ａ／Ｄ変換器４においてディジタル信号に変換された信号を信号処理装置１７はそのままメモリー５に蓄積する。
【００１４】
図６に示す干渉除去装置は信号処理装置１７内に構築されている。本装置は、ハードウェアあるいはソフトウェアで実現される。この干渉除去装置は、メモリー５の内容を読み出して、干渉信号のレプリカ１２２を作成する。そして、作成したレプリカを信号処理装置１７内の引き算回路１２６にてメモリーの内容から差し引く。これにより、強い干渉波となる近傍基地局の影響を大幅に削減する。干渉除去によりＳ／Ｉを高めてから計算される遅延プロファイルを使って、端末位置が推定される。
【００１５】
メモリー５内のデータには、パイロット信号、ページング信号、通話信号の３種類が含まれる。それぞれは全て別の拡散コードにより拡散されている。ブロック１０１〜１０６は個別の拡散コードにより逆拡散するブロックである。ブロック１０７はパイロット信号の逆拡散結果から伝搬路推定を行う。得られた推定結果を用い、ブロック１０９〜１１３は復調を行う。復調の結果、ＢＰＳＫ通信であれば、Ｉチャネルの符号の正負によって送信符号を判定することができる。ここでは伝搬路符号化も復調することで復調後の信号の符号誤り率を改善が考えられるが、一部の情報は秘匿のためにスクランブルもかかっているため、一般にはデスクランブル情報を持たない端末では伝搬路符号化を復調することはできない。しかし、信号自体のＳ／Ｉがよいため、伝搬路符号化を復調しなくても送信されている符号自体はほぼ正確に推定することが可能である。得られた仮判定結果はブロック１１４〜１１９にて再び同じ拡散符号で再拡散され、ブロック１２０にて加算されて特定基地局から送信された信号が再現される。ブロック１２３は伝搬路推定から遅延プロファイルを作成するブロックである。遅延プロファイルは、逆拡散で使用する拡散符号の系列の位相をずらしながらパイロット信号を逆拡散することで得られる。ここで遅延プロファイルは複素量である。ブロック１２０で加算された送信信号に、ブロック１２３で得られた遅延プロファイルをかけることで、受信信号のレプリカ１２２が完成する。このとき、レプリカを作成した基地局に関する遅延プロファイルは、ブロック１２３で既に作成されている。ブロック１２６は、メモリー５に蓄積された信号からレプリカ１２２を引き算する干渉除去機能である。干渉波が除去された信号は再びメモリー５に上書きされる。他の基地局からの信号についても同様に順次強いものから除去しながら遅延プロファイルを作成することでＳ／Ｉを高めてプロファイル作成が可能となる。強い信号が全て除去された後には、ブロック１０１、１０７、１２３を使って遅延プロファイルの作成のみを行い、受信電力が弱い基地局に関する遅延プロファイルを作成することができる。このようにして複数の基地局に関する遅延プロファイルが作成されれば、ＣＰＵ１６は各遅延プロファイルからパスを特定し、相対遅延差を計算する。またＣＰＵ１６は計算した相対遅延差と予め知っている基地局の位置情報から端末の位置を計算する。
【００１６】
以上の処理から遠近問題が発生している場合でも、複数の基地局からの信号を高いＳ／Ｉで受信することが可能となり、各基地局までの相対距離差が測定できる。よって端末の位置検出が可能となり、課題は解決される。
【００１７】
上記の方法は、Ｓ／Ｉを向上させるために、干渉電力（Ｉ）を抑圧する方法であるが、他方、信号電力（Ｓ）を向上させる方法がある。通常の通信と異なり、端末位置を測定する場合には即時性が要求されない。そこで、メモリー５に蓄積した全てのデータを使ってパス検出に必要な遅延プロファイルを作成することを考える。
【００１８】
信号電力（Ｓ）の向上に関し、原因と対策についてより詳細に説明する。
【００１９】
端末では、図１のＡＧＣ３が、受信した信号のレベルに応じてＲＦ部の利得を調整し、ダイナミックレンジを広く確保している。ＡＧＣは、例えば受信電力を累積平均して受信レベルを求め、その値が一定になるようにフィードバックしてＲＦ部の利得を調整している。ＡＧＣにより、Ａ／Ｄ変換器はビット幅を低減することができる。
【００２０】
ところで測位の場合には、最低３局以上の基地局が必要である。第１の基地局よりも遠い第２の基地局からの希望信号のＳ／Ｉを向上させるには、同相加算回数を大きくする方法が有効である。同相加算回数が大きいほど、遠くの基地局が観測できるようになり、観測基地局数の増加から、端末位置を推定する際の推定精度が向上する。
【００２１】
同相加算を行うには、２つの補償が必要である。１つは端末と基地局間で局部発信器の周波数偏差による受信信号の位相回転の補償である。もう１つは端末の位置や周囲環境が変化するために発生するフェージングによる位相回転の補償である。
【００２２】
端末−基地局間の同期ずれによる位相回転は、主に端末の周波数合わせ性能が十分でないことが原因で発生する。基地局間では、それぞれの基地局がＧＰＳ同期しているために十分な精度が確保されており、周波数偏差は非常に小さい。このため、複数ある基地局のうち、１つの基地局に同期ができれば、他の基地局の位相回転も同期することができる。同期する候補としては、最も信号品質の良い、すなわち相関演算で最大ピークとなった基地局信号の位相回転を基準にすることが有効である。
【００２３】
特定基地局のパイロット信号を逆拡散すると図７の様になる。雑音成分が残留するが、その位相回転を１次式で近似することで位相の平均回転数を推定することができる。推定された回転数で逆回転をする回転演算により、受信データに含まれる端末−基地局間キャリヤ周波数差はほぼ取り除かれる。補償後のデータは位相が揃っているため、ただ単純に加算するだけで同相加算演算が完了する。これによって受信信号のＳ／Ｉを向上させることができる。
【００２４】
１次近似を精度よく求めるためには、同相加算を行う全データを使って近似式を求めることである。このため端末は測位ジョブが開始されると、一定期間について受信データをメモリー５に蓄える。（メモリーの容量は、例えばチップレート１．２２８８ＭＨｚのＣＤＭＡセルラ方式で、４倍のオーバーサンプルを考え、１サンプルあたりＩ，Ｑ合わせて８ビットの情報がある場合には、１秒間のデータの収集に約５ＭＢの容量が必要である。）  端末は最も信号レベルの高い信号の受信タイミングを既に知っている。何故なら端末は通信を行うために周囲の基地局を常にサーチし、信号レベルが最も高い基地局を追い続けているからである。したがって端末はメモリー５に蓄えられた情報に対して、既知である最も信号レベルが高い基地局からの信号タイミングに合わせてパイロット信号を抽出する逆拡散演算を行い、伝搬路の推定を行う。ここでの逆拡散は０．１〜１０ｍｓ程度の比較的拡散率の低いものを実施し、複数の値を得る。図７は、実測で０．８秒間の逆拡散結果の位相を示した図である。図からわかるように得られた伝搬路推定結果は干渉電力や雑音電力の影響が大きく残る。しかしながら干渉＋雑音成分は、統計処理により取り除くことができる。これは、干渉＋雑音成分が広周波数帯域の信号であるのに対し、信号成分の回転ベクトルは周波数の低い信号であるため、平均操作により低周波数の成分を抽出すれば、その回転量が推定できることを利用している。図７から、回転量の１次式を推定することで、基地局と端末の周波数偏差は大幅に抑えられることが分かる。基地局と端末間の周波数偏差が抑えられたら、各基地局からの信号の受信タイミング付近において逆拡散演算を行うことで、各基地局からの信号の遅延プロファイルを算出することができる。
【００２５】
上記実施例では全ての測定データを一旦メモリーに蓄積する方法に関して説明したが、以下に示す別の方法も本発明に含まれる。
【００２６】
位相回転の推定結果が１次式で近似しているため、端末と基地局の周波数偏差が時間的にあまり変化しない場合には、過去の値を使って、端末に対する基地局信号の位相回転偏差を推定することが可能である。ただし基地局と端末の位相差は時間によって変わるため、最近の値を使う必要がある。受信品質が良好な特定の基地局のみを観測すればよいため、位相回転推定に必要となる逆拡散結果を蓄積するためのメモリーは大容量を必要としない。基地局からの報知情報から近傍に存在する基地局の送信タイミングは既知であるから、拡散符号系列をその近傍の基地局タイミング付近に合わせて逆拡散することで、遅延プロファイルを作成することができる。メモリーには逆拡散後の信号を蓄積することでメモリーの容量は大幅に削減できる。逆拡散結果は、先に予測されている位相回転量の補正を行って同相加算することで、上記と同等の同相加算の効果を得ることができる。よって大幅にメモリー量を削減することができる。
【００２７】
位相の回転がフェージングに起因するものである場合、複数の基地局の位相回転は独立に変動する。２つの基地局に関して位相回転推定を行い、その回転量の相関をとる。相関が例えば０．７以上であれば、互い基地局の位相回転量がほぼ同じであることを示しており、フェージングは発生していないと考えられる。逆に相関が例えば０．７以下になるようであれば、その長さの位相補正ができないことを示しており、より短い時間での同相加算しかできないことを示唆している。したがって同相加算回数が相関を閾値として定めることができる。
【００２８】
以上説明した手順で得られた遅延プロファイルより、従来例で説明した方法により端末の位置が計算される。
【００２９】
本実施例では端末だけで位置測定が可能であり、ＧＰＳ等の付加装置も不要であることから、端末が測位操作を行う際に必要となる消費電力を削減することが可能である。基地局の位置に関する情報は、ネットワークに繋がるサーバから得る方法や、基地局から送信される報知チャネルで周囲の基地局の位置を送信することで端末は知ることができる。一旦蓄積した基地局位置情報は、保持することで、サーバへのアクセスや報知情報の受信回数を減らすことができるため、端末の低消費電力化に有効である。
【００３０】
本実施例の説明ではＣＤＭＡ方式を例に挙げて説明したが、これに限るものではない。送信されている情報が既知であればＴＤＭＡやＦＤＭＡ方式においても本方法を適用することが可能である。ＴＤＭＡやＦＤＭＡの場合、各基地局からの信号が到来するタイミングや周波数が異なるため、端末では信号の含まれるタイミングや周波数を全て走査して得られた情報を端末位置管理局に伝送する必要がある。
【００３１】
以上に挙げた実施例では、位置測定用の端末に関して送信タイミングを測定し、端末位置を推定する方法について記載した。他方、本発明は位置測定用の端末だけでなく、基地局送信タイミング測定装置に対しても有効である。基地局送信タイミング測定装置は、各基地局の送信タイミングを測定し、測位をする際に含まれる各基地局の送信タイミング誤差を補正する場合に用いられる。基地局送信タイミング測定装置はセルラ回線、あるいはネットワークにつながっている。従って測位用の端末は、セルラ回線及びネットワークを介して基地局送信タイミング測定装置に送信タイミングの誤差を問い合わせ、正確な位置測定が可能となる。しかしながら基地局送信タイミング測定装置の設置可能な場所は、設置コストから限定されるため、１局の設置でできるだけ遠くの基地局まで測定したい。このような時には、端末の場合と同様に、遠近問題が発生するため、干渉除去や同相加算回数を長くするための工夫が必要である。本発明は、基地局送信タイミング測定装置に適用すれば遠近問題を回避し、遠くの基地局の送信タイミングまで測定することが可能である。
【００３２】
【発明の効果】
本発明によれば、セルラ通信システムを用いて端末の測位を行う場合に、メモリーに受信情報を貯めて推定精度を上げるために必要となる様々な処理を行うことができる。
【図面の簡単な説明】
【図１】本発明からなる第１の実施例の端末構成を示す図。
【図２】従来例を説明する図。
【図３】端末で受信した信号を相関演算した結果を示す図。
【図４】特定の基地局付近に端末がある場合の受信信号スペクトラムを示す図。
【図５】２つの基地局を例に挙げて基地局付近ではＳ／Ｉが低くなることを説明する図。
【図６】本発明からなる第３の実施例でレプリカ作成のための機能ブロック図。
【図７】特定基地局の逆拡散結果の位相回転量を示す図。
【符号の説明】
１．．．アンテナ、２．．．無線部（ＲＦ部）、３．．．ＡＧＣ、４．．．ＡＤ変換器、５．．．メモリー、６、７．．．逆拡散器、８．．．復調器、９．．．ＡＦＣ、１０．．．変調器、１１．．．拡散器、１２．．．ＤＡ変換器、１３．．．音声コーデック、１４．．．スピーカ、１５．．．マイク、１６．．．ＣＰＵ、１７．．．ＤＳＰ、２０．．．ＧＰＳ衛星、２１、２２、２３．．．基地局、２４．．．端末、２５、２６、２７．．．基地局から送信される信号、４０、４１、４２、４３．．．相関演算結果のピーク、４４．．．最短伝搬距離信号の差（相対遅延差）、６０．．．雑音電力、６１．．．最も近い基地局からの信号電力、６２、６３．．．他の基地局からの信号電力、６４．．．トータルの受信電力、７０．．．トータルの受信電力、７１．．．雑音電力、７２．．．基地局Ａからの信号電力、７３．．．基地局Ｂからの信号電力、７４．．．電子化雑音電力、１００．．．受信信号、１０１、１０２、１０３、１０４、１０５、１０６．．．逆拡散部、１０７．．．伝送路推定部、１０８、１０９、１１０、１１１、１１２、１１３．．．復調部、１１４、１１５、１１６、１１７、１１８、１１９．．．再拡散部、１２０．．．加算器、１２１．．．位相回転部、１２２．．．レプリカ信号、１２３．．．遅延プロファイル作成部、１２４．．．遅延プロファイル、１２６．．．引き算回路。
【図１】
0003573039.tif 000002
【図２】
0003573039.tif 000003
【図３】
0003573039.tif 000004
【図４】
0003573039.tif 000005
【図５】
0003573039.tif 000006
【図６】
0003573039.tif 000007
【図７】
0003573039.tif 000008