JP3566833B2

JP3566833B2 - ディジタル信号変調方法及び復調方法

Info

Publication number: JP3566833B2
Application number: JP18151797A
Authority: JP
Inventors: 和彦府川; 博鈴木; 純一京田
Original assignee: NTT Docomo Inc
Current assignee: NTT Docomo Inc
Priority date: 1997-07-07
Filing date: 1997-07-07
Publication date: 2004-09-15
Anticipated expiration: 2017-07-07
Also published as: JPH1127330A

Description

【０００１】
【発明の属する技術分野】
この発明は移動通信などの無線通信における変動する伝送路において信頼性の高い信号伝送が可能な変調方法、その復調方法を提供するものである。
【０００２】
【従来の技術】
通信用伝送路においては伝送特性、例えば送信機から受信機までのインパルスレスポンスが時間的に大きく変動することがある。マイクロ波無線伝送、移動通信などはその代表的な例である。このような伝送路では比較的雑音レベルの高い受信状態において、希望信号は時間的に変動する波形歪を伴って受信され、さらに変動する同一チャネル干渉が重畳されて伝送特性が大幅に劣化することがある。そのため、これらの原因による劣化を抑えつつ、信頼性の高い伝送方式を実現する必要がある。
従来の解決法と欠点
このように変動する伝送路において適応アルゴリズムを用いた最適受信機が利用されてきた。
【０００３】
変動する歪みが存在する伝送路では、従来、適応受信機として適応等化器が利用されてきた。この適応等化器はその構造から、線形等化器と非線形等化器に分類される。
まず、線形等化器について述べる。図９は、線形等化器の構成を示している。以下では、変調された受信信号を複素表示して動作を説明する。時点ｉにおける等化器への入力信号ｘ（ｉ）の実数成分は受信信号の同相成分の振幅、虚数成分は受信信号の直交成分の振幅を表わしている。入力信号ｘ（ｉ）が入力端子１１からトランスバーサルフィルタ１２へ入力され、トランスバーサルフィルタ１２のｊ番目のタップ係数ｗｊ（ｉ）を制御することにより歪みが除去されて判定回路１３に入力される。判定された信号は出力端子１４から出力される。判定回路１３の入力端と出力端の両信号は差分演算回路１５に入力されて誤差信号ｅ（ｉ）が計算される。誤差信号は制御部１６へ送られる。制御部１６では、誤差信号ｅ（ｉ）と入力信号ｘ（ｉ）から、タップ係数ｗｊが更新される。線形等化器の動作は、例えば、佐藤洋一，“線形等化理論”，丸善などに詳しく説明されている。
【０００４】
トランスバーサルフィルタ１２のタップ数をＭとして、トランスバーサルフィルタ１２の各タップ係数ｗ１（ｉ）からｗＭ（ｉ）を並べた伝送路伝達係数Ｗ（ｉ）と、入力信号ｘ（ｉ）からｘ（ｉ−Ｍ＋１）を並べた入力信号ベクトルＸ（ｉ）を用いると、伝送路伝達係数Ｗ（ｉ−１）は、入力信号ベクトルＸ（ｉ）と誤差信号ｅ（ｉ）を用いて、Ｗ（ｉ）に更新される。理想的に動作する場合には、最小２乗法の解となっており、入力信号ベクトルＸ（ｉ）から求められる自己相関行列Ｒ（ｉ）と、入力信号ベクトルＸ（ｉ）と判定信号との相互相関ベクトルＶ（ｉ）とを用いて、伝送路伝達係数はＷ（ｉ）＝Ｒ^−１（ｉ）Ｖ（ｉ）のように表わされる。これは、正規方程式：Ｒ（ｉ）Ｗ（ｉ）＝Ｖ（ｉ）の解である。しかしながら、自己相関行列の逆行列Ｒ^−１（ｉ）はＲ（ｉ）が正則でなければ存在しないという制限があった。そのため、自己相関行列のランクがＭより小さくなるような受信波、たとえば低雑音条件において同一符号が続く変調波を受信した場合などでは、逆行列が発散し、その結果として伝送路伝達係数が発散するという欠点があった。
【０００５】
前記正規方程式の解を逐次的に求める適応アルゴリズムも知られている。例えば、カルマンフィルタ，ＲＬＳ，ＬＭＳなどが代表的な例である。適応アルゴリズムについてはＨａｙｋｉｎ，“ＡｄａｐｔｉｖｅＦｉｌｔｅｒＴｈｅｏｒｙ”，Ｐｒｅｎｔｉｃｅ−Ｈａｌｌ，１９９１が詳しい。これらの方法においても自己相関行列が非正則のときには解の発散が起きる。
【０００６】
自己相関行列が非正則のときでも発散しない解が得られる方法として、一般逆行列を用いる方法が知られている。一般逆行列については、Ａ．Ａｌｂｅｒｔ，“ＲｅｇｒｅｓｓｉｏｎａｎｄｔｈｅＭｏｏｒＰｅｎｒｏｓｅＰｓｅｕｄｏｉｎｖｅｒｓｅ ”，ＡｃａｄｅｍｉｃＰｒｅｓｓ，１９７２、が詳しい。一般逆行列を用いると、自己相関行列が非正則のときにも解の発散が抑制されているが、その解Ｗ（ｉ）はタップ伝送路伝達係数のノルム‖Ｗ（ｉ）‖が最小となる最小ノルム解となり、時点ｉが推移すると、各時点ごとに最小ノルム解が異なるのでかならずしもその解が徐々に本来の解に収束するものではなかった。
【０００７】
自己相関行列が非正則のときにも徐々に本来の解へ漸近する方法としては、Ｍｏｏｒ−Ｐｅｎｒｏｓｅ一般逆行列を用いる直交射影法が知られている。この方法は、尾関，梅田，“アフィン部分空間への直交射影を用いたフィルタ・アルゴリズムとその諸特性”，電子情報通信学会論文誌，ｖｏｌ．Ｊ６７−Ａ，Ｎｏ．２，ｐｐ．１２６−１３２，１９８４年２月、に示されている。しかしながら、伝送特性が変動する伝送路では自己相関行列も変化していくので、Ｍｏｏｒ−Ｐｅｎｒｏｓｅ一般逆行列を次々と更新する必要があり、この更新には大量の演算を必要とするため実際にリアルタイム条件で使用することは困難であった。
【０００８】
以上説明したように線形等化器の適応性に関する改良が行われてきたが、非最小位相となる歪み条件では十分な等化効果が得られないという欠点があった。
非最小位相となる歪み条件でも十分な等化効果が得られるように非線形等化器が用いられてきた。非線形等化器の構成例を図１０に示す。この構成は適応最尤系列推定器としても知られている。入力信号ｘ（ｉ）が入力端子１１から差分回路２１へ入力され、トランスバーサルフィルタ１２から出力されるレプリカ信号ｙｍ（ｉ）との差分演算が行われて、誤差信号ｅｍ（ｉ）が得られる。誤差信号ｅｍ（ｉ）の絶対値の２乗が２乗器２２でとられて最尤系列推定部２３に入力され、信号が推定される。最尤系列推定部２３からは、符号系列候補｛ａｍ｝が出力される。符号系列候補｛ａｍ｝に対応して、変調部２４からは変調波候補ｓｍが出力される。変調波候補ｓｍはトランスバーサルフィルタ１２に入力され、レプリカ信号ｙｍ（ｉ）が形成される。最尤系列推定部２３から出力される候補は各状態に対して状態遷移の数だけあるので、同一の入力信号ｘ（ｉ）に対して、状態数と状態遷移数に対応した上述の演算を行う。制御部２５は、誤差信号ｅｍ（ｉ）と変調波候補ｓｍを用いて、トランスバーサルフィルタ１２のタップ係数を更新する。非線形等化器の動作は、例えば、府川，鈴木，“移動無線における逐次最小２乗形適応最尤系列推定（ＲＬＳ−ＭＬＳＥ）の特性”，電子情報通信学会論文誌（Ｂ−ＩＩ），ｖｏｌ．Ｊ７５−Ｂ−ＩＩ，Ｎｏ．８，ｐｐ．５３５−５４６，１９９２年８月、などに詳しく説明されている。
【０００９】
トランスバーサルフィルタ１２のタップ数をＭとして、トランスバーサルフィルタの各タップ係数ｗｍ，１（ｉ）からｗｍ，Ｍ（ｉ）を並べた伝送路伝達係数Ｗｍ（ｉ）と、変調波候補ｓｍ（ｉ）からｓｍ（ｉ−Ｍ＋１）を並べた変調波候補ベクトルＳｍ（ｉ）を用いて、伝送路伝達係数Ｗｍ（ｉ−１）は、変調波候補ベクトルＳｍ（ｉ）と誤差信号ｅｍ（ｉ）を用いて、Ｗｍ（ｉ）に更新される。理想的に動作する場合には、最小２乗法の解となっており、変調波候補ベクトルＳｍ（ｉ）が作る自己相関行列Ｒｍ（ｉ）と、変調波候補ベクトルＳｍ（ｉ）と入力信号ｘ（ｉ）との相互相関ベクトルＶｍ（ｉ）とを用いて、伝送路伝達係数はＷｍ（ｉ）＝Ｒｍ^−１（ｉ）Ｖｍ（ｉ）のように表わされる。しかしながら、自己相関行列の逆行列Ｒｍ^−１（ｉ）はＲｍ（ｉ）が正則でなければならないという制限があった。そのため、自己相関行列のランクがＭより小さくなるような変調波候補ベクトル、たとえば同一符号が続く変調波候補などでは、伝送路伝達係数が発散するという欠点があった。
【００１０】
このような同一符号が続く変調波候補において伝送路伝達係数が発散する問題への対処策として、以下の方法が知られている。
（１）府川，鈴木，“Ｂｌｉｎｄｉｎｔｅｒｆｅｒｅｎｃｅｃａｎｃｅｌｌｉｎｇｅｑｕａｌｉｚｅｒｆｏｒｍｏｂｉｌｅｒａｄｉｏｃｏｍｍｕｎｉｃａｔｉｏｎｓ ”，電子情報通信学会論文誌（Ｂ−ＩＩ），ｖｏｌ．Ｅ７７−Ｂ−ＩＩ，Ｎｏ．５，ｐｐ．５８０−５８８，１９９４年５月では、あらかじめすべての符号系列に対する自己相関行列のアンサンブル平均をとり、その平均値に固定する方法がとられている。
（２）古谷，後川，伊佐，佐藤，“ブラインドビタビ等化方式の一提案”，１９９１年春季電子情報通信学会全国大会では、自己相関行列が正則でないときには、アップデートしない方法がとられている。
（３）ＧｒａｈａｍＣ．Ｇｏｏｄｗｉｎ，ＫａｗａｉＳａｎｇＳｉｎ “ＡｄａｐｔｉｖｅＦｉｌｔｅｒｉｎｇ，Ｐｒｅｄｉｃｔｉｏｎ，ａｎｄＣｏｎｔｒｏｌ ”，ＰＲＥＮＴＩＣＥ−ＨＡＬＬ，１９８４年では、自己相関行列を求める前の入力ベクトルに雑音成分を混入する方法が示されている。
（４）鈴木，府川“一般逆行列によるキャリヤ推定を用いたブラインド形最尤推定法”，電子情報通信学会技法，ＲＣＳ９５−５８，１９９５年７月では、一般逆行列を用いてアップデートする方法が示されている（特願平７−１５１６０）。
【００１１】
以上の方法はいずれも復調側で対処する方法であるが、同一符号がかなり長く続く場合などには動作が不安定なるという欠点があった。
【００１２】
【発明が解決しようとする課題】
この発明は、適応等化器などで最尤系列推定を行う信号処理において、すべてのディジタル変調波系列候補に対して、伝送路伝達係数が正確に抽出できる変調方法とその復調方法を提供することを課題としている。
【００１３】
【課題を解決するための手段】
この変調方法と復調方法の構成は次のようなものである。変調方法では（１）ディジタル変調波の複素振幅を特定のコードによってシンボル毎に変化させる変調を行なう。復調方法では（２）該ディジタル変調波の受信信号から該特定コードのタイミングを抽出するタイミング再生過程、（３）該受信信号に含まれる受信ディジタル変調波成分に該特定コードを複素乗積して得られる受信複素振幅成分を用いて適応的に伝送路伝達係数を推定するキャリア成分推定過程、（４）該受信複素振幅成分を用いて受信ディジタル変調波のレプリカを生成して該ディジタル変調波のシンボル系列を推定する最尤系列推定過程、（５）最尤系列推定過程とキャリア成分推定過程が一体となり、再帰的にディジタル変調波シンボル系列の推定と伝送路伝達係数の推定を行う復調過程を有する。
作用
この発明における基本的な作用は次のようなものである。（１）変調方法ではディジタル変調波の複素振幅を特定のコードによってシンボル毎に変化させる。復調方法においては（２）タイミング再生過程ではディジタル変調波の受信信号から該特定コードのタイミングを抽出する。（３）キャリア成分推定過程では該受信信号に含まれる受信ディジタル変調波成分に該特定コードを複素乗積して得られる受信複素振幅成分を用いて適応的に伝送路伝達係数を推定する。（４）最尤系列推定過程では、該受信複素振幅成分を用いて受信ディジタル変調波のレプリカを生成して該ディジタル変調波のシンボル系列を推定する。（５）復調過程では、最尤系列推定過程で推定されたディジタル変調波シンボル系列を用いて、キャリア成分推定過程によって伝送路伝達係数の更新を行い、その更新された伝送路伝達係数を用いて再帰的に最尤系列推定によってディジタル変調波シンボル系列の推定を行う。
【００１４】
この発明は、従来の変調方法、復調方法と以下の点が異なる。
従来は同一符号が連続しているディジタル変調波候補系列において伝送路伝達係数が発散する問題に対し、復調側で対処を行っているが、この発明では変調側から同一のディジタル変調波が連続して送られないようにディジタル変調波の複素振幅を特定コードによってシンボル毎に変化させることにより、伝送路伝達係数の発散を回避している。
【００１５】
【発明の実施の形態】
実施例１
この発明の基本構成を図１に示す。図１Ａにおいて入力端子３１にディジタル変調波ｂ（ｉ）が変調手段３２に入力されると、ｂ（ｉ）の複素振幅を変化させるために特定コードＣ（ｉ）＝Ａ（ｉ）ｅｘｐ（ｊｄ（ｉ））が複素乗積される。このようにディジタル変調波の複素振幅を変化させて、復調側での伝送路伝達係数の推定に用いられる自己相関行列Ｒｍ（ｉ）を、すべてのディジタル変調波候補系列において正則となるようにする。このようにして得られた送信複素振幅ｓ（ｉ）で矩形パルスを矩形パルス発生器３３で生成し、つまり所定のシンボル長の矩形波として搬送波に乗せて送信機３４で送信する。
【００１６】
受信機３５からは受信した信号をベースバンド周波数に落とした信号ｘ（ｉ）が出力され、復調手段３６に入力される。図１Ｂに復調手段３６の原理図を示す。最尤系列推定手段３７では時刻ｉにおいて各状態に設定されている伝送路伝達係数Ｗｍ（ｉ−１）と入力信号ｘ（ｉ）から、受信複素振幅系列Ｓｍ（ｉ）を推定する。この受信複素振幅系列Ｓｍ（ｉ）の自己相関行列を演算部３８で求め、更にその逆行列を演算部３９で求める。更に入力信号ｘ（ｉ）と受信複素振幅系列Ｓｍ（ｉ）との相互相関行列Ｖｍ（ｉ）を相互相関部４１で求める。キャリア位相推定手段４２では、Ｓｍ（ｉ）を用いて作られる自己相関行列の逆行列Ｒｍ^−１（ｉ）を用いて、伝送路伝達係数を以下のように更新する。
【００１７】
Ｗｍ（ｉ）＝Ｒｍ^−１（ｉ）Ｖｍ（ｉ）
受信複素振幅系列Ｓｍ（ｉ）から作られるＲｍ（ｉ）は、すべてのＳｍ（ｉ）において正則となるように変調手段３２によってディジタル変調波の複素振幅を特定のコードによって変化させているため発散する事はなく、したがって伝送路伝達係数が発散することはない。
【００１８】
各状態ごとの状態遷移に対応する符号系列は通常固定されているので、受信複素振幅系列Ｓｍ（ｉ）も状態遷移ごとに固定されている。すなわちＳｍ（ｉ）から作られる自己相関行列Ｒｍ（ｉ）はＲｍとみなすことができる。したがって、キャリア位相推定手段４２ではＲｍからあらかじめ逆行列Ｒｍ^−１を計算しておき、メモリに蓄積しておくことが可能である。この性質は先に説明した線形等化器における直交射影法とは大きく異なる点である。
実施例２
上記の特定コードＣ（ｉ）として、Ａ（ｉ）＝１とし、さらにφ（ｉ）を２値、例えば、φ_１，φ_２としても、自己相関行列を正則にすることができる。このような簡単な特定コードを用いることにより変調手段３２を簡略化できるメリットがある。
実施例３
ＢＰＳＫについて、２値特定コードを利用する場合について詳しく述べる。図２Ａに示すように伝送されるシンボル系列ｂ（ｉ）に対して、偶数奇数のタイムスロットに分類しφ（２ｋ）＝φ_１＝０，φ（２ｋ＋１）＝φ_２＝θとし、偶数番目のシンボルでは図２Ｂに示すように０、奇数番目のシンボルでは図２Ｃに示すようにθ＝４５°の２値の特定コードによってディジタル変調波の複素振幅を変化させる。このようにして得られる受信複素振幅系列Ｓｍ（ｉ）から自己相関行列の逆行列Ｒｍ^−１（ｉ）を作り、その行列式を計算すると図３に示すようになる。図３のパラメータは送信符号列の状態である。伝送路伝達係数の推定誤差はＲｍ^−１（ｉ）の大きさに比例することが知られているので、Ｒｍ^−１（ｉ）の大きさにほぼ比例する行列式ｄｅｔ（Ｒｍ^−１（ｉ））の値を示している。ゆえにすべての受信複素振幅系列Ｓｍ（ｉ）において、それから作られるＲｍ^−１（ｉ）の大きさが小さくなるようにθを設定する。ＢＰＳＫの場合は図３からθ＝π／４あるいは３π／４とすればよいことが分かる。
【００１９】
このように、ＢＰＳＫでは、２値位相シンボル系列として０とπ／４の位相を持った複素シンボルを交互に並べた系列を用いることにより、逆行列の大きさを小さくすることができる。したがって、適応制御における係数の推定値の精度を上げることができる。このときのビット誤り率（ＢＥＲ）特性を図４に示す。同図は差動符号化と同期検波の両方の特性を示している。どちらもθが４５°と１３５°においてＢＥＲが最小になっている。この場合の特定コードは２値位相シンボル系列であって、一般的には位相をモジュロπ／２で表わしたとき、その２値位相差がπ／４である複素シンボル系列である。
実施例４
ＱＰＳＫについて、２値特定コードを利用する場合について詳しく述べる。例えば図５Ａ，Ｂに示すようにφ（２ｋ）＝φ_１＝０，φ（２ｋ＋１）＝φ_２＝θ＝２２．５°とし、偶数番目のシンボルでは０、奇数番目のシンボルではθの２値の特定コードによってディジタル変調波の複素振幅を変化させる。このようにして得られる受信複素振幅系列Ｓｍ（ｉ）から自己相関行列の逆行列Ｒｍ^−１（ｉ）を作り、その行列式を計算すると図６のようになる。伝送路伝達係数の推定誤差はＲｍ^−１（ｉ）の大きさに比例することが知られているので、すべての受信複素振幅系列Ｓｍ（ｉ）において、それから作られるＲｍ^−１（ｉ）の大きさが小さくなるようにθを設定する。ＱＰＳＫの場合は図６からθ＝π／８とすればよいことが分かる。
【００２０】
このように、ＱＰＳＫでは、２値位相シンボル系列として０とπ／８の位相を持った複素シンボルを交互に並べた系列を用いることにより、逆行列の大きさを小さくすることができる。したがって、適応制御における係数の推定値の精度を上げることができる。このときの誤り率特性を図７に示す。この場合の特定コードも２値位相シンボル系列であって一般的には位相をモジュロπ／４で表わすと、その２値位相の差がπ／８の複素シンボル系列である。
実施例５
上述したようにこの発明では、例えば偶数奇数のタイムスロットのように正しいタイミングが再生されないと、正しい結果が得られない場合がある。そこでタイミング再生手段としては図８に示すようなものを用いる。
【００２１】
まず送信側では図８Ａに示すようにタイミング再生用信号Ｓ_Ｔ（ｉ）と特定コードとの複素乗算を行いタイミング再生用送信信号系列を送出する。Ｓ_Ｔ（ｉ）は通常のデータ系列の中に一定間隔で挿入される。一方、受信側ではタイミング再生用信号Ｓ_Ｔ（ｉ）が既知であり、既知の特定コードとの乗算を乗積４５で行ってタイミング再生用受信信号系列を生成する。次にこの信号系列をトランスバーサルフィルタ１２のタップ係数として設定する。これにより、受信信号系列にタイミング再生用送信信号系列と同じ系列があるとトランスバーサルフィルタ１２の出力にピークが発生する。このピークをピーク検出器４６で検出して再生タイミング信号を得ることができる。実際にはこのピーク位置を位相同期ループ（ＰＬＬ）で平滑化し、誤検出が起きないようにする。このタイミング再生手段が図１Ａ中のタイミング再生手段５１として用いられ、この再生タイミングで復調手段３６が動作される。また復調手段３６では特定コードについて、レプリカを生成し、この特定コードレプリカを最尤系列推定で得られた符号系列候補に乗積して、レプリカ信号ｙｍ（ｉ）とすればよい。
【００２２】
上述した非線形等化器だけでなくブラインド形非線形等化器、あるいは非線形干渉キャンセラ、また単純な同期検波器あるいはダイバーシチ構成の同様な受信機などでもこの発明が適用できる。ＣＤＭＡにおける逆拡散後の信号をＲＡＫＥ合成する受信機においても有効である。また上述では特定コードとして、シンボルごとに２値位相の一方を交互にとる系列としたが、シンボルごとに２値振幅の一方を交互にとるようにしてもよい。つまり、そのような２値振幅シンボル系列をもち、その２値振幅の差が各種の値のものについてこれを横軸とした図３、図６に示す特性を調べ、１／ｄｅｔ〔Ｒｍ（ｉ）〕が常に小さい値をとる２値振幅を選定すればよい。特定コードとして、位相と振幅を共にシンボルごとに２値の一方をとる系列としてもよい。
【００２３】
【発明の効果】
この発明によればディジタル変調波の複素振幅を特定のコードによってシンボル毎に変化させることにより、伝送路伝達係数の推定に用いる自己相関行列をすべてのディジタル変調波候補系列に対して正則にすることができるため、係数の推定を安定に行うことができるとともに、伝送路変動に対する追従性が良い動作が期待される。状態遷移候補ごとに用意する逆行列は、各状態遷移候補で決まっているので計算してメモリに蓄積しておけばよく、計算処理量は逆行列を計算する一般の適応信号処理と比べると少なくなる。また、最尤系列推定を用いるので信号検波性能が優れている。
【００２４】
フェージング変動が速い移動通信、移動衛星通信、また、大容量伝送のために高精度の復調回路を必要とするマイクロ波などの固定無線伝送の分野に適している。これらの分野における適応等化器、適応干渉キャンセラ、ダイバーシチなどの適応受信機を容易に製作できる。
【図面の簡単な説明】
【図１】（ａ）はこの発明の変調方法、復調方法を適用したシステムの基本機能構成を示す図、（ｂ）は図（ａ）中の復調手段３６の変調波候補ｍにおける復調原理の機能的構成を示す図である。
【図２】Ａは複素振幅Ｓ（ｉ）のＢＰＳＫパターンの例を示す図、Ｂはその偶数タイムスロットの複素振幅、Ｃは奇数タイムスロットの複素振幅である。
【図３】ＢＰＳＫの場合の１／ｄｅｔ〔Ｒｍ（ｉ）〕と位相θとの関係例を示す図。
【図４】ＢＰＳＫにおけるビット誤り率と位相θとの関係を示す図。
【図５】ＱＰＳＫの複素振幅値を示す図。
【図６】ＱＰＳＫにおける１／ｄｅｔ〔Ｒｍ（ｉ）〕とθとの関係を示す図。
【図７】ＱＰＳＫにおけるビット誤り率とθとの関係を示す図。
【図８】Ａは既知タイミング再生用信号に特定コードを乗積する送信側の構成図、Ｂは既知タイミング再生信号によるタイミング信号再生手段の構成を示す図である。
【図９】従来の線形等化器の機能構成を示す図。
【図１０】従来の非線形等化器の機能構成を示す図。

