(57)【特許請求の範囲】
【請求項１】
ＢＰＳＫ復調されたディジタル信号に含まれる「０」および「１」の少なくとも一方の数を検出する検出手段と、
前記検出手段による検出結果に応じて前記ＢＰＳＫ復調されたディジタル信号の復調前の信号に振幅変調された信号が含まれていたか否かを判別する判別手段と、を備え、
前記判別手段は、
前記検出手段により検出された「０」の数、又は「１」の数が所定値以上の場合、前記復調前の信号に前記振幅変調された信号が含まれていたと判別する、データ識別装置。
【請求項２】
前記判別手段は、
前記検出手段により検出された「０」の数、及び、「１」の数のいずれもが前記所定値未満の場合、前記復調前の信号はＢＰＳＫ変調信号であったと判別する、請求項１記載のデータ識別装置。
【請求項３】
前記振幅変調された信号は、前記ＢＰＳＫ変調信号に直交することを特徴とする、請求項２記載のデータ識別装置。
【請求項４】
前記振幅変調された信号はＡＲＩ信号であることを特徴とする、請求項１記載のデータ識別装置。
【請求項５】
前記振幅変調された信号はＡＲＩ信号であり、前記ＢＰＳＫ変調信号はＲＤＳ信号であることを特徴とする、請求項２記載のデータ識別装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明はデータ識別装置に関し、特にたとえば、ＢＰＳＫ変調信号の復調信号に所定データが含まれるか否かを識別するデータ識別装置に関する。
【０００２】
【従来の技術】
ＢＰＳＫ変調信号を復調するためのコスタスループの一例が、「スペクトラム拡散通信方式」（R.C.Dixon 著、立野敏他訳、ジャテック出版）のp.１６６〜p.１６８に記載されている。このコスタスループ１を図１６を参照して説明すると、ＢＰＳＫ変調信号は、乗算器２ａでＶＣＯ(Voltage Controlled Oscillator)３からの発振周波数信号と乗算されるとともに、乗算器２ｂで、９０°移相回路４によって移相された発振周波数信号と乗算される。乗算器２ａおよび２ｂからの出力は、その後低域通過フィルタ（ＬＰＦ）５ａおよび５ｂで高周波成分を除去され、ＬＰＦ５ａからの出力がＢＰＳＫ復調信号として出力されるとともに、ＬＰＦ５ａおよび５ｂの出力が乗算器２ｃで乗算される。そして、乗算器２ｃからの出力がループフィルタ６を経てＶＣＯ３に与えられ、これによってＶＣＯ３の発振周波数が制御される。
【０００３】
ＢＰＳＫ変調信号は±Ａ cos（ωｔ）で表され、ＶＣＯ３の出力は cos（ωｔ＋Δ）と表される。なおΔは移相誤差である。これによって、乗算器２ａの出力は±Ａ（ cosΔ＋ cos（２ωｔ＋Δ））／２となり、これをＬＰＦ５ａに通すと、高周波成分が除去され、±Ａ cosΔ／２が残る。この成分は、Δ→０となることによって±Ａ／２となり、復調データとなる。一方、ＬＰＦ５ｂの出力は、上述と同様に考えて± sinΔ／２となる。ＬＰＦ５ａおよび５ｂの出力を乗算器２ｃで乗算すると、乗算結果は sin（２Δ）／４となる。したがって、Δが小さいときは sin（２Δ）／４≒Δ／２と近似され、ＶＣＯ３をΔ→０となるように制御することができる。なお、ループフィルタ６には制御系の定数が設定される。
【０００４】
【発明が解決しようとする課題】
このようなコスタスループ１はアナログ回路で構成されることが多いため、集積回路化することが困難であった。また、アナログ回路であるため、煩雑な回路調整が必要となるとともに、回路が温度変化の影響を受け易いという問題点もあった。
【０００５】
それゆえに、この発明の主たる目的は、ディジタル回路で構成されたコスタスループに適したデータ識別回路を提供することである。
【０００６】
【課題を解決するための手段】
本発明は、ＢＰＳＫ復調されたディジタル信号に含まれる「０」および「１」の少なくとも一方の数を検出する検出手段と、前記検出手段による検出結果に応じて前記ＢＰＳＫ復調されたディジタル信号の復調前の信号に振幅変調された信号が含まれていたか否かを判別する判別手段と、を備え、前記判別手段は、前記検出手段により検出された「０」の数、又は「１」の数が所定値以上の場合、前記復調前の信号に前記振幅変調された信号が含まれていたと判別する、データ識別装置である。
【０００７】
本発明では、たとえばＲＤＳ(Radio Data System) 信号のＢＰＳＫ変調信号の復調されたＲＤＳ信号の所定長に含まれる「０」および「１」が、たとえば第１カウンタおよび第２カウンタでカウントされ、第１カウンタおよび第２カウンタのいずれか一方のカウント値が所定値より大きいかどうか判別回路で判別される。そして、所定値より大きければ、ＲＤＳ信号にたとえばＡＲＩ(Autofahrer Rundfunk Imformations)信号が含まれるとして、判別回路からたとえばハイレベル信号が出力される。
【０００８】
【発明の実施の形態】
図１を参照して、本発明の実施の形態のコスタスループ１０は、コンパレータ１２を含む。このコンパレータ１２によってＲＤＳ信号のＢＰＳＫ変調信号が基準電圧と比較されることによって、ＢＰＳＫ変調信号が２値化される。２値化されたＢＰＳＫ変調信号は、Ｄ−ＦＦ回路１４ａおよび１４ｂに与えられる。Ｄ−ＦＦ回路１４ａに与えられたＢＰＳＫ変調信号は、ＶＣＯ１６から出力された発振周波数信号によってラッチされ、Ｄ−ＦＦ回路１４ｂに与えられたＢＰＳＫ変調信号は、ＶＣＯ１６から出力されかつ９０°移相器１８で移相された発振周波数信号によってラッチされる。Ｄ−ＦＦ回路１４ａからの出力信号は、復調信号（ＲＤＳ信号）として出力されるとともに位相比較器２０に与えられ、Ｄ−ＦＦ回路１４ｂからの出力は、位相比較器２０に与えられる。位相比較器２０にはまた、ＶＣＯ１６から出力された発振周波数信号，９０°移相器１８から出力された発振周波数信号および２値化されたＢＰＳＫ変調信号が与えられる。これらの信号に基づいて位相比較器２０から位相比較結果が出力され、これがループフィルタ２２を経てＶＣＯ１６に与えられ、これによってＶＣＯ１６の発振周波数が制御される。
【０００９】
なお、この実施の形態では、ＶＣＯ１６から第１発振周波数信号を取り出し、９０°移相器１８から第２発振周波数信号を取り出すようにした。しかしながら、１つのＶＣＯからの発振周波数信号をたとえば分周器に与え、この分周器から第１発振周波数信号および第２発振周波数信号を取り出すようにしてもよい。
【００１０】
位相比較器２０の構成を図２に示す。位相比較器２０は、図５（Ｃ）に示す高速パルスによってインクリメントされかつ図５（Ｂ）に示す９０°移相器１８からの発振周波数信号によってリセットされるカウンタ２４ａを含み、このカウンタ２４ａのカウント値がＤ−ＦＦ回路２６ａに与えられる。Ｄ−ＦＦ回路２６ａにはまた、図５（Ａ）に示すＢＰＳＫ変調信号のエッジのうち９０°移相器１８から発振周波数信号が入力された直後のエッジのみを選択するエッジ選択器２８からの選択信号が与えられ、これによってカウント値がラッチされる。したがって、ラッチされたカウント値は発振周波数信号とその発振周波数信号直後のＢＰＳＫ変調信号のエッジとの間の高速パルス数（この場合は「２」）となる。なお、発振周波数信号の入力直後のエッジのみを選択することとしたのは、ＢＰＳＫ変調信号の周波数は発振周波数の２倍であり、カウント値のラッチに用いることができないエッジが存在するからである。Ｄ−ＦＦ回路２６ａでラッチされたカウント値は加算器２９で定数０．５と加算され、加算結果が選択回路３０に与えられる。なお、加算器２８で定数０．５を加算するのは、ＢＰＳＫ変調信号のエッジと高速パルスとは同期しておらず、エッジは平均的に高速パルスと高速パルスとの真中に存在すると考えられるからである。
【００１１】
位相比較器２０はまた、図６（Ｃ）に示す高速パルスによってインクリメントされかつ図６（Ａ）に示すＢＰＳＫ変調信号によってリセットされるカウンタ２４ｂを含み、このカウンタ２４ｂのカウント値が、Ｄ−ＦＦ回路２６ｂで図６（Ｂ）に示す９０°移相器１８からの発振周波数信号によってラッチされる。したがって、ラッチされたカウント値は、ＢＰＳＫ変調信号のエッジと当該エッジ直後の９０°移相器１８からの発振周波数信号との間の高速パルス数（この場合は「３」）となる。そして、減算器３２で、ラッチされたカウント値から０．５が減算され、減算された値の符号が符号反転回路３４で反転され、反転された値が選択回路３０に与えられる。なお、減算器３２で０．５を減算するのは上述と同じ理由による。また、符号を反転させるのは、上述の図５（Ａ）〜（Ｃ）の場合と比べて位相制御を逆にする必要があるからである。
【００１２】
一方、ＥＸＯＲ回路３６にはラッチ回路１４ａおよび１４ｂからの出力信号が与えられ、排他的論理和をとられる。したがって、図７（Ａ）および（Ｃ）からわかるように発振周波数信号の位相がＢＰＳＫ変調信号よりも進んでいるときは、ＥＸＯＲ回路３６からは「０」が出力され、図７（Ａ）および（Ｄ）からわかるように発振周波数の位相がＢＰＳＫ変調信号よりも遅れているときは、ＥＸＯＲ回路３６からは「１」が出力される。そしてＥＸＯＲ回路３６の出力が、Ｄ−ＦＦ回路２６ｃでＶＣＯ１６からの発振周波数信号によってラッチされる。Ｄ−ＦＦ回路２６ｃの出力はその後選択回路３０に与えられ、これによって加算器２８および符号反転器３４から選択回路３０に与えられた信号のいずれかが選択される。すなわち、Ｄ−ＦＦ回路３８からの出力が「０」であれば加算器２８からの出力信号が出力され、Ｄ−ＦＦ回路３８からの出力が「１」であれば符号反転器３４からの出力信号が出力される。選択回路３０からの出力信号はその後、ＶＣＯ１６からの発振周波数信号によってラッチされるＤ−ＦＦ回路２６ｄを経て位相比較結果として出力される。
【００１３】
ループフィルタ２２の構成を図３に示す。入力された位相比較結果は加算器３８ａおよび割算器４０に与えられる。割算器４０では位相比較結果が１／Ｒ２され、割算器４０からの出力が加算器３８ｂに与えられる。加算器３８ｂにはまた、ＶＣＯ１６からの発振周波数信号によって加算器３８ｂからの加算結果をラッチするＤ−ＦＦ回路４２の出力が与えられ、これによって割算器４０の出力が積分される。Ｄ−ＦＦ回路４２からの出力はその後加算器３８ａで位相比較結果と加算され、加算器３８ａから加算結果の符号のみが出力される。
【００１４】
このループフィルタ２２をアナログループフィルタより類推して説明する。アナログループフィルタの伝達関数を
【００１５】
【数１】
0003702281.tif 000002
【００１６】
とする。ループフィルタ２２はディジタルフィルタであるため、数１を双１次変換すると、伝達関数は数２で表される。
【００１７】
【数２】
0003702281.tif 000003
【００１８】
なお、Ｔはサンプリング周期である。ここで、分子の１＋Ｚ-1は、或るサンプリング時刻の値とその次のサンプリング時刻の値の足し算を意味するが、その値の差は少ないと考え、１＋Ｚ-1≒２とする。その結果、数２は数３で表される。
【００１９】
【数３】
0003702281.tif 000004
【００２０】
ここで、Ｒ２／Ｒ１＝１およびＴ／Ｒ１Ｃ＝１／Ｒ２となるように定数を決めると、１／（１−Ｚ-1）は入力を順次積分するディジタル積分となるので、ループフィルタ２２では入力と入力をＲ２で割ったもののディジタル積分の出力とが加算される。
【００２１】
ＶＣＯ１６の構成を図４に示す。選択回路４４には定数１９０１９および１８９８１が与えられ、ループフィルタ２２からの出力によってこの定数のいずれか一方が選択される。すなわち、ループフィルタから出力される符号が＋であれば定数１９０１９が選択され、符号が−であれば定数１８９８１が選択される。一方、定数２５０は加算器４６に与えられ、ここでＤ−ＦＦ回路４８からの出力と定数２５０とが加算される。選択回路４４による選択結果および加算回路４６による加算結果は演算器５０に与えられ、演算器５０で両者が比較される。そして、選択結果＞加算結果であれば端子Ａから加算結果がそのまま出力されるとともに端子Ｂからローレベル信号が出力され、選択結果≦加算結果であれば、端子Ａから加算結果−選択結果が出力されるとともに端子Ｂからハイレベル信号が出力される。なお、Ｄ−ＦＦ回路４８は高速パルスによって演算器５０からの出力をラッチする。
【００２２】
ＶＣＯ１６の動作原理を説明する。定数ＮをモジュロＮ／ｆで順次加算することとすると、桁上げが起きる周波数はｆとなる。したがって、モジュロ（Ｎ／ｆ＋Δ）およびモジュロ（Ｎ／ｆ−Δ）を用意し、入力に従ってどちらかを選択するようにすると、１／（１／ｆ＋Δ／Ｎ）≒ｆ−Δ／Ｎおよび１／（１／ｆ−Δ／Ｎ）≒ｆ＋Δ／Ｎの２種類の周波数を得ることができる。位相を進ませたい場合はｆ＋Δ／Ｎを選択し、遅らせたい場合はｆ−Δ／Ｎを選択すればよい。
【００２３】
ＶＣＯ１６の動作を、データ転送速度を１ｂｐｓ、キャリア周波数を２４Ｈｚ、高速パルスを１８２４Ｈｚとして説明する。高速パルスの周波数は１８２４Ｈｚなので、選択回路４４で定数１９０１９が選択されたときは、発振周波数は１８２４×２５０／１９０１９≒２３．９７６Ｈｚとなり、選択回路で定数１８９８１が選択されたときは、発振周波数は１８２４×２５０／１８９８１≒２４．０２４Ｈｚとなる。したがって、位相を進ませるときは定数１８９８１を選択し、位相を遅らせるときは１９０１９を選択すればよい。なお、９０°移相器１８は遅延回路で構成され、１／（２４×４）秒の遅延で２４Ｈｚに対する９０°移相となる。
【００２４】
このように、ＶＣＯ１６の発振周波数を制御することによって、図８（Ａ）に示す２値化ＢＰＳＫ変調信号のエッジと図８（Ｂ）に示す発振周波数信号の９０°移相成分との位相タイミングが合えば、ＶＣＯ１６からの発振周波数信号の位相タイミングはＢＰＳＫ変調信号のキャリアの半波の中央と一致する。したがって、図８（Ａ）に示すＢＰＳＫ変調信号をＤ−ＦＦ回路１４ａに与え、ＶＣＯ１６からの発振周波数信号でラッチすれば、図８（Ｃ）に示すように正しいデータが復調される。
【００２５】
この実施形態によれば、コスタスループ１０をディジタル回路によって構成することができるため、コスタスループ１０を集積回路化することができ、煩雑な回路調整が不要となるとともに温度変化の影響を抑えることができる。さらに、従来のアナログによるコスタスループのようにＬＰＦを設ける必要がないため、部品点数を少なくすることができる。
【００２６】
なお、この実施形態では、図３に示す割算器４０の値を１／Ｒ２として説明したが、この割算値が１／２M （Ｍは整数）であれば、図９に示すように割算器４０の出力をＭＳＢよりｍビットずらして入力すればよい。また、図３において加算器３８ｂの出力の絶対値が大きくなりすぎると、オーバフローする可能性が出てくる。したがって、オーバフローをなくすためには、図１０に示すように加算器３８ｂの出力の絶対値がオーバフローしたことを検出するオーバフロー検出回路５２およびゲート５４を設け、オーバフロー検出回路５２でオーバフローが検出されたときは、発振周波数信号がゲート５４を通過しないようにすればよい。
【００２７】
さらにまた、位相比較器２０の構成を簡略化するために、位相比較器２０を図１１に示すようにＥＸＯＲ回路３６，Ｄ−ＦＦ回路２６ｃおよびＤ−ＦＦ回路２６ｃの出力「０」を＋１に変換し、「１」を−１に変換する変換回路５６によって構成し、ＥＸＯＲ回路３６に上述と同様にラッチ回路１４ａおよび１４ｂからの出力を与え、変換回路５６の出力に基づいてＶＣＯ１６を制御するようにしてもよい。すなわち、Ｄ−ＦＦ回路２６ｃの出力が「０」であれば、発振周波数信号の位相が進んでいるとして、変換回路５６の出力「＋１」に基づいて位相を遅らせ、Ｄ−ＦＦ回路２６ｃの出力が「−１」であれば、発振周波数信号の位相が遅れているとして、変換回路５６の出力「−１」に基づいて位相を進ませればよい。
【００２８】
ＡＲＩ信号はＲＤＳ信号の３倍の振幅をもち、ＲＤＳ信号とは直交の関係にある。したがって、両信号は数４で表される。
【００２９】
【数４】
0003702281.tif 000005
【００３０】
ここで±Ａ cos（ωｔ）がＲＤＳ信号であり、３Ａ sin（ωｔ）がＡＲＩ信号である。このような式で表される信号を図１６に示す従来のコスタスループ１に与えると、ＶＣＯ３の出力は sin（ωｔ）（または− sin（ωｔ））となり、乗算器２ｂの出力がデータとなる。したがって、ＡＲＩ信号がないときは乗算器２ａの出力からデータをとり、ＡＲＩ信号があるときは乗算器２ｂの出力からデータをとる必要があり、そのためにＡＲＩ信号の有無の判定が必要となる。そして、ＡＲＩ信号の有無の判定にあたっては、従来は乗算器２ａの出力レベルによって判定していた。すなわち、ＡＲＩ信号がないときは乗算器２ａの出力は±Ａであり、ＡＲＩ信号がないときは乗算器２ａの出力は±３Ａであるため、このレベルによってＡＲＩ信号の有無を判定していた。
【００３１】
しかし、図１に示す実施形態のコスタスループ１０では、ＢＰＳＫ変調信号をコンパレータ１２によって２値化するため、Ｄ−ＦＦ回路１４ａの出力レベルが一定になってしまい、レベルによってＡＲＩ信号の有無を判定することはできなかった。図１２に示すこの実施の形態例のデータ識別装置６０は、このような問題点を解決するためのものである。
【００３２】
図１に示すコスタスループ１０から出力された図１３（Ａ）に示す復調信号は、クロック再生回路６２に与えられ、これによってクロック再生回路６２から復調信号に同期した図１３（Ｂ）に示すようなクロックが再生される。クロック再生回路６２から再生されたクロックはパルス発生回路６４に与えられ、パルス発生回路６４からはクロックが４回与えられる毎に図１３（Ｃ）に示すパルスが出力される。クロック再生回路６２からのクロックとパルス発生回路６４からのパルスとは、カウンタ６６ａおよび６６ｂに与えられ、カウンタ６６ａおよび６６ｂは、クロックによってインクリメントされるとともに、パルスによってリセットされる。また、復調信号はカウントイネーブル信号としてカウンタ６６ａおよび６６ｂに与えられる。そして、カウンタ６６ａでは復調信号が「１」であるときカウントがイネーブルされ、カウンタ６６ｂでは復調信号が「０」であるときカウントがイネーブルされる。
【００３３】
カウンタ６６ａおよび６６ｂから出力されたカウント値は、Ｄ−ＦＦ回路６８ａおよび６８ｂでパルス発生回路６４からのパルスによってラッチされる。したがって、Ｄ−ＦＦ回路６８ａおよび６８ｂでそれぞれラッチされたカウント値は、パルスとパルスとの間における「１」および「０」の数となる。Ｄ−ＦＦ回路６８ａおよび６８ｂでラッチされたカウント値はそれぞれ判定器７０ａおよび７０ｂに与えられ、判定器７０ａおよび７０ｂで「１」および「０」の数が所定値以上であるかどうか判定される。そして、所定値以上であればハイレベルの信号が出力され、ＯＲ回路７２に与えられる。したがって、「１」および「０」のいずれか一方の数が所定値より大きければ、ＯＲ回路７２からハイレベル信号が出力され、「１」および「０」の数のいずれもが所定値よりも小さければＯＲ回路７２からはローレベル信号が出力される。
【００３４】
このように、復調信号に含まれる「１」および「０」の数に応じてデータを識別することとしたのは、復調信号にＡＲＩ信号が含まれていなければ復調信号の「０」と「１」の数は同じであるが、復調信号にＡＲＩ信号が含まれるときは、その部分の復調信号は「０」か「１」だけとなるからである。ただし、実際にはデータの誤りがあるため、復調信号に含まれる「１」または「０」の数が所定値以上であればＡＲＩ信号が含まれていると判定する。
【００３５】
図１４を参照して、他の実施の形態例のデータ識別装置６０は、シフトレジスタ７４を含む。このシフトレジスタ７４はクロック再生回路６２から再生されたクロックに従ってデータを取り込むとともに取り込んだデータをシフトさせる。論理回路７６は、シフトレジスタ内の「０」および「１」の数に応じてハイレベルまたはローレベルの信号を出力する。すなわち、図１５に示すようにシフトレジスタ７４内のデータに応じた信号、つまりデータが「１」のときハイレベルとなりデータが「０」のときローレベルとなる信号が、ＡＮＤ回路７８およびＮＡＮＤ回路８０に与えられ、ＡＮＤ回路７８およびＮＡＮＤ回路８０の出力信号がＯＲ回路８２を経て出力される。したがって、シフトレジスタ７４内の数字が全て「１」であるかまたは全て「０」であるとき論理回路７６からハイレベル信号が出力され、１つでも異なる数字が含まれていれば論理回路７６からローレベル信号が出力される。
【００３６】
これらの実施形態によれば、ディジタル回路だけでＡＲＩ信号の有無が判定できるので、ディジタル方式のコスタスループに適用でき、またデータ識別装置を集積回路化することができる。なお、この実施形態では、ＲＤＳ信号およびＡＲＩ信号を用いて説明したが、この発明はこの場合に限らず、信号に含まれる「１」および「０」の数によって識別されるあらゆる信号に適用できることはもちろんである。
【００３７】
【発明の効果】
本発明によれば、復調信号の所定長に含まれる「０」および「１」の数によってデータが識別されるので、復調された信号に含まれるデータを識別することができる。
【図面の簡単な説明】
【図１】本発明の実施形態を説明するブロック図である。
【図２】図１装置の一部を示すブロック図である。
【図３】図１装置の他の一部を示すブロック図である。
【図４】図１装置のその他の一部を示すブロック図である。
【図５】図１装置の動作の一部を示すタイミング図である。
【図６】図１装置の動作の一部を示すタイミング図である。
【図７】（Ａ）はＢＰＳＫ変調信号を示す波形図であり、（Ｂ）〜（Ｄ）は発振周波数信号を示す波形図である。
【図８】（Ａ）はＢＰＳＫ変調信号を示す波形図であり、（Ｂ）は発振周波数信号を示す波形図であり、（Ｃ）は復調信号を示す波形図である。
【図９】本発明の実施形態の他の一部を示すブロック図である。
【図１０】本発明の実施形態の他の一部を示すブロック図である。
【図１１】本発明の実施形態の他の一部を示すブロック図である。
【図１２】本発明の実施形態のデータ識別装置を示すブロック図である。
【図１３】（Ａ）はＢＰＳＫ変調信号を示す波形図であり、（Ｂ）は再生クロックを示す波形図であり、（Ｃ）はパルスを示す波形図である。
【図１４】本発明の実施形態の他のデータ識別装置を示すブロック図である。
【図１５】図１２装置の一部を示すブロック図である。
【図１６】従来技術を示すブロック図である。
【符号の説明】
１０  …コスタスループ
１２  …コンパレータ
１４ａ，１４ｂ，２６ａ〜２６ｄ，４２，４８，６８ａ，６８ｂ  …Ｄ−ＦＦ回路
１６  …ＶＣＯ
１８  …９０°移相器
２０  …位相比較器
２２  …ループフィルタ
６０  …データ識別装置
６２  …クロック再生回路
６４  …パルス発生回路
６６ａ，６６ｂ  …カウンタ
【図１】
0003702281.tif 000006
【図２】
0003702281.tif 000007
【図３】
0003702281.tif 000008
【図４】
0003702281.tif 000009
【図５】
0003702281.tif 000010
【図６】
0003702281.tif 000011
【図７】
0003702281.tif 000012
【図８】
0003702281.tif 000013
【図９】
0003702281.tif 000014
【図１０】
0003702281.tif 000015
【図１１】
0003702281.tif 000016
【図１２】
0003702281.tif 000017
【図１３】
0003702281.tif 000018
【図１４】
0003702281.tif 000019
【図１５】
0003702281.tif 000020
【図１６】
0003702281.tif 000021