(57)【要約】【課題】本発明は、高調波歪み分離方法、非線形特性確定方法、装置及びシステムを提供する。【解決手段】かかる高調波歪み分離方法は、測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波と生成高調波との間の位相差を確定し；及び、前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離することを含む。【選択図】図３
2017139762.tif 000002
【特許請求の範囲】
【請求項１】
送信機又は受信機又は通信システムに配置される高調波歪み分離装置であって、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波と生成高調波との間の位相差を確定するための確定ユニット；及び
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離するための分離ユニットを含む、装置。
【請求項２】
請求項1に記載の装置であって、
前記確定ユニットは、
【数１】
2017139762.tif 000054
を用いて、第二次高調波に対応する固有高調波と生成高調波との間の位相差Δφ2を確定するための第一確定モジュールを含み、
R2は、第二次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L2は、第二次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c2は、第二次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される、装置。
【請求項３】
請求項2に記載の装置であって、
前記第一確定モジュールは、
0乃至2πの範囲から複数の角度Δφ2を選択するための第一選択モジュール；
各Δφ2を用いて、数１の式の中の
（外１）
2017139762.tif 000055
及び
（外２）
2017139762.tif 000056
に対して線形フィッティングを行うための第一線形フィッティングモジュール；及び
フィッティングの相関係数が最大であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択するための第二選択モジュールを含む、装置。
【請求項４】
請求項2に記載の装置であって、
前記第一確定モジュールは、
0乃至2πの範囲から複数の角度Δφ2を選択するための第三選択モジュール；
各Δφ2を用いて、数１の式の中の
（外３）
2017139762.tif 000057
及び
（外４）
2017139762.tif 000058
に対して除算を行うための第一除算モジュール；及び
各組の商の分散が最小であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択するための第四選択モジュールを含む、装置。
【請求項５】
請求項2に記載の装置であって、
前記分離ユニットは、
【数２】
2017139762.tif 000059
を用いて、前記測量待ちシステムにより生成された第二次高調波のパワーH2を分離するための第一分離モジュールを含む、装置。
【請求項６】
請求項1に記載の装置であって、
前記確定ユニットは、
【数３】
2017139762.tif 000060
を用いて、第三次高調波に対応する固有高調波と生成高調波との間の位相差Δφ3を確定するための第二確定モジュールを含み、
R3は、第三次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L3は、第三次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c3は、第三次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される、装置。
【請求項７】
請求項6に記載の装置であって、
前記第二確定モジュールは、
0乃至2πの範囲から複数の角度Δφ3を選択するための第五選択モジュール；
各Δφ3を用いて、数３の式の中の
（外５）
2017139762.tif 000061
及び
（外６）
2017139762.tif 000062
に対して線形フィッティングを行うための第二線形フィッティングモジュール；及び
フィッティングの相関係数が最大であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択するための第六選択モジュールを含む、装置。
【請求項８】
請求項6に記載の装置であって、
前記第二確定モジュールは、
0乃至2πの範囲から複数の角度Δφ3を選択するための第七選択モジュール；
各Δφ3を用いて、数３の式の中の
（外７）
2017139762.tif 000063
及び
（外８）
2017139762.tif 000064
に対して除算を行うための第二除算モジュール；及び
各組の商の分散が最小であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択するための第八選択モジュールを含む、装置。
【請求項９】
請求項6に記載の装置であって、
前記分離ユニットは、
【数４】
2017139762.tif 000065
を用いて、前記測量待ちシステムにより生成された第三次高調波のパワーH3を分離するための第二分離モジュールを含む、装置。
【請求項１０】
送信機又は受信機又は通信システムに配置される非線形特性確定装置であって、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波と生成高調波との間の位相差を確定するための第一確定ユニット；
前記位相差を用いて、測量待ちシステムにより生成された高調波のパワーを分離するための分離ユニット；及び
分離された、前記測量待ちシステムにより生成された高調波のパワーを用いて、前記測量待ちシステムの非線形特性を確定するための第二確定ユニットを含む、装置。
【請求項１１】
高調波歪み分離方法であって、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波と生成高調波との間の位相差を確定し；及び
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離することを含む、方法。
【請求項１２】
請求項11に記載の方法であって、
前記固有高調波と生成高調波との間の位相差の確定は、
【数５】
2017139762.tif 000066
を用いて、第二次高調波に対応する固有高調波と生成高調波との間の位相差Δφ2を確定することを含み、
R2は、第二次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L2は、第二次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c2は、第二次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される、方法。
【請求項１３】
請求項12に記載の方法であって、
前記第二次高調波に対応する固有高調波と生成高調波との間の位相差Δφ2の確定は、
0乃至2πの範囲から複数の角度Δφ2を選択し；
各Δφ2を用いて、数５の式の中の
（外９）
2017139762.tif 000067
及び
（外１０）
2017139762.tif 000068
に対して線形フィッティングを行い；及び
フィッティングの相関係数が最大であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択することを含む、方法。
【請求項１４】
請求項12に記載の方法であって、
前記第二次高調波に対応する固有高調波と生成高調波との間の位相差Δφ2の確定は、
0乃至2πの範囲から複数の角度Δφ2を選択し；
各Δφ2を用いて、数５の式の中の
（外１１）
2017139762.tif 000069
及び
（外１２）
2017139762.tif 000070
に対して除算を行い；及び
各組の商の分散が最小であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択することを含む、方法。
【請求項１５】
請求項12に記載の方法であって、
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離することは、
【数６】
2017139762.tif 000071
を用いて、前記測量待ちシステムにより生成された第二次高調波のパワーH2を分離することを含む、方法。
【請求項１６】
請求項11に記載の方法であって、
前記固有高調波と生成高調波との間の位相差の確定は、
【数７】
2017139762.tif 000072
を用いて、第三次高調波に対応する固有高調波と生成高調波との間の位相差Δφ3を確定すすることを含み、
R3は、第三次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L3は、第三次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c3は、第三次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される、方法。
【請求項１７】
請求項16に記載の方法であって、
前記第三次高調波に対応する固有高調波と生成高調波との間の位相差Δφ3の確定は、
0乃至2πの範囲から複数の角度Δφ3を選択し；
各Δφ3を用いて、数７の式の中の
（外１３）
2017139762.tif 000073
及び
（外１４）
2017139762.tif 000074
に対して線形フィッティングを行い；及び
フィッティングの相関係数が最大であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択することを含む、方法。
【請求項１８】
請求項16に記載の方法であって、
前記第三次高調波に対応する固有高調波と生成高調波との間の位相差Δφ3の確定は、
0乃至2πの範囲から複数の角度Δφ3を選択し；
各Δφ3を用いて、数７の式の中の
（外１５）
2017139762.tif 000075
及び
（外１６）
2017139762.tif 000076
に対して除算を行い；及び
各組の商の分散が最小であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択することを含む、方法。
【請求項１９】
請求項16に記載の方法であって、
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離することは、
【数８】
2017139762.tif 000077
を用いて、前記測量待ちシステムにより生成された第三次高調波のパワーH3を分離することを含む、方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、通信技術分野に関し、特に、高調波歪み（harmonic  distortion）分離方法、非線形特性確定方法、装置及びシステムに関する。
【背景技術】
【０００２】
通信システムの伝送レートが絶えず増加するにつれて、システムにおける非線形特性は、パフォーマンスを制限する主な要因の１つになるため、重要視されている。システムの非線形特性に対処する各種の方法では、１つの必要な前提ステップは、非線形特性を測量及び推定することである。
【０００３】
今のところ、よく用いられる方法では、先ず、測量待ちシステム又は測量待ち素子（部品）の非線形高調波を測量し、その後、高調波のパワーを直接利用して非線形特性（素子のTHD（Total  harmonics  distortion、全高調波歪み）指標）を表し、或いは、非線形をモデル化し、高調波のパワーに基づいてモデルの係数を計算し、そして、モデルを用いて非線形特性を確定する。その利点は、測量方法が簡単であり、適用性が広いことにある。しかし、これらの方法には、１つの問題、即ち、高調波を測量する時に、測量待ちシステム又は測量待ち素子に入力されるテスト信号自身が比較的大きい高調波成分を有すれば（これは、測量時に不可避である）、測得された高調波のパワーも比較的大きい誤差を有するという問題が存在する。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
発明者は、本発明の実現過程において、次のようなこと、即ち、従来の方法では、固有高調波（テスト信号自身のもの）及び生成高調波（測量待ちシステム又は測量待ち素子により生成されたもの）の位相関係を知らないので、このような誤差を無視し、後続操作の許容度を上げる方法を採用することでこのような測量誤差による影響に対処するだけであるということを発見した。その欠点は、明らかであるように、後続操作の効率及びパフォーマンスを犠牲にすることにある。
【０００５】
上述の問題を解決するために、本発明の実施例は、高調波歪み分離方法、非線形特性確定方法、装置及びシステムを提供する。
【課題を解決するための手段】
【０００６】
本実施例の第一側面によれば、高調波歪み分離方法が提供され、そのうち、該方法は、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数（memoryless  nonlinear  transfer  function）の複数組の入力パワー、出力パワー及び基本波振幅（fundamental  amplitude）を用いて、固有高調波及び生成高調波の位相差を確定し；及び
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離することを含む。
【０００７】
本実施例の第二側面によれば、非線形特性確定方法が提供され、そのうち、該方法は、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定し；
前記位相差を用いて、測量待ちシステムにより生成された高調波のパワーを分離し；及び
分離された、前記測量待ちシステムにより生成された高調波のパワーを用いて、前記測量待ちシステムの非線形特性を確定することを含む。
【０００８】
本実施例の第三側面によれば、送信機又は受信機又は通信システムに配置される高調波歪み分離装置が提供され、そのうち、該装置は、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定するための確定ユニット；及び
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離するための分離ユニットを含む。
【０００９】
本実施例の第四側面によれば、送信機又は受信機又は通信システムに配置される非線形特性確定装置が提供され、そのうち、該装置は、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定するための第一確定ユニット；
前記位相差を用いて、測量待ちシステムにより生成された高調波のパワーを分離するための分離ユニット；及び
分離された、前記測量待ちシステムにより生成された高調波のパワーを用いて、前記測量待ちシステムの非線形特性を確定するための第二確定ユニットを含む。
【００１０】
本実施例の第五側面によれば、送信機が提供され、前記送信機は、前述の第三側面又は第四側面に記載の装置を含む。
【００１１】
本実施例の第六側面によれば、受信機が提供され、前記受信機は、前述の第三側面又は第四側面に記載の装置を含む。
【００１２】
本実施例の第七側面によれば、通信システムが提供され、そのうち、前記通信システムは、送信機及び受信機を含み、且つ前記通信システムは、さらに、前述の第三側面又は第四側面に記載の装置を含む。
【００１３】
本発明の有益な効果は、本発明の実施例により、テスト信号自身の固有高調波を除去し、測量待ちシステム又は測量待ち素子により生成された高調波のパワーを分離することにより、非線形測量の正確度及びシステムのパフォーマンスを向上させることができるということにある。
【図面の簡単な説明】
【００１４】
【図１】本出願の原理を示す図である。
【図２】本出願の測量待ちシステムの非線形モデルを示す図である。
【図３】本実施例における高調波歪み分離方法のフローチャートである。
【図４】測量待ちシステムの線形フィルタレスポンスの測量を示す図である。
【図５】パフォーマンスの比較を示す図である。
【図６】本実施例における非線形特性確定方法のフローチャートである。
【図７】本実施例における高調波歪み分離装置の１つの実施方式を示す図である。
【図８】図7に示す装置における第一確定モジュールの１つの実施方式を示す図である。
【図９】図8に示す実施方式のハードウェア構成を示す図である。
【図１０】図7に示す装置における第一確定モジュールの他の実施方式を示す図である。
【図１１】図10に示す実施方式のハードウェア構成を示す図である。
【図１２】図7に示す装置における第一分離モジュールのハードウェア構成を示す図である。
【図１３】図7に示す装置における第二確定モジュールの１つの実施方式を示す図である。
【図１４】図7に示す装置における第二確定モジュールの他の実施方式を示す図である。
【図１５】本実施例における非線形特性確定装置の１つの実施方式を示す図である。
【図１６】本実施例における受信機の１つの実施方式を示す図である。
【図１７】本実施例における受信機の他の実施方式を示す図である。
【図１８】本実施例における送信機の１つの実施方式を示す図である。
【図１９】本実施例における送信機の他の実施方式を示す図である。
【図２０】本実施例における通信システムの１つの実施方式を示す図である。
【発明を実施するための形態】
【００１５】
本発明は、高調波歪み分離方法、非線形特性確定方法、装置及びシステムを提供し、設定されている非線形モデルに基づいて、モデルの係数が定数であるという仮定を用いて、先ず、固有高調波及び生成高調波の位相差を確定し、そして、この位相差を用いて、高調波の分離を行い、測量待ちシステム又は測量待ち素子により生成された高調波のパワーを計算する。これにより、非線形測量の正確度及びシステムのパフォーマンスを向上させることができる。ここで、測量待ちシステム又は測量待ち素子は、例えば、通信システムにおける増幅器、送信機における変調器などであり、便宜のため、測量待ちシステム及び測量待ち素子を統一的に測量待ちシステムと称する。
【００１６】
以下、図面を参照しながら本発明の各種の実施方式について説明する。なお、これらの実施方式は、例示に過ぎず、本発明を限定するものでない。
【００１７】
図1は、本実施例の原理を示す図である。図1に示すように、{T1、T2、T3、…}、{R1、R2、R3、…}、及び{H1、H2、H3、…}は、信号の基本波、第二次高調波、第三次高調波などのパワーを表す。例えば、{T1、T2、T3、…}は、テスト信号（入力信号）の基本波のパワー、第二次高調波のパワー、第三次高調波のパワーなどを表し、{R1、R2、R3、…}は、出力信号の基本波のパワー、第二次高調波のパワー、第三次高調波のパワーなどを表し、{H1、H2、H3、…}は、高調波分離後で測量待ちシステムにより生成された信号の基本波のパワー、第二次高調波のパワー、第三次高調波のパワーなどを表す。従来技術では、測量待ちシステムの非線形モデルの出力パワー{R1、R2、R3、…}（図1における点線に示すように）のみ用いて、該測量待ちシステムの非線形モデルの係数を確定することで、該測量待ちシステムの非線形特性を得るが、本出願では、測量待ちシステムの非線形モデルの入力パワー及び出力パワーを用いて、該測量待ちシステムにより生成された高調波のパワーを分離することで、該測量待ちシステムの非線形モデルの係数を確定し、これにより、該測量待ちシステムの非線形特性を得る。
【００１８】
図2は、本実施例における測量待ちシステムの非線形モデルを示す図である。本実施例では、Wienerモデルを例とするが、これに限定されず、他の類型の非線形モデル、例えば、Hammersteinモデルなどを用いても良い。図2に示すように、図1と同様に、{T1、T2、T3、…}、{L1、L2、L3、…}、及び{R1、R2、R3、…}は、信号の基本波、第二次高調波、第三次高調波などのパワーを表す。
【実施例１】
【００１９】
本実施例は、高調波歪み分離方法を提供する。図3は、該方法を示す図である。図3に示すように、該方法は、次のようなステップを含む。
【００２０】
ステップ301：測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定し；
ステップ302：前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離する。
【００２１】
本実施例では、先ず、測量待ちシステムの線形フィルタレスポンスを測量しても良く、例えば、図4に示すように、複数の周波数点の正弦（サイン）信号（テスト信号）を送信することにより、該測量待ちシステムの線形フィルタレスポンスを測量しても良いが、本実施例では、具体的な測量方法について限定しない。
【００２２】
その後、テスト信号のパワーを変更し、複数組の測量待ちシステムの入力高調波のパワー（入力パワーと略称する）及び出力高調波のパワー（出力パワーと略称する）、例えば、図2に示す{T1、T2、T3、…}及び{R1、R2、R3、…}を測量する。ここで、テスト信号は、例えば、ノイズのような信号であるが、本実施例では、具体的な測量方法について限定しない。また、一組の測量結果とは、１つのテスト信号のパワーの場合に対応して得られた{T1、T2、T3、…}及び{R1、R2、R3、…}を指し、複数組の測量結果とは、複数のテスト信号のパワーの場合に対応して得られた{T1、T2、T3、…}及び{R1、R2、R3、…}を指す。
【００２３】
続いて、設定されている測量待ちシステムの非線形モデルに基づいて、測得された前記入力パワー、出力パワー、及び前述の線形フィルタレスポンスを用いて、該非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー及び出力パワー、例えば、図2に示す{L1、L2、L3、…}及び{R1、R2、R3、…}を取得し、そのうち、{L1、L2、L3、…}は、{T1、T2、T3、…}が測量により得された線形フィルタを通過した後に得られたものであり、{R1、R2、R3、…}は、該測量待ちシステムの非線形モデルにより直接得られたものである。
【００２４】
本実施例では、基本波について、出力基本波R1を用いて、基本波パワーH1を計算及び分離することができ、例えば、H1=R1である。また、分離された基本波パワーH1を用いて、基本波振幅Aを計算により得ることもでき、例えば、
【数１】
2017139762.tif 000003
である。
【００２５】
本実施例では、該測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー、及び基本波振幅を得たら、固有高調波及び生成高調波の位相差を確定することができる。本実施例では、具体的な確定方法について限定しない。以下、幾つかの実施方式を以て説明を行うが、当業者が理解すべきは、測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー、及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定し得る任意の方法は、本発明の保護の範囲に属するということである。
【００２６】
本実施例では、第二次高調波について、各組のメモリレス非線形伝達関数の入力パワーL2、出力パワーR2及び基本波振幅Aを用いて、固有高調波及び生成高調波の位相差Δφ2を確定することができる。例えば、次の式で上述の位相差Δφ2を確定する。即ち、
【数２】
2017139762.tif 000004
である。
【００２７】
ここで、R2は、第二次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L2は、第二次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c2は、第二次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される。
【００２８】
本実施例では、位相差Δφ2の確定とは、上の式の中の
（外１）
2017139762.tif 000005
及び
（外２）
2017139762.tif 000006
が、できるだけ、最も良い線形関係を有するように、Δφ2の値を選択することを指す。
【００２９】
１つの実施方式では、先ず、0乃至2πの範囲から複数の角度Δφ2を選択し、その後、各Δφ2を用いて、上の式の中の
（外３）
2017139762.tif 000007
及び
（外４）
2017139762.tif 000008
に対して線形（直線）フィッティングを行い、そして、フィッティングの相関係数が最大であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択する。
【００３０】
他の実施方式では、先ず、0乃至2πの範囲から複数の角度Δφ2を選択し、その後、各Δφ2を用いて、
（外５）
2017139762.tif 000009
及び
（外６）
2017139762.tif 000010
に対して除算を行い、そして、各組の商の分散（variance）が最小であるようにさせるΔφ2を、第二次高調波に対する前記位相差として選択する。
【００３１】
本実施方式では、所定のステップ長さ（step  length）を以て上述の角度を選択しても良く、ランダムに選択しても良く、又は、他のポリシーに基づいて選択しても良いが、本実施例では、これについて限定しない。本実施方式により、第二次高調波に対応する固有高調波及び生成高調波間の位相差Δφ2を取得し、これにより、測量待ちシステムにより生成された第二次高調波のパワーH2を分離することができる。
【００３２】
１つの実施方式では、次の式で該測量待ちシステムにより生成された第二次高調波のパワーH2を分離しても良い。即ち、
【数３】
2017139762.tif 000011
である。
【００３３】
該式の中の各項目の意味は、上述と同様であるため、ここでは、その詳細な説明を省略する。
【００３４】
本実施例では、第三次高調波について、その方法は、第二次高調波の場合と同様であり、即ち、各組のメモリレス非線形伝達関数の入力パワーL3、出力パワーR3及び基本波振幅Aを用いて、固有高調波及び生成高調波の位相差Δφ3を確定することができる。例えば、次の式で上述の位相差Δφ3を確定する。即ち、
【数４】
2017139762.tif 000012
である。
【００３５】
ここで、R3は、第三次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L3は、第三次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c3は、第三次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される。
【００３６】
本実施例では、第二次高調波の方法と同様に、１つの実施方式では、0乃至2πの範囲から複数の角度Δφ3を選択し、その後、各Δφ3を用いて、上の式の中の
（外７）
2017139762.tif 000013
及び
（外８）
2017139762.tif 000014
に対して線形フィッティングを行い、そして、フィッティングの相関係数が最大であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択しても良い。
【００３７】
他の実施方式では、0乃至2πの範囲から複数の角度Δφ3を選択し、その後、各Δφ3を用いて、
（外９）
2017139762.tif 000015
及び
（外１０）
2017139762.tif 000016
に対して除算を行い、そして、各組の商の分散が最小であるようにさせるΔφ3を、固有高調波及び第三次高調波間の前記位相差として選択しても良い。
【００３８】
上述のように、本実施例では、この２種類の方法に限定されない。また、上述のように、本実施方式では、所定のステップ長さを以て上述の角度を選択しても良く、ランダムに選択しても良く、又は、他のポリシーに従って選択しても良いが、本実施例では、これについて限定しない。
【００３９】
本実施例では、第三次高調波に対応する固有高調波及び生成高調波間の位相差Δφ3を得たら、測量待ちシステムにより生成された第三次高調波のパワーH3を分離することができる。例えば、次の式で、該測量待ちシステムにより生成された第三次高調波のパワーH3を分離しても良い。即ち、
【数５】
2017139762.tif 000017
である。
【００４０】
該式の中の各項目の意味は、上述と同様であるため、ここでは、その詳しい説明を省略する。
【００４１】
以上、第二次高調波及び第三次高調波を以て本実施例の高調波分離について説明したが、本実施例の方法は、これに限定されない。具体的な実施過程では、ニーズに応じて、類似した方法を用いて、より高い次数の高調波の高調波分離を行い、分離後の、測量待ちシステムにより生成された高調波のパワー{H1、H2、H3、…}を得ても良い。これにより、これらの高調波のパワーを用いて測量待ち素子の非線形特性を取得し、対処することができる。
【００４２】
図5は、測量待ちシステムの第二次高調波のパワーの測量のパフォーマンス検証を示す図である。図5に示すように、本出願の方法により測量された高調波のパワーは、より正確なものである。
【００４３】
本実施例に係る方法により、測量待ちシステムにより生成された高調波のパワーを分離し、これにより、非線形測量が比較的高い正確度を有するようにさせることができる。
【実施例２】
【００４４】
本発明の実施例は、さらに、非線形特性確定方法を提供する。図6は、該方法を示す図である。図6に示すように、該方法は、次のようなステップを含む。
【００４５】
ステップ601：測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定し；
ステップ602：前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離し；
ステップ603：分離された、前記測量待ちシステムにより生成された高調波のパワーを用いて、前記測量待ちシステムの非線形特性を確定する。
【００４６】
本実施例では、ステップ601及びステップ602の具体的な実施方式は、実施例1のステップ301及びステップ302と同じであるため、その内容は、ここに合併され、ここでは、重複説明が省略される。
【００４７】
本実施例では、分離された、測量待ちシステムにより生成された高調波のパワーを用いて、該測量待ちシステムの非線形特性を確定することができる。本実施例では、具体的な確定方法について限定せず、実施可能な任意の方法を用いても良い。
【００４８】
本実施例に係る方法により、分離された、測量待ちシステムにより生成された高調波のパワーを用いて、非線形特性の測量を行うことで、比較的高い正確度を得ることができる。
【実施例３】
【００４９】
本実施例は、高調波歪み分離装置を提供する。該装置は、送信機又は受信機又は通信システムに配置され得る。該装置が問題を解決する原理は、実施例1の方法と同様であるため、その具体的な実施は、実施例1の方法の実施を参照することができ、内容が同じである重複説明は、省略される。
【００５０】
図7は、該装置を示す図である。図7に示すように、該装置700は、確定ユニット701及び分離ユニット702を含む。
【００５１】
本実施例では、該確定ユニット701は、測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定する。該分離ユニット702は、前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離する。
【００５２】
１つの実施方式では、図7に示すように、該確定ユニット701は、第一確定モジュール7011を含み、該第一確定モジュール7011は、次のような式で、第二次高調波に対応する固有高調波及び生成高調波間の位相差Δφ2を確定することができる。即ち、
【数６】
2017139762.tif 000018
である。
【００５３】
そのうち、R2は、第二次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L2は、第二次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c2は、第二次高調波に対応する非線形モデルの係数であり、かつ定数と仮定される。
【００５４】
該第一確定モジュール7011の１つの実施方式では、図8に示すように、該第一確定モジュール7011は、第一選択モジュール801、第一線形フィッティングモジュール802、及び第二選択モジュール803を含む。該実施方式では、該第一選択モジュール801は、0乃至2πの範囲から複数の角度Δφ2を選択することができるが、本実施例では、具体的な選択方式について限定しない。該第一線形フィッティングモジュール802は、各Δφ2を用いて、上の式の中の
（外１１）
2017139762.tif 000019
及び
（外１２）
2017139762.tif 000020
に対して線形フィッティングを行うことができる。該第二選択モジュール803は、フィッティングの相関係数が最大であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択することができる。具体的には、上述と同様であるため、ここでは、その詳しい説明を省略する。
【００５５】
図9は、該実施方式の第一確定モジュール7011のハードウェア構成を示す図である。本実施例では、これに限定されない。該実施方式は、ソフトウェアにより実現することもできる。
【００５６】
該第一確定モジュール7011の他の実施方式では、図10に示すように、該第一確定モジュール7011は、第三選択モジュール1001、第一除算モジュール1002、及び第四選択モジュール1003を含む。該実施方式では、該第三選択モジュール1001は、0乃至2πの範囲から複数の角度Δφ2を選択することができるが、本実施例では、具体的な選択方式について限定しない。該第一除算モジュール1002は、各Δφ2を用いて、
（外１３）
2017139762.tif 000021
及び
（外１４）
2017139762.tif 000022
に対して除算を行うことができる。該第四選択モジュール1003は、各組の商の分散が最小であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択することができる。
【００５７】
図11は、該実施方式の第一確定モジュール7011のハードウェア構成を示す図である。本実施例では、これについて限定しない。該実施方式は、ソフトウェアにより実現することもできる。
【００５８】
本実施方式では、図7に示すように、該分離ユニット702は、第一分離モジュール7021を含み、それは、次のような式で、前記測量待ちシステムにより生成された第二次高調波のパワーH2を分離する。即ち、
【数７】
2017139762.tif 000023
である。
【００５９】
そのうち、各項目の意味は、上述と同様である。
【００６０】
図12は、該第一分離モジュール7021のハードウェア構成を示す図である。本実施例では、これについて限定しない。該実施方式は、ソフトウェアにより実現することもできる。
【００６１】
他の実施方式では、図7に示すように、該確定ユニット701は、さらに、第二確定モジュール7012を含んでも良く、該第二確定モジュール7012は、次のような式で、固有高調波及び第三次高調波間の位相差Δφ3を確定することができる。即ち、
【数８】
2017139762.tif 000024
である。
【００６２】
そのうち、R3は、第三次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L3は、第三次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c3は、第三次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される。
【００６３】
該第二確定モジュール7012の１つの実施方式では、第一確定モジュール7011の第１の実施方式と同様であり、図13に示すように、該第二確定モジュール7012は、第五選択モジュール1301、第二線形フィッティングモジュール1302、及び第六選択モジュール1303を含む。該第五選択モジュール1301は、0乃至2πの範囲から複数の角度Δφ3を選択し、本実施例では、具体的な選択方式について限定しない。第二線形フィッティングモジュール1302は、各Δφ3を用いて、上の式の中の
（外１５）
2017139762.tif 000025
及び
（外１６）
2017139762.tif 000026
に対して線形フィッティングを行う。第六選択モジュール1303は、フィッティングの相関係数が最大であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択する。
【００６４】
該第二確定モジュール7012の他の実施方式では、第一確定モジュール7011の第２の実施方式と同様であり、図14に示すように、該第二確定モジュール7012は、第七選択モジュール1401、第二除算モジュール1402、及び第八選択モジュール1403を含む。該第七選択モジュール1401は、0乃至2πの範囲から複数の角度Δφ3を選択し、本実施例では、具体的な選択方式について限定しない。該第二除算モジュール1402は、各Δφ3を用いて、
（外１７）
2017139762.tif 000027
及び
（外１８）
2017139762.tif 000028
に対して除算を行う。該第八選択モジュール1403は、各組の商の分散が最小であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択する。
【００６５】
本実施方式では、図7に示すように、該分離ユニット702は、さらに、第二分離モジュール7022を含み、それは、次のような式で、前記測量待ちシステムにより生成された第三次高調波のパワーH3を分離する。即ち、
【数９】
2017139762.tif 000029
である。
【００６６】
そのうち、各項目の意味は、上述と同様である。
【００６７】
本実施方式では、該第二確定モジュール7012及び該第二分離モジュール7022のハードウェアの実現方式は、前の実施方式と同様であるため、ここでは、その詳しい説明を省略する。
【００６８】
本実施例に係る装置により、測量待ちシステムにより生成された高調波のパワーを分離し、これにより、非線形測量が比較的高い正確度を有するようにさせることができる。
【実施例４】
【００６９】
本実施例は、非線形特性確定装置を提供する。該装置は、送信機又は受信機又は通信システムに配置され得る。該装置が問題を解決する原理は、実施例2の方法と同様であるため、その具体的な実施は、実施例2の方法の実施を参照することができ、内容が同じである重複説明は、省略される。
【００７０】
図15は、該装置を示す図である。図15に示すように、該装置1500は、第一確定ユニット1501、分離ユニット1502、及び第二確定ユニット1503を含む。
【００７１】
本実施例では、該第一確定ユニット1501は、測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定する。該分離ユニット1502は、前記位相差を用いて、測量待ちシステムにより生成された高調波のパワーを分離する。該第二確定ユニット1503は、分離された、前記測量待ちシステムにより生成された高調波のパワーを用いて、前記測量待ちシステムの非線形特性を確定する。
【００７２】
本実施例では、該第一確定ユニット1501及び該分離ユニット1502は、それぞれ、実施例3の確定ユニット701及び分離ユニット702により実現することができるため、その内容は、ここに合併され、ここでは、その詳しい説明を省略する。
【００７３】
本実施例に係る装置により、分離された、測量待ちシステムにより生成された高調波のパワーを用いて、非線形特性の測量を行うことで、比較的高い正確度を得ることができる。
【実施例５】
【００７４】
本実施例は、さらに、受信機を提供する、図16は、該受信機を示す図である。図16に示すように、該受信機1600は、実施例3に記載の高調波歪み分離装置又は実施例4に記載の非線形特性確定装置を含む。実施例3及び実施例4において該高調波歪み分離装置700及び該非線形特性確定装置1500について既に詳細に説明しているので、その内容は、ここに合併され、ここでは、その詳しい説明を省略する。
【００７５】
図17は、本発明の実施例における受信機のシステム構成を示す図である。図17に示すように、受信機1700は、フロントエンドを含み、フロントエンドは、入力された光信号を２つの偏光状態（polarization  state）上でのベースバンド信号に変換するためのものであり、本発明の実施例では、該２つの偏光状態は、H偏光状態及びV偏光状態を含む。
【００７６】
図17に示すように、該フロントエンドは、ローカル発振レーザー1710、光周波数混合器（Optical  90deg  hybrid）1701、光電検出器（O/E）1702、1704、1706及び1708、ADC  1703、1705、1707及び1709、分散補償器（dispersion  compensator）1711、イコライザー1712、及び高調波歪み分離装置又は非線形特性確定装置1713を含み、そのうち、高調波歪み分離装置又は非線形特性確定装置1713の構造及び機能は、実施例3又は実施例4と同じであるため、ここでは、この詳細な説明を省略する。ローカル発振レーザー1710は、ローカル光源を提供するために用いられる。光信号は、光周波数混合器（Optical  90deg  hybrid）1701、光電検出器（O/E）1702及び1704、ADC  1703及び1705を経て、１つの偏光状態上でのベースバンド信号に変換される。また、該光信号は、光周波数混合器（Optical  90deg  hybrid）1701、光電検出器（O/E）1706及び1708、ADC  707及び1709を経て、もう１つの偏光状態上でのベースバンド信号に変換される。その具体的なプロセスは、従来技術と同様であるため、ここでは、その詳しい説明を省略する。
【００７７】
また、周波数オフセット及び位相ノイズがOSNRの推定に影響すれば、受信機1700には、周波数オフセット補償器及び位相ノイズ補償器（図示せず）が含まれても良い。
【００７８】
本実施例に係る受信機により、測量待ちシステムにより生成された高調波のパワーを分離し、該高調波のパワーを用いて非線形特性の測量を行うことで、比較的高い正確度を得ることができる。
【実施例６】
【００７９】
本実施例は、さらに、送信機を提供する。図18は、本実施例における送信機の構成を示す図である。図18に示すように、該送信機1800は、実施例3に記載の高調波歪み分離装置又は実施例4に記載の非線形特性確定装置を含む。実施例3及び実施例4において該高調波歪み分離装置700及び該非線形特性確定装置1500について既に詳細に説明しているため、その内容は、ここに合併され、ここでは、その詳しい説明を省略する。
【００８０】
図19は、本実施例における送信機のシステム構成を示す図である。図19に示すように、送信機1900は、信号生成器1901、信号設置ユニット1902、DACユニット1903、及び光変調器1904を含む。
【００８１】
そのうち、信号生成器1901は、送信データに基づいてデジタル信号を生成し；信号設置ユニット1902は、生成されたデジタル信号においてパイロット信号を設置し；DACユニット1903は、該デジタル信号に対してDACを行い；光変調ユニット1904は、該DACユニット1903の変換後の信号を変調信号として、光に対して変調を行う。
【００８２】
また、該送信機は、さらに、選択ユニット1905を含んでも良い。光変調ユニット1904は、選択ユニット1905により選択されたコードワードに基づいて、信号の変調を行い、これにより、異なる偏光状態上でのパイロット信号の周波数が異なるようにさせる。
【００８３】
また、該送信機は、さらに、高調波歪み分離装置又は非線形特性確定装置1906を含んでも良く、それは、測量待ちシステムに対して高調波歪み分離又は非線形特性推定を行うために用いられる。
【００８４】
本実施例に係る受信機により、測量待ちシステムにより生成された高調波のパワーを分離し、該高調波のパワーを用いて非線形特性の測量を行うことで、比較的高い正確度を得ることができる。
【実施例７】
【００８５】
本発明の実施例は、さらに、通信システムを提供する。図20は、本実施例における通信システムの構造を示す図である。図20に示すように、通信システム2000は、送信機2001、伝送リンク2002及び受信機2003を含む。そのうち、送信機2001の構造及び機能は、実施例6と同じであり、受信機2003の構造及び機能は、実施例5と同じであるため、ここでは、その詳しい説明を省略する。伝送リンク2002は、従来の伝送リンクの構造及び機能を有しても良いので、本発明の実施例では、伝送リンクの構造及び機能について限定しない。
【００８６】
本発明の実施例によって提供される通信システムにより、測量待ちシステムにより生成された高調波のパワーを分離し、該高調波のパワーを用いて非線形特性の測量を行うことで、比較的高い正確度を得ることができる。
【００８７】
本発明の実施例は、さらに、コンピュータ可読プログラムを提供し、そのうち、送信機又は受信機又は通信システム中で前記プログラムを実行する時に、前記プログラムは、コンピュータに、前記送信機又は受信機又は通信システム中で実施例1又は実施例2に記載の方法を実行させる。
【００８８】
本発明の実施例は、さらに、コンピュータ可読プログラムを記憶したコンピュータ読み取り可能な記憶媒体を提供し、そのうち、前記コンピュータ可読プログラムは、コンピュータに、送信機又は受信機又は通信システム中で、実施例1又は実施例2に記載の方法を実行させる。
【００８９】
また、本発明の実施例による装置及び方法は、ソフトウェアにより実現されても良く、ハードウェアにより実現されてもよく、ハードウェア及びソフトウェアの組み合わせにより実現されても良い。また、本発明はこのようなコンピュータ可読プログラムにも関し、即ち、前記プログラムは、ロジック部品により実行される時に、前記ロジック部品に、上述の装置又は構成要素を実現させることができ、又は、前記ロジック部品に、上述の方法又はそのステップを実現させることができる。さらに、本発明は上述のプログラムを記憶するための記憶媒体、例えば、ハードディスク、磁気ディスク、光ディスク、DVD、flashメモリなどにも関する。
【００９０】
また、上述の実施例を含む実施形態に関し、更に以下の付記を開示する。
【００９１】
（付記1）
送信機又は受信機又は通信システムに配置される高調波歪み分離装置であって、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定するための確定ユニット；及び
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離するための分離ユニットを含む、装置。
【００９２】
（付記2）
付記1に記載の装置であって、
前記確定ユニットは、
【数１０】
2017139762.tif 000030
を用いて、第二次高調波に対応する固有高調波及び生成高調波間の位相差Δφ2を確定するための第一確定モジュールを含み、
そのうち、R2は、第二次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L2は、第二次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c2は、第二次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される、装置。
【００９３】
（付記3）
付記2に記載の装置であって、
前記第一確定モジュールは、
0乃至2πの範囲から複数の角度Δφ2を選択するための第一選択モジュール；
各Δφ2を用いて、上の式の中の
（外１９）
2017139762.tif 000031
及び
（外２０）
2017139762.tif 000032
に対して線形フィッティングを行うための第一線形フィッティングモジュール；及び
フィッティングの相関係数が最大であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択するための第二選択モジュールを含む、装置。
【００９４】
（付記4）
付記2に記載の装置であって、
前記第一確定モジュールは、
0乃至2πの範囲から複数の角度Δφ2を選択するための第三選択モジュール；
各Δφ2を用いて、上の式の中の
（外２１）
2017139762.tif 000033
及び
（外２２）
2017139762.tif 000034
に対して除算を行うための第一除算モジュール；及び
各組の商の分散が最小であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択するための第四選択モジュールを含む、装置。
【００９５】
（付記5）
付記2に記載の装置であって、
前記分離ユニットは、
【数１１】
2017139762.tif 000035
を用いて、前記測量待ちシステムにより生成された第二次高調波のパワーH2を分離するための第一分離モジュールを含む、装置。
【００９６】
（付記6）
付記1に記載の装置であって、
前記確定ユニットは、
【数１２】
2017139762.tif 000036
を用いて、第三次高調波に対応する固有高調波及び生成高調波間の位相差Δφ3を確定するための第二確定モジュールを含み、
そのうち、R3は、第三次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L3は、第三次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c3は、第三次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される、装置。
【００９７】
（付記7）
付記6に記載の装置であって、
前記第二確定モジュールは、
0乃至2πの範囲から複数の角度Δφ3を選択するための第五選択モジュール；
各Δφ3を用いて、上の式の中の
（外２３）
2017139762.tif 000037
及び
（外２４）
2017139762.tif 000038
に対して線形フィッティングを行うための第二線形フィッティングモジュール；及び
フィッティングの相関係数が最大であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択するための第六選択モジュールを含む、装置。
【００９８】
（付記8）
付記6に記載の装置であって、
前記第二確定モジュールは、
0乃至2πの範囲から複数の角度Δφ3を選択するための第七選択モジュール；
各Δφ3を用いて、上の式の中の
（外２５）
2017139762.tif 000039
及び
（外２６）
2017139762.tif 000040
に対して除算を行うための第二除算モジュール；及び
各組の商の分散が最小であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択するための第八選択モジュールを含む、装置。
【００９９】
（付記9）
付記6に記載の装置であって、
前記分離ユニットは、
【数１３】
2017139762.tif 000041
を用いて、前記測量待ちシステムにより生成された第三次高調波のパワーH3を分離するための第二分離モジュールを含む、装置。
【０１００】
（付記10）
送信機又は受信機又は通信システムに配置される非線形特性確定装置であって、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定するための第一確定ユニット；
前記位相差を用いて、測量待ちシステムにより生成された高調波のパワーを分離するための分離ユニット；及び
分離された前記測量待ちシステムにより生成された高調波のパワーを用いて、前記測量待ちシステムの非線形特性を確定するための第二確定ユニットを含む、装置。
【０１０１】
（付記11）
高調波歪み分離方法であって、
測量待ちシステムの非線形モデルのメモリレス非線形伝達関数の複数組の入力パワー、出力パワー及び基本波振幅を用いて、固有高調波及び生成高調波の位相差を確定し；及び
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離することを含む、方法。
【０１０２】
（付記12）
付記11に記載の方法であって、
前記固有高調波及び生成高調波の位相差の確定は、
【数１４】
2017139762.tif 000042
を用いて、第二次高調波に対応する固有高調波及び生成高調波間の位相差Δφ2を確定することを含み、
そのうち、R2は、第二次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L2は、第二次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c2は、第二次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される、方法。
【０１０３】
（付記13）
付記12に記載の方法であって、
前記第二次高調波に対応する固有高調波及び生成高調波間の位相差Δφ2の確定は、
0乃至2πの範囲から複数の角度Δφ2を選択し；
各Δφ2を用いて、上の式の中の
（外２７）
2017139762.tif 000043
及び
（外２８）
2017139762.tif 000044
に対して線形フィッティングを行い；及び
フィッティングの相関係数が最大であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択することを含む、方法。
【０１０４】
（付記14）
付記12に記載の方法であって、
前記第二次高調波に対応する固有高調波及び生成高調波間の位相差Δφ2の確定は、
0乃至2πの範囲から複数の角度Δφ2を選択し；
各Δφ2を用いて、上の式の中の
（外２９）
2017139762.tif 000045
及び
（外３０）
2017139762.tif 000046
に対して除算を行い；及び
各組の商の分散が最小であるようにさせるΔφ2を、第二次高調波に対応する前記位相差として選択することを含む、方法。
【０１０５】
（付記15）
付記12に記載の方法であって、
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離することは、
【数１５】
2017139762.tif 000047
を用いて、前記測量待ちシステムにより生成された第二次高調波のパワーH2を分離することを含む、方法。
【０１０６】
（付記16）
付記11に記載の方法であって、
前記固有高調波及び生成高調波の位相差の確定は、
【数１６】
2017139762.tif 000048
を用いて、第三次高調波に対応する固有高調波及び生成高調波間の位相差Δφ3を確定すすることを含み、
そのうち、R3は、第三次高調波に対応する前記メモリレス非線形伝達関数の出力パワーであり、L3は、第三次高調波に対応する前記メモリレス非線形伝達関数の入力パワーであり、Aは、基本波振幅であり、c3は、第三次高調波に対応する非線形モデルの係数であり、且つ定数と仮定される、方法。
【０１０７】
（付記17）
付記16に記載の方法であって、
前記第三次高調波に対応する固有高調波及び生成高調波間の位相差Δφ3の確定は、
0乃至2πの範囲から複数の角度Δφ3を選択し；
各Δφ3を用いて、上の式の中の
（外３１）
2017139762.tif 000049
及び
（外３２）
2017139762.tif 000050
に対して線形フィッティングを行い；及び
フィッティングの相関係数が最大であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択することを含む、方法。
【０１０８】
（付記18）
付記16に記載の方法であって、
前記第三次高調波に対応する固有高調波及び生成高調波間の位相差Δφ3の確定は、
0乃至2πの範囲から複数の角度Δφ3を選択し；
各Δφ3を用いて、上の式の中の
（外３３）
2017139762.tif 000051
及び
（外３４）
2017139762.tif 000052
に対して除算を行い；及び
各組の商の分散が最小であるようにさせるΔφ3を、第三次高調波に対応する前記位相差として選択することを含む、方法。
【０１０９】
（付記19）
付記16に記載の方法であって、
前記位相差を用いて、前記測量待ちシステムにより生成された高調波のパワーを分離することは、
【数１７】
2017139762.tif 000053
を用いて、前記測量待ちシステムにより生成された第三次高調波のパワーH3を分離することを含む、方法。
【０１１０】
以上、本発明の好ましい実施形態を説明したが、本発明はこの実施形態に限定されず、本発明の趣旨を離脱しない限り、本発明に対するあらゆる変更は本発明の技術的範囲に属する。
【図１】
2017139762.tif 000078
【図２】
2017139762.tif 000079
【図３】
2017139762.tif 000080
【図４】
2017139762.tif 000081
【図５】
2017139762.tif 000082
【図６】
2017139762.tif 000083
【図７】
2017139762.tif 000084
【図８】
2017139762.tif 000085
【図９】
2017139762.tif 000086
【図１０】
2017139762.tif 000087
【図１１】
2017139762.tif 000088
【図１２】
2017139762.tif 000089
【図１３】
2017139762.tif 000090
【図１４】
2017139762.tif 000091
【図１５】
2017139762.tif 000092
【図１６】
2017139762.tif 000093
【図１７】
2017139762.tif 000094
【図１８】
2017139762.tif 000095
【図１９】
2017139762.tif 000096
【図２０】
2017139762.tif 000097