Claims

シンボル系列によりディジタル変調を行い、そのディジタル変調波を搬送波に乗せて送信する変調方法において、
前記ディジタル変調波の複素振幅系列Ｓｍ（ｉ）から作られる自己相関行列がすべてのＳｍ（ｉ）において正則となるように、特定コードによりシンボル毎に前記ディジタル変調波の搬送波複素振幅を変化させることを特徴とするディジタル信号変調方法。
請求項１の変調方法において、
前記特定コードは２値位相シンボル系列であることを特徴とする変調方法。
請求項２の変調方法において、
前記ディジタル変調波は２値ＰＳＫ波であり、かつ前記２値位相シンボル系列は位相をモジュロπ／２で表したときその差がπ／４となる２つの位相を持った複素シンボルを交互に並べた系列であることを特徴とする変調方法。
請求項２の変調方法において、
前記ディジタル変調波は４値ＰＳＫ波であり、かつ前記２値位相シンボル系列は位相をモジュロπ／４で表したときその差がπ／８となる２つの位相を持った複素シンボルを交互に並べた系列であることを特徴とする変調方法。
請求項１乃至４の何れかの変調方法において、
送信シンボル系列に既知のタイミング再生用信号を挿入することを特徴とする変調方法。
ディジタル変調波の複素振幅系列Ｓｍ（ｉ）から作られる自己相関行列がすべてのＳｍ（ｉ）において正則となるように、特定コードによりシンボル毎にディジタル変調波の搬送波複素振幅が変化されたディジタル変調波の復調方法であって、
受信されたディジタル変調波から上記特定コードのタイミングを抽出し、
上記ディジタル変調波に上記特定コードを複素乗積して受信複素振幅成分を得、
その受信複素振幅成分を用いて適応的に伝送路伝達係数を推定し、
上記受信複素振幅成分を用いて受信ディジタル変調波のレプリカを生成し、
上記ディジタル変調波の送信シンボル系列を最尤系列推定により推定し、
その推定されたディジタル変調波シンボル系列を用いて上記伝送路伝達係数の更新を行い、その更新された伝送路伝達係数を用いて再帰的に最尤系列推定によってディジタル変調波シンボル系列の推定を行うディジタル信号復調方法。
請求項６の復調方法において、
上記タイミングの抽出を、既知の該タイミング再生用信号と既知の特定コードを用いて相関をとることにより行うことを特徴とするディジタル信号復調方法。