(57)【要約】  本発明は中継機受信チャネルの逆チャネル推定装置およびその方法に関するもので、歪み推定信号が挿入されている受信信号から送信機と中継機間の伝送路によって設定されたチャネルを推定して、推定されたチャネルから中継機受信チャネルの逆チャネルの情報を有する時間領域でのフィルタ係数を安全性と因果性が保障されるように生成するための、中継機受信チャネルの逆チャネル推定装置およびその方法を提供するものである。このために、本発明は逆チャネル推定装置において、約束された信号が挿入された受信信号から伝送路によって設定されたチャネルを推定するためのチャネル推定手段と、前記チャネル推定手段で推定したチャネルから最小位相チャネルの逆チャネルを生成して時間領域のフィルタ係数に変換するための逆最小位相フィルタ係数生成手段と、前記チャネル推定手段で推定したチャネルと前記最小位相チャネルを用いて逆グループ遅延成分を有する時間領域のフィルタ係数を生成するための逆グループ遅延フィルタ係数生成手段と、前記逆最小位相フィルタ係数生成手段で変換されたフィルタ係数と前記逆グループ遅延フィルタ係数生成手段で生成したフィルタ係数を結合して時間領域で受信チャネルの逆チャネルを生成するための結合手段とを備える。
2010538566.tif 000002
【特許請求の範囲】
【請求項１】
逆チャネル推定装置において、
約束された信号が挿入された受信信号から伝送路によって設定されたチャネルを推定するためのチャネル推定手段と、
前記チャネル推定手段で推定したチャネルから最小位相チャネルの逆チャネルを生成して時間領域のフィルタ係数に変換するための逆最小位相フィルタ係数生成手段と、
前記チャネル推定手段で推定したチャネルと前記最小位相チャネルを用いて逆グループ遅延成分を有する時間領域のフィルタ係数を生成するための逆グループ遅延フィルタ係数生成手段と、
前記逆最小位相フィルタ係数生成手段で変換されたフィルタ係数と前記逆グループ遅延フィルタ係数生成手段で生成したフィルタ係数を結合して時間領域で受信チャネルの逆チャネルを生成するための結合手段と、
を備えることを特徴とする逆チャネル推定装置。
【請求項２】
前記逆最小位相フィルタ係数生成手段は、
前記チャネル推定手段で推定したチャネルからｚ−平面上のポールゼロが単位円内に存在する最小位相チャネルを生成するための最小位相チャネル生成手段と、
前記最小位相チャネル生成手段で生成した最小位相チャネルの逆チャネルを生成するための逆チャネル生成手段と、
前記逆チャネル生成手段で生成した周波数領域の最小位相チャネルの逆チャネルを時間領域フィルタ係数に変換するための第１逆高速フーリエ変換手段と、
を備えることを特徴とする請求項１に記載の逆チャネル推定装置。
【請求項３】
前記最小位相チャネル生成手段は、
前記チャネル推定手段で推定したチャネルをログ変換するためのログ変換手段と、
前記ログ変換手段で変換したチャネルを時間領域の係数に変換するための第２逆高速フーリエ変換手段と、
周波数不変の線形フィルタ係数を生成するための線形フィルタ係数生成手段と、
前記第２逆高速フーリエ変換手段で変換された係数と前記線形フィルタ係数生成手段で生成した周波数不変線形フィルタ係数を乗算するための乗算手段と、
前記乗算手段からの出力係数を周波数領域の係数に変換するための高速フーリエ変換手段と、
前記高速フーリエ変換手段で変換された係数を指数変換して最小位相チャネルを生成するための指数変換手段と、
を備えることを特徴とする請求項２に記載の逆チャネル推定装置。
【請求項４】
前記逆グループ遅延フィルタ係数生成手段は、
前記チャネル推定手段で推定したチャネルを前記最小位相チャネルで分けてチャネルのオールパス成分を生成するためのオールパス成分生成手段と、
前記オールパス成分生成手段で生成したチャネルのオールパス成分の位相情報からグループ遅延成分を計算するためのグループ遅延計算手段と、
前記グループ遅延計算手段で計算したグループ遅延成分の逆に該当する逆グループ遅延成分を有する時間領域のフィルタ係数を生成するための逆グループ遅延係数生成手段と、
を備えることを特徴とする請求項１に記載の逆チャネル推定装置。
【請求項５】
前記チャネル推定手段は、
受信チャネルの歪みを推定するのに用いられる基準信号を保存している基準信号保存手段と、
前記受信信号から歪み推定信号を抽出するための歪み推定信号抽出手段と、
前記歪み推定信号抽出手段で抽出した歪み推定信号と前記基準信号保存手段に保存されている基準信号とを比較して前記受信チャネルの歪みを推定するためのチャネル歪み推定手段と、
を備えることを特徴とする請求項４に記載の逆チャネル推定装置。
【請求項６】
前記約束された信号は、
パイロット信号またはトレーニングシーケンスであることを特徴とする請求項４に記載の逆チャネル推定装置。
【請求項７】
逆チャネル推定方法において、
約束された信号が挿入された受信信号から伝送路によって設定されたチャネルを推定するチャネル推定ステップと、
前記推定したチャネルから最小位相チャネルの逆チャネルを生成して時間領域のフィルタ係数に変換する逆最小位相フィルタ係数生成ステップと、
前記推定したチャネルと前記最小位相チャネルを用いて逆グループ遅延成分を有する時間領域のフィルタ係数を生成する逆グループ遅延フィルタ係数生成ステップと、
前記変換したフィルタ係数と前記生成したフィルタ係数を結合して時間領域で受信チャネルの逆チャネルを生成する逆チャネル生成ステップと、
を含むことを特徴とする逆チャネル推定方法。
【請求項８】
前記逆最小位相フィルタ係数生成ステップは、
前記推定したチャネルからｚ−平面上のポールゼロが単位円内に存在する最小位相チャネルを生成する最小位相チャネル生成ステップと、
前記生成した最小位相チャネルの逆チャネルを生成するステップと、
前記生成した周波数領域の最小位相チャネルの逆チャネルを時間領域フィルタ係数に変換するステップと、
を含むことを特徴とする請求項７に記載の逆チャネル推定方法。
【請求項９】
前記最小位相チャネル生成ステップは、
前記推定したチャネルをログ変換するステップと、
前記ログ変換したチャネルを時間領域の係数に変換するステップと、
周波数不変の線形フィルタ係数を生成するステップと、
前記変換した時間領域の係数と前記生成した周波数不変線形フィルタ係数を乗算するステップと、
前記乗算した結果を周波数領域の係数に変換するステップと、
前記変換した周波数領域の係数を指数変換して最小位相チャネルを生成するステップと、
を含むことを特徴とする請求項８に記載の逆チャネル推定方法。
【請求項１０】
前記逆グループ遅延フィルタ係数生成ステップは、
前記推定したチャネルを前記最小位相チャネルで分けてチャネルのオールパス成分を生成するステップと、
前記生成したチャネルのオールパス成分の位相情報からグループ遅延成分を計算するステップと、
前記計算したグループ遅延成分の逆に該当する逆グループ遅延成分を有する時間領域のフィルタ係数を生成するステップと、
を含むことを特徴とする請求項７に記載の逆チャネル推定方法。
【請求項１１】
前記チャネル推定ステップは、
受信チャネルの歪みを推定するのに用いられる基準信号を保存しているステップと、
前記受信信号から歪み推定信号を抽出するステップと、
前記抽出した歪み推定信号と前記保存している基準信号とを比較して前記受信チャネルの歪みを推定するステップと、
を含むことを特徴とする請求項１０に記載の逆チャネル推定方法。
【請求項１２】
前記約束された信号は、
パイロット信号またはトレーニングシーケンスであることを特徴とする請求項１０に記載の逆チャネル推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、中継機受信チャネルの逆チャネル推定装置およびその方法に関し、より詳細には、主送信機または他の中継機から伝送された信号から、中継機受信チャネルの逆チャネル情報を有する時間領域でのフィルタ係数を安全性（ｓｔａｂｉｌｉｔｙ）と因果性（ｃａｕｓａｌｉｔｙ）が保障されるように生成するための、中継機受信チャネルの逆チャネル推定装置およびその方法に関する。
【背景技術】
【０００２】
本発明は、情報通信部および情報通信研究振興院のＩＴ新成長動力核心技術開発事業の一環として行われた研究から導出されたものである［課題管理番号：２００７−Ｓ−００６−０１、課題名：ＯＦＤＭ伝送方式での同一チャネルＲＦ中継技術］。
【０００３】
中継機は、主送信機の信号が弱く受信される地域に設置され、受信の不安を解消し主送信機信号の伝送領域を広げる役割を果たす。
【０００４】
一般的に、主送信機と中継機の間にはチャネル歪みが発生し、中継機は、主送信機と中継機の間の受信チャネルの歪みを補償して再伝送する。このために中継機は、受信チャネルを推定し、その逆チャネルの情報を有する時間領域でのフィルタ係数を生成し、前記生成したフィルタ係数を用いて受信信号をフィルタリングすることによって受信チャネルの歪みを補償する。
【０００５】
以下、図１を参照し、従来の中継機受信チャネルの逆チャネル推定方法に対して説明する。
【０００６】
図１は、一般的な周波数領域での中継機受信チャネルとその逆チャネルに対する一例示図である。
【０００７】
図１に示したように、推定した受信チャネル（Ｈ）１１からその逆チャネル（ＨINV）１２を生成する従来の方法は、除算を用いた直接的な逆変換で下記の数１のとおりである。
【０００８】
【数１】
2010538566.tif 000003
【０００９】
ここで、Ｎは、推定したチャネルのサンプル（ｓａｍｐｌｅ）個数を表し、Ｔは転置（ｔｒａｎｓｐｏｓｅ）を表す。
【００１０】
このような除算を用いた逆変換方法が安全性と因果性を保障するためには、推定した受信チャネルに対するｚ−平面上のポールゼロが図２に示したように単位円内に存在しなければならない。
【００１１】
しかし、実際の環境で推定した受信チャネルのｚ−平面上のポールゼロの一部は図３に示したように単位円の外に存在する。
【発明の概要】
【発明が解決しようとする課題】
【００１２】
前記のような従来技術は推定した受信チャネルに対するｚ−平面上のポールゼロ（ｐｏｌｅ−ｚｅｒｏ）が単位円の外に存在するため、安全性と因果性を保障できないという問題があり、このような問題を解決しようとすることが本発明の課題である。
【００１３】
したがって、本発明は約束された信号（パイロット信号あるいはトレーニングシーケンス）が挿入されて伝送された信号から送信機と中継機間の伝送路によって設定されたチャネルを推定し、推定されたチャネルから中継機受信チャネルの逆チャネルの情報を有する時間領域でのフィルタ係数を安全性と因果性が保障されるように生成するための、中継機受信チャネルの逆チャネル推定装置およびその方法を提供することにその目的がある。
【００１４】
本発明の目的は以上で言及した目的で制限されず、言及されなかった本発明の他の目的および長所は、下記の説明によって理解することができ、本発明の実施形態によって分明に知り得るものである。また、本発明の目的および長所は、特許請求の範囲に示した手段およびその組合せによって実現される可能性があることを容易に知り得るものである。
【課題を解決するための手段】
【００１５】
前記目的を達成するための本発明の装置は、逆チャネル推定装置において、約束された信号が挿入された受信信号から伝送路によって設定されたチャネルを推定するためのチャネル推定手段と、前記チャネル推定手段で推定したチャネルから最小位相チャネルの逆チャネルを生成して時間領域のフィルタ係数に変換するための逆最小位相フィルタ係数生成手段と、前記チャネル推定手段で推定したチャネルと前記最小位相チャネルを用いて逆グループ遅延成分を有する時間領域のフィルタ係数を生成するための逆グループ遅延フィルタ係数生成手段と、前記逆最小位相フィルタ係数生成手段で変換されたフィルタ係数と前記逆グループ遅延フィルタ係数生成手段で生成したフィルタ係数を結合し、時間領域で受信チャネルの逆チャネルを生成するための結合手段とを備える。
【００１６】
一方、前記目的を達成するための本発明の方法は、逆チャネル推定方法において、約束された信号が挿入された受信信号から伝送路によって設定されたチャネルを推定するチャネル推定ステップと、前記推定したチャネルから最小位相チャネルの逆チャネルを生成して時間領域のフィルタ係数に変換する逆最小位相フィルタ係数生成ステップと、前記推定したチャネルと前記最小位相チャネルを用いて逆グループ遅延成分を有する時間領域のフィルタ係数を生成する逆グループ遅延フィルタ係数生成ステップと、前記変換したフィルタ係数と前記生成したフィルタ係数を結合して時間領域で受信チャネルの逆チャネルを生成する逆チャネル生成ステップとを含む。
【発明の効果】
【００１７】
前記のような本発明は、約束された信号（パイロット信号あるいはトレーニングシーケンス）が挿入されて伝送された信号から送信機と中継機間の伝送路によって設定されたチャネルを推定し、前記推定したチャネルから中継機受信チャネルの逆チャネルの情報を有する時間領域でのフィルタ係数を逆最小位相フィルタ係数（ｈinv、min）と逆グループ遅延フィルタ係数（ｈinv、ap）に分けて生成することによって、安全性と因果性を保障することができるという効果がある。
【図面の簡単な説明】
【００１８】
【図１】一般的な周波数領域での中継機受信チャネルとその逆チャネルに対する一例示図である。
【図２】単位円内に存在する中継機受信チャネルのｚ−平面上のポールゼロに対する一例示図である。
【図３】単位円の外に存在する中継機受信チャネルのｚ−平面上のポールゼロに対する一例示図である。
【図４】本発明による中継機受信チャネルの逆チャネル推定装置に対する一実施形態図である。
【図５】本発明による中継機受信チャネルの逆チャネル推定装置のチャネル推定部に対する一実施形態の詳細図である。
【図６】本発明による中継機受信チャネルの逆チャネル推定装置の最小位相チャネル生成部に対する一実施形態の詳細図である。
【図７】本発明による中継機受信チャネルの逆チャネル推定方法に対する一実施形態のフローチャートである。
【発明を実施するための形態】
【００１９】
上述した目的、特徴および長所は添付された図を参照して詳細に後述されている詳細な説明によってより明確になるものであり、それに応じて本発明が属する技術分野で通常の知識を有する者が本発明の技術的思想を容易に実施できるだろう。また、本発明を説明するにおいて本発明と関連した公知技術に対する具体的な説明が本発明の要旨を不必要にぼかし得ると判断される場合にその詳細な説明を省略するものとする。以下、添付された図を参照して本発明による好ましい実施形態を詳細に説明することにする。
【００２０】
図４は、本発明による中継機受信チャネルの逆チャネル推定装置に対する一実施形態図である。
【００２１】
2010538566.tif 000004
【００２２】
ここで、前記逆最小位相フィルタ係数生成部４２は、最小位相チャネル生成部４２１、逆チャネル生成部４２２、および逆高速フーリエ変換部（ＩＦＦＴ）４２３を備える。
【００２３】
また、前記逆グループ遅延フィルタ係数生成部４３は、除算部４３１、グループ遅延計算部４３２、および逆グループ遅延係数生成部４３３を備える。
【００２４】
以下、本発明による中継機受信チャネルの逆チャネル推定装置の各構成要素に対して詳細に説明する。
【００２５】
まず、チャネル推定部４１は、約束された信号（パイロット信号あるいはトレーニングシーケンス）が挿入された受信信号から送信機と中継機間の伝送路によって設定されたチャネル（Ｈ）を推定する。この時、前記推定したチャネル（Ｈ）は一般的に下記の［数２］のように、最小位相成分とオールパス（ａｌｌ  ｐａｓｓ）成分を含む。
【００２６】
【数２】
2010538566.tif 000005
【００２７】
この時、前記［数２］を簡略化すれば下記の［数３］の通りである。
【００２８】
【数３】
2010538566.tif 000006
【００２９】
ここで、ＨMIN（ｉ）はチャネルの最小位相成分を表し、ＨAP（ｉ）はチャネルのオールパス成分を表す。
【００３０】
その後、最小位相チャネル生成部４２１は、前記推定したチャネル（Ｈ）からｚ−平面上のポールゼロが単位円内に存在する最小位相チャネル（ＨMIN）を生成する。
【００３１】
その後、逆チャネル生成部４２２は、下記の［数４］を用いて前記生成した最小位相チャネル（ＨMIN）から逆チャネル（ＨINV、MIN）を生成する。
【００３２】
【数４】
2010538566.tif 000007
【００３３】
ここで、前記生成した最小位相チャネル（ＨMIN）は、ｚ−平面上のポールゼロが単位円内に存在するため、すなわち安全性と因果性が保障されるため除算を利用した逆変換が可能である。
【００３４】
その後、逆高速フーリエ変換部（ＩＦＦＴ）４２３は、前記周波数領域の最小位相チャネルの逆チャネル（ＨINV、MIN）を時間領域のフィルタ係数（ｈinv、min）に変換する。
【００３５】
一方、除算部４３１は、下記の［数５］を用いて前記推定したチャネル（Ｈ）を前記生成した最小位相チャネル（ＨMIN）に分けてチャネルのオールパス成分（ＨAP）を生成する。
【００３６】
【数５】
2010538566.tif 000008
【００３７】
そして、グループ遅延計算部４３２は、下記の［数６］を用いて前記生成したチャネルのオールパス成分（ＨAP）の位相情報から微分（ｄｉｆｆｅｒｅｎｔｉａｔｉｏｎ）によってグループ遅延成分（ＧAP）を計算する。
【００３８】
【数６】
2010538566.tif 000009
【００３９】
2010538566.tif 000010
【００４０】
2010538566.tif 000011
【００４１】
その後、結合部４４は、下記の［数７］を用いて逆高速フーリエ変換部４２３で逆高速フーリエ変換されたフィルタ係数（ｈinv、min）と、逆グループ遅延係数生成部４３３で生成した逆グループ遅延成分を有する時間領域のフィルタ係数（ｈinv、ap）を結合して時間領域での受信チャネルの逆チャネルを生成する。
【００４２】
【数７】
2010538566.tif 000012
【００４３】
ここで、＊は畳み込み（ｃｏｎｖｏｌｕｔｉｏｎ）演算を意味する。
【００４４】
図５は、本発明による中継機受信チャネルの逆チャネル推定装置のチャネル推定部に対する一実施形態の詳細図である。
【００４５】
図５に示したように、本発明による中継機受信チャネルの逆チャネル推定装置のチャネル推定部は、受信チャネルの歪み推定に用いられる基準信号（Ｘ）を保存している基準信号保存部５０、受信信号から送／受信端の間に約束された信号（Ｙ）を抽出するための約束された信号抽出部５１、および前記約束された信号抽出部５１で抽出した約束された信号（Ｙ）と前記基準信号保存部５０に保存されている基準信号（Ｘ）とを比較して受信チャネルの歪みを推定するためのチャネル歪み推定部５２を備える。
【００４６】
ここで、前記チャネル歪み推定部５２は、下記の［数８］を用いて受信チャネルの歪みを推定する。
【００４７】
【数８】
2010538566.tif 000013
【００４８】
2010538566.tif 000014
【００４９】
図６は、本発明による中継機受信チャネルの逆チャネル推定装置の最小位相チャネル生成部に対する一実施形態の詳細図である。
【００５０】
図６に示したように、本発明による中継機受信チャネルの逆チャネル推定装置の最小位相チャネル生成部は、ログ関数部６１、逆高速フーリエ変換部（ＩＦＦＴ）６２、周波数不変線形フィルタ係数生成部６３、乗算部６４、高速フーリエ変換部（ＦＦＴ）６５、および指数関数部６６を備える。
【００５１】
以下、前記各構成要素に対してもう少し詳細に説明する。
【００５２】
まず、ログ関数部（ログ変換部）６１は、下記の［数９］を用いてチャネル推定部４１で推定したチャネル（Ｈ）をログ変換する。
【００５３】
【数９】
2010538566.tif 000015
【００５４】
2010538566.tif 000016
【００５５】
その後、逆高速フーリエ変換部６２は、前記ログ変換されたチャネル（ＨL）を時間領域の係数（ｃx）に変換する。
【００５６】
そして、周波数不変（ｆｒｅｑｕｅｎｃｙ  ｉｎｖａｒｉａｎｔ）線形フィルタ係数生成部６３は、下記の［数１０］のような周波数不変の線形フィルタ係数を生成する。
【００５７】
【数１０】
2010538566.tif 000017
【００５８】
その後、乗算部６４は、前記逆高速フーリエ変換部６２で逆高速フーリエ変換された係数（ｃx）と周波数不変（ｆｒｅｑｕｅｎｃｙ  ｉｎｖａｒｉａｎｔ）線形フィルタ係数生成部６３で生成した周波数不変線形フィルタ係数（ｌmin）を下記の［数１１］のように乗算する。
【００５９】
【数１１】
2010538566.tif 000018
【００６０】
2010538566.tif 000019
【００６１】
その後、高速フーリエ変換部（ＦＦＴ）６５は、前記乗算部６４からの出力係数（ｘmin）を周波数領域の係数（Ｘmin）に変換する。
【００６２】
その後、指数関数部（指数変換部）６６は高速フーリエ変換された係数（Ｘmin）を下記の［数１２］のように指数変換して最小位相チャネル（ＨMIN）を生成する。
【００６３】
【数１２】
2010538566.tif 000020
【００６４】
図７は、本発明による中継機受信チャネルの逆チャネル推定方法に対する一実施形態のフローチャートである。
【００６５】
まず、約束された信号（パイロット信号あるいはトレーニングシーケンス）が挿入されて伝送された信号から送信機と中継機間の伝送路によって設定されたチャネル（Ｈ）を推定する（７０１）。
【００６６】
その後、前記推定したチャネル（Ｈ）からｚ−平面上のポールゼロが単位円内に存在する最小位相チャネル（ＨMIN）の逆チャネル（ＨINV、MIN）を生成して時間領域のフィルタ係数（ｈinv、min）に変換する（７０２）。
【００６７】
すなわち、前記推定したチャネル（Ｈ）からｚ−平面上のポールゼロが単位円内に存在する最小位相チャネル（ＨMIN）を生成して、前記生成した最小位相チャネル（ＨMIN）から逆チャネル（ＨINV、MIN）を生成した後前記生成した周波数領域の最小位相チャネルの逆チャネル（ＨINV、MIN）を時間領域のフィルタ係数（ｈinv、min）に変換する。
【００６８】
2010538566.tif 000021
【００６９】
2010538566.tif 000022
【００７０】
その後、前記変換されたフィルタ係数（ｈinv、min）と前記生成したフィルタ係数（ｈinv、ap）を結合して時間領域で受信チャネルの逆チャネルを生成する（７０４）。
【００７１】
このような本発明による中継機受信チャネルの逆チャネル推定装置は、放送（ＡＴＳＣ、ＤＶＢ、ＤＭＢ、ＩＳＤＢ−Ｔなど）および通信（ＷｉｂｒｏおよびＣＤＭＡなど）に適しているが、これに限定されるものではなく、一般的な中継機が必要な環境ではどこにでも適用が可能である。
【００７２】
一方、前述したことのような本発明の方法はコンピュータプログラムで作成が可能である。そして前記プログラムを構成するコードおよびコードセグメントは当該分野のコンピュータプログラマーによって容易に推論されることができる。また、前記作成されたプログラムはコンピュータが読み込み可能な記録媒体（情報保存媒体）に保存され、コンピュータによって判読されて実行されることによって本発明の方法を具現する。そして前記記録媒体はコンピュータが判読可能なすべての形態の記録媒体を含む。
【００７３】
以上で説明した本発明は、本発明が属する技術分野で通常の知識を有する者にあって本発明の技術的思想に外れない範囲内で色々な置換、変形および変更が可能なため前述した実施形態および添付された図面によって限定されるものではない。
【符号の説明】
【００７４】
１１  受信チャネル
１２  逆チャネル
４１  チャネル推定部
４２  逆最小位相フィルタ係数生成部
４２１  最小位相チャネル生成部
４２２  逆チャネル生成部
４２３  逆高速フーリエ変換部
４３  逆グループ遅延フィルタ係数生成部
４３１  除算部
４３２  グループ遅延計算部
４３３  逆グループ遅延係数生成部
５０  基準信号保存部
５１  約束された信号抽出部
５２  チャネル歪み推定部
６１  ログ関数部
６２  逆高速フーリエ変換部（ＩＦＦＴ）
６３  周波数不変線形フィルタ係数生成部
６４  乗算部
６５  高速フーリエ変換部（ＦＦＴ）
６６  指数関数部
【図１】
2010538566.tif 000023
【図２】
2010538566.tif 000024
【図３】
2010538566.tif 000025
【図４】
2010538566.tif 000026
【図５】
2010538566.tif 000027
【図６】
2010538566.tif 000028
【図７】
2010538566.tif 000029
【国際調査報告】
2010538566.tif 000030
2010538566.tif 000031