JPH0131807B2 - - Google Patents

Info

Publication number
JPH0131807B2
JPH0131807B2 JP58023948A JP2394883A JPH0131807B2 JP H0131807 B2 JPH0131807 B2 JP H0131807B2 JP 58023948 A JP58023948 A JP 58023948A JP 2394883 A JP2394883 A JP 2394883A JP H0131807 B2 JPH0131807 B2 JP H0131807B2
Authority
JP
Japan
Prior art keywords
signal
phase
signals
timing
orthogonal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired
Application number
JP58023948A
Other languages
Japanese (ja)
Other versions
JPS59149439A (en
Inventor
Botaro Hirosaki
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
Nippon Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Electric Co Ltd filed Critical Nippon Electric Co Ltd
Priority to JP58023948A priority Critical patent/JPS59149439A/en
Publication of JPS59149439A publication Critical patent/JPS59149439A/en
Publication of JPH0131807B2 publication Critical patent/JPH0131807B2/ja
Granted legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J11/00Orthogonal multiplex systems, e.g. using WALSH codes

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Description

【発明の詳細な説明】 本発明は直交多重された複数個の直交振幅変調
(以下QAMと略称する)信号を伝送路を介して
受信し、所定の標本化クロツクにて標本化し、ア
ナログ−デジタル(以下ADと略称する)変換、
波操作および離散フーリエ変換操作(以下
DFTと略称する)等により複数個の原基底帯域
信号を復調する直交多重信号の受信系におけるキ
ヤリア・タイミング結合制御方式に関する。
Detailed Description of the Invention The present invention receives a plurality of orthogonally multiplexed amplitude modulation (hereinafter abbreviated as QAM) signals via a transmission line, samples them at a predetermined sampling clock, and converts them from analog to digital. (hereinafter abbreviated as AD) conversion,
Wave operations and discrete Fourier transform operations (below
This invention relates to a carrier timing combination control method in a receiving system for orthogonal multiplexed signals that demodulates a plurality of original baseband signals using DFT (abbreviated as DFT) or the like.

複数個のQAM信号を直交多重して送信し、受
信側にて各ベースバンド信号を復調する直交多重
伝送方式は高能率データ伝送を可能にし、しかも
そのデジタル信号処理過程にDFTを導入すれば
送受信装置が、著しく簡単化できることが、既に
知られている。(例えば、文献として1981年7月
発行の刊行物「IEEE TRANSACTIONS ON
COMMUNICATIONS、Vol.COM−29、No.7」
の第982頁−第989頁記載の論文“AN
ORTHOGONALLY MULTIPLEXED QAM
SYSTEM USING THE DISCRETE
FOURIER TRANSFORM”および文献2とし
て1981年登録の米国特許第4300229号記載の
“TRANSMITTER AND RECEIVER FOR
AN ORTHOGONALLT MULTIPLEXED
QAM SIGNAL OF A SAMPLING RATE
N TIMES THAT OF PAM SIGNALS、
COMPRISING、AN TRANSFORM
PROCESSOR”を参照されたい)。
The orthogonal multiplex transmission method, in which multiple QAM signals are orthogonally multiplexed and transmitted, and each baseband signal is demodulated on the receiving side, enables highly efficient data transmission, and if DFT is introduced into the digital signal processing process, transmission and reception will be possible. It is already known that the device can be significantly simplified. (For example, as a reference, the publication "IEEE TRANSACTIONS ON
COMMUNICATIONS, Vol.COM−29, No.7”
The paper “AN” on pages 982-989 of
ORTHOGONALLY MULTIPLEXED QAM
SYSTEM USING THE DISCRETE
FOURIER TRANSFORM” and “TRANSMITTER AND RECEIVER FOR” described in US Patent No. 4300229, registered in 1981 as Document 2
AN ORTHOGONALLT MULTIPLEXED
QAM SIGNAL OF A SAMPLING RATE
N TIMES THAT OF PAM SIGNALS,
COMPRISING, AN TRANSFORM
(Please refer to “PROCESSOR”).

また、復調された複素ベースバンド信号に対す
る2次元自動等化器も既に提案されており、この
2次元自動等化器を用いれば、伝送路にて被つた
回線歪を完全に等化することができ、更に、静的
なキヤリア位相ずれ、タイミング位相ずれをも同
時に抑圧してしまうことが知られている(例え
ば、文献3として1980年1月発行の刊行物
「IEEE TRANSACTIONS ON
COMMUNICATIONS、Vol.COM−28、No.1」
の第73頁−第83頁記載の論文“AN ANALYSIS
OF AUTOMATIC EQUALIZERS FOR
ORTHOGONALLY MULTIPLEXED QAM
SYSTEMS”を参照されたい)。
Additionally, a two-dimensional automatic equalizer for demodulated complex baseband signals has already been proposed, and if this two-dimensional automatic equalizer is used, it is possible to completely equalize the line distortion experienced in the transmission path. Furthermore, it is known that static carrier phase deviations and timing phase deviations are also suppressed at the same time (for example, in the publication "IEEE TRANSACTIONS ON
COMMUNICATIONS, Vol.COM−28, No.1”
The paper “AN ANALYSIS” on pages 73-83 of
OF AUTOMATIC EQUALIZERS FOR
ORTHOGONALLY MULTIPLEXED QAM
SYSTEMS).

しかし乍ら、前記の自動等化器のみでは、謂る
動的なキヤリア位相変動、タイミング位相変動に
充分な精度で追従し得ない。そうした動的位相変
動を抑圧する方法として送信信号と共にパイロツ
トを送出し、受信側ではこのパイロツトの位相変
動を検出して複素ベースバンド信号に適当な補償
用位相回転を施す方法が、既に知られている(例
えば、文献4として昭和55年特許願4382号を参照
されたい)。
However, the automatic equalizer alone cannot follow so-called dynamic carrier phase fluctuations and timing phase fluctuations with sufficient accuracy. As a method for suppressing such dynamic phase fluctuations, a method is already known in which a pilot is sent along with the transmission signal, and the receiving side detects the phase fluctuation of this pilot and applies appropriate compensation phase rotation to the complex baseband signal. (For example, please refer to Patent Application No. 4382 filed in 1982 as Document 4).

しかし、上記の方法ではパイロツトを伝送する
為の余分な帯域を要することとなり、特に伝送効
率の高さを重視する場合には、これが欠点となる
可能性があつた。
However, the above method requires an extra band for pilot transmission, which may be a drawback, especially when high transmission efficiency is important.

本発明の目的は、以上の観点に立脚し、余分な
パイロツトを使用する事なく、動的位相変動を抑
圧し得るキヤリア・タイミング結合制御方式を提
供することに有る。
An object of the present invention is to provide a carrier timing coupling control system that can suppress dynamic phase fluctuations without using an extra pilot, based on the above-mentioned viewpoints.

即ち、本発明は直交多重されたQAM信号が全
てスタガQAM(同相、直交間にシンボル周期の
丁度半分の遅延差があるQAM)である点、およ
び全てのQAM信号のシンボルタイミングが同期
している点に着目し、これらの性質を巧みに組み
合わせた最尤受信系を構成することにより動的位
相変動を精度良く抑圧し得るキヤリア・タイミン
グ結合制御方式を提供するものである。
That is, the present invention is characterized in that all orthogonally multiplexed QAM signals are staggered QAM (QAM with a delay difference of exactly half the symbol period between in-phase and orthogonal signals), and that the symbol timings of all QAM signals are synchronized. By focusing on these points and configuring a maximum likelihood receiving system that skillfully combines these properties, we provide a carrier timing coupling control system that can suppress dynamic phase fluctuations with high accuracy.

以下、図面を用いて、本発明の原理及びその構
成を説明する。
Hereinafter, the principle of the present invention and its configuration will be explained using the drawings.

第1図は、二次元変調信号の一般的な生成過程
を示すブロツク図であり、1,2は各々同相デー
タ直交データの入力される入力端であり、3,4
は各々同相側ベースバンド整形フイルタ、直交側
ベースバンド整形フイルタであり、5,6は乗算
器、7は加算器を表わしている。8は出力端であ
り、この出力端を介して二次元変調信号が伝送路
へ送出される。第1図において、入力端1,2に
入力される同相データ{ak}、直交データ{bk
のクロツク周期をT秒とし、3で示される同相側
ベースバンド整形フイルタの単位応答をp(t)、
4で示される直交側ベースバンド整形フイルタの
単位応答をq(t)とする。また、変調キヤリア
角周波数をωc(ラジアン/秒)とすれば、出力端
8に得られる二次元変調信号s(t)は、次のよ
うに表現される。
FIG. 1 is a block diagram showing a general process of generating a two-dimensional modulation signal, in which 1 and 2 are input terminals into which in-phase data and quadrature data are input, respectively;
are an in-phase side baseband shaping filter and an orthogonal side baseband shaping filter, respectively, 5 and 6 represent multipliers, and 7 represents an adder. 8 is an output end, and a two-dimensional modulated signal is sent out to the transmission path via this output end. In FIG. 1, in-phase data {a k } and orthogonal data {b k } are input to input terminals 1 and 2.
The clock period of is T seconds, and the unit response of the in-phase baseband shaping filter shown by 3 is p(t),
Let q(t) be the unit response of the orthogonal side baseband shaping filter shown by 4. Further, if the modulation carrier angular frequency is ω c (radian/second), the two-dimensional modulation signal s(t) obtained at the output end 8 is expressed as follows.

s(t)= 〓k {akp(t−kT)cosωct +bkq(t−kT)sin−ωct} ここで、上記の二次元変調信号s(t)に対す
る受信操作について考えてみる。簡単のため、伝
送路における回線歪は無視し得るものとし、専ら
伝送遅延t′およびキヤリアオフセツト△ωを被る
ものとする。この時、受信信号r(t)は、s
(t)が伝送遅延t′とキヤリアオフセツト△ωを
受けて得られる信号ro(t)と白色雑音n(t)と
の和として表わされる。ここで、ro(t)は、ω′c
=ωc+△ωとして、 ro(t)= 〓k {akp(t−t′−kT)cosωc′(t−t′) +bkq(t−t′−kT)sinωc′(t−t′)} と表わされる。
s(t)=〓 k {a k p(t-kT) cosω c t +b k q(t-kT) sin-ω c t} Here, regarding the reception operation for the above two-dimensional modulated signal s(t) I'll think about it. For the sake of simplicity, it is assumed that line distortion in the transmission path can be ignored, and that only transmission delay t' and carrier offset Δω are suffered. At this time, the received signal r(t) is s
(t) is expressed as the sum of a signal ro(t) obtained by receiving a transmission delay t' and a carrier offset Δω, and a white noise n(t). Here, ro(t) is ω′ c
As =ω c +△ω, ro(t)= 〓 k {a k p(t−t′−kT) cosω c ′(t−t′) +b k q(t−t′−kT) sinω c ′ (t-t')}.

t′=mT+τ(但し、|τ|T/2) Θ=△ωt−ω′ct′ とすれば、結局ro(t)は、次のように表現され
る。
If t'=mT+τ (where |τ|T/2) Θ=Δωt−ω' c t', ro(t) can be expressed as follows.

ro(t)= 〓k {akp(t−τ−kT)cos(ωct+Θ) +bkp(t−τ−kT)cos(ωct+Θ} 即ち、ro(t)は受信側では未知量であるタイ
ミングずれτおよびキヤリア位相ずれΘを受けて
いることになる。実際の受信信号r(t)は、ro
(t)が、雑音n(t)に埋ずもれたもの、即ち、
r(t)=ro(t)+n(t)であるが、受信側では
この信号r(t)から元のデータ系列{ak}{bk
を復元せねばならない。この復元に際しては、ま
ず、タイミングずれτ、キヤリア位相ずれΘを補
償しなければならないが、その方法として従来、
最尤推定法が知られている。以下に最尤推定法の
基本的な考え方を示す。
ro(t)= 〓 k {a k p(t-τ-kT) cos(ω c t+Θ) +b k p(t-τ-kT) cos(ω c t+Θ} In other words, ro(t) is This means that the timing shift τ and carrier phase shift Θ are unknown quantities.The actual received signal r(t) is ro
(t) is buried in noise n(t), that is,
r(t)=ro(t)+n(t), but on the receiving side, from this signal r(t), the original data sequence {a k }{b k }
must be restored. When performing this restoration, it is first necessary to compensate for the timing deviation τ and the carrier phase deviation Θ.
Maximum likelihood estimation method is known. The basic idea of maximum likelihood estimation method is shown below.

いま、{ak}、{bk}、τ、Θの推定側を{a^k}、
{b^k}、τ^、Θ^とする。これらの推定値を用いて、
次の基準信号u(t)を設定する。
Now, let the estimated sides of {a k }, {b k }, τ, and Θ be {a^ k },
Let {b^ k }, τ^, Θ^. Using these estimates,
Set the next reference signal u(t).

u(t)= 〓k {a^kp(t−τ−kT)cos(ωct+Θ^) +b^kq(t−τ^−kT)sin(ωct+Θ^)} 最尤推定の立場では、“送信側にてu(t)が送
られたと仮定した時に受信信号r(t)を得る事
後確率が最大となるような推定値を以て最も確か
らしい推定値”と判断する。この推定過程におい
て、データ{ak}、{bk}が正確に知られていると
いう前提の上で推定を行なう方法を判定帰還形
(以後DA形と略称する)と呼び、{ak}、{bk}に
ついては、その統計的性質しか判つていないとい
う前提の上で推定を行なう方法を判定無帰還形
(以後NDA形と略称する)と呼ぶ。一般にτ、Θ
の微小変動に対する制御系の定常パフオーマンス
はDA形の方が優れているが、τ、Θの変動が大
きい時は、データ判定結果に誤りを生じ易く、そ
のため制御系のパフオーマンスも著しく劣化す
る。また、制御系の初期引込時には正しい判定結
果が得られていないため、DA形では引込特性が
劣化する。即ちτ、Θに対する高いトラツキング
能力を期待するならば、NDA形最尤推定法を用
いるのが望ましい。
u(t)= 〓 k {a^ k p(t-τ-kT) cos(ω c t+Θ^) +b^ k q(t-τ^-kT) sin(ω c t+Θ^)} Maximum likelihood estimation From our point of view, the most probable estimate is the one that maximizes the a posteriori probability of obtaining the received signal r(t) when it is assumed that u(t) is sent on the transmitting side. In this estimation process, the method in which estimation is performed on the assumption that the data {a k }, { b k } are accurately known is called the decision feedback type (hereinafter abbreviated as DA type). , {b k }, a method that performs estimation on the premise that only its statistical properties are known is called a decision-free feedback method (hereinafter abbreviated as NDA method). Generally τ, Θ
The DA type is better in steady-state performance of the control system against small fluctuations in , but when the fluctuations in τ and Θ are large, errors tend to occur in the data judgment results, and as a result, the performance of the control system deteriorates significantly. In addition, since correct judgment results are not obtained during the initial pull-in of the control system, the pull-in characteristics deteriorate in the DA type. That is, if high tracking ability for τ and Θ is expected, it is desirable to use the NDA type maximum likelihood estimation method.

前記二次元変調信号に対するNDA形尤度関数
L(τ^、Θ^)は、n(t)を白色ガウス雑音として
次式で定義される。
The NDA type likelihood function L(τ^, Θ^) for the two-dimensional modulated signal is defined by the following equation, where n(t) is white Gaussian noise.

但し、σ2はn(t)の分散であり、E{・}は、
{ak}{bk}に関する平均化操作を表わしている。
However, σ 2 is the variance of n(t), and E{・} is
It represents an averaging operation on {a k }{b k }.

L(τ^、Θ^)は確率密度関数であり、必ず正の値
をとり、しかもτ^、Θ^の最適点にて唯一の最大値
を持つので、L(τ^、Θ^)の対数をとつたもの、即
ち、対数尤度関数Ω(τ^、Θ^)も唯一の最大値を
有する。従つて、τ^、Θ^に関する制御系は、Ω
(τ^、Θ^)を最大にすべく、τ^、Θ^を制御すればよ
い事が判る。ここで、{ak}、{bk}を、ガウス分
布をする一種の雑音とみなせば、良く知られてい
るようにΩ(τ^、Θ^)は、次式で表わされる。
L(τ^, Θ^) is a probability density function, always takes a positive value, and has a unique maximum value at the optimal point of τ^, Θ^, so L(τ^, Θ^) The logarithm, ie, the log-likelihood function Ω(τ^, Θ^), also has a unique maximum value. Therefore, the control system for τ^ and Θ^ is Ω
It turns out that it is sufficient to control τ^ and Θ^ in order to maximize (τ^, Θ^). Here, if {a k } and {b k } are regarded as a type of noise having a Gaussian distribution, Ω(τ^, Θ^) is expressed by the following equation, as is well known.

Ω(τ^、Θ^)=定数× 〓m v2 1(mT+τ^) ・ 〓m v2 2(nT+τ^) 但し、v1(t)=∫T0w1(s)・p(s−t)ds v2(t)=∫T0w1(s)q(s−t)ds w1(t)=r(t)・cos(ωct+Θ^) w2(t)=r(t)×sin(ωct+Θ^) 即ち、v1(t)は、同相側復調信号w1(t)を、
同相側整形フイルタに対する整合フイルタに通し
て得られる出力であり、v2(t)は直交側復調信
号w2(t)を直交側整形フイルタに対する整合フ
イルタに通して得られる出力である。
Ω (τ^, Θ^) = constant × 〓 m v 2 1 (mT + τ^) ・ 〓 m v 2 2 (nT + τ^) However, v 1 (t) = ∫ T0 w 1 (s)・p(s- t)ds v 2 (t)=∫ T0 w 1 (s)q(s-t)ds w 1 (t)=r(t)・cos(ω c t+Θ^) w 2 (t)=r(t ) × sin (ω c t + Θ^) In other words, v 1 (t) is the in-phase side demodulated signal w 1 (t),
v 2 (t) is an output obtained by passing the orthogonal side demodulated signal w 2 (t) through a matching filter for the orthogonal side shaping filter.

τ^、Θ^に対する制御は次のアルゴリズムにて実
行される。即ち、現在の推定値τ^k、Θ^kより次の
推定値τ^k+1、Θ^k+1を得るには、K〓、K〓をループ
ゲインとして、 とする。ここで、∂Ω/∂τ^、∂Ω/∂Θ^を計算す
ると ∂Ω/∂τ定数× 〓m {v1(mT+τ^)v〓(mT+τ^) +v2(mT+τ^)v〓2(mT+τ^)} ∂Ω/∂Θ^定数× 〓m {−v1(mT+τ^) +y2(mT+τ^)+v2(mT+τ^)y1(mT+τ^)} {担し、y1(t)=∫T0w1(s)q(s-t)ds y2(t)=∫T0w2(s)p(s-t)ds} となる。
Control over τ^ and Θ^ is performed using the following algorithm. That is, to obtain the next estimated values τ^ k +1 and Θ^ k+1 from the current estimated values τ^ k and Θ^ k, set K〓 and K〓 as loop gains, shall be. Here, when calculating ∂Ω/∂τ^ and ∂Ω/∂Θ^, ∂Ω/∂τ constant × 〓 m {v 1 (mT+τ^)v〓 (mT+τ^) +v 2 (mT+τ^)v〓 2 (mT+τ^)} ∂Ω/∂Θ^ constant× 〓 m {−v 1 (mT+τ^) +y 2 (mT+τ^)+v 2 (mT+τ^)y 1 (mT+τ^)} {carry, y 1 (t )=∫ T0 w 1 (s)q(st)ds y 2 (t)=∫ T0 w 2 (s)p(st)ds}.

通常のQAM信号では、p(t)=q(t)であ
るからy1(t)=v1(t)、y2(t)=v2(t)となり

∂Ω/∂Θ^が恒等的に0となつてしまう。
In a normal QAM signal, p(t)=q(t), so y 1 (t)=v 1 (t), y 2 (t)=v 2 (t),
∂Ω/∂Θ^ becomes equal to 0.

即ち、良く知られているように通常のQAM信
号に対しては、NDA形のキヤリア位相制御を施
すことができない。これに対し、謂るスタガ
QAMにおいては、q(t)=p(t−t0)(但し、
t0は同相−直交間の時間ずれ)であり、τ^、Θ^い
ずれに対してもNDA形の位相制御をかけること
ができ、特にt0=T/2の時、その制御能力が最
大となる。
That is, as is well known, NDA type carrier phase control cannot be applied to normal QAM signals. In contrast, the so-called staga
In QAM, q(t)=p(t-t 0 ) (however,
t 0 is the time difference between in-phase and quadrature), and NDA type phase control can be applied to both τ^ and Θ^, and especially when t 0 = T/2, the control ability is maximum. becomes.

第2図は、t0=T/2のスタガQAM信号に対
するタイミング・キヤリア結合制御系の構成を示
すブロツク図である。即ち、第2図において、入
力端11に入力された受信信号は、乗算器12,
14において、電圧制御発振器17の出力および
その90゜位相シフトしたものとして得られる同相
キヤリア直交キヤリアにより各々逆変調を受け
る。
FIG. 2 is a block diagram showing the configuration of a timing/carrier combination control system for a staggered QAM signal with t 0 =T/2. That is, in FIG. 2, the received signal input to the input terminal 11 is transmitted to the multiplier 12,
At 14, the output signal is inversely modulated by the output of the voltage controlled oscillator 17 and an in-phase carrier and a quadrature carrier obtained by phase-shifting the output by 90 degrees.

同相側出力におよび直交側出力は、更に整合フ
イルタ15および18に各々供給される。整合フ
イルタ15,18の出力は、乗算器19にて相互
の相関がとられ、この相関値により前記電圧制御
発振器17の発振位相が制御される。また、タイ
ミング制御に関しては、タイミング位相を定める
電圧制御発振器25に対し、整合フイルタ15の
出力をT/2秒遅延させた信号と、その微分信号
との相関および整合フイルタ18の出力信号とそ
の微分信号との相関とを加算器24にて加算して
得られた制御信号を供給することにより実行され
ている。サンプラー26,27は、こうして得ら
れたタイミングクロツクにより、同相信号、直交
信号を各々標本化するものであり、それらの出力
は出力端28,29を介して後続の自動等化器等
へ供給される。
The in-phase side output and the quadrature side output are further supplied to matching filters 15 and 18, respectively. The outputs of the matching filters 15 and 18 are correlated with each other in a multiplier 19, and the oscillation phase of the voltage controlled oscillator 17 is controlled by this correlation value. Regarding timing control, for the voltage controlled oscillator 25 that determines the timing phase, the output of the matched filter 15 is delayed by T/2 seconds and its differential signal correlates with the output signal of the matched filter 18 and its differential signal. This is executed by adding the signal and the correlation in an adder 24 and supplying the obtained control signal. The samplers 26 and 27 sample the in-phase signal and the quadrature signal, respectively, using the timing clocks obtained in this way, and their outputs are sent to a subsequent automatic equalizer, etc. via output terminals 28 and 29. Supplied.

第3図は、第2図にて示した制御系を全てデイ
ジタル処理にて実現するための構成を示したブロ
ツク図であり、端子31に入力された受信信号は
サンプラ32にて電圧制御発振器46の出力とし
て得られる標本化クロツクで標本化され、乗算器
33,34に供給される。乗算器33,34は、
入力信号に対して前記標本化クロツクにて読出専
用メモリ35を駆動して得られる余弦サンプル値
系列および正弦サンプル値系列を各々掛け合わ
せ、同相信号および直交信号を得る。こうして得
られた同相信号および直交信号は位相回転回路3
8にて位相回転を受けるが、この位相回転量は乗
算器40の出力として得られる同相−直交間の相
関値によつて制御される。位相回転回路38の出
力のうち同相信号は、遅延回路39にてT/2秒
の遅延を受けた後、出力端47に至る。これに対
し直交信号はそのまま出力端48に至る。電圧制
御発振器46に対する制御信号は加算器45より
供給されるがこの加算器45はT/2秒遅延を受
けた同相信号とその微分信号との相関値および直
交信号とその微分信号との相関値との加算をおこ
なう。
FIG. 3 is a block diagram showing a configuration for realizing the control system shown in FIG. 2 entirely through digital processing. The signal is sampled by the sampling clock obtained as the output of the signal, and is supplied to multipliers 33 and 34. The multipliers 33 and 34 are
The input signal is multiplied by a cosine sample value series and a sine sample value series obtained by driving the read-only memory 35 with the sampling clock to obtain an in-phase signal and a quadrature signal. The in-phase signal and quadrature signal thus obtained are sent to the phase rotation circuit 3.
8, the amount of phase rotation is controlled by the in-phase-quadrature correlation value obtained as the output of the multiplier 40. The in-phase signal among the outputs of the phase rotation circuit 38 is delayed by T/2 seconds in the delay circuit 39 and then reaches the output terminal 47 . On the other hand, the orthogonal signal reaches the output terminal 48 as it is. A control signal for the voltage controlled oscillator 46 is supplied from an adder 45, and this adder 45 calculates the correlation value between the in-phase signal delayed by T/2 seconds and its differential signal, and the correlation between the quadrature signal and its differential signal. Performs addition with the value.

本発明によるキヤリア・タイミング結合制御方
式の適用対象は、直交QAM伝送系であるが、こ
の直交QAM伝送系は基本的には前記のスタガ
QAM信号を複数個互いにスペクトラム重なりを
許容しつつ周波数多重して伝送する方式である。
The carrier timing coupling control method according to the present invention is applied to an orthogonal QAM transmission system, but this orthogonal QAM transmission system is basically the staggered
This is a method in which multiple QAM signals are frequency-multiplexed and transmitted while allowing spectrum overlap.

第4図に直交QAM伝送系の原理的な構成図を
示し、第5図に伝送路に送出される直交QAM信
号のスペクトラムの例を示す。
FIG. 4 shows the basic configuration of the orthogonal QAM transmission system, and FIG. 5 shows an example of the spectrum of the orthogonal QAM signal sent to the transmission path.

第4図において、端子51,52には、各々第
1チヤネル同相データ{a(1) k}、直交データ{b(1) k
がクロツク周期T秒にて入力される。同相データ
{a(1) k}は遅延回路59にてT/2秒の遅延を受け
た後、成形フイルタ63を通つて乗算器71に至
る。一方、直交データ{b(1) k}は、成形フイルタ
64を通つて乗算器72に至る。成形フイルタ6
3および64は同一の周波数伝達特性G(ω)を
有しており、このG(ω)はいわゆるルートナイ
キスト形に設定されている。G(ω)のインパル
ス応答をg(t)とすれば、成形フイルタ63の
出力信号u(1) R(t)および成形フイルタ64の出力
信号u(1) I(t)は、各々 u(1) R(t)= 〓k a(1) k・g(t−T/2−kT) ……(1) u(1) I(t)=Σb(1) k・g(t−kT) ……(2) と表わされる。乗算器71、乗算器72、加算器
79はキヤリア角周波数ω1の直交振幅変調器を
構成しており、加算器79の出力には第1QAM
信号q(1)(t)は、 q(1)(t)=u(1) R(t)cosω1t+u(1) I(t)sinω1t =Re{〔u(1) R(t)+ju(1) I(t)〕e-j1t
……(3) と表わされる。但しRe{・}は、実数部のみをと
る演算を表わしている。全く同様にして加算器8
0,81,82の出力には各々第2QAM信号q(2)
(t)、第3QAM信号q(3)(t)、第4QAM信号q(4)
(t)が得られる。これらの信号は次式で表わさ
れる。
In FIG. 4, terminals 51 and 52 have first channel in-phase data {a (1) k } and quadrature data {b (1) k }, respectively.
is input with a clock cycle of T seconds. The in-phase data {a (1) k } is delayed by T/2 seconds in the delay circuit 59 and then passes through the shaping filter 63 to the multiplier 71 . On the other hand, the orthogonal data {b (1) k } passes through the shaping filter 64 and reaches the multiplier 72 . Molded filter 6
3 and 64 have the same frequency transfer characteristic G(ω), and this G(ω) is set to a so-called root Nyquist shape. If the impulse response of G(ω) is g(t), the output signal u (1) R (t) of the shaping filter 63 and the output signal u (1) I (t) of the shaping filter 64 are respectively u ( 1) R (t) = 〓 k a (1) k・g (t-T/2-kT) ...(1) u (1) I (t)=Σb (1) k・g (t-kT ) ...(2) The multiplier 71, the multiplier 72, and the adder 79 constitute a quadrature amplitude modulator with a carrier angular frequency ω 1 , and the output of the adder 79 has the first QAM
The signal q (1) (t) is q (1) (t) = u (1) R (t) cosω 1 t + u ( 1 ) I (t) sinω 1 t = Re )+ju (1) I (t)〕e -j1t }
...(3) It is expressed as. However, Re{·} represents an operation that takes only the real part. Adder 8 in exactly the same way
The outputs of 0, 81, and 82 each have the second QAM signal q (2)
(t), 3rd QAM signal q (3) (t), 4th QAM signal q (4)
(t) is obtained. These signals are expressed by the following equations.

q(2)(t)=Re{〔u(2) R(t)+ju(2) I(t)〕e-j
2t
……(4) q(3)(t)=Re{〔u(3) R(t)+ju(3) I(t)〕e-j
3t
……(5) q(4)(t)=Re{〔u(4) R(t)+ju(4) I(t)〕e-j
4t
……(6) 但し、 u(2) R(t)= 〓k a(2) kg(t−kT) ……(7) u(2) I(t)= 〓k b(2) kg(t−T/2−kT) ……(8) u(3) R(t)= 〓k a(3) kg(t−T/2−kT) ……(9) u(3) I(t)= 〓k b(3) kg(t−kT) ……(10) u(4) R(t)= 〓k a(4) kg(t−kT) ……(11) u(4) I(t)= 〓k b(4) kg(t−T/2−kT) ……(12) であり、{a(2) k}、{b(2) k}は、各々端子53,54
に入力される第2チヤネル同相データおよび直交
データを、{a(3) k}、{b(3) k}は、端子55,56に
入力される第3チヤネル同相データおよび直交デ
ータを{a(4) k}、{b(4) k}は各々端子57,58に入
力される第4チヤネル同相データおよび直交デー
タを各々示している。ここで各キヤリア角周波数
ω1、ω2、ω3、ω4の間には、ωk−ωk-1=2π/Tの
関係が成立しており、更にT/2秒遅延回路5
9,60,61,62は各々同相側、直交側、同
側、直交側と交互に挿入されている。これによ
り第1乃至第4QAM信号は、周波数軸上で直交
多重配置されることになる。この直交多重信号を
s(t)とすると、s(t)は加算回路83の出力
として得られ、図のように4チヤネル多重の例で
は、 s(t)=4n=1 g(n)(t) ……(13) と表わされる。また、直交多重信号s(t)のス
ペクトラムは、第5図に示す如くなる。即ち、第
5図において、110,111,112,113
は、各々第1、第2、第3、第4QAM信号のス
ペクトラムを表わしており、隣接チヤネル間に
は、必ずスペクトラム重なりが生ずる。直交多重
伝送系の特徴は、こうした隣接スペクトラム重な
りがあつても受信側にて原ベースバンドデータを
誤まり無く復元できる点にあり、このため、高能
率データ伝送系が実現できる。
q (2) (t) = Re {[u (2) R (t) + ju (2) I (t)] e -j
2t
...(4) q (3) (t) = Re {[u (3) R (t) + ju (3) I (t)] e -j
3t
...(5) q (4) (t) = Re {[u (4) R (t) + ju (4) I (t)] e -j
4t
……(6) However, u (2) R (t)= 〓 k a (2) k g(t−kT) ……(7) u (2) I (t)= 〓 k b (2) k g(t-T/2-kT) ……(8) u (3) R (t)= 〓 k a (3) k g(t-T/2-kT) ……(9) u (3) I (t)= 〓 k b (3) k g (t-kT) ……(10) u (4) R (t)= 〓 k a (4) k g (t-kT) ……(11) u (4) I (t) = 〓 k b (4) k g (t-T/2-kT) ...(12), and {a (2) k } and {b (2) k } are , terminals 53 and 54, respectively.
{a (3) k } and {b (3) k } represent the second channel in-phase data and quadrature data input to the terminals 55 and 56, respectively. (4) k } and {b (4) k } respectively indicate fourth channel in-phase data and quadrature data input to terminals 57 and 58, respectively. Here, the relationship ω k −ω k-1 = 2π/T is established between each carrier angular frequency ω 1 , ω 2 , ω 3 , ω 4 , and furthermore, the T/2 second delay circuit 5
9, 60, 61, and 62 are inserted alternately on the in-phase side, the orthogonal side, the same side, and the orthogonal side, respectively. As a result, the first to fourth QAM signals are orthogonally multiplexed on the frequency axis. Letting this orthogonal multiplexed signal be s(t), s(t) is obtained as the output of the adder circuit 83, and in the example of 4-channel multiplexing as shown in the figure, s(t)= 4n=1 g (n ) (t) ...(13) Further, the spectrum of the orthogonal multiplexed signal s(t) is as shown in FIG. That is, in FIG. 5, 110, 111, 112, 113
represent the spectra of the first, second, third, and fourth QAM signals, respectively, and spectrum overlap always occurs between adjacent channels. A feature of the orthogonal multiplex transmission system is that even if such adjacent spectrums overlap, the original baseband data can be restored without error on the receiving side, and thus a highly efficient data transmission system can be realized.

さて、第4図において、受信側では送信側で施
された演算の逆演算がおこなわれる。例えば、乗
算器84,85および成形フイルタ92,93
は、第1QAM信号から同相ベースバンド信号v(1) R
(t)、直交ベースバンド信号v(1) I(t)を復元する
ために用いられる。ここで成形フイルタ92,9
3の周波数伝達特性は、逆信側の成形フイルタ6
3,64等と全く同一のルートナイキスト特性を
有している。従つて、v(1) R(t)、v(1) I(t)は、伝
送路での雑音が無い場合、各々次式で与えられ
る。(但しω1は2π/Tより大としている。) v(1) R(t)=(u(1) R(τ)cosΘ1−u(1) I(t)sin
Θ1)*g(t)+(u(2) R(t)cos(Θ1−2π/Tt
) −u(2) I(t)sin(Θ1−2π/Tt)*g(t
)……(14) v(1) I(t)=(u(1) R(t)sinΘ1+u(1) I(t)cos
Θ1)*g(t)+(u(2) R(τ)sin(Θ1−2π/Tt
) +u(2) I(t)cos(Θ1−2π/Tt))*g(
t)……(15) ここで、Θ1は受信側の第1キヤリアの位相ず
れであり、*は畳み込み演算を表わしている。
(14)、(15)式に(1)、(2)、(7)、(8)式を代入すると
、 v(1) R(t)= 〓k a(1) kh0(t−T/2−kT)cosΘ1− 〓k b(1) kh0(t−kT)sinΘ1+ 〓k a(2) k(hc(t−kT)cosΘ1 hs(t−kT)sinΘ1)+ 〓 〓k b(2) k(hc(t−T/2−kT)sinΘ1−hs(t−T/2
−kT)cosΘ1……(16) v(1) I(t)= 〓k a(1) kh0(t−T/2−kT)sinΘ1− 〓k b(1) kh0(t−kT)cosΘ1+ 〓k a(2) k(hc(t−kT)sinΘ1 −hs(t−kT)cosΘ1)− 〓 〓k b(2) k(hc(t−T/2−kT)cosΘ1−hs(t−T/2
−kT)sinΘ1)……(17) を得る。但し、h0(t)、hs(t)は以下に定義さ
れる関数である。
Now, in FIG. 4, on the receiving side, the inverse calculation of the calculation performed on the transmitting side is performed. For example, multipliers 84, 85 and shaping filters 92, 93
is the in-phase baseband signal v (1) R from the first QAM signal
(t), is used to restore the orthogonal baseband signal v (1) I (t). Here, the molded filters 92, 9
The frequency transfer characteristic of No. 3 is based on the shaping filter 6 on the reverse signal side.
It has exactly the same root Nyquist characteristics as 3, 64, etc. Therefore, v (1) R (t) and v (1) I (t) are each given by the following equations when there is no noise in the transmission path. (However, ω 1 is assumed to be larger than 2π/T.) v (1) R (t) = (u (1) R (τ) cosΘ 1 −u (1) I (t) sin
Θ 1 ) * g (t) + (u (2) R (t) cos (Θ 1 −2π/Tt
) −u (2) I (t) sin(Θ 1 −2π/Tt) * g(t
)...(14) v (1) I (t) = (u (1) R (t) sinΘ 1 +u (1) I (t) cos
Θ 1 ) * g (t) + (u (2) R (τ) sin (Θ 1 −2π/Tt
) +u (2) I (t) cos (Θ 1 −2π/Tt)) * g (
t)...(15) Here, Θ 1 is the phase shift of the first carrier on the receiving side, and * represents the convolution operation.
Substituting equations (1), (2), (7), and (8) into equations (14) and (15), v (1) R (t)= 〓 k a (1) k h 0 (t− T/2−kT) cosΘ 1 − 〓 k b (1) k h 0 (t−kT) sinΘ 1 + 〓 k a (2) k (h c (t−kT) cosΘ 1 h s (t−kT) sinΘ 1 )+ 〓 〓 k b (2) k (hc (t-T/2-kT) sinΘ 1 -h s (t-T/2
−kT) cosΘ 1 …(16) v (1) I (t)= 〓 k a (1) k h 0 (t−T/2−kT) sinΘ 1 − 〓 k b (1) k h 0 ( t-kT) cosΘ 1 + 〓 k a (2) k (h c (t-kT) sinΘ 1 -h s (t-kT) cosΘ 1 )- 〓 〓 k b (2) k (h c (t- T/2-kT) cosΘ1 - hs (t-T/2
−kT) sinΘ 1 )...(17) is obtained. However, h 0 (t) and h s (t) are functions defined below.

hp(t)=g(t)*g(t) hc(t)=g(t)cos2π/Tt*g(t) hs(t)=g(t)sin2π/Tt*g(t) g(t)はルートナイキスト特性を有するフイ
ルタのインパルス応答であるから、整数値mに対
し、 hp(mT)=1 m=0の時 0 m≠0の時 であり、しかも、任意の整数mに対し、 hs(mT)=0、hc(mT+T/2)=0 が成立することは、既に知られている。従つて、
(16)、(17)式において、Θ1=0とし、v(1) R(t)

v(2) R(t)を各々t=mT+T/2、t=mTのタイミ ング時刻にて標本化すれば、 v(1) R(mT+T/2)=a(1) n v(1) I(mT)=b(1) n となり、第1チヤネルの同相データ、直交データ
を正しく復元できることがわかる。第2乃至第4
チヤネルの同相データ、直交データも全く同様に
して正しく復元することができる。しかし乍ら、
一般に受信側では、キヤリア位相ずれ、タイミン
グずれが存在するため、このずれを制御する事が
必要となる。本発明はこうした直交QAM伝送系
において、各チヤネルがスタガQAMとなつてい
る点および全チヤネルが同期クロツクにて動作し
ている点に着目し、上記のキヤリア位相ずれ、タ
イミングずれを精度良く制御し得るキヤリア・タ
イミング結合制御方式を提供する。
h p (t)=g(t)*g(t) h c (t)=g(t) cos2π/Tt*g(t) h s (t)=g(t) sin2π/Tt*g(t ) Since g(t) is the impulse response of a filter having root Nyquist characteristics, for an integer value m, h p (mT) = 1 when m = 0, 0 when m≠0, and for any arbitrary It is already known that h s (mT)=0 and h c (mT+T/2)=0 hold for an integer m. Therefore,
In equations (16) and (17), Θ 1 =0 and v (1) R (t)
,
If v (2) R (t) is sampled at the timings of t=mT+T/2 and t=mT, then v (1) R (mT+T/2)=a (1) n v (1) I (mT)=b (1) n , and it can be seen that the in-phase data and orthogonal data of the first channel can be correctly restored. 2nd to 4th
Channel in-phase data and quadrature data can also be correctly restored in exactly the same manner. However,
Generally, on the receiving side, there is a carrier phase shift and a timing shift, so it is necessary to control these shifts. In this orthogonal QAM transmission system, the present invention focuses on the fact that each channel is a staggered QAM and that all channels operate with a synchronous clock, and precisely controls the carrier phase shift and timing shift described above. A carrier timing coupling control method is provided.

いま、第4図において、受信側の真のタイミン
グずれをτ*とし、第1乃至第4QAMチヤネルに
対する真のキヤリア位相ずれを各々Θ* 1、Θ* 2
Θ* 3、Θ* 4とする。第kチヤネルの送信複素エンベ
ローブ、受信複素エンベローブを各々βk(t)、γk
(t)とすれば、 γk(t)=βk(t+τ*)ej*k+βk-1(t+τ*
)ej0(t+*)+jk*+βk+1(t+τ*)ej0(t+*)
+j
*k +nk(t) k=1、2、3、4 ……(19) と表わされる。但し、β0(t)=β5(t)=0とす
る。
Now, in FIG. 4, the true timing shift on the receiving side is τ * , and the true carrier phase shifts for the first to fourth QAM channels are respectively Θ * 1 , Θ * 2 ,
Let Θ * 3 and Θ * 4 . The transmission complex envelope and reception complex envelope of the k-th channel are β k (t) and γ k , respectively.
(t), then γ k (t)=β k (t+τ * )e j*kk-1 (t+τ *
)e j0(t+*)+jk*k+1 (t+τ * )e j0(t+*)
+j
*k +n k (t) k=1, 2, 3, 4 (19). However, β 0 (t)=β 5 (t)=0.

また、ω0=2π/Tであり、nk(t)は第kチヤ
ネルの複素ベースバンド雑音である。(19)式の
第2項、第3項は第1項の所望信号とは独立なチ
ヤネル間干渉であるから、これらを近似的にガウ
ス雑音と考え、γk(t)を次のように書き表わす。
Also, ω 0 =2π/T, and n k (t) is the complex baseband noise of the k-th channel. Since the second and third terms in equation (19) are inter-channel interference independent of the desired signal in the first term, they can be approximately considered Gaussian noise, and γ k (t) can be expressed as follows: Write down.

γk(t)=βk(t+τ*)ej*k+Nk(t)……(20
) 但し、Nk(t)は、(19)式中、第1項以外を
全て加算したものである。この系における尤度関
数をL(τ、Θ1、Θ2、Θ3、Θ4)とすれば、L
(τ、Θ1、Θ2、Θ3、Θ4) ={送信複素エンベローブ集合β1(t+τ)ej1
β2(t+τ)ej3、β4(t+τ)ej4を送つたとし
た時に受信複素エンベローブ集合γ1(t)、γ2
(t)、γ3(t)、γ4(t)を受ける事後確率} と定義される。ここで雑音N1(t)、N2(t)、N3
(t)、N4(t)は、互いに無相関なガウス雑音と
みなせるから、 L(τ、Θ1、Θ2、Θ3、Θ4) =4k=1 {βk(t+τ)ejk を送つたとしたときにγk(t)を受ける事後確率}
となり、その対数をとつた対数尤度関数をΩ
(τ、Θ1、Θ2、Θ3、Θ4)は、Ω(τ、Θ1、Θ2
Θ3、Θ4)=Σ log{βk(t+τ)ejkを送つた時

γk(t)を受ける事後確率}と表わされる。
γ k (t)=β k (t+τ * )e j*k +N k (t)……(20
) However, N k (t) is the sum of all terms other than the first term in equation (19). If the likelihood function in this system is L (τ, Θ 1 , Θ 2 , Θ 3 , Θ 4 ), then L
(τ, Θ 1 , Θ 2 , Θ 3 , Θ 4 ) = {transmission complex envelope set β 1 (t+τ)e j1 ,
When β 2 (t+τ)e j3 , β 4 (t+τ)e j4 are sent, the received complex envelope set γ 1 (t), γ 2
(t), γ 3 (t), γ 4 (t)}. Here, the noise N 1 (t), N 2 (t), N 3
(t) and N 4 (t) can be regarded as mutually uncorrelated Gaussian noise, so L(τ, Θ 1 , Θ 2 , Θ 3 , Θ 4 ) = 4k=1k (t+τ)e j 〓 Posterior probability of receiving γ k (t) when sending k }
The log likelihood function obtained by taking the logarithm is Ω
(τ, Θ 1 , Θ 2 , Θ 3 , Θ 4 ) is Ω(τ, Θ 1 , Θ 2 ,
Θ 3 , Θ 4 )=Σ log {β k (t+τ) e j 〓 posterior probability of receiving γ k (t) when sending k }.

Ω(τ、Θ1、Θ2、Θ3、Θ4)をデータに関し平
均化したものをとすれば、 =4k=1 k (τ、Θk) ……(21) 但し、kはkか奇数の時、 (τ、Θk)=定数×Σv2 1k(mT+T/2+τ)Σv2 2k(nT+τ) ……(22) でありkが偶数の時、 (τ、Θk)=定数× 〓m v2 1k(mT+τ)〓ov2 2k(nT+T/2+τ)……(23) であり、ここでv1k(t)、v2k(t)は、
各々 v1k(t)∫Re{γk(s)e-jk}g(s−t)ds v2k(t)∫lm{γk(s)e-jk}g(s−t)ds であり、これに相当する信号は、第4図の出力端
101〜108にて得られる。即ち、出力端10
1,102には、第1チヤネル復調ベースバンド
信号の同相分v1,1(t)および直交分v2,1(t)が得
られ、出力端103,104にはv1,2(t)、v2,2
(t)が、出力端105,106には、v1,3(t)、
v2,3が、出力端107,108にはv1,4(t)、v2,4
(t)が各々得られる。
If Ω (τ, Θ 1 , Θ 2 , Θ 3 , Θ 4 ) is averaged over the data, then = 4k=1 k (τ, Θ k ) ...(21) However, k is When k is an odd number, (τ, Θ k ) = constant × Σv 2 1 , k (mT + T / 2 + τ) Σv 2 2 , k (nT + τ) ... (22) When k is an even number, (τ, Θ k ) = constant × 〓 m v 2 1 , k (mT+τ)〓 o v 2 2 , k (nT+T/2+τ)...(23) where v 1 , k (t) , v 2 , k (t) is
respectively v 1 , k (t)∫Re{γ k (s)e -jk }g(s-t)ds v 2 , k (t)∫lm{γ k (s)e -jk } g(s-t)ds, and signals corresponding to this are obtained at output terminals 101-108 in FIG. That is, the output end 10
1,102, the in-phase component v 1,1 (t) and quadrature component v 2,1 (t) of the first channel demodulated baseband signal are obtained, and the output terminals 103, 104 have v 1,2 (t ), v 2,2
(t), but the output terminals 105 and 106 have v 1,3 (t),
v 2,3 , output terminals 107, 108 have v 1,4 (t), v 2,4
(t) are obtained respectively.

上記(21)式よりタイミングずれτ、キヤリア
位相ずれΘ1、Θ2、Θ3、Θ4に対する以下のような
制御方式が得られる。
From the above equation (21), the following control method for the timing deviation τ and the carrier phase deviations Θ 1 , Θ 2 , Θ 3 , and Θ 4 can be obtained.

即ち、現在の推定値τ、Θ(i) 1、Θ(i) 2、Θ(i) 3、Θ
(i) 4
ら次の推定値τ(i+1)、Θ1 (i+1)、Θ2 (i+1)、Θ3 (i+1)

Θ4 (i+1)を得るには、K〓、K〓をループゲインとし
て、 とすればよい。直交QAM伝送系では各チヤネル
がスタガQAMであるため、既に述べた如く、タ
イミング、キヤリア位相に関しNDA制御が可能
となる。更に、タイミングに関しては、全チヤネ
ルが同期関係にあるため、各チヤネルから得られ
たタイミング制御信号を全て加算したもので、共
通制御をかければよい事が判る。
That is, the current estimates τ, Θ (i) 1 , Θ (i) 2 , Θ (i) 3 , Θ
(i) Estimates from 4 to τ (i+1) , Θ 1 (i+1) , Θ 2 (i+1) , Θ 3 (i+1)
,
To obtain Θ 4 (i+1) , let K〓, K〓 be the loop gains, And it is sufficient. In the orthogonal QAM transmission system, each channel is staggered QAM, so as mentioned above, NDA control is possible with respect to timing and carrier phase. Furthermore, regarding timing, since all channels are in a synchronous relationship, it is understood that common control can be applied by adding up all the timing control signals obtained from each channel.

第6図は、以上の原理に基いた本発明の一実施
例を示すブロツク図である。即ち、第6図は第4
図の受信部処理を全てデジタル化したものに、本
発明によるキヤリア・タイミング結合制御方式を
適用した例であつて、121,122,123は
各々多重分離回路、ポリフエーズ回路、逆オフセ
ツトフーリエ変換回路を表わし、逆オフセツトフ
ーリエ変換回路123の4つの出力は、第4図の
出力端101乃至108に得られる第1チヤネル
複素ベースバンド信号乃至第4チヤネル複素ベー
スハンド信号のT/2秒標本化系列となつてい
る。なお、多重分離回路、ポリフエーズ回路、逆
オフセツトフーリエ変換回路の詳細な内容につい
ては、文献5として昭和55年特許願28740号を参
照されたい。
FIG. 6 is a block diagram showing an embodiment of the present invention based on the above principle. That is, Fig. 6 is the 4th
The figure shows an example in which the carrier timing coupling control method according to the present invention is applied to a system in which all of the receiving section processing has been digitized, and 121, 122, and 123 are demultiplexing circuits, polyphase circuits, and inverse offset Fourier transform circuits, respectively. The four outputs of the inverse offset Fourier transform circuit 123 are T/2 seconds sampling of the first channel complex baseband signal to the fourth channel complex base hand signal obtained at the output terminals 101 to 108 in FIG. It has become a series. For details of the demultiplexing circuit, the polyphase circuit, and the inverse offset Fourier transform circuit, please refer to Patent Application No. 28740 filed in 1987 as Document 5.

第6図において、標本化回路138は入力端1
20に得られた受信信号を周波数fsなる標本化ク
ロツクにて標本化している。ここでfsと各チヤネ
ルのボー速度1/Tとは、通常、簡単な整数比関
係を保つており、従つて、いわゆるタイミング時
刻は、この標本化クロツクを基準にして定められ
る。従つて、この系ではタイミング制御信号は、
標本化クロツクを発生するVCO137を制御す
ることになる。124乃至127は、第1チヤネ
ル乃至第4チヤネルに対する位相回転回路であ
り、128乃至131は各チヤネルに対する位相
差検出回路である。124乃至127の位相回転
回路は、例えば第9図に示す如く構成される。即
ち、第9図において160および161は各々復
調ベースバンド信号の同相分および直交分の入力
される入力端であり、162,163,164,
165は乗算器であり、166は加算器であり、
167は減算器であり、168は正弦波系列、余
弦波系列の書き込まれたROMであり、169は
1サンプル遅延回路であり、170は加算器であ
り172,173は位相回転を施されたベースバ
ンド信号の同相分および直交分の出力される出力
端である。いま、制御端子171にφなる位相制
御信号が入力されたとすると、ROM168はφ
に相当する2進符号にてアドレスのアクセスがな
される。ROM168にはφの入力に対し、cosφ
とsinφとを出力するよう内部データが書込まれて
いる。ROM168の第1の出力として得られた
cosφは乗算器162,165にて各々同相信号、
直交信号に乗ぜられる。同様にして第2の出力で
あるsinφは乗算器163,164にて各々同相信
号、直交信号に乗ぜられる。更に、乗算器162
の出力と乗算器163の出力とは、加算器166
にてお互いの和がとられ、この和信号を新たな同
相信号として出力端172に出力する。乗算器1
64の出力と乗算器165の出力とは減算器16
7にて互いの差がとられ、この差信号を新たな直
交信号として出力端173に出力する。この操作
により入力複素信号は(−φ)の位相回転を受け
て出力されることになる。
In FIG. 6, the sampling circuit 138 is connected to the input terminal 1.
The received signal obtained at 20 is sampled by a sampling clock having a frequency fs . Here, f s and the baud rate 1/T of each channel usually maintain a simple integer ratio relationship, and therefore, the so-called timing time is determined based on this sampling clock. Therefore, in this system, the timing control signal is
It will control the VCO 137 which generates the sampling clock. 124 to 127 are phase rotation circuits for the first to fourth channels, and 128 to 131 are phase difference detection circuits for each channel. The phase rotation circuits 124 to 127 are configured as shown in FIG. 9, for example. That is, in FIG. 9, 160 and 161 are input terminals into which the in-phase and quadrature components of the demodulated baseband signal are input, respectively, and 162, 163, 164,
165 is a multiplier, 166 is an adder,
167 is a subtracter, 168 is a ROM in which sine wave series and cosine wave series are written, 169 is a 1 sample delay circuit, 170 is an adder, and 172 and 173 are bases subjected to phase rotation. This is the output terminal where the in-phase and quadrature components of the band signal are output. Now, if a phase control signal φ is input to the control terminal 171, the ROM 168
Addresses are accessed using binary codes corresponding to . In response to the input of φ, the ROM168 has cosφ
Internal data is written to output sinφ. Obtained as the first output of ROM168
cosφ is an in-phase signal at multipliers 162 and 165, respectively;
Multiplied by orthogonal signal. Similarly, the second output sinφ is multiplied by the in-phase signal and the quadrature signal by multipliers 163 and 164, respectively. Furthermore, the multiplier 162
The output of the adder 166 and the output of the multiplier 163 are
The sum signal is summed at the output terminal 172, and this sum signal is outputted to the output terminal 172 as a new in-phase signal. Multiplier 1
64 and the output of the multiplier 165 are the subtracter 16
7, the difference between them is taken, and this difference signal is outputted to the output terminal 173 as a new orthogonal signal. By this operation, the input complex signal undergoes a phase rotation of (-φ) and is output.

位相制御信号は第6図の位相差検出回路12
8,129,130,131の出力として得られ
るが、これら位相差検出回路は、例えば第7図の
ように構成される。即ち、第7図において、14
0,141は複素ベースバンド信号の同相分、直
交分の入力される入力端であり、142,143
はT/2秒遅延回路であり、144,145は乗
算器であり、146は減算器であり、147は位
相制御信号の出力される出力端である。入力端1
40,141のいずれに同相分を入力させるか
は、送信側でのT/2秒遅延の付与の仕方に依存
する。即ち、第1、第3チヤネルの如く、送信側
にて同相分がT/2秒遅れているチヤネルに対し
ては、入力端140に同相分、入力端141に直
交分という割振りになる。第2、第4チヤネルの
如く、送信側にて直交分がT/2秒遅れているチ
ヤネルに対しては、その逆の割振りを行なう。こ
のようにして出力端147には、二つの入力の相
関値と二つのT/2秒遅延入力の相関値との差成
分が位相制御信号として出力される。
The phase control signal is output from the phase difference detection circuit 12 in FIG.
These phase difference detection circuits are configured as shown in FIG. 7, for example. That is, in FIG. 7, 14
0,141 are input terminals into which the in-phase and quadrature components of the complex baseband signal are input, and 142,143
is a T/2 second delay circuit, 144 and 145 are multipliers, 146 is a subtracter, and 147 is an output terminal from which a phase control signal is output. Input end 1
40 or 141 to which the in-phase component is input depends on how the T/2 second delay is applied on the transmitting side. That is, for channels such as the first and third channels in which the in-phase component is delayed by T/2 seconds on the transmitting side, the in-phase component is allocated to the input end 140 and the orthogonal component is allocated to the input end 141. For channels such as the second and fourth channels where the orthogonal portion is delayed by T/2 seconds on the transmitting side, the opposite allocation is performed. In this way, the difference component between the correlation values of the two inputs and the correlation values of the two T/2 seconds delayed inputs is outputted to the output terminal 147 as a phase control signal.

以上の如くキヤリア位相ずれΘ1、Θ2、Θ3、Θ4
は各々位相回転回路124,125,126,1
27にて補正される。この制御ループは、その内
部に1サンプル程度の遅延しか含まないために、
キヤリア周波数にオフセツトがあつたとしても、
それに対するトラツキング特性は良好である。
As described above, the carrier phase deviations Θ 1 , Θ 2 , Θ 3 , Θ 4
are phase rotation circuits 124, 125, 126, 1, respectively.
Corrected in step 27. This control loop only includes a delay of about 1 sample, so
Even if there is an offset in the carrier frequency,
The tracking characteristics for this are good.

次にタイミングずれτは、第6図の132,1
33,134,135にて示されるタイミングず
れ検出回路で、各々独立に検出され、これら検出
出力は全て加算回路136にて加算され、VCO
137を制御することになる。各タイミングずれ
検出回路は、例えば、第8図の如く構成される。
即ち、第8図において、150,151は位相補
正された複素ベースバンド信号の同相分、直交分
の入力される入力端であり、152,153は、
離散値系として実現される微分回路であり、15
4,155は乗算器であり、156は加算器であ
り、158はT/2秒遅延回路であり157は出
力端である。第6図における132,134の如
く送信側にて同相信号が直交信号に対しT/2秒
遅延している奇数番目のチヤネルに対するタイミ
ングずれ検出回路においては、入力端150に同
相信号、入力端151に直交信号が各々入力され
る。
Next, the timing deviation τ is 132,1 in FIG.
Timing deviation detection circuits 33, 134, and 135 detect each independently, and all of these detection outputs are added together in an adder circuit 136, and the VCO
137 will be controlled. Each timing deviation detection circuit is configured as shown in FIG. 8, for example.
That is, in FIG. 8, 150 and 151 are input terminals into which the in-phase and quadrature components of the phase-corrected complex baseband signal are input, and 152 and 153 are as follows.
It is a differential circuit realized as a discrete value system, and 15
4, 155 is a multiplier, 156 is an adder, 158 is a T/2 second delay circuit, and 157 is an output terminal. In the timing shift detection circuit for odd-numbered channels such as 132 and 134 in FIG. Orthogonal signals are input to the ends 151, respectively.

逆に、第6図における133,135の如く送
信側にて直交信号が同相信号に対しT/2秒遅延
している偶数番目チヤネルに対するタイミングず
れ検出回路においては、入力端150に直交信
号、入力端151に同相信号が各々入力される。
Conversely, in a timing shift detection circuit for even-numbered channels such as 133 and 135 in FIG. In-phase signals are input to input terminals 151, respectively.

第8図において、乗算器154は入力端150
に入力された信号とその微分信号との相関をと
り、乗算器153は入力端151に入力された信
号をT/2秒遅延して得られる信号とその微分信
号との相関をとり、これ等の相関値は加算器15
6にて加算された後、端子157にタイミングず
れ信号として出力される。このタイミングずれ信
号は各チヤネル間で互いに独立に混入して来る雑
音分を除いて全てのチヤネルで共通であるから、
第6図の加算回路の出力には、信号対雑音比の高
いタイミング制御信号が得られる。
In FIG. 8, multiplier 154 is connected to input terminal 150.
The multiplier 153 calculates the correlation between the signal input to the input terminal 151 and its differential signal and the signal obtained by delaying the signal input to the input terminal 151 by T/2 seconds. The correlation value of
After being added in step 6, the signals are outputted to a terminal 157 as a timing shift signal. This timing shift signal is common to all channels except for noise that mixes independently between each channel, so
A timing control signal with a high signal-to-noise ratio is obtained at the output of the adder circuit of FIG.

第6図に示した本発明によるキヤリア・タイミ
ング結合制御方式の一実施例は、伝送帯域内での
群遅延歪振幅がT秒に比し充分小さい場合には、
そのまま適用できる。しかし乍ら、通常の搬送帯
域では、特に帯域端近傍にて著しく大きな群遅延
歪を生ずる事がある。このような回線に対しては
タイミング制御信号を全チヤネル共通化する第6
図の実施例は使用できない。しかし、こうした回
線に対しても、帯域の中央近傍にスペクトラム配
置された数チヤネルのみでタイミング制御信号の
共通化を図る事により、キヤリア・タイミングに
対するトラツキングを可能にする事ができる。
In one embodiment of the carrier timing coupling control method according to the present invention shown in FIG. 6, when the group delay distortion amplitude within the transmission band is sufficiently small compared to T seconds,
It can be applied as is. However, in a normal carrier band, extremely large group delay distortion may occur, especially near the band edges. For such lines, there is a sixth method that makes the timing control signal common to all channels.
The illustrated embodiment cannot be used. However, even for such lines, tracking of carrier timing can be made possible by sharing timing control signals among only a few channels whose spectrum is arranged near the center of the band.

第10図は、上記の観点より導かれた本発明に
よるキヤリア・タイミング結合制御方式の一般的
な一実施例を示すブロツク図であつて、181,
182,183,184は、第6図の標本化回路
138、多重分離回路121、ポリフエーズ回路
122、逆オフセツトフーリエ変換回路123に
各々対応している。また、185乃至188は第
6図の124乃至127に対応する第1チヤネル
乃至第4チヤネル用位相回転回路である。更に、
189,190は第6図の129,130に対応
した位相差検出回路であり、191,192は第
6図の133,134に対応したタイミングずれ
検出回路である。第10図に示した本発明の一般
的な一実施例では、第2チヤネル、第3チヤネル
より得られたタイミングずれ信号のみを、加算回
路193にて共通化してVCO194に対する制
御信号を得ている。これによつて第2チヤネルと
第3チヤネルとは各々タイミングずれ、キヤリア
ずれ共に制御されてしまう。これに対し第1チヤ
ネルと第4チヤネルについては、位相差検出回
路、タイミングずれ検出回路を備えていない、し
かし乍ら、少なく共タイミングずれの時間変動量
は、回線歪に関わり無く全チヤネルに共通であ
り、この共通なずれ分は、VCO194、標本化
回路181等で成る制御ループにより補償され
る。また、第1チヤネルおよび第4チヤネルには
位相回転回路185,188が設けられており、
各々の位相回転量は、位相差検出回路189にて
制御されているため、キヤリア位相ずれの時間変
動分についても補償される。即ち、加算回路19
3にて共通化されない第1および第4チヤネルに
ついても、その出力複素信号にはタイミング・キ
ヤリア位相の時間的な変動分は存在せず静的な位
相ずれのみが残留することになる。こうした静的
なタイミングずれ、キヤリア位相ずれは、文献3
に示された2次元自動等化器にて吸収し得ること
が知られている。
FIG. 10 is a block diagram showing a general embodiment of the carrier timing coupling control method according to the present invention derived from the above viewpoint, and includes 181,
182, 183, and 184 correspond to the sampling circuit 138, demultiplexing circuit 121, polyphase circuit 122, and inverse offset Fourier transform circuit 123 in FIG. 6, respectively. Further, 185 to 188 are phase rotation circuits for the first to fourth channels corresponding to 124 to 127 in FIG. Furthermore,
189 and 190 are phase difference detection circuits corresponding to 129 and 130 in FIG. 6, and 191 and 192 are timing shift detection circuits corresponding to 133 and 134 in FIG. In a general embodiment of the present invention shown in FIG. 10, only the timing deviation signals obtained from the second channel and the third channel are shared in an adder circuit 193 to obtain a control signal for the VCO 194. . As a result, both the timing shift and the carrier shift of the second channel and the third channel are controlled. On the other hand, the first channel and the fourth channel are not equipped with a phase difference detection circuit or a timing shift detection circuit, however, the amount of time variation of the timing shift is common to all channels regardless of line distortion. This common deviation is compensated by a control loop made up of the VCO 194, sampling circuit 181, and the like. Further, phase rotation circuits 185 and 188 are provided in the first channel and the fourth channel,
Since each phase rotation amount is controlled by the phase difference detection circuit 189, the time variation of the carrier phase shift is also compensated for. That is, the addition circuit 19
Regarding the first and fourth channels that are not shared in the third channel, there is no temporal variation in timing/carrier phase in the output complex signal, and only a static phase shift remains. These static timing deviations and carrier phase deviations are
It is known that this can be absorbed by the two-dimensional automatic equalizer shown in .

以上述べた如く、本発明によればトラツキング
特性の優れたキヤリア・タイミング結合制御方式
を得る事ができる。
As described above, according to the present invention, a carrier timing coupling control system with excellent tracking characteristics can be obtained.

尚、第6図の実施例は第10図の一般的な実施
例の特例と見做し得る事は明らかである。
It is clear that the embodiment shown in FIG. 6 can be regarded as a special case of the general embodiment shown in FIG.

【図面の簡単な説明】[Brief explanation of drawings]

第1図は二次元変調信号の一般的な生成過程を
示すブロツク図であり、3,4は各々同相側ベー
スバンド整形フイルタ、直交側ベースバンド整形
フイルタである。第2図は同相一直交間にT/2
秒の遅延差を有するスタガQAM信号に対するキ
ヤリア・タイミング結合制御系の構成を示すブロ
ツク図であり、12,14,19,21,22は
乗算器であり、13は90゜位相シフタ、15,1
8は整合フイルタ、16はT/2秒遅延回路、1
7,25はVCO、20,23は微分回路、24
は加算器である。第3図は第2図のタイミング・
キヤリア結合制御系を全てデジタル処理にて実現
するための構成を示したブロツク図であり、32
はサンプラー、35はROM、38は位相回転回
路である。 第4図は直交QAM伝送系の原理的な構成図で
あり、第5図はこの伝送系にて伝送路に送出され
る直交QAM信号のスペクトラム配置例を示した
図である。第6図は本発明によるキヤリア・タイ
ミング結合制御方式の具体的な一実施例を示すブ
ロツク図であり、121,122,123は、
各々多重分離回路、ポリフエーズ回路、逆オフセ
ツトフーリエ変換回路、124乃至127は位相
回転回路、128乃至131は位相差検出回路、
132乃至135はタイミングずれ検出回路、1
36は加算回路であり、137はVCO、138
は標本化回路である。第7図は位相差検出回路の
具体的な構成例を示したブロツク図であり、14
2,143はT/2秒遅延回路、144,145
は乗算器、146は減算器である。第8図はタイ
ミングずれ検出回路の具体的な構成例を示したブ
ロツク図であり、152,153は微分回路、1
54,155は乗算器、156は加算器、158
はT/2秒遅延回路である。第9図は位相回転回
路の具体的な構成例を示したブロツク図であつ
て、162乃至165は乗算器、166は加算
器、167は減算器、168はROM、169は
1サンプル遅延回路、170は加算器である。 第10図は本発明になるキヤリア・タイミング
結合制御方式の一般的な一実施例を示すブロツク
図であつて、185乃至188は位相回転回路、
189,190は位差検出回路、191,192
はタイミングずれ検出回路、193は加算回路、
194はVCOである。
FIG. 1 is a block diagram showing a general process of generating a two-dimensional modulated signal, and 3 and 4 are an in-phase baseband shaping filter and an orthogonal baseband shaping filter, respectively. Figure 2 shows T/2 between in-phase and orthogonal
FIG. 2 is a block diagram showing the configuration of a carrier timing coupling control system for staggered QAM signals having a delay difference of seconds, in which 12, 14, 19, 21, and 22 are multipliers, 13 is a 90° phase shifter, 15, 1
8 is a matching filter, 16 is a T/2 second delay circuit, 1
7, 25 are VCOs, 20, 23 are differential circuits, 24
is an adder. Figure 3 shows the timing and timing of Figure 2.
32 is a block diagram showing a configuration for realizing a carrier coupling control system entirely through digital processing;
35 is a ROM, and 38 is a phase rotation circuit. FIG. 4 is a diagram showing the basic configuration of an orthogonal QAM transmission system, and FIG. 5 is a diagram showing an example of spectrum arrangement of orthogonal QAM signals sent to a transmission line in this transmission system. FIG. 6 is a block diagram showing a specific embodiment of the carrier timing coupling control system according to the present invention, and 121, 122, 123 are
Each includes a demultiplexing circuit, a polyphase circuit, an inverse offset Fourier transform circuit, 124 to 127 are phase rotation circuits, 128 to 131 are phase difference detection circuits,
132 to 135 are timing deviation detection circuits;
36 is an adder circuit, 137 is a VCO, 138
is the sampling circuit. FIG. 7 is a block diagram showing a specific example of the configuration of the phase difference detection circuit.
2,143 is T/2 second delay circuit, 144,145
is a multiplier, and 146 is a subtracter. FIG. 8 is a block diagram showing a specific example of the configuration of the timing deviation detection circuit, in which 152 and 153 are differentiating circuits;
54, 155 are multipliers, 156 are adders, 158
is a T/2 second delay circuit. FIG. 9 is a block diagram showing a specific configuration example of the phase rotation circuit, in which 162 to 165 are multipliers, 166 is an adder, 167 is a subtracter, 168 is a ROM, 169 is a one-sample delay circuit, 170 is an adder. FIG. 10 is a block diagram showing a general embodiment of the carrier timing coupling control system according to the present invention, in which 185 to 188 are phase rotation circuits;
189, 190 are phase difference detection circuits, 191, 192
193 is a timing shift detection circuit, 193 is an addition circuit,
194 is a VCO.

Claims (1)

【特許請求の範囲】[Claims] 1 第1乃至第N番目の直交振幅変調信号が、周
波数領域にて直交多重された信号を伝送路を介し
て受信し、これを標本化クロツクにて標本化しア
ナログ−デジタル変換、波操作および離散フー
リエ変換により第1乃至第N番目の複素ベースバ
ンド信号に復調する直交多重信号の受信系におい
て該複数個の複素ベースバンド信号に対し、各々
独立に位相回転を与える位相回転操作を施し、予
め定められた第k1乃至第km番目の複素ベースバ
ンド信号の前記位相回転後の出力およびそれらの
微分信号より得られる。m個のタイミング制御信
号を全て加算して得られた信号にて前記標本化ク
ロツクの位相を制御すると共に前記第k1乃至第
km番目の複素ベースバンド信号に対する位相回
転量は、各複素ベースバンド信号の実数部と虚数
部との相関値によつて各々制御され、第k1乃至第
km番目以外の複素ベースバンド信号に対する位
相回転量は、前記相関値のうちの少なくとも1つ
によつて制御されることを特徴とするキヤリアタ
イミング結合制御方式。
1 The first to N-th orthogonal amplitude modulation signals receive orthogonally multiplexed signals in the frequency domain via a transmission line, are sampled by a sampling clock, and are used for analog-to-digital conversion, wave manipulation, and discrete In the receiving system of the orthogonal multiplexed signal which is demodulated into the first to Nth complex baseband signals by Fourier transform, a phase rotation operation is performed to independently give a phase rotation to the plurality of complex baseband signals, and a predetermined is obtained from the phase-rotated outputs of the k 1 -km th complex baseband signals and their differential signals. The phase of the sampling clock is controlled by the signal obtained by adding all m timing control signals, and
The amount of phase rotation for the kmth complex baseband signal is controlled by the correlation value between the real part and the imaginary part of each complex baseband signal, and
A carrier timing coupling control method, wherein an amount of phase rotation for a complex baseband signal other than the km-th signal is controlled by at least one of the correlation values.
JP58023948A 1983-02-16 1983-02-16 Control system for carrier timing coupling Granted JPS59149439A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP58023948A JPS59149439A (en) 1983-02-16 1983-02-16 Control system for carrier timing coupling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP58023948A JPS59149439A (en) 1983-02-16 1983-02-16 Control system for carrier timing coupling

Publications (2)

Publication Number Publication Date
JPS59149439A JPS59149439A (en) 1984-08-27
JPH0131807B2 true JPH0131807B2 (en) 1989-06-28

Family

ID=12124766

Family Applications (1)

Application Number Title Priority Date Filing Date
JP58023948A Granted JPS59149439A (en) 1983-02-16 1983-02-16 Control system for carrier timing coupling

Country Status (1)

Country Link
JP (1) JPS59149439A (en)

Also Published As

Publication number Publication date
JPS59149439A (en) 1984-08-27

Similar Documents

Publication Publication Date Title
US5787123A (en) Receiver for orthogonal frequency division multiplexed signals
KR100377257B1 (en) Method and apparatus for fine frequency synchronization in multi-carrier demodulation systems
JP4149044B2 (en) Method and circuit apparatus for correcting phase and / or frequency error of digital multi-carrier signal
JP2526931B2 (en) PSK signal demodulator
US6243423B1 (en) Receiver, transmitter-receiver, and communication method
US6720824B2 (en) Demodulation method and apparatus
JP4331405B2 (en) How to form a training sequence
US5867532A (en) Data reception apparatus, data transmission apparatus and method thereof
KR19990043408A (en) Method and apparatus for simple frequency acquisition of orthogonal split-band systems
JP3058870B1 (en) AFC circuit
JPH0131821B2 (en)
JP2000059448A (en) Method and system for detecting cpm-modulated information symbol
KR100341200B1 (en) Quadrature frequency division multiplexing demodulator
EP0614582A1 (en) Digital demodulator for preamble-less burst communications
JPH0131808B2 (en)
JPH0131807B2 (en)
JP3495962B2 (en) Distortion estimation device
KR20070095138A (en) Up link signal receiving apparatus and method of multiuser detection using successive interference cancellation in wireless ofdma system
JP2000201193A (en) Demodulating device and method and communication system including the device
JP3148090B2 (en) OFDM signal synchronous demodulator
US6075826A (en) Method and apparatus for obtaining initial carrier and symbol phase estimates for use in synchronizing transmitting data
JPH0141054B2 (en)
JP2005086552A (en) Mfsk receiving system
JPH08274693A (en) Mlse equalizer and demodulator using the same
JP2004096290A (en) Amplitude and phase correction method