(57)【要約】システム及び方法は、３個のパリティチェックビットで伝送されるデータの各ｎ−３ビットを符号化してｎビットのブロック（ｎ−３個の情報ビットに、ｎ個の情報ビットに対応付けられた３個のパリティビットを加えたもの）を生成するためのパリティビット符号器（２０２）を含む。ｎビットのブロックの各々は、光波長に変調されて受信機に伝送される３個のＱＡＭシンボルにグレイマッピングされる。サイクルスリップを訂正するのに最大事後確率（ＭＡＰ）復号器が受信機で使用される。【選択図】図２
2016533127.tif 000002
【特許請求の範囲】
【請求項１】
ｎ−３ビットの連続するブロックを３個のパリティビットで符号化してｎビットの連続するブロックを供給するように構成されたパリティビット符号器（２０２）、
前記パリティビット符号器（２０２）に結合され、前記ｎビットのブロックの各々を３個の対応の直交振幅変調（ＱＡＭ）シンボルにマッピングするように構成されたグレイマッパ（２０４）、
前記グレイマッパ（２０４）に結合され、前記グレイマッパの出力に応じて光信号を変調して、前記３個の対応のＱＡＭシンボルを含む変調光信号を供給するように構成された変調器（２０６）、
前記変調光信号を受信し、前記光信号を表す電気信号を供給するための検出器（３０２）、及び
前記電気信号に応じてｎ−３ビットの前記ブロックを表すデマッパ出力を供給するように構成されたデマッパ（３０４、３０４ａ）
を備えたシステム。
【請求項２】
Ｐ１、Ｐ２及びＰ３をそれぞれ前記３個のパリティビットの第１、第２及び第３のものとし、ｘ１、ｘ２、・・・ｘｎ−３をそれぞれ前記ｎ−３ビットとして、前記３個のパリティビットが、以下の式：
【数１】
2016533127.tif 000008
【数２】
2016533127.tif 000009
【数３】
2016533127.tif 000010
及び
【数４】
2016533127.tif 000011
に従って計算される、請求項１に記載のシステム。
【請求項３】
前記パリティビット符号器（２０２）が以下の：
Ｐ１が最初の２ｎ／３ビットのいずれかに位置決めされること、
Ｐ２が２ｎ／３＋１番目のビットと５ｎ／６番目のビットとの間のいずれかに位置決めされること、及び
Ｐ３が２ｎ／３＋１番目のビットと、前記Ｐ２の位置とは異なるｎ番目のビットとの間のいずれかに位置決めされること  
に従って前記ｎビット内に前記３個のパリティビットを位置決めするように構成された請求項２に記載のシステム。
【請求項４】
前記デマッパ（３０４、３０４ａ）が、前記３個の対応のＱＡＭシンボルのうちの第１のもののハミング重みに基づいて、前記３個の対応のＱＡＭシンボルのうちの第２のものを選択することによって前記デマッパ出力を供給するように構成された請求項１に記載のシステム。
【請求項５】
Ｍを前記３個の対応のＱＡＭシンボルの各シンボルによって表されるビットの数として、前記デマッパ（３０４、３０４ａ）が、前記３個の対応のＱＡＭシンボルのうちの前記第１及び第２のものの｛ｌｏｇ２（Ｍ２）｝／２個の最大有意ビットハミング重みに基づいて前記３個の対応のＱＡＭシンボルのうちの第３のものを選択することによって前記デマッパ出力を供給するように構成された請求項４に記載のシステム。
【請求項６】
前記デマッパ（３０４、３０４ａ）がさらに、前記３個のパリティビットによって示されるパリティを用いてサイクルスリップの訂正をもたらすように構成された請求項１に記載のシステム。
【請求項７】
シリアル入力データストリームを受信し、前記シリアル入力データストリームを複数の多重分離データストリームに多重分離するように構成されたデマルチプレクサ（７０１）、
複数の前方誤り訂正（ＦＥＣ）符号器（７０２−１・・・７０２−（ｎ−３））であって、該ＦＥＣ符号器の各々が、ＦＥＣ符号を用いて前記多重分離データストリームの対応のものを符号化し、対応のＦＥＣ符号化出力を供給するように構成された、複数のＦＥＣ符号器（７０２−１・・・７０２−（ｎ−３））、及び
前記複数のＦＥＣ符号器（７０２−１・・・７０２−（ｎ−３））に結合され、前記ｎ−３ビットの連続するブロックを含むインタリーブ出力を前記パリティビット符号器（２０２）に供給するように構成されたインタリーバ（７０４）
をさらに備えた請求項１に記載のシステム。
【請求項８】
前記デマッパに結合され、前記デマッパ出力に応じて複数のデインタリーブ出力を供給するように構成されたデインタリーバ（８０２）、及び
複数のＦＥＣ復号器（８０４−１・・・８０４−（ｎ−３））であって、該ＦＥＣ復号器（８０４−１・・・８０４−（ｎ−３））の各々が、前記デインタリーブ出力のうちの対応のものを受信し、対応のＦＥＣ復号出力信号を供給するように構成された複数のＦＥＣ復号器（８０４−１・・・８０４−（ｎ−３））
をさらに備えた請求項１に記載のシステム。
【請求項９】
３個の対応の直交振幅変調（ＱＡＭ）信号にマッピングされるｎ−３個の情報ビット及び３個のパリティビットからなるｎビットのブロックを表す変調光信号を受信し、前記光信号を表す電気信号を供給するための検出器（３０２）、及び
前記電気信号に応じて、前記ｎ−３個の情報ビットを表すデマッパ出力を供給するように構成されたデマッパ（３０４、３０４ａ）
を備えた光信号受信機。
【請求項１０】
前記デマッパ（３０４、３０４ａ）が、前記３個の対応のＱＡＭシンボルのうちの第１のもののハミング重みに基づいて、前記３個の対応のＱＡＭシンボルのうちの第２のものを選択することによって前記デマッパ出力を供給するように構成された請求項９に記載の光信号受信機。
【請求項１１】
Ｍを前記３個の対応のＱＡＭシンボルの各シンボルによって表されるビットの数として、前記デマッパ（３０４、３０４ａ）が、前記３個の対応のＱＡＭシンボルのうちの前記第１及び第２のものの｛ｌｏｇ２（Ｍ２）｝／２個の最大有意ビットハミング重みに基づいて前記３個の対応のＱＡＭシンボルのうちの第３のものを選択することによって前記デマッパ出力を供給するように構成された請求項１０に記載の光信号受信機。
【請求項１２】
前記デマッパ（３０４、３０４ａ）がさらに、前記３個のパリティビットによって示されるパリティを用いてサイクルスリップの訂正をもたらすように構成された請求項９に記載の光信号受信機。
【請求項１３】
前記デマッパ（３０４、３０４ａ）に結合され、前記デマッパ出力に応じて複数のデインタリーブ出力を供給するように構成されたデインタリーバ（８０２）、及び
複数のＦＥＣ復号器（８０４−１・・・８０４−（ｎ−３））であって、該ＦＥＣ復号器（８０４−１・・・８０４−（ｎ−３））の各々が、前記デインタリーブ出力のうちの対応のものを受信し、対応のＦＥＣ復号出力信号を供給するように構成された複数のＦＥＣ復号器（８０４−１・・・８０４−（ｎ−３））
をさらに備えた請求項９に記載の光信号受信機。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、光信号データ検出に関し、より詳細には、カスケード単一パリティチェック符号化を用いるシステム及び方法に関する。
【背景技術】
【０００２】
波長分割多重（ＷＤＭ）光通信システムでは、多数の異なる光キャリア波長がデータによって個別に変調されて変調光信号を生成する。変調光信号は統合信号に合波され、光伝送経路を介して受信機に伝送される。受信機は、データを検出及び復調する。
【０００３】
光通信システムで使用され得る変調の一種として、位相シフトキーイング（ＰＳＫ）がある。様々なバリエーションのＰＳＫによると、光波長の位相又は位相遷移が、１以上のビットを符号化するシンボルを表すように、光波長の位相を変調することによってデータが伝送される。バイナリ位相シフトキーイング（ＢＰＳＫ）変調法では、例えば、２位相が用いられてシンボルあたり１ビットを表す。直交位相シフトキーイング（ＱＰＳＫ）変調法では、４位相が用いられてシンボルあたり２ビットを符号化する。他の位相シフトキーイング形式として、差動位相シフトキーイング（ＤＰＳＫ）形式、並びにリターントゥゼロＤＰＳＫ（ＲＺ−ＤＰＳＫ）及び偏波多重ＱＰＳＫ（ＰＤＭ−ＱＰＳＫ）のようなＰＳＫ及びＤＰＳＫ形式の変形がある。
【０００４】
複数の情報ビットが単一の伝送シンボルに符号化されることになるＱＰＳＫなどの変調形式は、一般にマルチレベル変調形式ともいわれる。マルチレベル変調技術は、例えば、伝送速度を増加させ、チャネル間隔を減少させ、それによりＷＤＭシステムの各チャネルのスペクトル効率（ＳＥ）を増加させることを可能とするのに使用されてきた。１つのスペクトル効率の良いマルチレベル変調形式は直交振幅変調（ＱＡＭ）である。ＱＡＭ信号では、例えば、シンボルあたり複数ビットを符号化するのに、位相シフトキーイング及び振幅シフトキーイングの組合せを用いて情報が変調される。シンボルあたりＭビットを符号化するのにＭ２−ＱＡＭ信号を用いることができる。例えば、シンボルあたり４ビットを符号化するのに１６−ＱＡＭ変調形式を用いることができる。ＰＳＫ変調法（例えば、ＢＰＳＫ及びＱＰＳＫ）は、ＱＡＭの１つのレベル（例えば、それぞれ２ＱＡＭ及び４ＱＡＭ）ともいえる。
【０００５】
より高い次数のＱＡＭ法は、異なるＳＥを用いて複数のデータレートが同一帯域幅（又はシンボルレート）にわたって達成され得るような変動レートのトランスポンダ技術を実現する際に有用である。複数のＳＥは、設計における冗長性の量を変動させることによってＱＡＭを用いて実現でき、例えば、冗長性が高いほどＳＥは低くなる。冗長性は、その手法において用いられる前方誤り訂正符号（ＦＥＣ）のオーバーヘッド（ＯＨ）を変更することから、又は主に符号化された変調の目的のために用いられる符号化の簡素な追加レイヤを導入することによって生ずる。符号化された変調は、シンボルレート、コンステレーショングリッド及びＦＥＣオーバーヘッドが固定されている場合に特に有用である。
【０００６】
例えば、ＱＡＭ法を用いる位相変調光通信システムでは、受信機は、コヒーレント検出、例えばホモダイン又はヘテロダイン検出を用いて変調光信号を検出するコヒーレント受信機となり得る。受信機に関してここで使用される場合の用語「コヒーレント」とは、受信信号を復調するためのローカル発振器（ＬＯ）を含む受信機のことをいう。受信信号を処理して復調データを供給するようなシステムでは、デジタル信号処理（ＤＳＰ）を行うことができる。受信信号のデジタル信号処理は速度及び柔軟性を与え、波長分散、偏波モード分散などのような光伝送経路に関連する非線形性の補正を含む種々の機能を実行するのに使用され得る。
【０００７】
位相変調システムのためのコヒーレント検出法は、絶対位相検出を利用することがある。絶対位相検出は、推定位相に基づいて受信データストリームにおける各ビットの値に関する決定、例えばソフト決定を行うことを伴う。残念ながら、Ｍ２−ＱＡＭ、例えばＱＰＳＫ及び１６ＱＡＭについては、信号コンステレーションは角度π／２の位相回転の下で不変である。推定位相を判定するのに使用されるキャリア位相推定器は、角度θと角度θ＋π／２の間を区別できない。結果として、推定されたキャリア位相は、現在の安定動作点から、隣接安定動作点の誘引力の領域に押し込まれることになり、これにより信号コンステレーションがπ／２だけ事実上回転する。この現象をサイクルスリップという。サイクルスリップによって、サイクルスリップのイベント後に大量の決定エラーが発生し得る。サイクルスリップの影響は、情報シンボルの差動復号化を用いることによって実際のスリップの継続時間に制限され得る。しかし、差動復号化は、絶対位相検出の約２倍のビットエラーレートを有することになる。
【０００８】
サイクルスリップを訂正するための１つの手法は、既知の情報シンボルにパイロットシンボルを導入することである。受信パイロットシンボルと既知の情報シンボルとの間の位相差を計算することによってパイロットシンボルのキャリア位相は明確に推定され得るので、パイロットシンボルは位相の曖昧さをなくす。しかし、パイロットシンボルのオーバーヘッドによってシンボルレートが大きくなり、感度ペナルティがもたらされる。これに対処するために、パイロットシンボルが長い期間で挿入され得る。一般に、サイクルスリップが検出されてキャリア位相の基準が訂正されるまでに、パイロットシンボル間のシンボル数の約半分の期間を要することになる。訂正を行うのに要する時間は、検出データにおけるバーストエラーをもたらすことになる。
【０００９】
添付図面を参照して、発明が例示として記載される。
【図面の簡単な説明】
【００１０】
【図１】図１は、本開示によるシステムの一例の実施形態のブロック図である。
【図２】図２は、本開示による送信機の一例の実施形態のブロック図である。
【図３】図３は、本開示による受信機の一例の実施形態のブロック図である。
【図４】図４は、本開示による一例の送信機の変調出力を模式的に示す。
【図５】図５は、本開示によるシステムにおける９０度位相エラーに関連するシンボルを示すように相互接続されたコンステレーション点を有する１６−ＱＡＭ信号の一例の実施形態のコンステレーション図である。
【図６Ａ】図６Ａは、ＱＰＳＫを用いる本開示による集合分割を模式的に示すコンステレーション図を含む。
【図６Ｂ】図６Ｂは、１６−ＱＡＭを用いる本開示による集合分割を模式的に示すコンステレーション図を含む。
【図６Ｃ】図６Ｃは、６４−ＱＡＭを用いる本開示による集合分割を模式的に示すコンステレーション図を含む。
【図７】図７は、本開示による送信機の他の例の実施形態のブロック図である。
【図８】図８は、本開示による受信機の他の例の実施形態のブロック図である。
【図９】図９は、ビットインタリーブ符号変調及び反復復号化を用いる本開示による一例の手法の性能を示すＢＥＲ対ＳＮＲのプロットである。
【図１０】図１０は、本開示による方法の一例を示すフローチャートである。
【発明を実施するための形態】
【００１１】
一般に、本開示によるシステムは、送信されるデータのｎ−３ビットの各ブロックに対して３個のパリティビットを挿入してｎビットのブロック（ｎ−３個の情報ビットに、ｎ個の情報ビットに対応付けられた３個のパリティビットを加えたもの）を生成する。ｎビットのブロックの各々は、光キャリア上に変調されて受信機に伝送される３個の対応のＱＡＭシンボルにグレイマッピングされる。受信機は、最大事後確率（ＭＡＰ）を用いてシンボルを検出し、サイクルスリップを訂正する。本開示によるシステムを、カスケード単一パリティチェック（ＣＳＰＣ）符号化を用いて説明することができる。
【００１２】
一実施形態では、例えばｎビットの各ブロックに対応付けられたシンボルは、正しいパリティ及び検出シンボルに対する最小ユークリッド距離を有するシンボルを選択することによって復号され得る。シンボル決定を行う際に検出器は正しいパリティを用いるので、検出器におけるエラー信号はサイクルスリップとともに増加して、サイクルスリップを自動的に検出及び訂正し、検出器が９０度（π／２）位相エラーに収束することを防止する。
【００１３】
ここで使用する「グレイマップ」又は「グレイマッピング」とは、隣接コードワードが１ビットだけ異なるようにビットの連続集合の各々に符号が割り当てられ、追加のビットをデータストリームに付加することを伴わない公知のグレイマッピング法のことをいう（すなわち、グレイマッピングはオーバーヘッドを有しない）。ここで使用される用語「結合される」とは、それによって１つのシステム要素により搬送される信号が、「結合される」要素に与えられるような任意の接続、結合、リンクなどをいう。そのような「結合される」装置、又は信号及び装置は、必ずしも相互に直接接続されず、そのような信号を操作又は修正する中間構成要素又は装置によって分離されていてもよい。
【００１４】
図１は、本開示によるＷＤＭ伝送システム１００の一例の実施形態の簡略化されたブロック図である。伝送システムは、複数の光チャネルを光情報経路１０２上で送信端末１０４から、遠隔位置にある１以上の受信端末１０６に伝送するように作用する。例示のシステム１００は、チャネルを送信機から５０００ｋｍ以上離れた受信機に伝送するために構成された長距離海底システムであってもよい。例示の実施形態は光システムの関連で記載され、長距離ＷＤＭ光学システムとの関係において有用であるが、ここに記載される広範な概念は、他のタイプの信号を送受信する他の通信システムにおいても実施され得る。
【００１５】
当業者であれば、システム１００は説明の容易化のために大幅に簡略化されたポイントトゥポイントシステムとして図示されていることを認識するはずである。例えば、送信端末１０４及び受信端末１０６は、もちろん、双方とも送受信機として構成されてもよく、それによって各々が送信機能及び受信機能の双方を実行するように構成され得る。ただし、説明の容易化のために、端末は、ここでは送信機能又は受信機能に関してのみ図示及び記載される。開示によるシステム及び方法が多種多様なネットワーク要素及び構成に取り入れられ得ることが理解されるはずである。ここで説明する例示の実施形態は、限定ではなく、例示のみとして提供される。
【００１６】
説明する例示の実施形態では、複数の送信機ＴＸ１、ＴＸ２・・・ＴＸＮの各々は対応の入力ポート１０８−１、１０８−２・・・１０８−Ｎでデータ信号を受信し、データ信号を対応の波長λ１、λ２・・・λＮで送信する。送信機ＴＸ１、ＴＸ２・・・ＴＸＮの１以上は、本開示による変形ＢＩＣＭ−ＩＤ法を用いて、データを対応の波長で変調するように構成され得る。送信機は、もちろん、説明の容易化のために大幅に簡略化された形式で示されている。当業者であれば、データ信号を所望の振幅及び変調でその対応の波長で送信するために構成された電気的及び光学的部品を各送信機が含み得ることを認識するはずである。
【００１７】
伝送波長又はチャネルは、複数の経路１１０−１、１１０−２・・・１１０−Ｎでそれぞれ搬送される。データチャネルは、マルチプレクサ又はコンバイナ１１２によって光経路１０２上で統合信号に合波される。光情報経路１０２は、光ファイバ導波路、光増幅器、光学フィルタ、分散補償モジュール並びに他の能動及び受動部品を含み得る。
【００１８】
統合信号は、１以上の遠隔受信端末１０６で受信される。デマルチプレクサ１１４は、波長λ１、λ２・・・λＮの伝送チャネルを、対応の受信機ＲＸ１、ＲＸ２・・・ＲＸＮに結合された対応の経路１１６−１、１１６−２・・・１１６−Ｎ上に分離する。受信機ＲＸ１、ＲＸ２・・・ＲＸＮの１以上は、本開示による変形ＢＩＣＭ−ＩＤ法に関連する反復復号化を用いて伝送信号を復調するように構成され、対応の出力データ信号を対応の出力経路１１８−１、１１８−２、１１８−３、１１８−Ｎ上に供給する。
【００１９】
図２は、本開示による一例の送信機２００の簡略化されたブロック図である。説明する例示の実施形態２００は、パリティビット符号器２０２、グレイマッパ２０４、並びに符号化及び変調された出力をキャリア波長λＮ上に供給するための連続波レーザ２０８の出力を変調する変調器２０６を含む。
【００２０】
パリティビット符号器２０２は、３個の対応のパリティビットで入力経路１１８−Ｎ上に供給されたデータストリームのｎ−３個の情報ビットの各ブロックを符号化するように構成され得る。パリティビット符号器２０２の出力は、ｎビットの連続するブロック、すなわち、ｎ−３個の情報ビットに３個のパリティビットを加えたものを含む。より詳細を以下に示すような態様で、パリティビット符号器２０２によって付加されたパリティビットの各々は、ｎビットの異なる対応の部分のパリティを特定する。
【００２１】
周知のように、パリティビットは偶数パリティビット又は奇数パリティビットとなり得る。偶数パリティを用いる場合、パリティビットは、パリティビットに対応付けられた情報ビット（パリティビットを含まず）における１の数が奇数の場合には１の値にセットされる。パリティビットに対応付けられた情報ビットにおける１の数が既に偶数の場合には、偶数パリティビットは０にセットされる。奇数パリティを用いる場合、パリティビットは、パリティビットに対応付けられた情報ビット（パリティビットを含まず）における１の数が偶数の場合に１にセットされる。パリティビットに対応付けられた情報ビットにおける１の数が既に奇数の場合には、奇数パリティビットは０にセットされる。パリティビット符号器２０２において入力データのビットを偶数又は奇数パリティビットで符号化するためのハードウェア及びソフトウェア構成は、当業者には周知である。
【００２２】
パリティビット符号器の符号化出力は、グレイマッパ２０４に結合される。グレイマッパ２０４は、ｎビットの各ブロックを３個の対応のＱＡＭシンボルにマッピングするように構成される。ｎビットの各ブロックに対応付けられた３個のＱＡＭシンボルは、変調器２０６を用いて連続波レーザ２０８の光キャリア波長λＮ上に変調される。変調器２０６は、任意の周知の変調方法を用いて複数のＱＡＭシンボルをキャリア波長λＮ上に変調することができる。変調器２０６の符号化、マッピング及び変調された出力は、ＷＤＭシステムにおけるマルチプレクサ１１２（図１）に結合され得る。
【００２３】
図３は、本開示による一例の受信機３００の簡略化されたブロック図である。説明する例示の実施形態３００は、光信号検出器３０２及びデマッパ３０４を含む。検出器３０２は周知のコヒーレント受信機、例えば、偏波ダイバシティコヒーレント受信機を含むことができ、その周知のコヒーレント受信機は、光キャリア波長λＮ上で信号を受信し、変調器２０６（図２）によって光キャリア波長λＮ上に変調されたＱＡＭシンボルを表す１以上の対応の電気的出力（例えば、偏波多重変調フォーマットにおいて各偏波に対応付けられた出力）に光信号を変換するように構成されている。
【００２４】
デマッパ３０４は、デジタル信号処理（ＤＳＰ）回路３０８の一部として構成され得る。一般に、ＤＳＰは、例えば、直接に及び／又はソフトウェアインストラクションの制御下で、特定のインストラクションシーケンスを実行するために構成された１以上の特定用途向け集積回路（ＡＳＩＣ）及び／又は特殊用途プロセッサを用いた信号処理を伴う。検出器すなわちコヒーレント受信機を内蔵する受信機、及びコヒーレント受信機のデジタル出力を処理するためのキャリア位相推定を用いるＤＳＰ回路の一例は、米国特許第８２９５７１３号に記載され、その教示が参照としてここに取り込まれる。
【００２５】
図２及び３の双方を参照すると、ＤＳＰ回路３０８は、検出器３０２の出力を処理し、入力１０８−Ｎで送信機２００に供給されたデータを再生する出力を供給することができる。デマッパ３０４は、検出器３０２の電気的出力を受信し、キャリア位相推定機能を用いて、グレイマッパ２０４によって適用されたマッピングを逆処理し、パリティビット符号器２０２によって適用されたパリティビットを除去する。デマッパの出力は、入力１０８−Ｎで送信機２００に供給されたデータのｎ−３ビットの連続するブロックを表すデマッピング出力である。
【００２６】
デマッピングは、例えば、最大事後確率（ＭＡＰ）検出器を用いて実行され、受信機の出力からの先験的尤度比（ＬＬＲ）帰還に応じて反復的に実行され得る。本開示によるシステムにおけるデマッパ３０４は、パリティビット符号器２０２によって適用されるパリティビットで示されるパリティを用いてサイクルスリップの訂正をもたらすことができる。ある実施形態では、例えば、デマッパ３０４は、パリティビットを用いて、サイクルスリップの訂正を自動的にもたらすようにデマッピングを実行することができる。
【００２７】
Ｍを正の整数として、所与のＭ２−ＱＡＭについて、各シンボルはｌｏｇ２（Ｍ２）個のビットによって表される。本開示によるＣＳＰＣ法は、一度に３個のＱＡＭシンボルを用いる。例えば、１６−ＱＡＭについて、３シンボルＣＳＰＣブロック毎の合計ビット数は、３×ｌｏｇ２（１６）＝１２となる。ｎ＝３×ｌｏｇ２（Ｍ２）を符号ビット数として、３シンボルＣＳＰＣブロック符号は、（ｎ−３）／ｎの符号レートを有する。
【００２８】
例えば、図４は、本開示による１２ビットＣＳＰＣ１６−ＱＡＭ（ここではＣＳＰＣ（９，１２）ともいう）の出力４００を模式的に示すものであり、パリティビット符号器２０２は、３個のパリティビットを各連続する９個の情報ビットに加えて１２ビットのブロックをグレイマッパ２０４に供給する。グレイマッパ２０４は、１２ビットの連続ブロックを３個の対応のシンボルにマッピングする。図５は、グレイマッパ２０４による１６−ＱＡＭ信号のグレイマッピングを示すコンステレーション図５０２である。周知の態様において、図５のコンステレーション図は、シンボルに対応付けられたビット（コードワード）とともに、各ＱＡＭシンボルの振幅及び位相を示す複数のコンステレーション点を示す。説明の簡略化及び容易化のために、ここに記載される例示の実施形態はＣＳＰＣ（９，１２）を用いるものとする。ただし、当業者であれば、本開示によるシステムは、任意のＭ２−ＱＡＭコンステレーションに対して一般化され得ることを認識するはずである。
【００２９】
パリティビット符号器２０２によるデータストリームに挿入されたパリティビットＰ１、Ｐ２及びＰ３の各々は、グレイマッパ２０４に供給されるＣＳＰＣブロックにおける対応のｎビットの集合のパリティ（情報ビットｘ１、ｘ２・・・ｘｎ−３及びパリティビットＰ１、Ｐ２及びＰ３）を示す。特に、
１．パリティビットＰ１は、
【数１】
2016533127.tif 000003
（式１）
として計算されるので、２ｎ／３ビットの合計のパリティは奇数となる。
２．パリティビットＰ２は、
【数２】
2016533127.tif 000004
及び
【数３】
2016533127.tif 000005
（式２）
として計算されるので、ｎ／２ビットの合計のパリティは偶数となる。
３．パリティビットＰ３は、
【数４】
2016533127.tif 000006
（式３）
として計算されるので、全ｎビットのパリティは奇数となる。パリティビットＰ１、Ｐ２及びＰ３はそれらの対応のビットの集合についての特定の（奇数又は偶数の）パリティを示すものとして記載されるが、当業者であれば、パリティビット符号器がパリティビットＰ１、Ｐ２及びＰ３について逆のパリティ（奇数又は偶数）も使用できることを認識するはずである。
【００３０】
上式１−３に記載されるように計算された３個のパリティビットＰ１、Ｐ２及びＰ３は、以下のように、ｎビットのＣＳＰＣブロック内での既知の位置にパリティビット符号器２０２によって挿入される。
１．Ｐ１は、最初の２ｎ／３ビットのいずれかに位置決めされる。
２．Ｐ２は、２ｎ／３＋１番目のビットと５ｎ／６番目のビットとの間のいずれかに位置決めされる。
３．Ｐ３は、２ｎ／３＋１番目のビットと、Ｐ２の位置とは異なるｎ番目のビットとの間のいずれかに位置決めされる。
パリティビットをｎビットのＣＳＰＣブロック内でこれらの位置に配置すると、ＣＳＰＣブロックに対応付けられたＱＡＭシンボルにおける集合分割が確立する。本開示によるＣＳＰＣを用いる集合分割によって、コンステレーション点間の増加する最小ユークリッド距離の集合から受信器において変調シンボルが選択可能となる。これによって、低い信号対雑音比（ＳＮＲ）又は同等に低いＱファクタでの検出が可能となる。
【００３１】
図６Ａ、６Ｂ及び６Ｃは、例えば、それぞれＱＰＳＫ、１６−ＱＡＭ及び６４−ＱＡＭで使用されるＣＳＰＣについての本開示による集合分割を示す。集合分割処理は、図６Ａ〜６Ｃに示す実施形態に共通し、本開示による任意のＣＳＰＣ符号化Ｍ２−ＱＡＭに適用可能である。図示するように、本開示によるＣＳＰＣ符号化Ｍ２−ＱＡＭにおける集合分割は、ｄを第１のシンボルにおけるコンステレーション点間の最小ユークリッド距離として、第２のシンボルにおけるコンステレーション点間のユークリッド距離を
【数５】
2016533127.tif 000007
に増加させる。第３のシンボルにおける点間のユークリッド距離は２ｄに増加する。
【００３２】
図６Ａ、６Ｂ及び６Ｃを参照すると、例えば、シーケンスにおいて第１のシンボル（Ｓ１）があるとすると、本開示によるＣＳＰＣの結果として、第２のシンボル（Ｓ２）は、Ｓ１のハミング重みによって決定される部分集合から選択され得る。当業者であれば、ビットのストリングのハミング重みはストリングにおける「１」の数であることを認識する。例えば、Ｓ１のハミング重みが偶数である場合には、Ｓ２は図６Ａ〜６Ｃにおける左側に示す部分集合から選択され、Ｓ１のハミング重みが奇数である場合には、Ｓ２は図６Ａ〜６Ｃにおける右側に示す部分集合から選択される。
【００３３】
第３のシンボル（Ｓ３）がＳ１及びＳ２によって決定される。Ｓ３は、Ｓ１及びＳ２の｛ｌｏｇ２（Ｍ２）｝／２個の最大有効ビット（ＭＳＢ）（ＩＳ１＾ＩＳ２という）のハミング重みによって決定される部分集合から選択されることになる。例えば、ＩＳ１＾ＩＳ２のハミング重みが偶数である場合には、Ｓ３は図６Ａ〜６Ｃにおける右側に示す部分集合から選択され、ＩＳ１＾ＩＳ２のハミング重みが奇数である場合には、Ｓ３は図６Ａ〜６Ｃで左側に示す部分集合から選択される。
【００３４】
受信機３００において、デマッパ３０４は、グレイマッピングされたＱＡＭではマッピングされたビットに対応付けられたパリティが各９０度位相回転後に変化することを利用して、サイクルスリップを検出及び訂正することができる。この特性は、コンステレーション点を相互接続する四角形５０４、５０６、５０８及び５１０によって図５に示される。各四角形の角は、隣接角部のコンステレーション点から９０度となるコンステレーション点に位置決めされる。任意のコンステレーション点について、そこから±９０度となるコンステレーション点は異なるパリティを有する。例えば、ビット（１１１０）に対応付けられた１６−ＱＡＭシンボルは、９０度位相回転後にビット（００１１）又はビット（１００１）に対応付けられた１６−ＱＡＭシンボルに変化することになる。
【００３５】
この規則は、Ｍをシンボルによって符号化されたビットの数として、グレイマッピングでの任意のＭ２−ＱＡＭについて当てはまる。特に、Ｍ２−ＱＡＭについて、シンボルｘ＋ｙｉが、ｌｏｇ２（Ｍ）個のバイナリビットを実部及び虚部の双方において独立してＭパルス振幅変調（ＰＡＭ）にマッピングすることによって生成される。ｘ及びｙを実数としてグレイマッピングすると、ＰＡＭ信号ｙは、
Ｐａｒｉｔｙ（ｙ）≠Ｐａｒｉｔｙ（−ｙ）
のプロパティを有する。これは、グレイマッピングを行えば、反対の符号を有する２つの最も近いシンボルは１ビット差のみを有するのでこれらの２つのシンボルのパリティが異なるためである。ＱＡＭシンボルのパリティは、
Ｐａｒｉｔｙ（ｘ＋ｙｉ）＝Ｐａｒｉｔｙ（ｘ）  ＸＯＲ  Ｐａｒｉｔｙ（ｙ）
Ｐａｒｉｔｙ（ｘ−ｙｉ）＝Ｐａｒｉｔｙ（ｘ）  ＸＯＲ  Ｐａｒｉｔｙ（−ｙ）
となる。
したがって、
Ｐａｒｉｔｙ（ｘ＋ｙｉ）≠Ｐａｒｉｔｙ（ｘ−ｙｉ）
となる。
【００３６】
一方、Ｍ２−ＱＡＭシンボルの実部及び虚部の入替えによってはパリティは変化しない。
Ｐａｒｉｔｙ（ｘ−ｙｉ）＝Ｐａｒｉｔｙ（−ｙ＋ｘｉ）
したがって、シンボル（ｘ＋ｙｉ）ｅｊπ／２＝−ｙ＋ｘｉのπ／２位相回転によってパリティが変化する。
Ｐａｒｉｔｙ（ｘ＋ｙｉ）≠Ｐａｒｉｔｙ（−ｙ＋ｘｉ）
グレイマッピングされるＱＡＭシンボルに対応付けられたビットのパリティは９０度位相回転毎に変化する。そのため、シンボルシーケンスエラーの確率が９０度位相回転によって劇的に増加することになるので、任意の９０度位相回転、すなわち、サイクルスリップは、受信機においてＭＡＰ復号器によって検出及び訂正可能となる。
【００３７】
例えば、本開示によるＣＳＰＣ（９，１２）の実施形態において、デマッパ３０４は、μを適応フィルタのステップサイズとして、
ｅ（ｍ）＝ＭＡＰ（ｕ（ｍ−１）、ｕ（ｍ）、ｕ（ｍ＋１））−ｕ（ｍ）
ｈｘ（ｋ）＝ｈｘ（ｋ）＋μｅ（ｍ）ｘ＊（ｍ−ｋ）
ｈｙ（ｋ）＝ｈｙ（ｋ）＋μｅ（ｍ）ｙ＊（ｍ−ｋ）
としてＭＡＰ復号器からエラー信号ｅ（ｍ）によって適応的に更新される線形歪及びキャリア位相オフセットを等化するのに、マルチタップリニアフィルタｈｘ（ｋ）、ｈｙ（ｋ）、ｋ＝−Ｎ、−Ｎ＋１・・・Ｎを含み得る。ＭＡＰ復号器は、フィルタリングされた信号ｕ（ｎ）からの最小ユークリッド距離を有するパリティビットＰ１、Ｐ２及びＰ３（９個の情報ビットで合計２ｎ−３＝５１２個の可能な対がある）によって示される正しいパリティを有する３個の１６ＱＡＭシンボルの集合を決定することができる。一実施形態では、例えば、ＭＡＰ復号器は、ｎ＝１０のユークリッド距離計算しか必要としないＷａｎｇｅｒのアルゴリズムによって効率的に実施可能である。
【００３８】
エラー信号ｅ（ｍ）は正しいパリティでＭＡＰ復号器によって駆動されるので、エラー信号は９０度（π／２）位相エラーで大きくなる。言い換えると、デマッパ３０４における位相追跡のための適応等化器の解決手段は、決して９０度位相エラーには収束し得ず、それにより隣接象限へのサイクルスリップを防止することができる。さらに、本開示によるＣＳＰＣから得られるシンボルシーケンス間での増加したユークリッド距離に起因して、ＭＡＰ復号器の出力は、低いＱファクタでさえも向上した信頼性及び高い訂正能力を有する。これによって、シンボルエラーの確率が２５％よりも大きくなる３．５ｄＢの低いＱファクタ値においてさえも、９０度位相の不確定さを回避しつつ等化が可能となる。したがって、さらに、エラーのない伝送が低いＳＮＲで達成され、より高いネット符号化利得（ＮＣＧ）、例えば１２．８ｄＢが、所与の同じＦＥＣ符号で達成され、変動レートのトランスポンダ技術を適応させる適度なスペクトル効率が達成可能となる。
【００３９】
本開示によるＣＳＰＣ符号化システムは様々な構成で実施され得る。図７は、ビットインタリーブ符号化変調（ＢＩＣＭ）構成でＣＳＰＣ符号化を実施するように構成された本開示による一例の送信機７００の簡略化されたブロック図である。説明する例示の実施形態７００は、デマルチプレクサ７０１、複数のＦＥＣ符号器７０２−（ｎ−３）、７０２−（ｎ−２）・・・７０２−１、インタリーバ７０４、図２との関係で説明したようなパリティビット符号器２０２、及び変調器２０６を駆動するためのグレイマッパ２０４を含む。
【００４０】
デマルチプレクサ７０１は、シリアル入力データストリームを経路１０８−Ｎ上で受信し、入力データストリームをｎ−３個の個別の並列多重分離データストリームに多重分離するための周知の構成を採ることができる。ｎ−３個のデータストリームの各々は、ＦＥＣ符号器７０２−（ｎ−３）、７０２−（ｎ−２）・・・７０２−１の１つに対応付けられたものに結合される。ＦＥＣ符号器７０２−（ｎ−３）、７０２−（ｎ−２）・・・７０２−１の各々は、それにより対応のＦＥＣ符号で受信されたデータストリームを符号化して対応のＦＥＣ符号化出力を供給するように構成され得る。
【００４１】
多数のＦＥＣ符号が知られており、各々は、どのように符号が生成されてその結果それらがどのように機能するかに関して異なるプロパティを有する。周知の誤り訂正符号の例は、線形及び巡回ハミング符号、巡回ＢＣＨ（Ｂｏｓｅ−Ｃｈａｕｄｈｕｒｉ−Ｈｏｃｑｕｅｎｇｈｅｍ）符号、畳込み（ビタビ）符号、巡回Ｇｏｌａｙ及びＦｉｒｅ符号、ターボ畳込み及び積符号（ＴＣＣ、ＴＰＣ）並びに低密度パリティチェック符号（ＬＤＰＣ）を含む。ＬＤＰＣ符号は、その低い対応のオーバーヘッドのために、ＦＥＣ符号器７０２−（ｎ−３）、７０２−（ｎ−２）・・・７０２−１について特に有用な符号となる。ＦＥＣ符号器７０２−（ｎ−３）、７０２−（ｎ−２）・・・７０２−１及び対応の復号器８０２−（ｎ−３）、８０２−（ｎ−２）・・・８０２−１（図８）において種々のエラー訂正符号を実施するためのハードウェア及びソフトウェア構成は当業者には周知である。
【００４２】
ＦＥＣ符号器７０２−（ｎ−３）、７０２−（ｎ−２）・・・７０２−１の各々の符号化出力は、インタリーバ７０４に結合される。様々なインタリーバ構成が知られている。説明する実施形態では、インタリーバ７０４は、ＦＥＣ符号器７０２−（ｎ−３）、７０２−（ｎ−２）・・・７０２−１からｎ−３個の出力コードワードを受信し、ｎ−３ビットの並列インタリーブ出力をパリティビット符号器２０２に供給する。
【００４３】
上述のように、パリティビット符号器２０２は、３個の対応のパリティビットでｎ−３ビットの各ブロックを符号化するように構成され得る。パリティビット符号器２０２の出力は、ｎビット、すなわち、ｎ−３個の情報ビットに３個のパリティビットを加えた連続するブロックを含む。パリティビット符号器２０２の符号化出力は、グレイマッパ２０４に結合される。グレイマッパ２０４は、ｎビットの各ブロックを３個の対応のＱＡＭシンボルにマッピングするように構成される。ｎビットの各ブロックに対応付けられたＱＡＭシンボルは、シンボルを連続波レーザ２０８の光キャリア波長λＮ上に変調するために変調器２０６に供給される。
【００４４】
図８は、図７に示す送信機７００を用いて変調された信号を受信するための、本開示による一例の受信機８００の簡略化されたブロック図である。説明する例示の実施形態８００は、検出器３０２、デマッパ３０４ａ、デインタリーバ８０２、複数のＦＥＣ復号器８０４−（ｎ−３）、８０４−（ｎ−２）・・・８０４−１、マルチプレクサ８０６及びインタリーバ８０８を含む。
【００４５】
上述したように、検出器３０２は、キャリア波長λＮ上に変調された光信号を受信し、光信号を電気信号に変換するように構成される。デマッパ３０４ａは、検出器３０２の出力を受信し、データのマッピングを送信機においてグレイマッパ２０４によって与えられた変調形式に戻し、例えば上述したようにＭＡＰ復号器を用いてパリティビット符号器によって適用されるパリティビットを除去する。図示するように、デマッピングは、受信機の出力、すなわち、説明する実施形態におけるマルチプレクサ８０６の出力からの先験的ＬＬＲ帰還に応じて反復的に実行され得る。デマッパ３０４ａは、ｎ−３ビットのブロックを含むデマッパ出力を供給し、送信機７００におけるインタリーバ７０４の出力を再生する。
【００４６】
デマッパ３０４ａのデマッパ出力は、送信機７００において対応のインタリーバ７０４によって実行されるビットインタリーブを戻してｎ−３個の対応のデインタリーブ出力をＦＥＣ復号器８０４−（ｎ−３）、８０４−（ｎ−２）・・・８０４−１に供給するデインタリーバ８０２に結合される。ＦＥＣ復号器８０４−（ｎ−３）、８０４−（ｎ−２）・・・８０４−１の各々は、送信機７００において適用されるＦＥＣ符号を用いてそれにより受信されたデータストリームを復号し、対応のＦＥＣ復号出力を供給するように構成され得る。ＦＥＣ復号器８０４−（ｎ−３）、８０４−（ｎ−２）・・・８０４−１の各々のｎ−３個のＦＥＣ復号出力は、既知のマルチプレクサ８０６に結合される。マルチプレクサ８０６は、ＦＥＣ符号の復号器８０４−（ｎ−３）、８０４−（ｎ−２）・・・８０４−１の出力を多重化して、送信機７００で経路１１８−Ｎ上にデータを再生するシリアルデジタルデータ出力を生成する。
【００４７】
マルチプレクサ８０６のデータ出力は、インタリーバ８０８を介してデマッパ３０４ａにフィードバックされて、そこへの入力を復号する際にデマッパ３０４ａによって使用される先験的ＬＬＲ情報を供給する。インタリーバ８０８は、基本的に、デインタリーバ８０２によって実行されるデインタリーブを戻す。
【００４８】
４．８ｂ／ｓＨｚのＳＥとなる第１のＦＥＣ符号として２０％ＬＤＰＣ符号を用いる反復復号（ＩＤ）の実施形態での上述のＣＳＰＣ（９，１２）−ＢＩＣＭの性能を図９との関係で説明する。図９は、ラインビットエラーレート（ＢＥＲ）の測定プロット７０２及びＭＡＰ検出器出力におけるＢＥＲのプロット７０４、並びに受信機の出力におけるビットエラーレート（ＢＥＲ）対受信機の入力におけるＳＮＲ／ビット（ｄＢ）のプロット７０６、７０８、７１０、７１２及び７１４を含む。各プロット７０６、７０８、７１０、７１２及び７１４は、受信機における復号の異なる反復に対応付けられる。図示するプロットは、本開示による変形ＣＳＰＣ−ＢＩＣＭ−ＩＤ法によって累進的に低いＳＮＲで反復した等化が可能となることを示す。
【００４９】
図１０は、本開示による方法１０００を示すフローチャートである。処理１００２は、３個のパリティビットで信号のｎ−３ビットの連続するブロックを符号化してｎビットの連続するブロックを供給することを含む。ｎビットの連続するブロックは、光キャリア波長で変調されて変調光信号を供給する（１００６）３個の対応のＱＡＭシンボルにグレイマッピングされる（１００４）。変調光信号が検出されて電気信号を供給し（１００８）、３個の対応のＱＡＭシンボルが電気信号からデマッピングされてｎ−３ビットのブロックを表すデマッパ出力を供給する（１０１０）。
【００５０】
図１０は実施形態による種々の処理を示すが、図１０に図示する処理の全てが必ずしも他の実施形態に対して必要とはならないことが理解されるべきである。もちろん、本開示の他の実施形態において、図１０に図示する処理及び／又はここに記載される他の処理が、いずれの図面にも具体的示されていないがそれでも本開示に完全に従う態様で組み合わされてもよいことが、ここに充分に考慮される。したがって、ある図面に厳密に示されていない構成及び／又は処理に関する請求項は、本開示の範囲及び内容のうちにあるものとみなされる。
【００５１】
開示の一形態によると、ｎ−３ビットの連続するブロックを３個のパリティビットで符号化してｎビットの連続するブロックを供給するように構成されたパリティビット符号器、パリティビット符号器に結合され、ｎビットのブロックの各々を３個の対応の直交振幅変調（ＱＡＭ）シンボルにマッピングするように構成されたグレイマッパ、グレイマッパに結合され、グレイマッパの出力に応じて光信号を変調して、３個の対応のＱＡＭシンボルを含む変調光信号を供給するように構成された変調器、変調光信号を受信し、光信号を表す電気信号を供給するための検出器、及び電気信号に応じてｎ−３ビットのブロックを表すデマッパ出力を供給するように構成されたデマッパを含むシステムが提供される。
【００５２】
開示の他の形態によると、３個の対応の直交振幅変調（ＱＡＭ）信号にマッピングされるｎ−３個の情報ビット及び３個のパリティビットからなるｎビットのブロックを表す変調光信号を受信し、光信号を表す電気信号を供給するための検出器、及び電気信号に応じて、ｎ−３個の情報ビットを表すデマッパ出力を供給するように構成されたデマッパを含む光信号受信機が提供される。
【００５３】
開示の他の形態によると、ｎ−３ビットの連続するブロックを３個のパリティビットで符号化してｎビットの連続するブロックを供給すること、ｎビットのブロックの連続する３個の対応の直交振幅変調（ＱＡＭ）シンボルにグレイマッピングすること、３個の対応のＱＡＭシンボルを光キャリア波長で変調して変調光信号を供給すること、変調光信号を検出して電気信号を供給すること、及び電気信号から３個の対応のＱＡＭシンボルをデマッピングしてｎ−３ビットのブロックを表すデマッパ出力を供給することを含む方法が提供される。
【００５４】
ここに記載される方法の実施形態は、プロセッサ及び／又は他のプログラム可能な装置を用いて実施できる。この目的のため、ここに記載される方法は、１以上のプロセッサによって実行されたときに当該方法を実行するようにそこに記憶されたインストラクションを有する有形のコンピュータ可読記憶媒体で実施できる。したがって、例えば、送信機及び／又は受信機は、ここに記載される処理を実行する（例えば、ファームウェア又はソフトウェアでの）インストラクションを記憶する記憶媒体（不図示）を含み得る。記憶媒体は、任意のタイプの非一時的な有形媒体、例えば、フロッピーディスク、光学ディスク、コンパクトディスク読出し専用メモリ（ＣＤ−ＲＯＭ）、再書込み可能コンパクトディスク（ＣＤ−ＲＷ）及び電磁光学ディスクを含む任意のタイプのディスク、読出し専用メモリ（ＲＯＭ）、ダイナミック及びスタティックＲＡＭなどのランダムアクセスメモリ（ＲＡＭ）、消去可能プログラマブル読出し専用メモリ（ＥＰＲＯＭ）、電気的に消去可能なプログラマブル読出し専用メモリ（ＥＥＰＲＯＭ）、フラッシュメモリ、磁気若しくは光学カードなどの半導体装置、又は電子的インストラクションを記憶するのに適した任意のタイプの媒体を含み得る。
【００５５】
ここでのいずれのブロック図も、開示の原理を具現する例示的回路の概念的な図を表していることが当業者には分かるはずである。同様に、いずれのフローチャート、フロー図、状態遷移図、疑似符号などはコンピュータ可読媒体において実質的に表され、コンピュータ又はプロセッサによって、そのようなコンピュータ又はプロセッサが明示的に示されているか否かにかかわらず、そのように実行され得る種々の処理を表すことが分かるはずである。ソフトウェアモジュール、又は単にソフトウェアであることが示唆されるモジュールは、フローチャートの要素又は処理ステップ及び／若しくはテキストの記述の実行を示す他の要素の任意の組合せとしてここに示され得る。そのようなモジュールは、明示的に又は暗に示されるハードウェアによって実行され得る。
【００５６】
任意の機能的ブロックを含む図面に示す種々の要素の機能は、専用ハードウェア又は適切なソフトウェアとの関連でソフトウェアを実行することができるハードウェアの使用によって与えられる。プロセッサによって与えられる場合、機能は、単一の専用プロセッサによって、単一の共有プロセッサによって、又はそのいくつかは共有され得る複数の個別のプロセッサによって与えられる。さらに、用語「プロセッサ」又は「コントローラ」の明示的な使用は、ソフトウェアを実行することができるハードウェアを排他的に意味するものと解釈されるべきではなく、暗に、限定ではなく、デジタル信号プロセッサ（ＤＳＰ）ハードウェア、ネットワークプロセッサ、特定用途向け集積回路（ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、ソフトウェアを記憶するための読出し専用メモリ（ＲＯＭ）、ランダムアクセスメモリ（ＲＡＭ）及び不揮発性ストレージを含み得る。他の従来的な及び／又はカスタムのハードウェアも含まれ得る。
【００５７】
いずれかの実施形態において使用される「回路」は、例えば、単独で又は任意の組合せにおいて、ハードウェア回路、プログラム可能な回路、状態機械回路、及び／又はプログラム可能な回路によって実行されるインストラクションを記憶するファームウェアを備え得る。少なくとも１つの実施形態において、送信機及び受信機は１以上の集積回路を備え得る。「集積回路」は、例えば、ただし限定するわけではなく、半導体集積回路チップなどの、デジタル、アナログ若しくは混合信号用の半導体装置及び／又は小型電子装置であればよい。
【図１】
2016533127.tif 000012
【図２】
2016533127.tif 000013
【図３】
2016533127.tif 000014
【図４】
2016533127.tif 000015
【図５】
2016533127.tif 000016
【図６Ａ】
2016533127.tif 000017
【図６Ｂ】
2016533127.tif 000018
【図６Ｃ】
2016533127.tif 000019
【図７】
2016533127.tif 000020
【図８】
2016533127.tif 000021
【図１０】
2016533127.tif 000022
【図９】
2016533127.tif 000023
【国際調査報告】
2016533127.tif 000024
2016533127.tif 000025
2016533127.tif 000026