JPH03503347A - 多重通路歪みを補正するシステム - Google Patents

多重通路歪みを補正するシステム

Info

Publication number
JPH03503347A
JPH03503347A JP1502655A JP50265589A JPH03503347A JP H03503347 A JPH03503347 A JP H03503347A JP 1502655 A JP1502655 A JP 1502655A JP 50265589 A JP50265589 A JP 50265589A JP H03503347 A JPH03503347 A JP H03503347A
Authority
JP
Japan
Prior art keywords
signal
umbrella
zero
book
delay
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.)
Granted
Application number
JP1502655A
Other languages
English (en)
Other versions
JP3051984B2 (ja
Inventor
チヤオ ツアイ‐ホング
キヤンベル エドワード リチヤード ザ サード
ルイス ヘンリイ ガートン
Original Assignee
アールシーエー トムソン ライセンシング コーポレーシヨン
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 アールシーエー トムソン ライセンシング コーポレーシヨン filed Critical アールシーエー トムソン ライセンシング コーポレーシヨン
Publication of JPH03503347A publication Critical patent/JPH03503347A/ja
Application granted granted Critical
Publication of JP3051984B2 publication Critical patent/JP3051984B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/14Picture signal circuitry for video frequency region
    • H04N5/21Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03012Arrangements for removing intersymbol interference operating in the time domain
    • H04L25/03019Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
    • H04L25/03038Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception with a non-recursive structure
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/14Picture signal circuitry for video frequency region
    • H04N5/21Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
    • H04N5/211Ghost signal cancellation

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Multimedia (AREA)
  • Power Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Picture Signal Circuits (AREA)
  • Color Television Systems (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるため要約のデータは記録されません。

Description

【発明の詳細な説明】 適応型テレビジョンゴースト消去システム本発明は、ゴーストの無い入力信号に 応答してゴーストのある信号を発生するチャネル・モデルを最初に発生するテレ ビジョンゴースト信号消去システムに関する。
それから、このモデルに基づいて、このシステムはゴーストのある信号からゴー スト信号成分を実質的に取り除くゴースト消去フィルタを形成する。
テレビジョン受信は多重通路歪み、すなわち不所望の多重信号の受信に長いこと 悩まされてきた。建築物および他の大きな物体から反射されるか、または十分に 終端されていないケーブル回路網から生じるこれらの不所望の信号は、直接のテ レビジョン信号が遅延された変形信号きして現われ、通常、再生画像のゴースト 信号と呼ばれる。
このゴースト信号は、直接信号とゴースト信号間の信号経路の長さに関係する関 数として直接信号から遅延している。成る受像機の位置から別の受像機の位置ま でのこの関係の不規則性は、ゴースト搬送波信号の位相が直接信号の位相と任意 の関係を有することを示す。直接信号からゴースト信号を完全に取り除くために は、ゴースト信号の振幅、その遅延および直接のテレビジョン信号の搬送波位相 に対するゴースト信号の搬送波位相を考察することが必要である。
ゴースト信号の相対位相の重要性を理解するためには、テレビジョン信号そのも のについて更に知ることが有用である。NTSC規格の下では、テレビジョン信 号は残留側波帯形式で伝送される。ベースバンド信号の比較的低い周波数成分( 0−1、25M Hz )は両側波帯変調(1)SM)され、一方、高い周波数 成分(1,25−4゜75 M Hz )は単側波帯変調(SSM)される。こ の信号のD S M部分の2つの側波帯域の直角成分は互いに相殺するから、D SMのビデオ信号の直角成分は実質的に零である。しかしながら、この信号のS SM部分の直角成分は零でなく、ゴースト信号として、変調されたビデオ信号の 同相の部分に干渉することがある。
解析すると、変調されたビデオ信号v (t)の同相成分および直角成分は(1 )式で定義される複素ベースパンI・と等価な次式により表される。
v (t)=v 1  (t)+jvq (t)     (1)ここで、jは −1の平方根に相当する複素量であり、Vl (1)およびv、(t)は、信号 v (t)が、例えば、画像搬送波信号に関してそれぞれ同相であり、また直角 位相である振動信号を使用して中間周波(IF)信号から同期復調されると得ら れるベースバンド信号である。
信号v (t)が多重通路の伝送チャネルに供給されると、ゴースト歪み信号が 発生する。
第1図は直接信号およびゴースト信号の相対位相の重要性を示す。例えば、直接 信号が波形10で表される2Tパルスである場合、直接の搬送波信号とゴースト の搬送波信号間の相対位相角度かそれぞれOo、900.180°、または−9 0’  (2700)であるなら、ゴースト信号は波形1.0,1.2.14ま たは16で表される。
さらに、直接信号経路とゴースト信号経路間の関係がランダムであるがら、任意 の中間波形も起こり得る。
直接信号とゴースト信号の相対振幅および相対位相情報は、テレビジョン信号を 同相(I)成分および直角(Q)成分に復調することにより決定することができ る。
■成分はテレビジョン信号の画像搬送波と同相であり、Q成分は画像搬送波に対 して90’位相シフトされている信号と同相である。これらの成分は複素平面に おけるテレビジョン信号を表し、複素平面においてI成分およびQ成分はそれぞ れ実軸および虚軸の座標に対応する。
ビデオ信号の同相成分および直角成分をそれぞれ実成分および虚成分と呼ぶ約束 がこの出願を通じて使われる。
1977年5月に発行されたコンシューマ エレクトロニクスに関するアイ・イ ー・イー・イー・トランザクション(IEEE T「ansactions o n Consumer Elect「onics)の第175頁−第181頁に 掲載されたエイチ、セディック(H,Thedick)氏による「テレビジョン 放送のための適応型多重通路等化」という題目の論文に述べられているように、 ゴースト信号を発生する伝送経路は、直接信号が減衰定数aにより振幅が減少さ れ、また、時間期間τたけ遅延されてゴースト信号を形成するフィードフォワー ド・システムとしてモデル化される。多重通路のゴースト信号を含んでいるビデ オ信号は(2)式で表される。
ここで、a、は複素係数、τ、は1番目の信号経路の相対時間遅延である。この 式のn (t)の項は受信雑音である。直接信号経路を含め全部でに個の信号経 路(すなわち、K−1個のゴースト信号)がある。複素係数a。
の大きさおよび位相はそれぞれi番目の信号経路の相対減衰定数および搬送波位 相である。i=1を直接信号経路とする。理想的には、多重通路の問題がないと 、a1=1、τ =0、またi−2,・・・Kについてはa■+ 一〇である。
(2)式で表されるゴーストのある信号は、元の信号にノイズが加わったものの 多様な変形信号の加重合計にすぎない。ゴーストを除去する構成を持たない受像 機では、r (t)の実部または大きさのどちらがが復号される。これがテレビ ジョン表示において重なり合った画像、すなわちゴーストを発生する。
ビデオ信号の同相成分および直角成分に作用するゴースト消去システムが提案さ れている。このようなシステムの一例は「適応型テレビジョンのゴースト除去シ ステム」という名称の米国特許第4,703,357号明細書に示されている。
この特許に示されるシステムは、受信されるビデオ信号からゴースト信号を適応 的に取り除くために複素フィルタ係数を有するディジタルIIRフィルタを使用 する。このIIRフィルタは、最初は直接信号に対してゴースト信号の時間遅延 を決定するために相関器として動作し、その後ゴースト消去フィルタとして動作 するよう制御信号により条件付けられる。この相関器によって決定される遅延値 はFIRフィルタのタップを形成する可変遅延要素中に予めセットされ、またタ ップ係数は予め定められる値に設定される。次いでゴースト消去システムはトレ ーニング期間の間濾波されたビデオ信号をモニターし、またトレーニング期間に おけるフィルタの出力のゴースト信号を最小にする方向に係数値を変える。
前述の特許に記載されたシステムは一定のサンプリング周波数でゴーストのある ビデオ信号をサンプリングするアナログ・ディジタル変換器(ADC)を含んで いるディジタル・システムである。したがって、ゴースト消去フィルタにより処 理される信号は離散的サンプリング点てのみ有効である。これは、ゴースト信号 がザンプリング期間の整数倍でない時間量だけ主信号に対して遅延している場合 問題になる。この場合、フィルタにより消去されるゴースト信号は受信ビデオ信 号に混入する信号と同じではないかも知れず、従って、ゴースト信号のかなりの アーティファクト(a+1ilacl)が処理済みのビデオ信号中に残っている 。
より簡単なIIRフィルタを使用し、サンプリング期間の整数倍でない時間遅延 を有するゴースト信号を補正することのできるディジタルのゴースト消去システ ムが提供されるならば好都合である。
発明の概要 本発明は、トレーニング期間の間、受信されたゴーストのあるビデオ信号を解析 する制御処理回路を含み、ゴーストのあるビデオ信号を発生した伝送チャネルを 記述するチャネル・モデルを定式化するゴースト消去システムにおいて具体化さ れる。次いで制御処理回路はこのチャネル・モデルを複素11Rフイルタ用に1 組の係数値に変換する。本発明の第1の実施例では、データ処理回路により発生 される1組の係数値は、サンプリング期間の整数倍てない相対時間遅延を有する ゴースト信号を有効に消去するために複数の時間順次のサンプルの中から遅延サ ンプルを補間するようIIRフィルタを条件付ける。本発明の第2の実施例では 、このフィルタはサンプリング期間の整数倍でない相対時間遅延を有するゴース ト信号を有効に消去するためにIIRフィルタの少なくとも1つのタップに結合 されるサンプル値補間回路を含んでいる。
図面の簡単な説明 第1図は、前述したように、本発明が動作する環境を説明するのに自相な直接信 号に対する任意のゴースト信号の波形図である。
第2図は、本発明の一実施例を含んでいるテレビジョン受像機の信号処理部のブ ロック図である。
第3図および第4図は、第2図に示す回路の動作を説明するのに有用な信号を示 す振幅対時間の波形図である。
第5図は、第2図に示すマイクロプロセッサの動作を説明するのに有用なフロー チャートである。
第6図は、第2図に示す回路で使用するのに適したゴースト除去フィルタのブロ ック図である。
第7図は、第6図に示すゴースト除去フィルタに使用するのに適した複素乗算器 のブロック図である。
第8図は、本発明の他の実施例を示す第6図に示したゴースト除去フィルタの1 部のブロック図である。
第9図は、第8図に示す本発明の他の実施例で使用するのに適した補間器のブロ ック図である。
第10図は、第9図に示すROMのプログラミングを示す値のテーブル図である 。
先に引用した工rチ・セディック(H,TI+edick)氏による論文で説明 されているように、テレビジョン信号に多重通路歪みを加える伝送チャネルは有 限インパルス応答(FIR)フィルタとしてモデル化される。セディック氏の論 文で更に述べられているように、この伝送チャネルにより導入されるゴースト信 号は、伝送チャネルの伝達関数の逆関数である伝達関数を有する無限インパルス 応答(IIR)フィルタにより効果的に消去される。
後述する本発明の実施例では、調整可能な遅延値および調整可能な係数値を有す る複素11Rフイルタはゴースト消去フィルタとして使用される。時間遅延に関 する値および係数値はNTSC方式のテレビジョン信号と共に伝送されるトレー ニング信号からマイクロプロセッサにより決定される。これらの値は如何なるゴ ースト消去動作よりも111にゴーストの混入した信号を解析することにより決 定される。それから遅延値および係数はIIRフィルタに供給され、受信された ビデオ信号からゴースト信号を消去するようにIIRフィルタを条件付ける。
次に、第2図に示すゴースト消去システムを簡単に説明する。その後ゴースト消 去フィルタに関する遅延値および係数値が決定される方法を詳細に説明し、それ からゴースト消去フィルタそれ自体について更に詳細に説明する。
各図において、幅広の矢印は多ビットの並列ディジタル信号を伝達するバスを表 し、また単線の矢印はアナログ信号もしくは単一ビットのディジタル信号を伝達 する結線を表す。各装置の処理速度に依り、補償用遅延が信号経路の成るものに おいて必要とされる。ディジタル信号処理回路の設計分野の当業者はこの種の遅 延が個々のシステムのどこで必要とされるかを知っている。
第2図を参照すると、テレビジョン受像機の信号処理部が示されている。無線周 波(RF)信号はアンテナ208によって受信され、チューナおよびIF回路2 10にflj給される。回路210は、例えば、普通のテレビジョン・チューナ 、中間周波CIF)フィルタおよび増幅器を含んでいる。この実施例では、IF フィルタの通過帯域は変調された音声信号を含んでいることが望ましい。
回路210により発生されるIF倍信号ベースバンドの複合ビデオ信号C■を発 生する普通の包路線検波器242に供給される。普通の同期分離回路244は信 号C■に応答し、複合ビデオ信号から複合同期信号C8を取り除く。また同期分 離回路244はビデオ信号の各水平線から色同期バースト信号成分を抽出するの に使用されるパーストゲート信号BGを発生する。
複合同期信号C8に応答する検出器246は垂直同期パルス期間に先行する最後 (6番目)の前置等化パルスを検出する。検出器246は複合ビデオ信号の各フ ィールドの6番目の前置等化パルスと実質的に一致する出力パルス信号VSを発 生する。以下に説明するように、このパルスは直接信号に対してゴースト信号の 相対遅延、相対振幅および相対位相を決定するために使用されるトレーニング信 号を位置決めするために使用される。
チューナおよびIF回路210により発生される信号は、第1の同期検波器22 0、画像搬送波抽出回路222および第2の同期検波器230に供給される。画 像搬送波抽出回路222は、直接のビデオIF信号の画像搬送波と位相および周 波数が整合している第1の基準信号を発生する。この第1の基準信号は第1の同 期検波器220および90°の位相シフト回路224に供給される。
位相シフト回路224は、第1の基準信号に対して直角位相である第2の基準信 号を発生する。この第2の基準信号は第2の同期検波器230に供給される。同 期検波器220および230はIF倍信号それぞれ同相成分および直角位相成分 に復調する。同相信号は、システムのクロック信号CKに応答してディジタル信 号r1を発生するアナログ・ディンタル変換器(ADC)232に供給される。
同様に、直角位相信号は、クロック信号CKに応答してディジタル信号rQを発 生するADC234に供給される。例えば、NTSC方式の色副搬送波周波数f  の4倍に実質的に等しい周波数4f を有するりCC ロック信号CKは以下に説明する位相ロックループ(PLL)260により発生 される。
信号r およびrQはゴースト除去フィルタ280お■ よびマイクロプロセッサ282に供給される。以下に説明するように、ゴースト 除去フィルタ280は複素サンプル中データIIRフィルタを含んでいる。フィ ルタ280は、マイクロプロセッサ282の制御の下で、ゴーストが混入してい る信号「 およびrQを濾波し、任意のゴース)・信号が実質的に排除された直 接信号の同相成分を近似する信号r ′を発生する。信号r ′は、デI インタル信号r ′を表わすアナログのベースバンド複合ビデオ信号を発生する ディジタル・アナログ変換器(DAC)286に供給される。
アナログのヘースバンド複合ビデオ信号は、同期分離回路244から供給される パーストゲート信号BGに応答し、複合ビデオ信号の各水平線から色同期バース ト成分を分離する普通のバースト分離器288に供給される。
分離されたバースト信号は、例えば、はぼ4f の共振周波数を有する共振クリ スタル261を含んでいる普通のP L L 260に供給される。PLL26 0はバースト信号により制御され4f のクロック信号CKを発生する。
DAC286からの複合ビデオ信号は、普通のビデオ信号処理回路290とイン ターキャリヤ音声IF増幅器および検波回路292にも供給される。ビデオ信号 処理回路290は、例えば、複合ビデオ信号からルミナンスおよびクロミナンス 信号成分を分離し、これらの成分を処理して表示装置(図示せず)に供給される 赤、緑および青の原色信号(R,GおよびB)を発生する回路を含んでいる。イ ンターキャリヤ音声回路292は、複合ビデオ信号から4.5MHzの音声搬送 波を分離する共振同調回路、4.5MHzのIF増幅器および音声信号を発生す るためのFM復調器を含んでいる。音声信号はスピーカ(図示せず)に供給され る音声信号を発生する音声信号処理回路294に供給される。
マイクロプロセッサ282は直接記憶アクセス(DMA)命令、標準算術命令お よび割り込み処理機能を有し現在市販されている多数のマイクロプロセッサの中 の任意の1ってよい。マイクロプロセッサ282はランダムアクセスメモリ(R AM)284および読出し専用メモリ(ROM)285に結合される。またマイ クロプロセッサ282は、現在選択されているチャネルを示すチューナおよびI F回路210からの信号SELを受は取り、6番目の等化パルス検出器246か ら供給される信号VSを受は取り、クロック信号CKを受は取るように結合され る。マイクロプロセッサ282は、以下に説明するようにバスADおよびDAT Aを介して種々の信号をゴースト除去フィルタ280に供給するように結合され る。
ゴースト除去動作は、信号SELが新しいチャネルからの信号が選択されたこと を示すとき開始する。この時、パルス信号■Sに応答するマイクロプロセッサ2 82は、6番目の等化パルスに続く期間の間に発生する768個のr およびr QのサンプルをRAM284に記憶する■ ためにDMA命令を実行する。これらの768個のサンプルは入来信号の1つの 水平ライン期間のほぼ85%を構成する。この期間には、垂直同期パルスの前縁 に先行する少なくとも少数のサンプル、垂直同期パルスの前縁を表すサンプルお よび第1−の切り込みに先行する垂直同期パルスの部分を表すサンプルが含まれ ている。この信号はゴースト除去フィルタ280に供給される遅延値および係数 値を決定するためにマイクロプロセッサ282によりトレーニング信号として使 用される。
このトレーニング信号の同相部分の波形が第3図に示されている。6番目の等化 パルスと垂直同期の前縁間の信号の部分は水平ライン期間の0.46倍(0,4 6H)の継続時間およびQ  IRE単位の公称振幅を有する。
垂直同期の前縁と垂直同期パルスの第1の切り込み間の信号の部分は0.43H の継続時間および−40IRE単位の公称振幅を有する。雑音が無い場合、信号 のこの部分の間−40IRE単位の振幅値からの如何なる偏差も垂直同期の前縁 が遅延され、減衰され、またおそらく位相シフトされた変形であるゴースト信号 の結果であるものと仮定する。
この信号はほぼ一定のタイミングおよび振幅特性を有するよう連邦通信委員会( FCC)の規定により制約を受ける。この信号かFCCの規定に近い許容値に保 持されるので、この信号のモデルは受像機の製造中ROM285に記憶され、次 いで受信された信号と比較し、また直接信号に対するゴースト信号のタイミング 、振幅および位相を決定するためにマイクロプロセッサ282により使用される 。このモデルはビデオ信号の同相成分および直角成分の両方を含んでいる。
以下に説明するように、ROM285に記憶されるトレーニング信号のサンプル 値は受信されたビデオ信号より更に高い実効サンプリング周波数を有することが ある。
第4図は垂直同期の前縁の同相成分Iおよび直角位相成分Qを示す波形図である 。第4図に示す各々の波形図を描く144個のサンプルはROM285に保持さ れる。
これは64f  の実効サンプリング周波数を与える。トレーニング信号は垂直 同期の前縁に続く領域において実質的に平坦であるので、トレーニング信号の残 りのサンプル、すなわち、垂直同期の前縁と第1の切り込み間の期間を描くサン プルは2つの波形図の最後のサンプルを複製することによって発生される。
前述したように、本発明の実施例において、マイクロプロセッサ282は、トレ ーニング期間の間に抽出されるサンプルを処理し、伝送チャネルの伝達関数を近 似するモデルの伝達関数を発生する。このモデルの伝達関数から得られる遅延値 および係数値を使用して、マイクロプロセッサ282は受信されたビデオ信号か らゴースト信号を消去するようにIIRゴースト消去フィルタ280をプログラ ミングする。
以下に示すチャネル・モデルを発生するために使用されるアルゴリズムの説明に おいて、ベクトル表記法の拡大使用が行われる。一般に、信号は、例えば、s  (t)のように小文字活字によって表され、信号の倍数サンプル値を含んでいる ベクトルは、例えば、Sのように下線を引いた小文字活字によって表され、サン プル値の行列は、例えば、Sのように大文字活字によって表される。
以下に説明するアルゴリズムは、ディジタル・エクイップメントψコーポレーシ ョン(口1g1tal Equipme++t C。
+po+alion)によって作られたVAXコンピュータで実行するためにF ORTRAN77プログラミング言語で実行されてきた。これらのプログラムの コピーはこの明細書の説明部分の終わりに含まれている。これらのプログラムは 、プログラミング・リストに含まれていないサブルーチンC3VDSに対する参 照を含んでいる。このサブルーチンは複素要素を有する行列の特異値分解を実行 する。これはL I NPACKライブラリーパッケージから一般大衆にも入手 可能であり、また応用工業数学協会(SIAM)から発行されているジェイ・ジ ェイ・トンガラ(1,J、  Donga+a)氏にょるrLINPAcKの利 用者用手引き」の C122頁−0129頁に記載されている。
これらのプログラムは、第4図を参照して先に説明したように垂直同期の前縁の 同相成分および直角成分を表し、ROM285に読み込まれる値のリストも含ん でいる。これらの値は後述する基準信号行列および基準信号ベクトルを発生する ために使用される。本発明の実施例において、付録に含まれているようなコンバ イリング・プログラムによって得られる目的コードは垂直同期の前縁を示すデー タと共にROM285に記憶される。
伝達関数を発生するアルゴリズムにおける第1のステることである。先に述べた ように、マイクロプロセッサ282は、6番目の等化パルス検出器246から発 生される信号■Sに応答し、各ADC232および234により供給される信号 r1 (t)およびr Q (t )の各々の768個のサンプルを取り入れる 。信号rx(t)およびr、(t)のサンプルを含んでいる複素ベクトルは−L である。ROM285に記憶される基準信号の16番目毎のサンプルを抽出する ことによって得られる基準トレーニング信号のサンプルを含んでいる複素ベクト ルは王で表される。したがって、垂直同期の前縁は9個の複素サンプル値で表さ れる。
アルゴリズムにおける次のステップはベクトルrおよびユを微分することである 。これはベクトル中の各成分について(i+1)番目の成分から1番目の成分を 引くか、あるいはベクトルLおよび互の各々に(3)式で定められる行列Aを掛 けることによって達成される。
1−1  +1.  OO・・・0 0110−1.+10・・・0 01 A =   I  OO−1+1−−−0 011・・・・・・・・・・・・・ ・・・・・1(3)10 0 0 0 −1+11 この結果得られる微分されたベクトルはそれぞれr′および1二で表される。こ の手順における次のステップはニー−および−シj−を相関させてベクトルr′  を発生す□ m ることである。このベクトルの成分は(4)式に従ってここで、s′ *はS′ の共役複素数を示す。r′□[n]の絶対値の大きさは垂直同期の前縁のゴース ト画像がサンプル時間n+4に存在するかどうかを示す。
r ’ m [n ]の大きさがしきい値以上ならば、マイクロプロセッサは指 数nの値に隣接するL′□の値を調べ、しきい値を越える1群の値の中に中央値 qを検出する。
サンプル指数q+4は予備のゴースト遅延値として使用される。このテストで使 用されるしきい値は受信されたビデオ信号の信号対雑音比(SNR)の関数であ る。しきい値は低レベルのゴースト信号でも検出できるように出来るだけ低いこ とが望ましいが、入力ビデオ信号中の雑音がゴースト信号として解釈される程低 くはない。本発明の実施例において、しきい値は差分雑音値の自乗平均の平方根 (RMS)の2倍に設定される。RMS雑音値は、例えば、幾つかのフィールド 期間に亘って垂直帰線消去期間の間ビデオ信号をモニターすることにより決定さ れる。
受信されたビデオ信号中に直接信号およびに一1個(但し、Kは整数)のゴース ト信号が存在していると仮定すると、前述のアルゴリズムによって与えられる推 定遅延値は(5)式で記述されるベクトルヱで表される。
ここで、?、は1番目の経路についての遅延推定値である。ゴーストのあるトレ ーニング期間はこれらの遅延値を使用してROM285から得られる基準同期値 を重み付けすることにより合成される。合成されたトレーニング期間は(6)式 で表される。
システム雑音のため、この誤差は零ではあり得ない。入についての(10)式の 最小2乗(LS)解は誤差ベクトルの最小ノルムを与える。この解は次式で表さ れる。
ここで、I−■は共役−転置演算子である。最小2乗誤差ベクトル且のノルム平 方として定義される最小2乗誤差εここで、IはNxNの恒等行列である。
この最小平方誤差およびベクトルa  についての最LS 小平方解は推定遅延ベクトルエの関数である。前述のエツジ検出アルゴリズムは ゴースト信号の相対遅延のおおよその推定値を与えることができるだけなので、 この最小平方誤差は遅延推定値の精確さを示すために使用することが出来る。す なわち、誤差が小さければ小さいほどする基準信号行列は合成係数ベクトルの最 適の推定値を計算するために(11)に代入される。検索手続きの間jの全体は サンプリング期間の整数倍であるものと制限されていないことに注意しなければ ならない。
ROM285に記憶されるオーバサンプリングされた垂直同期信号は、先に計算 された2′の各々に加えられる端数の遅延時間を決定するために本発明の実施例 においては使用される。後述のアルゴリズムにおいて、最小た遅延値より実際の 遅延に近いかどうかを決定するために使用される。
発明者らは、最小平方誤差を定める(12)式が行列ユニタリ行列、■はK x  Kのユニタリ行列、Aは(14)式で表わされるNxKの対角行列である。
1 λ  0 0  ・・・ 01 A−10λ   O・・・ 01     (14)1・・・ ・・・ ・・・  ・・・ ・・・110 00  ・・・  λに1 1・・・ ・・・ ・・・ ・・・ ・・・1+000001 ここで、λ  λ  ・・・λえの値は特異値分解演算1 ′   2′ により決定される。
UlがUの一番左のに列として定義されるならば、(12)式は次式のようにな る。
行列の特異値分解は1979年ジョン ウィリー・アンl’−サンズ(Iohn  Wiley & 5ons)から発行された[コンピュータのためのコンパク ト数値方法−線形代数および関数最小化」という題目のジエイ・シー・ナツシュ (J、  C,Na+h)による本の2.5章の第21頁−第22頁に記載され ている。
最小平方誤差を求めて?空間を検索するためのアルゴリズムは第5図のフローチ ャートに示される。
第5図のフローチャートにより記述されるアルゴリズムは一度に一回ゴースト遅 延値を調整する。各々の調整によって、この値の局部最小が検出されるまで最小 平方誤差は求められる。
ステップ512で、変数■は1の値に設定され、この変数は調整される遅延時間 の値を保持する。ステップ514で、ベクトルΔ−1のすべての値は1番目の成 分を除いて零に設定される。この成分はROM285に記憶された基準サンプル ・ベクトルのサンプル時間T /16(すなわち、1./(64f  ))に設 定される。
する基準データからのサンプルを使用することと同じである。関数ε(9)の値 が1番目のゴーストについてのゴースト遅延が増加するとき小さくなるならば、 ステップ518および520は関数ε(、、i)についての最小値に達するまで このゴースト遅延を増大し続ける。
小値に達するまでステップ522.524および526は1番目のゴーストにつ いての遅延を減らす。1番目のゴースト信号について最適の遅延値が得られると 、ステップ528はもっと多くのゴースト信号が最適化されるべきかどうかを決 定する。IがKより小さいと、Iはステップ530で増加され、次のゴースト遅 延値は最適化され、さもなければアルゴリズムはステップ532で終rする。
前述の処理手順を使用すると、合成遅延ベクトルおよび合成係数ベクトルの最適 推定値を得ることができる。
送チャネルについての伝達関数の最適推定値、すなわち、ho、t (t)は( 16)式で与えられる。
、および土。、tのそれぞれのi番目の成分である。
尤に述べたように、一旦、伝送チャネルの伝達関数が決まると、この伝達関数の 遅延値および係数値をIIRフィルタに挿入することによりゴースト消去フィル タを設計することができる。前述のアルゴリズムにおいて、ゴースト遅延時間の ベクトルヱの開始時間は固定されていない。6番目の等化パルス検出器246に おける変動に因り、マイクロプロセッサ282により抽出されるサンプルのタイ ミングはフィールドがらフィールドで著しく変動する。この変動を補償するため に、本発明の実施例で使用されるゴースト消去システムは、IIRフィルタにお ける遅延値を最大の大きさを有する伝達関数における項とそれに続くすべての項 間の相対遅延値に設定する。最大の大きさを有する項は垂直同期の前縁に対応し 、先に述べたように、それに続くすべての項は垂直同期の前縁のゴーストに対応 するものと仮定する。垂直同期の前縁に対して相対的となるように伝達関数)1      <t)pt の遅延値を調整するのに加え、マイクロプロセッサ282は、垂直同期の前縁の 遅延時間以上がまたは等しい遅延時間を有する伝達関数の項の係数値を垂直同期 の前縁についての係数値の逆数により調整する。例えば、L番目の係数値が垂直 同期の前縁に対応すると、(16)式%式% 垂直同期の前縁値の時間遅延値より小さい時間遅延値フィルタを合成する際に無 視されるから、これらの項に対応する観測ベクトル玉の値も無視される。これら の値を除去する効果はNの値を減少させ、先に述べたアルゴリズムについての計 算オーバヘッドを減少させる。サンプル数におけるこの減少は観測ベクトル玉に おいて垂直同期の前縁を位置決めし、例えば、10個のサンプル期間より多く垂 直同期の前縁に先行するサンプルを除去することによって達成される。これは前 述の行列演算セットのどれよりも先に行われる。
リング期間の整数倍で発生する場合、IIRフィルタの応がある。
本発明の実施例に使用されるゴースト除去フィルタ280が第6図に示される。
このゴースト除去フィルタにおいて、受信されるビデオ信号の同相成分子1およ び直角位相成分子Qは複素乗算器610の信号入力ボートに供給される。複素係 数a、。9.の逆数の実数部および虚数部は、アドレスバスADおよびデータバ スDATAを介してで−rクロプロセッサ282により複素乗算器610の係数 入力ボートに供給される。本発明の実施例において、マイクロプロセッサ282 はあたかもRAM284中に値を記憶しているかのようにフィルタ係数値を複素 乗算器中に記憶する。
第7図は本発明の実施例で使用するのに適した複素乗算器のブロック図である。
第7図において、バスDATAは8ビツトの並列人力、並列出力の2つのレジス タ704および706のそれぞれの入力ボートに結合される。
アドレスバスADは復号回路702の入力ボートに結合される。回路702は、 バスADを介して供給される1つの予め定められる値に応答し、レジスタ704 のための負荷入力信号をパルス化し、別の予め定められる値に応答し、レジスタ 706のための別のロード入力信号をパルス化する。マイクロプロセッサ282 がゴースト除去フィルタ280で使用される複素乗算器の1つに係数値を取り入 れるとき、マイクロプロセッサ282はバスDATAを介して係数の実数部およ びバスADを介してレジスタ704のためのロード信号をパルス化するアドレス 値を同時に供給する。次に、マイクロプロセッサはバスDATAを介して係数値 の虚数部およびバスADを介してレジスタ706のためのロード信号をパルス化 するアドレス値を供給する。
ビデオ信号I およびQ□の同相成分および直角成分には、乗算器710.71 2.716および718、減算器714および加算器720を含んでいる回路に よりて、同相および直角係数値【2およびQ2のそれぞれが掛けられる。乗算器 710は信号11および係数値I2の債を形成し、その結果を減算器714に供 給する。減算器714は乗算器710によって発生される積から、乗算器712 によって発生される信号Q1および係数Q2の積を引き算し、同相の出力値I3 を発生する。乗算器716は信号Q および係数値I2の積を形成し、この積を 加算器720の1つの人力、ボートに供給する。乗算器718は信号I に係数 値Q2を掛は算し、加算器720への第2の入力信号を発生する。加算器720 から発生される出力信号は直角信号Q3である。
乗算器610によって実行される乗算演算は(1−7)式の総和の各項の分母( すなわち、1/aLOPT)に対応する。これは、この係数が垂直同期の前縁に 対応するので、ビデオ信号の直接成分の配分であると考えることができる 第6図を参照すると、乗算器610にから供給される同相信号および直角信号は それぞれ減算器618および620の各披減数入カボートに供給される。減算器 618および620の減数人力ボートは加算器666および668のそれぞれか らゴースト補正信号を受は取るように結合される。減算器618および620か ら供給される出力信号は補正された同相信号r ′および直角信号rQ′である 。先に述べたように、第2図を参照すると、信号r ′はゴースト除去フィルタ 280の出力信号で■ ある。
ゴースト消去システムの説明を簡単にするために、ゴースト除去フィルタ280 はこの最初の説明の目的のために3つだけのタップを有するように制限される。
後述するように、このフィルタは多数のタップを有することか望ましい3、この 種のフィルタは第6図に示すフィルタの直接拡張としてディジタル信号処理回路 の設計分野の当業者によって容易に設計される。
プログラム可能J能な3つのタップ付き遅延線622はIIRのゴースト除去フ ィルタに対して3つの遅延された信号を供給する。遅延線622は直列接続され たN個の遅延要素の2つの群を含んでいる。これらの遅延要素の各々はサンプリ ング・クロック信号CKに結合され、約70nsの時間遅延(1/ (4f   ) )を与える。最初の2つ、624と626、および最後のもの628により 示される遅延要素の第1の群はビデオ信号同相成分を遅延する。ビデオ信号の直 角成分は最初の2つ、630と632、および最後のもの634により示される 遅延要素の第2の群によって遅延される。
同相および直角遅延チェーンの対応する遅延要素の出力ポートはクロスバ−切り 換え行列における共通のスイッチ素子に接続される。例えば、遅延要素624お よび630の出力ポートはスイッチ素子636に結合され、切り換え素子636 を介してスイッチ素子638および640に結合される。同様に、遅延要素62 6および632の出力ポートはスイッチ素子642.644および646に結合 され、遅延要素628および634の出力ポートはスイッチ素子648.650 および652に結合される。これらの3つの組のスイッチ素子の各々はクロスバ −切り換え行列の列を定める。行列の各行のスイッチは互いに結合され、また、 それぞれ異なる複素係数乗算器に結合される。スイッチ636.642および6 48は乗算器656に結合される。スイッチ638および6・14および650 は乗算器658に結合される。スイッチ640.646および652は乗算器6 60に結合される。複素乗算器656.658および660は第7図に関して前 述した複素乗算器610と同一である。
行列のスイッチ素子の各々は、スイッチ選択論理回路654によって供給される 行選択信号(R8I、R82またはR33)および列選択信号(C3I、C82 またはC33)に応答し、選択された遅延要素対によって供給される出力信号を 乗算器に結合する。本発明の実施例において、どの行でも1つのスイッチだけが 、またどの列でも1つのスイッチだけが所定の時間に付勢される。
スイッチ選択論理回路654はバスADおよびDATを介してマイクロプロセッ サ282から供給される値によって制御される。バスADから供給される値は、 3つの起こり得る行/列選択値の中のどれがバスDATを介してマイクロプロセ ッサ282から供給されているかを示す。これらの値はスイッチ選択論理回路に よって内部的に記憶され、遅延された信号を各乗算器に供給する指示されたスイ ッチ素子を付勢するために使用される。この種のプログラム可能なタップ付き遅 延線は「ゴースト消去システムに関するクロスバ−スイッチ行列を含んでいるサ ンプル・データ・濾波システム」という名称の米国特許第4,727,424号 明細書に記載されている。
前述したように、本発明の実施例で、(17)式の遅延値はサンプリング期間の 整数の項の中にあるものとする。従って、(17)式の総和の中の第1項以外の 各項は同相および直角遅延線上の1つのタップおよび係数乗算器656.658 および660の中の1つに対応する。
従って、マイクロプロセッサ282は(17)式の係数値のそれぞれ異なる成分 である同相(実数)成分および直角(虚数)成分をバスADとDATAを介して 乗算器656.658および660の各々に供給する。
乗算器658および660から供給される同相信号は加算器662で合計される 。この加算器から供給される合成信号は、乗算器656から供給される同相信号 と加算器666で合計される。前述したように、加算器666から供給される信 号は減算器618に供給されるゴースト補正信号の同相成分である。
乗算器658および660から供給される直角信号は加算器664によって合計 され、その結果生じる和は加算器668に供給される。加算器668は乗算器6 56によって発生される直角信号にこの値を加え、減算器620に供給されるゴ ースト消去信号の直角成分を発生する。セディック氏の文献に述べられているよ うに、第6図に示すようなIIRフィルタは、受信されたビデオ信号の各々のゴ ースト信号成分に対して遅延線の1つのタップおよび1つの係数乗算器かあると 、ビデオ信号中の多重通路の歪みを効果的に消去する。
全体のゴースト消去システムの説明を簡単化するため、先に仮定したように、各 ゴースト信号はサンプリング期間の整数倍である時間量だけ直接信号に対して遅 延されている。これが当てはまらない場合、すなわち、ゴースト信号がサンプリ ング期間の整数倍でない時間量だけ遅延される場合、有効なゴースト信号消去を 達成するために前述のシステムの変更が望ましい。一般に、IIRフィルタは、 間に入る遅延された信号が利用可能な遅延信号から補間され、IIRのタップと して使用されるならば、タップ付き遅延線から得られる2つの連続する遅延時間 の間に入る遅延時間を有するゴースト信号を消去するように変更される。補間の 2つの方法については以下に説明する。第1の方法は、クロスバ−行列により多 くのスイッチ素子の行およびより多くの複素乗算器を加える第6図に示す構成の 直接拡張を含んでいる。また、マイクロプロセッサ282は、端数のサンプル遅 延を取り囲む夕・ツブに対して追加の係数値を計算するようにプログラミングさ れる。第2の方法は、可変遅延線の各出力ボートに補間回路を挿入することであ る。この補間回路は、マイクロプロセッサ282によって発生される信号に応答 し、ゴースト信号と事実上同じ遅延を有する同相信号および直角信号を発生する 。
次は、前述の第1の方法に対する追加の係数値を得るためのアルゴリズムを構成 する計算の説明である。すべてのゴースト遅延τ を次式のようにする。
k、=INT(τ /T  )       (18)1                 1       Sここで、T はサンプリング期間(1/4 f  o)であり、INT(x)はXより大きくない最大の整数を与える。
さらに、k を含み、またi番目のゴースト信号遅延を取り囲むことがあり得る 遅延線タップに対応する遅延値の群の中の1つの遅延値をP9.が表わすものと する。
J 後述するアルゴリズムにおいて、i番目のゴースト消去信号は遅延線の連続する J、のタップから補間される。
等価チャネル応答を得るためにアルゴリズムで使用された係数値からこのアルゴ リズムの係数値を差分するため、この係数値はす1.なる記号で示される。さら に、このJ アルゴリズムで使用される複素行列式を簡単にするために、(17)式の値しは 1に等しく、また垂直同期の前縁はτ1=0で発生すると仮定する。
前述した仮定を(17)式に取り入れると、(19)式が得られる。
この式で、係数す1.の最適値は既知でない。
J 前述したように、観測ベクトルエは(20)式によってきっちりと近似される。
ここで、5optはアール。ptである(8)式によって示される基準信号行列 である。
このアルゴリズムの場合、基準信号行列の最後の行、1(N)だけを考察する必 すなわち、ベクトルi。pt 要がある。この制限下で、(20)式は次のようになる。
ここで、 (19)式によって定められる伝達関数から、(21)式と同様の(23)式が 得られる。
y[N]=i[N]与              (23)ここで、y [N ]はr [N]の近似であり、s [Nコは(24)式によって定められる基準 信号ベクトルである。
Is ((N−Pll) T5) 1 s [N]  =  l s ((N P21) T5) l          (24)また、係数ベクトル互は(25)の式によって定められる。
により係数値すを定めることを決定した。この解は(26)式で表わされる。
−T −1へT br、yusa−(h (i、*s ) )   (E (i*s   ) )  a    (26)−m−−opt  −opt ここでベクトルがランダムベクトルとして扱われると、*は共役演算子を表し、 Tは転置演算子を表し、E()は期待値演算子を表わす。時間指数Nは処理が定 常状態(s * s ” o p t )はi*とj−”opt間の相互分散行 列と等価である。(26)式の数値を求める計算において、ベクトルエおよび直 。2.は0平均ベクトルに変換され、ビデオ周波数範囲内で平坦なスペクトル密 度を有するものと仮定する。分散行列の解説を含んでいるランダム変数およびそ れらの特性の解説はジー・アール・クーパー(G、  R,Cooper)氏お よびシー・ディー・マクギレム(C,D、 McGille+n)氏による「信 号およびシステム解析の方法」という題名の本(1967年ウィンストンのリネ ハルトのホルト社発行)の第9章−第10章に見られる。
この方法によって決定される係数ベクトルSEは、遅延値が(27)式で定めら れるベクトルPによって決定されるIIRフィルタに対する係数値を設定するた めに使用される。
IP    + 1.1 +p2.11 P=lP      l                 (27)−2,2 発明台らは、ゴースト遅延時間τ、かに一、 T 8に実質的に等しい場合はJ  が1の値に設定されることを定めた。
しかしながら、τ、とに、T  間に著しい差がある場合、l       I     S 4の値にJ を設定すると満足な結果が生じる。この方法を使用して、端数のゴ ースト遅延前後の2つのサンプルは端数のゴースト遅延を有する等価サンプルを 得るように補間される。
前述のアルゴリズムはFORTRAN  77のコンピュータ言語で書かれたコ ンピュータ・プログラムで実行される。このプログラムは本出願の付録に含まれ る。
前述の捕間方法に代わるものは、第8図に示すように、プログラム可能な遅延線 622の出力ポートと係数乗算器656.658および660の信号入力ポート 間にサンプル値補間回路を挿入することである。第8図に示す補間器810.8 12および814の各々は、例えば、9図に示す型式のものである2つの補償線 形補間器を含んでいる。第9図に示す補償補間器は、補間されたサンプル値を発 生ずるように比例因数Fによって連続するサンプル値を平均化する線形補間器9 20を含んで0る。
サンプルを平均化することによって、補間器は低域通過フィルタの機能を実行す る。ビデオ・サンプルのこのフィルタリングは補間された信号中に望ましくない 位相および振幅誤差を発生する。これらの誤差を補償するために第9図に示す補 間器は補償フィルタ950を含んでいる。補間された信号に供給される補償量は 比例因数Fの値による。
第8図を参照すると、本発明の実施例では、比例因数はマイクロプロセッサ28 2によって補間器810.812および814に供給される。補間器の中の両方 の補間回路には同じ因数値Fか供給される。しかしながら、補間器810.81 −2および814の各々に供給される値Fは、この値かi番目のゴースト信号( すなわち、τi=τ1opt   Lopt)に関連する遅延値τ、の−τ                        l端数部であるから異なる。値F は、バスADを介して復号器930に供給される値が入力ポートに供給される値 をロードさせるレジスタ932にパルスを与えるように復号器を条件付ける一方 、バスDATAを介して補間器のレジスタ932の入力ポートに供給される。レ ジスタ932の値は、値(1−F)を線形補間器920に供給し、値Cを乗算器 936に供給するROM934のアドレス入力ボートに供給される。Cの値は加 算器960により捕間信号に加えられる補償信号の大きさを決定する。
第10図は、(1/8)T  の粗さを有する補間回路のためのROM934の プログラミングを示す図である。
第9図に示す補間器は「振幅および位相補償を有するディジタル遅延補間フィル タ」という名称の米国特許第11.694,414号明細書に詳細に述べられて いる。
伝送チャネルについての計算された伝達関数からゴースト消去フィルタを発生す るための前述のアルゴリズムに関連して、この補間器は次のように使用される。
マイクロプロセッサ282が(17)式からの遅延値をゴースト除去フィルタ2 80のプログラム可能な遅延線622にf41.給するとき、遅延線622と補 間器810.812または814の中の選択されたものとの間の遅延値を配分す る。プログラム可能な遅延622によって実現される遅延値τ、の量は((τ  /T  ) −1)xT8な+                   1        3る量の整数部に実質的に等しい。遅延値T の残り十T、の分数部 は、比例因数Fとして値τ、の端数部を供給することによって選択された補間器 で実現される。遅延の余分のサンプル期間は前述の第10図のテーブルの合計遅 延列に示すように、選択された補間器に固有のものである。
[111述したゴースト除去システムについてのシステム雑音の効果を補償する ために、伝送チャネルの幾つかの推定値(すなわち、(17)式)は幾つかのビ デオフィールド期間に亘って形成され、これらの推定値からの係数値および遅延 値はゴースト除去フィルタ280による使用のため最終的な係数値および遅延値 を得るために平均化される。
以上説明したゴースト消去システムはビデオ信号から多重通路の歪みを除去する ディジタル・システムであるが、本発明はアナ0グ要素を使用してすることもで き、また、他の形式の信号から多重通路の歪みを除去するために使用することも 考えられる。
付   録 先に説明した本発明の詳細な説明に関連する例示的コンピュータープログラムの リストを以下に示す。先に述べたように、これらの例示的なプログラムはデジタ ル二りイプメントコーポレーション(DigiNI EquipmentCo「 po「ation )により製造されるパックス(■^X)コンピュータで実行 するためにフォートラン(FORTIIAN) 77のプログラミング言語で書 かれている。本発明の原理の所定の応用例において無限インパルス応答フィルタ の係数−ラの要求を満たすようにこれらのプログラムに対して種々の変更が加え られることを理解されたい。
c d  for(de  ghosle+1common/WINDOW/N X、 NY、 NTcha+acle+ネ20  In injege+$2 II  I(910,10501,rx  1(910, 101050)inle+ nl  vs、  npx  vs  bgint ege+ 1xo11. 1yoll、  iysepinteger npx   ob intege「:xliy inlege+ ians inlege+npx  vsl、  npx  vs2inleget nr I   vs complex Hvs(910) integeInrxvsl inlegernl  av、  npx  av  bgc    comp lex+x  av[lO)complex sl(201 injege+ nsf、  n5flCO呻few sld(200) intege「n5ld、  n5ldl、  n5sldc dvl 5ta rt (define sampling +ate、 lilenamei) call 1nill  (240) call dvlmap(3) paint *、” paint $、’define 1nput file real pa「t ’call cMel  (1,4,1,0)川nl $、” print t、’define 1nput l山imagina「y pu t’cBll chdel  (2,4,L O)p目11凱″ p「int 本、’define outpu口ile real paI+’ call cMel  (3,4,1,G)paint !、’enfe+ s ummary l山name、  In ?’read (5,1)  In l   1ofmat  (a201 open  (10,file=ln//’、 hum’、 ttalus=’  unknown’)open (II、 file=fn//’ 、 daf ’、 s+alus=’ unknovi’1c  「ead  in vid eo data  fileprint t、’ −−−−−−−+eadin g 1npul 1ilei −一−−−−−’call +eadil  ( 1,4,0,o  +)p「int s、I 「eal pa口1nput d one ’type *、 cba+ (7) call +eadil (2,4,0,TX  1lprint 本、’ i maginary part 1nput done’type j cha「 (7) lr9e t、 cha「(7] c vertical 57nc l1nde+paint $、’−−−−− −−vertical 5ync 1inde+−−−−−−−’100 pr int *、’enter 1ine number、 begining p ixel、 to be viewed、’prinf札’  nl    v s、  npx    vs   bg  7’+pad f5.本) nl   vs、 npx  vs  bgcall 1nill  (240) ixoll =50 iyoH−100 iysep =300 c no、  of pixel+ to be observednpx   ob =91G Call wavefoIm (IX  1.II  i、 nl  Vl、  npx  vt  bg、npx ob。
&        1xoll、 1yoll、 1ysep)call mo vabs (ixol[、700)call anmode p「int *、’ 「ead 1G901 bottom wavef+om 、 I(topllo(botl ?’read  (5,*)  1ans 百(ians、 eq、 I) thencall 5can  read ( ix、  iy、  1xotl、  1yofl+1ysep)1se call +c+n  read (ix、  iy、  1xofl、  1 yoll)end:1 prinl  章、’ob+ene  signal、  Ilo  ?’「e ad  (5,?)  1ana 百(iant、ne、0) golo 100c deline ve山cal  observation 1nterval in 1ine nl  vs 200 print 転’enle+ begining of the ab setvalion 1nterval npx  ysl@?’ call +c+n  read (ix、  iy、  1xoll、  1 yoll)npx  vsl = it + npx  vs  bg −tp rint 札’bIgining at vs absetvalion 1n terval ’、npx  vs1p+inl t、’cbange、  I lo ?’tead  (5,11anS 百(iant、 ne、 0) goto 200300 p[int t、’ enter end of the ab+ervafiu 1nterval  npx  vi2 ?’call 5ctn  read (it、  iy 、  1xoll、  1yall)npx  vi2 = 1x + npx   vs  bg −1paint  車、’end  of  vs  ob servation  1nte 「val  ’、npx    o2p「i nt傘、’change、  Ilo ?’「ead (5,掌1 1ans 百(ians、 ne、 01 gofo 300c  +Qpy  VB   IXOinm  TK   I   II    1nrxvs=npx  v s2−npx  vsl+1print t、 ’ n+x vs’、 nrx   vsdoi=l、+++x  vs II vi(i) =cmplx (float(+x  +(npx  vs l+i−1,nl vs)l。
&           (float(+x  i (npx  vsl+i −1,nl  vs)))endd。
c active video monito「call 1nill  (2 401 paint 札’ −−−−−−−active video monito+  −−−−−−−’400   paint  傘、’enle+  1ine   numbe「 to  be  viewed、  nl    av、  @npx    ay   bg  ?’read (5,$1 nl  ay 、  npx  av  bgc 1xoli=5fl、  1yofl=10 0. 1−ysep=300.  npx  ob=910call 1nil l  (240) call wavetorm (rt  1.II  i、 nl  av、  npx  ay  bg、9]0.50,100,300)call mova bs  (5(1,7H)call  anmode print t、’ read top of bottom wavel+o m、  I (topllo(bat) ?’Iead  f5j)   1a ns 百(iant、 eq、 11 thencall +c+n  read ( ix、  iy、  5G、  4001else call +c+n  read (it、  i7. 50. 100)nd il paint !、’obse+ve anothe+ 1inc ol vid eo wayelo+m、 Ilo ?’+ead  (5,傘)  1ans 百(ians、 ne、 fil golo 400c  define  「 eference、  t ransform、  and  1ocale   vertical  5ync@ edges  Ioughly call 1ldell  (sl、  nsl、  n511.20..23 G、)call an++Iel  (am、  nam、  11call  amtlm (+y  vs、 nrx  vs、 am、 nam、  IX p、 nrxplcall amllm (sl、  nsl、  am、   naa  slp、  n5lp)call +xp match  slp  (「xpm、 n+xpm、 +xp、 n+xp、山、n山)ixoll =  50 iyoll =100 call 1nill (2401 call  vt    pH−山 −c  (+xpm、n+xpm、1xo ll、iyoll)Calしく「n −IM (lh、 1xo11.1yol l)call IM abs (+xpm、 n+xpm、 lh、 ltxp m、 n1rxp+++)call pause psep=3 write (10,3)  lh、  pS!II3   format ( /、’ tb+eshold = ’、 gll、 4e3.’ p  sep  =’、 151c 1ocate peaks in rx  vscall  peak  5earch (pk、npk、p  sep、lrxpm、n l+xpm、rxpm、nrxpBnsll)c 5earch min of  1east +qua+e lit in the vicini17 ol  pkQnwinl−10 nwin2 = 10 nals=npk Idsln+ =910 call 5lddell  (sld、  n5ld、  n5fdl、   ntsld、  20.、 230.)Call pk  ledge (le dge、nledg、nsl畦nrx of、pk、npk、nvinl、nv in2)call sl  mat (slm、n51m1.n51m2、ld slm、 ledge、nledge、 sf、nsl、n5fl)call  ledge  laug (taug、 nlaug、 n55g、  led ge、 nledge)call  +gg  del  (s師、n5ln+ l、n51m2.ldslm、  n55g、  tang、nlaug。
&         sld、  n5ld、  n5ldl、  n5sld lcall lsl  no  dc (a  Is、na   Is、  百 m、 ldslm、 n51m1. nslm2゜&+x  vs、nrx   vsl、  cond、 errot)W日1e  (6,41n   5ea rch、  errot、  condwhile (1G、41 i 5ea rch、 er「or、 cnnd4  10111111 (/、’n  s ea+ch=’、i5.’ c++or=’、gll、4e3゜&         cnnd :’、 gll、 4e3)write f6,5) (i、p k(D、a lt(i)、i=1.npklwrite (10,51(i、p k(i)、 a  lt(i)、i=I、npk)5   fo「mal t’   i = ’、 i5.’ pk(i) = ’、 i5゜&     ’  a  1fil ・(’、g11.4e3.’ 、 ’、g11.4e3.’) ’)do ipk = I、  npk pk  1p(ipk) :pk(ipk)endd。
do ipk =I、  npk pk  lp(ipk)・pk  tp(ipk)  + 1n  +ea+c h =n  5eatch + 1call pk  ledge (ledg e、 nledg、 nsl畦nrx ol、pk  lp、npk、nwin l、nwin2j call sl  mat (slm、n+1ml、n51m2.91f1.   ledge、nledge、  cl、nsl、n5llj call ledge  laug (faug、  nlxug、  n55 g、  ledge、  nledge)call sgg  del (sl m、n51m1.n51m2.ld+Im、  n55g、  tang、nt aB。
&          sl4  n5ld、  n5ld1.  n5sld )call lsl  no  dc (a  lt  +p、na  Is、   slm、910.nsl+nl、n51m2゜&              I!  vs、nrx  vcl、  cond  lp、  er「oTl p)write (6,4) n  tea「ch、  merlp、  co nd  lpW「山(10,41n−5earch、 errorlp、 co nd lp=+山(ら、5!  (i、  pk   1p(i)、  a    lt(i)、i=1.npk)write (10,5) (i、pk 1p (i)、 a lt(i)、i=l、npk)if (etrot  lp、   le、euu) then500   call copyl (a  Is   tp、 a  Is、 na  lt)pk(ipk)・pk  1p(i pkler「oI−eHllI   Ip cond = cond  Ip pi  lp(ipk)・pk  tp(ipk)  + 1n    5ea rch  =  n   5earch  +1call pk  lsdge  (ledge、 nledge、 n51m1. nrx  tsl。
&          pk  tp、npk、nv+nl、nv+n2)ca ll sl  mat (slm、nsfml、n5la+2.910. le dge、山dge、 sl、nsl、n5ll)call ledge lau g (laug、 nlaug、 n55g、 ledge、 nledge) Call sgg def (slm、n51m1.n51m2.Idslm、  n55g、 +aug、出■。
&        sld、 n5ld、 n5ld!、 n5sld)cal l lsl no dc (!  Is  lp、na  Is、 sim、9 ]0.nsfml、nsfm2゜&rx  vs、nrx  vsl、 con d  lp、  error 1p)W「由(6,41n  5earch、  error  fp、 cond  lpwrite (10,4) n 5e arch、 error  lp、 cond lpwrite (6,5)  (i、 llk tp(i)、 a lt(i)、i=I、npk)while  (1(1,5) fi、pk tp(il、 a Is(ト)、 i=l、  npk)百(errot  tp、 Ie、errot) gone 500p k  lp(ipk) :pk(ipk)1se pk  +pfipk) =pk  tp(ipk)  −1600pk  t p(ipk) =pk  lp(ipk)  −1n  +earch = n   5earch +1call pk  ledge (ledge、nle dge、n51m1. nrx  vsl。
&             pk、 npL nwinl、 nvin2)c all sl  mat (+Im、n51m1.n51m2,910. le dg!、nledge、 sf、++sl、nsfl)call lclIge   laug (laug、 nfaug、 n55g、 ledge、 nl edgelcall  sgg   del  (slm、n51m1.na1 m2.ldslm、  n55g、  taug、ntaug。
&         sld、  n1ld、  n+Id1.  nold) call lsf  no  dc (a  lt  ip、na  Is、  +Im、910.n+1ml、ns[m2゜&tx  vc、nrx  vsl 、 cond  lp、  carat  lp1w口le  (6,4)   n   sca+ch、  errorip、  cond   tpwrit e (10,4) n  +euch、 euor  jp、 cond  l pw目le  (6,5)  (i、  pk  tp(i)、  a  1t (i)、i−1,npk)!r山(10,5) (i、pk tp(i)、 a  lt(il、i=I、npk)百(e「rorlp、le、euor) th encall copyl (a  Is  lp、 a  lt、 na   It)pk(ipkl・pk  tp(ipk)eTIOI : e目a「1p con+1=cond   lp ndil it  (e「「o「   lp、le、erroi  golo  600p k  1p(ipk)・pk(1pk)ndil end+i。
c 5earch 5ub−It in laug。
call pkヨedge (ledge、山dg、 n51m1. nrx  vtl、 pk、 npk、 nvinl、 nwin2)call ledg e  laug (tang、  nfaug、  n+sg、  ledge 、  nledge)do ilaug = 1.  nlaugdlau−1 25 call copyl  t (taug、  laug  lp、  nla ug)do 1dlau :1. 4 c 5euch down to I/16 of l5dlau = dla u / 2゜ laug lp(itaugl = laug(ilaug) +dtauca ll  sgg  del  (slm、n51m1.n51m2.ldsfm 、  n55g、  laug  lp、nlaug。
&            if4  n+ld、  n5ld1.  ns+ 1d)call  l百−no  dc  (a   It  tp、na    Is、  slm、910.nsfml、n51m2゜l               Hvs、nrx  nl、  cond  jp、  euor   1p)n  5earch = n  5earch + 1white (6 ,4) n +ea+ch、 euor  Ip、 cond jpwhile  (10,4] n sea「ch、 er「orlp、 cond lpwt ile (6,21)  (i、taug  jp(il、 a  It  f p(i)、i=I、nlaug)while (10,21)  (i、lau g 1p(i)、 a Is 1p(i)、i=l、nlaug)21     1oonal (’  i =’、i5.’  laug(i)  =’、g1 1.4e3゜&        ’   a  lt(i)  ・ (’、g1 1.4e3.’  、  ’、g!1.4e3.’)’)if (euor   lp、lLe++o+) thentaug (iliug) = laug   tp (ilaug)etrot=eno+   1p cond = cond  1p call copyl (a  Is  lp、  a  Is、  na   Is)lse taug lp(ilaugj = laugfifaug) −dlauca ll +gg  del (slm、n51m1.n51m2.ldslm、   ngsg。
&            Iaug  lp、nlaug、  +ld、   n5ld、 n5ldl、  ns+Id)call lcl  no  dc  (a  Is  lp、na  Is、  slm、910.n51m1.n 51m2゜l              rx  vs、nrx  vsl、   cond −【p、  error  tp)n  Bea+ch = n   +earch + 1W「ile (6,4)  n  sea「ch、   errortp、  cond  jpwhile (10,41n  +e a+ch、  etfot  tp、  cond  tpwrite (6, 21)  (i、taug  1p(i)、  a  lt  tp(i)、i =I、nlaug)while (10,21)  住Laug  1p(i) 、  a  Is  1p(i)、 i=I、 ntaug)if (e「ro r  lp、It、morl thenlaug  (ilaug)  = l aug   lp  (自aug)c++or=c++o+   1p cond = cond   1p call copyl  (!  Is  lp、  a  Is、na  I t)ndil c+++l1l endd。
endd。
call +gg  del (slm、n51m1.nstm2.Idslm 、n+sg。
&           laug、 ntaug、 sld、n5ld、   n+ldl、 n5sld)call dc  lsl (d、 a  Is、   slm、910.n51m1.n51m2.  II  vs、nrx   vsl)write (6,6) v「ile (10,6) 6 1o+mat (//、’ −−−−−−−−−−−−−−−−−−−−一 −−−−−−−−−−−−−−−−−−−−−−−−−−−|’+/。
&        ′ 傘傘傘Hminimun  legsト1quare− sca+ch  result    傘傘*孝ネ1./)while (6, 4) n  +earch、  error、  condwrite (10 ,41n  5eaIch、  error、  condwrite (6, 22)  (i、pk(i)、taug(i)、a  lt(i)、i=1.n pk)write (10,22)  (i、pk(i)、laug(i)、   a  lt(i)、i=l、npk122  fo「mal Ci =’、i 5.’  pk(il =’、i5.’  Iaug(i) =’、g11.4 e3゜&     ’ a  lt(i) ・(’、g11.4e3.’ 、  ’、g11.4e3.’)’)c calculate filte「pa+a melerscalllaug   a  hg  4(hx、kx、kkx、 hy、ky、kkyjaug、a  Is、npklwhile  (6,9) wtite (10,9) 9  foImal (//、” −−−−ftu deghcsting p 「ocesso「−−−−−’)write  (6,7) d write  (10,7) d ?  to+mat (’ dc oflsel     d = (’、g1 1.4e3.’ 、 ’、g11.4e3.’ )’)write  (6,8 )  hx(1)write (10,8) lu(1)11  format  (’ amp−ph aius+ 1(1)=(’、g11.4e3.’ 、  ’、 gll、4e3.’ l’)nile (6,If)  (i、 ky (i)、 by(il、 i=1. kk7 )w「ite (10,II)   (i、ky(i)、by(i)、i=I、kky)It    fo 「ma l  ぐ   i  =  ’、  i5.’  ky(i)  : ’、   i5゜&     ’   hy(i) ・(’、g11.4e3.’ 、 ’ 、g11.4e3.’ )’)c deghosjing I 41irlds  !11 video daIx)c no ol oye「IIow/und e+Ilow、 nal pad on17if  (mod(山ne、  1 0)、 eQ、 01  paint  傘、’−−−  p+ocessiB   1ine  ’、 @1line do 1pixel =1,910 x=cmplx(float(rx +(ipixcl、1line))、Il ogj(rx  1(ipixel、1lins)))蟲     −d calldl山et (x、y、 hx、kx、kkx、 hy、に7.kk7 .5b)II  T   dg  (ipixel、 山ne)  =  1n l(teal(V))II   i    dg  (ipixel、  山n e’l  =  inl(aimag(y)1百(rx r  dg (ipi xel、山ne1. gl、 255) 1henn  ovll = n   ull + 11I   I   dg  (ipixel、 己1ne)   = 255elseif (rx r  dg fipixel、山ne)、  11. Ql thenn  udll−n  udll + 11I l   dg (ipixel、1linel 二〇ndil endd。
endd。
c output +csult summaryvtile (6,12)  n oyll、 n  udflnile (10,12) n  OF目、+ +ud月12  for+nal (/、’ −−−# of oye+Ilo w (255)/underllov(01in real …レー戟B &     /、’  n  ovll=’、i&、’  n  udll=’ 、18)close(10) c w+iHng deghosfed resultsprint *、’− −−−−−− w「iting deghosted +esulls −−− −−−−’call +uit!il (3,4,0,rx r  dg)do i=1.5 t7pe  j chir (7) endd。
c show vuelroms ofye山cxl 57nc(nl vs)  and active widen (nl it) 5■撃■モ撃■■ c before and gaie「deghosjing、  5tore  comparison in tile 11vtile (If、13)  nl  vslg    jovial  (’  ネ傘を傘傘参   ver tical  Bnt  H車掌ネ本ネ  fin  ’、i5./。
&    ipx    +eil in−−−−imB in      r eal out−−−imxg out’)do i = I、 910 while (II、14) i、rx +(i、nl vs)、【xi(i、 nl vs)。
&         rx  r  dg(i、nl  vI)、rx  i   dg(i、nl  ws)endd。
14  jovial  (lx、 i5.5x、 ill、 3x、 ill 、 7x、 i8.3L ill)while (11,15) 15   1+u+ul  Cl/、’ 傘傘傘傘 active  vide o  本市#Iネ   nl   xw=’。
l       i5.’npx  mV  hg=’、i5./。
&’   ipx    teal in−−−−imxg in        tell ouト−−−inlB out’)do i :1.910 ipX = i + npl  u  bg −1write (11,14)   ipx、  tx  +lpx、nl  av)、rx  1(ipx、n l  u)。
&        tx  +  dg(ipx、nl  av)、rx  +   dg(ipx、nl  av)endd。
close(II) call inl口(2401 call wavelo+m (rx  +、「1  i、nl  vs、   ]、 910. 50.450. 150)call waveform (I X  T  dg、 IX  i  dg、 nl  vs、 ]、 910. 50.100.150)Call pause call 1nill (240) call wavelo+m (4x  L[i、nl av、npx  ay  bg、910.5G、l[lO,300)call pau+e call 1nill  (240) call wavefo「m(+x−rdg、 IX i  dg、 nl a v、npx aマーbg、 910.50.100.300j call pause sl(Ill nd C*零ネ零零傘零零*参孝*$$$参*傘NN$**N傘孝傘*傘傘傘傘傘傘孝 ネ孝傘*傘傘*傘*傘傘ネ*傘ネ傘傘市*****傘傘零傘趨{$参市参傘ネ$ サブルーチンliog  a  tuY4 (hx、 kx、 kkx、 h7 . ky、口1. l*B、 s  Is、 nfsuOC#傘#零本ネ*ネ* 傘*ネ$参本零本本章傘ネ零*ネ*本本孝本本ネネ*ネ*本零零ネ傘傘ネネ*ネ *零零孝電傘ネ牟IIネ本$$参ネ*l本孝ネ参孝傘拳ネネ傘$ c generate comp!n5ajion l1lle+ coeff icients and deli7c  I+olIlidentified  channel ct+aracje+is+icc  for 2 pxt h@trial +oulinec cl: lpd for notatio nc  41aps for relative daIx7 which i s not infege+ times ojc   sampling i nlenalcomplexlu(11,hy(1)、  a  l5(1)t eal laug(1) inlege+ kx(1)、  kkx、 ky(Il、  kkycomp lex b(4) 「eat eu  lpd、  dlau2exteInal To  fla t it (nlaug、gl。2)  thenp「int 傘、’$1ttt  100 many palhs fat corrent implemenl ilion…本′lop ndil fu(1)・1. a  ls(1) 百(nlaugoge、 2) thenky  tp = 1nl((lau g(2)−1aug(1)) /Is )dlau2 =laug(21−1a ug(11−1+ $ ky  lpd(dtau2. It、 10. e− 3) 1henc smallest 5tep is O,015625ky (1)・ky  tp hym ・a It(2) / a ls(1)kky = 1 1se dlau2 = dlau2 + 1scall lpd (b、 4. et t  lpd、 dtau2. ls、 To  flat、 !、 15)k y (1)・kylp−1 by(II = a  1t(2) 柿(1) / a  1S(1)kV(2 +・k7  tp hy(21=a  1S(21tb(2)/ a  l+mky(3)・kyt p+1 by(3) = a  1t(2)傘b(3) / a  ls(1)ky(4 ) =ky  lp + 2 hy(4)  □  a    ls(2)  章b(4)  /  a     ls(1)hhy = 4 ndil ndil 「eturn nd c Id、to+ (litterdesign)c design degh oSling ime「c −−−−−−一−−−−−−−−−−data d elinilion+  −−−−−−−一−−−−−−−−−c Irefe 「ence vecto「and its 「elated datac  s l (s  +ele+ence)、 nsl (n  s  +ele+en ce)c  +Ip (3tel  pained、  n5lp (n  l   tel  p+imelc  +slp (「o  s  ref  p「 ime)c  coslp(coya+  s  ref  prime)   ncoslp(n  cova+  s  ref  prime)c     (colume scanning mat+1xlc 2.  t+anil o+malion mat口xc  am (i  matrix)v、  n am (n  a  mat「1x)c    (uSe onlVnone  xe+o entry in 「ow of the whole am)c  3.  ob+e「マafion datac  fx (received   vector)、  n+x (n  lee  vector)c  +x p (+eceived+  vectorp「ime)、  nTxpc   4.  match fille+ lh+e+holdingc   「xp m (lee  vector  prime  matched)、n+xp m(n  lee  rector  pr奄高■@ matched) c   lhl (lh+es!+nld  l)c   d+xpm(det ected  tec  prime  matched)。
nd+xpm(n  detected−rec  prime  match ed)c   bkr (block  +angel、 nbk+ (n   block −range)c  5. 1easl 5quare lit  1c  blsl (b Is 11. nblsl (n b  Is I) c   Ih1sl (lt++eshold  Is  I)、 dlsl  (detected  Is  I)c   tpls (lp  Is、   n1pls (n  lp  ls)(use bkr、  then gen erate new bkr Itom lh+esholding blsl )c  5. 1east +quare lit 2c   bls2 (b   Is  2)、 nbls2 (n  b  Is  2)c   1hl i2 (threshold  1lI)、 dl+2 (detected   l+  2)C(use bkr) C***傘傘*本**零tネ零本ネネネ傘傘傘**本**本拳京本*ネネUネ傘 *零***傘孝本*零孝零ネ孝*本本本本ネ*零ネ*サブルーチン+1del  (sl、 nsl、 null)C#本*ネ本ネU本I*孝孝零零零**孝零# #*零零傘*傘傘本**本孝本傘傘*本零車零ネ傘零**零零本*傘傘ネ*** 本ネネc  5ldel  (s   ref    defined):   define  「eference  vectorc +I −s  te letence vectorc nsl −no、  ol ’none t eto’  enf+ie+ in slcomplex slm integer nsl、  n5llc some 1nle+m1dial e signal +ele+udc (not Ycl implement ed)s1m=cmplx(−1,8649517,−3,6012861)s H2)zcfflplx(−1,2861736,−5,6591640)sl (31= cmplxf−0,77170418,−11,897106)tl (4) =cmplx(−6,4951768,−20,482315)sl( 51:cmplx(−19,099679,−25,241158)+I(61 ・cmplx(−32,411575,−20,7395491s+(7) □  cmplx(−39,163987,−Il、 189711)sl(8)  = cmplx(−39,742765,−5,1125401)sl(9)  = cmplx(−38,906753,−3,0546624)「eturn nd Cネ本**零零本#零傘傘本零傘零ネ参傘ネ本ネネ拳零本孝ネ本本本ネ孝章ネ本 本*ネ本本ネ牟本**ネ***牟*ネ参ネ傘*$*ネ*サブルーチンamdel  (am、 nam、  iam)c  rttt$tHt$ttt孝零#I孝 零傘傘傘孝孝傘孝零孝孝傘孝孝孝孝零零孝傘零を零季孝零孝零nnmnn傘孝零 c  aa+  (Iir)  definitionCII (a  mil +ix) −the ’none tcto’ put ol s vecto r Of 1  m5trixc  nzm  (n   s    matr ix)  −no、  ol  ’none  xeto’  amc iam  (i  a  mal+ixl −5election of dille+ ent Amcomplex  am(1) intege「nam、  目m il  (iam、B、])  thenelseif  lxa、eq、2]   then請1)  =  −0,5 1se paint t、’ enle+ 1lalll ?’「ead (54) n xm print 零、’enter un(naml ?’read (5,$)  (am(i)、  i=l、 namlndi1 nd C本市傘零参傘#傘本字本**零零零*零零零$零零零零孝零零孝孝零ネ零零孝 孝孝傘ネ*tuunu傘傘傘を中傘傘季季ネ孝サブルーチンamll+n (x 、 nx、  am、 nam、  xp、nxpIC零*傘傘本***ネネ* 孝*孝傘傘孝傘*傘傘**傘章傘ネ傘傘傘傘本零**電傘傘零ネネネ*傘零傘傘 傘孝零本零零零傘傘傘零本c amllm (a matrix l+ansl o+m) to L and generate xpc  x  −1npu t  vectorCIII  −dimension  01  Kc  a m  −a   matrix  vecto+cnam−n  a  i+a frixno、 alnonezeroa  mjlr:xvectorcom plex x(1)、  am(1)、  xp(1)integer nx、   nim、  nxpnxp = nx −naL11+ 1do i =  I、  nxp xp(i)・ 0゜ do j = l、 nxm n(i)  ・ xp(i)  +  xm(iDx(i+i−1)endd。
endd。
eturn nd C*傘零*傘ネ**ネ傘ネネ零拳傘***零ネ傘*ネ**零零傘ネ*ネ傘ネ傘零 孝孝#傘*ネ*傘傘孝傘ネ**傘零*傘孝ネ傘傘孝零零すブルーチンlhc   norm (xlnx、  th、  II、 IIIK)C#本本本ネ本章零 零I傘tネ零傘傘傘*零$$零零零章傘I*本*ネ*本零*幸ネ*零零本傘本* ネ*孝ネ*ネ*零傘孝傘零孝ネ*ネc Ibresholdinglbe no tm of x、  put detected IOc!口on in lx C百(x(i)tconjg(x(i)l > th、 It(k) = iC X −1nput vecto+ c nx −dimension 01 !c th −th「eshold c lx −detected K 1ocalion+c nlx −110 011ocalions detectedcomplex Im integer IIX、  1x(11,nlx百(cabs (1(i)) 、 gf、 thl thennlx = nlx + 1 1x(nlx)  = i ndil endd。
eturn en+1 il (nlx、ge、2) thenC#***零**本零零零傘*#I#零 本本零本傘ネ傘幸**車本傘傘*零傘傘傘**本本零傘傘***傘市**零零零 零本本孝***サブルーチンbk+ld (bk+、 nbk+、  II、  nlx、 n5lp)C傘傘傘零***零車掌零傘零傘傘傘傘章***傘傘ネ傘 中傘**零零零ネ零傘*傘**ネネ傘傘ネ*参傘零参拳本本零章零零傘本ネ零c bkrl+I (block  range  1inde+)c find  blocks f「om detected positionc if b( i)−lx(j) < nam、  eve[ything in betwe en 1x(i) and 1x(j)(belong  to  the   5une  blockc note  目tea+l−+qua+e−1it  does not con+ide+ +xpm that not in  lxc     it 5hould use +xpm and II、   di「ectly。
c bk+ (block−「ange) −integer array c onlain+ block bounda「ie@c                        in  two  consecutive  el ementsc nbk+ (n  block  +angel −no、   of blocksCII (location  x) −detecte d 1ocationsc nlx (n  1ocation  X) −n o of 1ocajion+ detectedc n1lp In  !   +el  prime) −minimum 5eparation +an ge to be isola狽■■ intege+ bk+(1)、  nbk+、  1x(11,nlx、   n5lpi1p−11(11 do i = 2.  III! clsll (lease  5quare  lit  11nbls =  nbls + 1 1blt(nblsl  = bk+fi傘2jbl+(nbls) = +x pm(bk+(目2))#5lp(1)1Ie do !l :1.  nbk do I2 :I、 nbk if (iabill−I2)、ge、nclp) thencovlp(il +(I2−1)本nbk)  = 0゜elseil (il、ge、I2)  thencovlp(jl+(I2−1)ネnbk) = +5lp(if−i 2+I)1se covlpfil+12−1)掌nbkl =conjg(+5lp(I2−j l+l))ndil endd。
endd。
e  5hould  be  replaced  b7  +outine s  th!+  use  cond 口inn  noA c 「1ter ion call m1vl  (covlp、  nbk、  covtp、  di p)do il ” 1. nbk bls(nnls+il) ;0゜ lbl+(nbl++ill  = bk+(i傘2−1)  +  if   −1do I2 ;l、 nbk bls(nbl山11 ;bls(nbls+il)&     + covt p(il+(I2−1)傘nbk)ネ+xpm (bk+ fi$2−11 + 12−1)endd。
endd。
nbls = nbls + nb& ndil endd。
Te1llln nd C#ネ*本*傘傘ネ**傘孝ネ*st*uu傘傘傘**傘傘孝傘孝ネ*ネネ*ネ ネ本ネ**ネネ傘参ネ孝孝ネ本傘傘ネ*ネ**零本サブルーチン1hrhdl  (1,lr、 nx、  th、 dlx、  ndlx)C傘**ネ傘**傘 **零孝*章傘傘傘孝*孝孝**章傘*ネ傘ネネ傘ネ傘孝*II#ネ*#傘傘掌 *ネ傘傘傘傘*孝ネネ傘孝ネ孝ネ*孝c lh+esholding L  p ut detected 1ocation in lxCil  conig (x(i))傘(x(i))  >  lh、  dlx(k)  = 1x( ilCI  −1nput  vecto+CII (locilion  I ) c nx −dimPnsinn 01 L  Itc  lh  −thre shold c dlx fdelCcl!d  1ocxlion  x) −detec ted ! Iocxlionsc nlx −no ol 1ocation s detectedcomplex x(1) integer 1x(1)、 nx、 dlx(1)、 ndlx+cil   th ndlx = Q do i = l、  nx if (+exl(conjg(x(i))孝x(i))、 gl、 thl  thenndlx = ndlx + 1 dlx(ndLr) = 1x(i) ndil end+i。
eturn nd c I*傘傘*傘ネ傘本市本市**ネ零**ネ本章零章市ネ*$傘傘傘*傘零章 零傘傘市ネ零孝ネ傘傘傘*傘傘傘傘傘傘市$ネ傘ネ傘孝*サブルーチンfchd el(Ich、  1chlC$@$*ネ孝傘傘#傘本傘本本本零本$零本本本 ネ孝孝*本本本傘零本零傘#章本本本*I*ネ本本傘ネネ*零ネ市ネネ傘*孝ネ 孝本c Ichdel (ieq channel  delin山o++)c     (no+malized v、 t、 1. coior 5ubca ++ie+)c  lch  (1+cq   channel)c  ich   (i   channel)  −channel  n。
integer ich 百(ich、 eq、 I) thenIch = 0゜ 1se Ich ・(43,25+ in6)/3.579545ndil eturn nd Cネ牟*傘#***掌I**孝ネネ傘ネネ傘零傘傘ネネネ傘傘孝傘零ネ孝孝零参 孝傘傘**ネ傘ネネネ傘傘参傘傘傘ネ零零零零傘零*ネネ傘l本孝参本拳 サブルーチンm  mulch  sap (rxpOIIn+xpm、  u p、 nrxp、 slp、 n5lp)C*傘*零零傘*孝#*傘傘零傘ネ* 零**ネ零傘孝傘傘傘傘*ネネ傘*傘ネ傘ネ$零孝*孝**ネ傘傘傘**ネ零傘 零零傘零*ネネ*$孝P傘ネ傘ネ傘 c txp  match−山to generate +xpm (rece ived X  prime  avtched)c +xpm (recei ved  x  prim!matchedlcn+xpm (n recei ved x  prime  matched)c ltp (+eceive d  I  prime)cn+xp (n  +eceived  I  p rime)c ltp (s  +el  prIme)c n5lp (n   +el  p「ime)complex +xpm(1)、  +xpm、   1lp(1)integeIn+xpm、  nrxp、  n5lpnlt pm = nrxp −n5lp + 1do i = 1.  n+xpm +xpmfi) :O do 1= I、  n5lp +xpm(i)  =  +xpm(i)  +  IU(山−1)本5lp( ilendd。
endd。
「eturn nd C#傘零本傘本*傘*傘中本本ネ*京*零*I*本t#本**ネ*I零章傘** *を章**孝*傘ネネ本孝*本ネ*本章**孝傘傘本本サブルーチン+xdel  (+L ntx、  lnlC零**零*掌*牢**零を傘tI*傘*本** 零章零零傘本零零零傘傘本***傘傘傘傘***零ネ本*本*ネネ本*本傘本本 本傘*c  udel  (「eceived   vector  deli n自1on)c read in II from file Inc )il e lormaf at txdet、dat(ls+ data c 2nd data c +x (+ecieved  vectorlc ntx (n  +ec ieved  vector)cha+acler*2tl In complex +x(1) Integer ntx open  (4,l1lt=In、5tatus=’ old’)read  (4,杓ntx do i = I、  ll11 「ead (4,1xr、  xi +x(i) = cmplx(xr、 xilendd。
close[4) eturn nd cl零章**零ネ傘章**零章市零*傘零本章****零零零零傘零**す零掌 ネ*車キ孝零零傘孝掌零孝孝5unup傘傘零サブルーチンplol  vec tor (px、 py、 nplCネネ傘***孝傘傘孝牢傘**零本*牢ネ 零*市孝傘傘**孝*本傘キネ傘*傘市傘傘零孝ネ零零傘傘傘傘孝零本本本本零 零傘傘傘本c plot (px、 p7) inlege+ numbe+  paintcp裏fplot  X) a…マ CpT (plot  xi arrayinteger px(1)、  p y(1)目+tege+  np call movai+5(pi(1)、 py(1))doi=2.np call drwabs(px(i)、 ay(il)endd。
eturn ND C*零本零**#キ**零牢ネ*ネ零*本ネ*本本本本*本**n*零零章ネ傘 本本ネ本ネ本傘*傘本本*傘本*本本*中本**本ネ*ネ*{ すブルーチンvl  pit  abs  c (vt、 nvl、  txo ll、  17off)C*ネ**拳零を本傘傘傘**傘傘傘本章***傘零本 零*傘牢ネ*孝傘零ネ零傘*市*傘市ネ本傘**孝零零*ネネ*ネ零零零零傘零 ネ零P* c plot abs value of complex vector r t with (ixoll、 1yoll) ollse■ c by using plotlO「oulinesc nvj −dime nsion of vector to be plottedcomplu  yl fl) integer nvl、  1xoll、  i7o11intege+   ix、  i7 1x −1xoll + 1 iy  −1yoll  + 1nt(cabi(vt(1)))call m ovabs  (ix、 17)call pnlabs (ir、 iy)百 (nvl、 ge、 2) thendo i = 2.  nvl ix = 1xoll + i iy −17o11 + 1nl(cabi(vf(i)))call drw abs  (ix、 1ylendd。
ndit 「eturn nd C傘傘ttom零ネI零ネ零*京*傘零ネI零本傘零章章傘*傘傘本*傘ネ孝傘 孝ネネ傘章零***本市*零ネ零傘**ネネサブルーチンvl  pHc (v l、 nvl、  ixo[1,1yoltlc **章本季**零傘本本零章 零本本*市零傘t*零零零*ネネ#傘*零ネ本傘本傘ネ車*本*本章*本零傘零 零宰*傘零本本零*c plot complex vector vf wi th (ixofl、 1yoll) oll+elc by using p loflOroutineSCvl  −vecto[ c nvl −dimension of veclo+ to beplot tedcomplex vj(1) intege「 nvl、   1xofl、   1yol[intege「 it、  1y ix = 1xolf + 1 iy  =  1yoll  +  int(+eal(vlml)p+int ?、’  i、ix、1y call movabs (ix、 1y)call pnjaba (it、  1y)if  (nvL、ge、2)  thendoi=2.nマ( ix = 1xoil + 1 iy= 1yall + 1nl(「eal(vl(i)))call d+w ab+ (ix、  1y)c+++Id。
en+1il ix = 1xot[+ 1 il =  1yoll  +  int (aimag(vf (1)))c all mowabs (ix、 171call pnlabs (ix、  1y)if  (nvf、ge、21  thendo i = 2.  nv l ix =  1xoll + i i7 = 1yoH+inl(aimag(vl(i)))call drwa bs (ix、  1y)endd。
ndil eturn nd C幸*傘市書オ傘*掌本*章*幸*孝本***幸*零零*孝零**孝*本本零傘 *ネ傘*孝傘傘零傘傘傘*本傘零ネ傘傘孝傘*傘ネサブルーチン+c+n  r ead (it、iy、  1xoll、  1YallC本##零tt$tH t*象U傘零零ネネ傘孝*本*ネ零零を京章傘ネネ零傘零零車傘ネ本零零ネ傘零 宰*本傘孝零*京傘傘零c read 5creen coordinate  f「om plot on  the 5creen b7 using pl ol撃n +outines (vbere (ixoll、 1yoll is 5creen plot  O山e(c ix、  iy −coo「dinate +elutnedin teger ix、  i7 integerixoll、  1yolt200  call scu+s+  (icha+、  ix、  1Y)call movabs (ixoll 、 700)call  anmode ix = it −1xoll iy = iy −1yol+ Wr口e (6,$)  ix、  1yprint t、’ tead an olhe+ coordinate、  Ilo ?’read  (5,ネ)    1ansi[(ian+、ne、0)  golo 200+elu+n nd c I#零零零傘掌を章車零*ネ#**I零傘t11本*零傘車*中零零傘車ネ 傘**零傘ネ傘本*tl傘*零零本寥零傘*傘サブルーチン+c+n  lhd  (th、  1xol[、1yofl)C傘*I零零#掌*を車$孝中*傘* ネ零京傘傘傘**本章***本ネ本本ネ本零宰ネ傘*本零零傘本*本本拳零本傘 本本傘本本本c find threshold loom plot on  the 5cteen by using plollO+ouline■ c where (ixoll、 1yoll) is +c+cen plo t oll+elc th −1h+eshold intege「1xofl、  1yollin+eger  ix、  1j 200  call +cu+s+ (icha+、  it、  1y)ca ll  movabs(ixoll、  )00)call  anmode write (6,$)  1x−ixoll、  1y−iyollpunt   章、’change  liteihold、  Ilo  ?’t e  a d  (5,章)  1ans百(iana、 ne、 0) goto  200th ” float(iy−iyolf)+elu+n +nd c  ****本*ネ零車零***中**零**#零零傘傘**本零零**** 掌本****零本傘本*傘零ネ本**ネ***市*ネ本ネ*Pネ零孝傘*** サブルーチンsl  mall (+fm、n51m1.n51m2.ledg e、nledge、sl、n5Lns目)C零**傘傘*#本ネ****I*本 *傘零市章車本傘本本*ネtネ**ネ本*傘*傘****零ネ傘*零*傘ネ傘* ネ本*零零*孝掌零傘{本傘ネネ掌 c generate +I  matrix loom sl and ed ge 1ocationsc ns[ml、  n51m2 (+Iimens ion ol slm)c   n51m1− lit dimension  of +Im、  length of each edge int。
c   n+1m2−2nd dimension of slm、 no、   of edgesc ledgeQ  (edge 1ocations)c  nledge (n  edges)c nclしn s  +et (ed ge 1ocation in s[)complex slm(1)、  a ll)intege「n51m1. n51m2. leage(1)、 nl edge、 n+1ns1m2−nledge do i = I、 nledge nfau = ledge(1) 百((nlau、 le、 nc印and、 (nslml、 le、 (ns l−nsll+n1au))) thendo i :1.  n51m1 s 1m (j+ (i−1)傘n+1m1)  = ill+n5ll−nt au)endd。
山eil ((nlau、 le、 n1ll)、 and、 (nilml、  gl、 (nsl−nsll+山u))l thendo i−1,(nsl −nill+n+aul+1m(i+(i−1)車n+1m1) =if(j+ n5ll−ntaulendd。
do J= (nsl−n+II+ntau+l)、 n+1m1s In ( i + (i−1)傘n51m1)  =  5l(nsl)endd。
elseil ((nlau、 gl、 n5ll)、 and、 (nslm l、 le、 (nsl−n+lI+n1au))) th■■ do i = 1. (nlau−nsll)slm(j+(i−1)零n51 m1)  =  +1)endd。
’do 1−(nlau−nsll+l)、 n+1m1S師(i+(i−1) $n+1m1l :5l(i−nlau+n5ll)endd。
山eil ((nlau、 gL n5fl)、 and、 (nslml、  gt、 (n+1−nsll+n1au))) thendo i =l、   (ntau−nsll)slm(i+(i−1)tnslml) □ 5ol) en+Id。
do J= (ntau−ncl+I)、 (nsl+ntau−nsll ) 11m(i+(j−1)tnslml)  ・5l(i−nlau+n百1)e ndd。
do j= (n+I+n1au−nsll+I)、  nslmlslm(i Ni−1)tnslml) =sl(nsl)endd。
ndil endd。
「eturn nd Cネ*#零*車章電市本*本*本*零ton零***本*市*本零市本本*本* *傘本ネ***本*零京本傘本本*本ネ参零本本本*傘**■{本*傘傘傘*傘 零 サブルーチンsf  mat  (slm、 n51m1. nctm2. l dsfm、 ledge、 nlle、sl。
c   nsl、 n5ll) C*………ネ零…………………傘傘……*本……*……本傘傘*傘傘本本#零* 1本本零本I零c generate sl  matrix iomsl a nd edge 1ocationsc +1m (s−「el  mat)c  n51m1. n17m2 (dimeuion at stm))c    nsfmI −1sl dimensIon of slm、  IeBht  of each edge int。
c   nsl+n2−2nd dimension of slm、  no 、  of edgesc ldtlm −leading dimensio n ofslm mat「ix stoIagec ledgeO(edge  1ocations)c nledge (−edges) c n5ll −n  s  tel (edge 1ocation in  sllcomplex +1m(ldalm、 IC+1mintegerns 1ml、 n51m2. ledge(1) 、 nledge、 n1ln+ +m2 = nledge do i = I、nledge nlau = ledge(i) if ((nlaa、 Ie、 n5l11. and、 (n+fm1. l e、 (n百−nslトn1auJ)l thendo i = l、  ns lml slm(i、 il = 5f(ifngll−nlau)endd。
elseil ((nlau、le、nsI)、 and、 (n+1ml、  gL (nsl−nsll+njau))) thendo j = I、   (nsI−nsll+ntau)slm(i、 il = sll+n5ll− ntau)endd。
do i −(nsl−nsll+n1au月)、  nslmlslml、  il  = sl (nsl)endd。
elseif  ((nlau、gt、n5ll)、and、(nil畦 le 、(niト1sll+n1au)))  1hendo i−1,(nlau− n+Il)+Imfi、 i) = 5ol endd。
do j= (nlau−nslHl、  n51m[slm(j、 i) =  ++(j−ntau+n5ll)endd。
elseif ((nlau、 gt、 n1f11. and、 (nsl畦 gt、 (nsl−nsll+n1tul)) thendo i = 1.   (ntau−nslllilm(j、 i) = ++(1) endd。
do i = (ntau−n+II+I)、  (nsf+njau−nsl l)sfm(i、 il = 5f(i−nlau+n5ll)endd。
da i = (nsl−nlaa−nsLl+1)、  n+1m1gfmQ 、 i) −5l(nil) endd。
ndil endd。
+el旧1 nd c t…………を京…章ネネ………………*……傘…………傘サブルーチンno   dc (x、 ldx、 nxl、 II!2)C…傘傘章……………傘… …ネ………傘…孝…零……ネ零…本ネc eliminate do olie l in each column of !c  xQ  −1nput   (nrl、nx21  matrix  會i1t+  leading  s torage  dime獅唐奄盾氏@ Idx c  Idr  −1eadiB  slragc  dimension   O+  1r  nxl、  nx2 −  づ1mtnsIon  ol   data  in  xcomplex x(ldx、 I) integeInxl、  nx2 complex dc  1p dc   lp  : dc    +p+x(j、1)endd。
da   lp  =  dc    lp/float (nil)do i  :1.  nxl x(j、 i) = x(i、 i) −dc  lp   ”endd。
endd。
+etu+n nd c *本章傘を本章**孝*零#傘*零tl零章零零**ネ零傘零*零孝*** 傘本零***零傘零孝孝傘孝ネ*傘章零零傘傘*本本**零P本本章傘零ネ傘本 零*傘 サブルーチンvl  plj  abs  cl (vl nvl、  m、1 xoll、  1YolllC#傘車ネ**零零零零零零零零章零*零傘零傘章 I傘零本零零本零零**零*零零*傘**零零ネ零傘傘ネ**本市零ネ*零*零 ネ**零■Pネ零零零零零本零零零 c plot ab+ value ol complex vecto+ v l with (ixoll、 17off) ollse■ c by using plollOroutinesc nvl −dime nsion Ql vecto[to be plottedc 1vtQ − 1ocalion of vlQ componentcomplex vlc I) integer nvl、  1xoH,1yoll、  lvt(1)int eger  ix、  i7 it−m(i!  + 1xall i7 = 1Yall call movabs (ix、 17)call pnlabs (ix、  iりix = 1vl(il + 1xol+i7 = 17off + 1 nl(yl(illconig(ylfi)))C!II d+vab+  ( ix、 iりendd。
「eturn C1章10傘傘傘傘傘傘I傘傘傘傘傘を傘傘傘傘傘章傘傘傘ネ市ネ傘傘**傘参 傘本傘*零市孝孝零*本零傘ネ本$零ネネネ孝***市零ネl傘零孝零傘傘零傘 零傘 サブルーチンpeak  sea「ch (pk、npk、p  sep、lr xpm、nlrlpm、「!pm、+oxpm、n5ll)C*孝零本本本#零 掌電傘傘京傘傘傘零本章#傘孝傘傘**傘市孝*ネ本傘書零ネ傘*$ネ*傘傘傘 傘*本傘ネ本ネ傘零掌*本ネ*零拳本ネl本零傘$ネ本本ネ*零孝 c sea+xh peak in abs(rxpmQ)  for 1oc ations 1ocated b7 IrxpmQc vhen 5epa「 ation (lrxpm(il−1xpm(i−1) >=(p  5ep) 。
c    5cpuale peik 5hould be conside+ cdCpko −pea目ocalion 1ocatedc npx −nu mbe「of peaksc p  sep −min、 5epa+alio n between山山o1 peak considerationc mp mQ −1ocation rx prime  mx+ched where  1bs(+xpa+0) > threshol■ c nl+xpm −number ol 1ocations in l+x pmOc rxpmQ −t  vaeclo+  prime  match ede nrxpm −data dimension of +xpminl ege+ pk(1)、 npk、 p  sep、  lrxpm(1)、  nl+xpm、 n+xpmcnmplex +xpmfl) integer pk  1p pk  jp = mpIl(i) do i =2.  nl+xpm i1 ((ixpm(i)山xpm(i−1))、ge、p 5epl the nIIpk = npk + 1 pk(npkl = pk  tp pk  jp = lrxpm(i) 山C目(cabs(+xpm(lrxpm(i)))、gt、ch(rxpH1 (pk 1p))) thenpk  lp :lrxpm(i) ndif end+i。
npk = npk + 1 pi(npk)・pk  lp mnsme peek locglions to the 1ocation s in +10do i = I、  npk pk(i) = PK(il+n5ll −Iendd。
+e1口rn nd C零傘…零…傘*傘掌傘傘*…傘…傘…………………傘傘ネネ傘車ネ……傘ネ… ネ$ネ…掌サブルーチンlsl  no  dc la  II、 na  I s、  +fa  hsla n51m1. n5ln2゜&                          Il、  nrxl、   cond、    error)C……孝傘傘零…ネネ幸…傘…………傘傘…傘……零……本… 零*零章零………零c  find  1east  sq■e  5olut ion  a   Is  01  tx  =  slm  t  ac a   Is  −1east 5quoe 5olutionc nm  Is  −dimension of 5olution in a  Is wecl orc slm   −reference signal mgHixc l dgla+  −lerding dimension of sto「age c  nslmI    n51m2 − 1sl  gnd  2nd  d imension  ol  dtli  in  slmc rx   −o b@erved dalic nrxl −slanting positio n of dalx in rx Io+ the legsl 1IIoir e f奄■ e cond  −condHion number loom csvdc  rouline。
c      il cond = 0.  (info =/ 0)、  i ll conditioned。
c mor −nore+ 5quare ol the euor yecl or of 1east 5qosre 1ftc I  eliminate  dc in slm and +xc2sla+=u*sv*v’ c3   a   l5=v*5y(−1)  零u’$+xlpc  4     error  =  noon((I  −UlネDI’)  *   t x   jp)complexa  l5(1)、 slm(ld+Im、I) 、  rx(1)inlegcr na  Is、  ldslm、  n51 m1.n51m2.  nnl「eal cond、  errot c nole: nsfm2 < 50.  n51m1 < 100con+ 山1 y  1p(50,511)、 wo「k(9101,5v(100)  、e(]θ0)、 rx  tp (910,I)integer ldv   lp、  job、  int。
complcx cfp ldv  lp =50 do i = 1. ufml II  1p(i、I) =+x(i+n+xl−l)cndd。
call no  dc (u  Ip、10(1,++sfm1.l1cal l csydc (slm、 Idslm、 n+1ml、 n51m2、sv 、 e、 +Im、ldi1m。
&            v    lp、  Idv    lp、  w o「に、 job、  1nfol百(inlo、 me、 O) thenp rint傘、’1llcond口1onedsin’cond = O 1se cond = cabs (sv (Ill/cabs (KV (ns1m2 11ndil c note: when 1nto、 ne、 0.  the main  p「ogram Ihould be 1nfo「med。
do ! :1. n+l+n2 vorkl) = 0゜ dok=1.n百m1 wo+k(i) :wo+k(i) + conig(slm(k、 i))$ +x  1p(k、 1)endd。
endd。
do i ・l、  n+1m2 a  l5(i)=a  l5(1+v  tp(i、jl*wo+k(i)/ 5y(i)endd。
endd。
c calculate no「m +qua+e of et+ordo i  = 1.  ns1m2 wo+k(i) = O do i = I、  ns1m1 wo+kii) = wo+kii) + conig(slm(i、 り)$ +x  1p(i、 1)endd。
endd。
do i ;1.  nsfml cpl = Q do i =1. ’n51m2 clp = ctp + +fm(i、 i)傘watkmendd。
clp = rx  1p(i÷l) −clpe+tor= e++o+ +  cabs(cfp孝conig(ctpl)endd。
retu「n nd c  tttt吹電孝す傘孝孝孝キ零傘傘ネネネ車車車拳傘ネ事ネネ傘ネ車車零 京t11零寥参掌*孝傘零傘零零傘牟傘孝**零本零零本本■■■鮪P*車ネ サブルーチンvl  pH2(vl、 v2. nv、  ixl、  iYl 、  it2. 1dy)c  $H本#本本零市事掌本本本本専零本本傘本本 本本卑キ*本*本本*キ本孝*本*を京市**零零零零零孝零*零ネネ傘章幸零 傘ネ拳F零零零傘孝傘傘孝 c plot 2 veclots using plollfl routi nesc vlO−「eal vector 1c v2Q −real ve clor 2c nv −numberof elements to be  piO口ed in each veciotc ixl −x−axis o riginc iyl −y−uis origin lot yeclor  yl Oc it2− y−axis origin lot vector  v20c idy −mat amp mapping on y axis  from eilhe+ veciotc real ylfl)、  v2m integer nv、  ixl、  iyl、  it2. 1dydo  i = 2. 910 百(vl(i)、gl、amx t) thenamx  r = vl (i l elSeil (vlfi)、11.amn  i thenamn  t =  vl (i) ndil if (v2(i)、gl、amx  i) thenamx   i = v 2(il elseil (v2(il、Il、amn  i) thenamn  i  = vHi) ndil endd。
a6百−1=  amX    I−amn    Iadit  i = a mx  i−amn  1adit = amaxl (adit  t、 a dit  1)it = ixl iy  =  iyl  +  1nl(idy傘(vll)−amn     +)/adit)C!II movabs(ix、 iりdo i =2. 9 10 1x = i −1+ 1xl i7  =  i71  +  1nl(idyt(71(i)−av+    +)/adil)call d+vabs(it、 17)endd。
it = lI] 17=i72+1nl(idy傘(v2(I)−amn    +)/xd百) call movabs (it、 17)do i =2. 910 it=i−1+i夏1 iy=iy2+int(idyl(v2(i)−amn   +)/adit) call d+wxbs(ir、 17)endd。
「etu「n nd C#傘ネI*I零市市傘傘傘傘傘傘傘ネ傘傘傘傘傘傘$#参零零傘零零傘傘ネ傘 $傘傘傘孝ネ傘孝*ネ市零孝$孝ネ傘傘傘傘零ネネ傘ネ傘傘P傘ネネネ*零傘 すブルーチンvl  pHcabs (vl、 nv、  ixl、  i71 . 1drlC*傘#を章傘ネ零傘傘キ**本傘ネ*ネ傘#零ネ零傘**傘*ネ 傘零傘傘*零零京章ネ傘傘ネ*傘傘本本車ネ本傘ネ市ネ傘零毒ネ*傘*ネ■Q傘 車傘本ネ傘 c  plot  mBnHude  ol  a  complex  ve ctor  using  plollO+oofinesc  vl Q   −complex  vecto「 1c  nv  −number  ol   elci+enls  to  plotted  in  each   vectorc  ixl  −x−axis  0riginc  i71  − 7−axis  OIigin  fo「 vector  vl。
c 1dy−maxamp mapping on y axiscomple x  vl (1) integer  nv、   ixl、   i71.   idyamx     t  =  cabs (yl (1))amn  r = cabs  (vl (It)dll  i  −2,91Q ■(cabs(vl(i))、gl、amx I) then1+111 1  = cabs (yl (i))elseil (cabs(wlm)、lLa mn−「)the++amn    +=cxbs(vl(il)ndil endd。
adit = amx  +−xmn  +ix = 1xl i7 = iYl + 1nt(id7(+cIbs(vl(1))−amn   +)/xdiりC11l moyxbs(ix、 17)do i = 2.  910 ix = i + ixl −1 i7 = i71 + 1nl(idytcabs(yl(il)−amn   +)/xdillcsll  d+wabs(it、17)endd。
+elu+n nd C#ネ車ネ傘零Iネ傘*傘傘傘傘ネ傘傘*ネネネネネネネ零零零本*章ネ*傘# #傘傘傘ネネネ傘ネ孝傘孝*参$ネネ傘*ネ傘傘ネ傘牟傘ネl零零ネネ零傘傘 すブル−チンvt  pH−cabal (TI、 nv、  ixl、  i dx、  i71. 1dy)C寥ネ傘傘ネ傘*ネ傘#孝傘傘ネ零零It*孝ネ 傘*ネ本傘*零傘*孝*零傘ネ傘#本本零傘零章傘*ネ零ネ傘零$傘*傘$ネネ 零本零零傘P掌ネ*ネ本本ネ c plol magnitude of a complex vector  using plollQ routinesc vlQ −complex  vector 1c nv −number of elemenl+ to  bc plotted in each vectorc ixl −I − atis o「1g1nc idx −5tep on x−hisc iYl  −F −axis origin lot vectorvlOCid7−1 1111 xmp mapping on 7 axiscomplexマ1( 1) integer ay、  ixl、  i71.  id7amx  r =  cabs (vl (1))imn  t = cabs (vl (1)1 do i =2. 9111 if (cabs(vl(i))、gl、amx  +) thenamx   t =cabs (yl (i))else目(cxbs(vl(il)、H, am  r)thenamn  r :cabs (vl (i))ndil endd。
iy  =  i71  +  1nl(idy?cabsfv団1)−amn    +)/adit)call movabs(ix、 1y)do i = 2. 910 ix =ixl + (i−1)$1dxi7 = iyl + inl(id ytcabs(vl(i))−amn  +)/adit)call dtwa b+ (ix、 1y)endd。
「etu「n nd cHネネ*零市ネ*傘**#市章*傘傘傘本*ネ傘*ネ零傘傘ネ孝傘孝*零傘市 傘傘*本*傘傘傘*ネ零零零零本本本本本傘ネ章零零零本拳■l零ネ零零零 サすルーチン+xp  coal  slp (+xpm、 nrxpm、   +xp、 n+xp、  +Ip、 n5lp)C傘傘傘傘傘傘*傘傘傘傘傘零 零傘傘*本傘傘傘傘掌傘傘傘*#*電を傘ネ**傘傘傘傘宇傘掌零ネ傘*ネ*傘 零傘傘傘傘*ネ孝零本零零零■棊■■Q本本本 c txp  coal  slp to gene「ate +xpm (r eceived  I  prime  correlaf■р■ c +xpm (「eceived  X  prime  matched) c nrxpm (n  received  X  p「ime  matc hed)c tip (received  I  prirIle)cn+x p (n  received  x  prime)c slp (s +e l  p+imel   (to be conjugated when c ouelating)c n1lp (n  +el  prime)comp lex npm(1)、  tip(1)、  1lp(1)intege「n rxpm、  n+xp、  n+lpn+xpm = n+xp −n5lp   + 1do i = 1.  nopm Txpm(i)・0゜ do j= I、  n5lp +xpm(i) = npm(i) + tip(i+i−1)tconig( slpl))endd。
endd。
eturn nd c    fo「’          4−28−86ct傘零本零傘*傘$ ***零零孝本ネ孝傘傘傘本ネ本I傘傘傘*傘本傘本本本傘*本ネ*掌**零本 傘傘掌*ネ零傘本本ネ*傘傘傘傘サブルーチンupdb+(x、 sb、 km x)CI零零市零中傘傘*****傘**傘拳を傘傘傘ネ*零傘*零傘章傘本ネ 本零零傘*本傘傘零傘本本本傘零本零零傘傘傘*傘零零**complex I u(1)、  by(1)integer口(])、 kkx、 ky(1)、  kkycomplex tp (p:x 百(kky、 ge、 I) thendOj・I、  kky ip = lp −by(i)*5bfky(i))endd。
e口dif 百 (kx fll、 eq、 0)  theny  : hx(1)  市  (p 1se y =hx(1) 電+b(kg(Illndil 百fkkt gl、 l) then do i = 2.  kkx 7 = F + hx(i)+5bfkx(i))endd。
endd call updbl(Ip、  sb、  max(kxfkkx )、   ky(ktr))1eturn nd c #*本#傘#ネ傘傘傘中#*#$零傘傘*$傘零*零傘傘*孝**市傘傘* *傘傘傘零傘傘*ネ傘*零**零傘*孝$傘零零傘孝零傘零P傘*孝零零$ サブルーチンpk  a −hxy (tu kx。kkx、 h7.に7.k k7. pk、 l  Is、npk)cH孝ネ孝ネ傘傘$$H$傘傘##I* 零傘$孝傘本**ネ孝孝ネ傘ネ*ネ孝*傘傘***傘傘ネ孝零拳*零傘ネ孝ネ$ 本本***孝参傘本l*傘ネ c cxlcola+e direcl 1ntense digilxl t iller hxQ、  h70c  loom pkQ and a  I@ Odirectlyc pkO−peak 1ocations in ter ms of posilons in II  ysQc      (1+a nsl+omed loom abs(+xpmO) pexks back  to rx  wsO)c a  lsQ −1east−square−1i t 5olo口ons using vindoved II  ysQc       and xssumirlgedge 1ocations pkQc  npk −no、 ol pkQ c  −−−on  return  −−−−c hxQ −numerat ors coellicienlsc kxQ −dela7 of hxQc kローno、 of coellicienl+ in hxQc h7Q − denomenalo+ coefficientsc k7Q −delB  ol h7Qc kk7− no、  ol coellicienjs in  hyOcomplex hx(I)、 hy(1)、 a  Is(1)in teger  kxm、kkx、  k7(1)、  kky、  pk(11 ,npkhx(11□ l/a  l5(1) ロー=O il  (npk、gl、I)  thenkk7 = npk  −1 do i = ’l、  npk ky(i−11・pk(i) −pk(1)hy(i−l) ・a  l5(i )/a  il(1)end+i。
ndil return nd C本II*零傘掌傘*ネ*本零零ネ*参零*傘傘傘ネ傘孝零傘傘ネ**#零零# 市本*ネ*傘傘*傘*孝孝***孝*傘本$ネ零傘ネ$本*l*ネ$傘傘ネ傘零 *ネ サブルーチンwavelor+a (iy!  +、ivf  i、 nl   yl、npx’ wl  bg、npx  plf。
&1xoff、  17o1f、  17scp)C#傘零傘ネネ*零傘傘ネ牟 *傘牟傘**掌章市ネ孝****傘零幸ネ傘**市ネ傘**傘ネ孝傘傘ネ傘本零 零*孝傘$$$零傘中本傘傘ネlネ傘ネ傘*ネ*** e plot video vawe+onn ol 1ine nl  vt  of vlQ tinted traIlnpx  vl@ bg c   ollset and npx  ob long on 5ctee n with (ixolf、 17o11) offse煤@jor Cjvl  I  O,(ixo目、 1Yall+17sep)  olls el lor ivl−10c   by uiing plollOrout inesc vl  +Q −「eal part ol video dat ac vl  iQ −iamg pail ol video dalgc  nl  rl −1ine lo be plovedcnpx  vl  b g−beginol pixel in l1nenl  vl fobepl olledc npx  pH−no、  of pixels to bc  plollcdc 17sep −5apualion between tw o plotsinlege+t2 iv[+ (910,11,ivl  i  (91(1,I)integer nl  vl、 npx  vl  bg 、 npx  pltinteger 1xoll、  17o11.  iy sepintegerix、  1y iX :1xoll + 1 iy=iyoll+17cep+jvl   +(npx   vl   bg 、nl    vl)call movabs (ix、 1ylcall p ntahs (it、 1Y)if (npx  plfge、2) the+ +do i :2.  npx  pltix = 1xofl +  i iy = 1Yall + 1yiep + ivt  +(npx  vl   bg+i−1,nl  vl)call d+vabs (ix、  1y) en+Id。
ndil ix = 1xoll + 1 iy = 1yoll + ivl  i (npx  vl  bg、 nl   vl)call  m0v2、t+t  fix、1y)call pnl abs  (it 1y)il (npx  plx4e、2)  thend o i =2.  npx  pltix =  1xoll +  i iy = 1yoll + jvl  i (npx  vl  bg+i−1 ,nl  vl)call dnab+  (ix、  1y)endd。
ndN +elu+n nd CI#零本傘本零零!#零本零I傘ネ**本**本I傘零*本**章ネ零京*本 **本#零零傘傘**章零ネ*本ネ零傘市*傘本本本本サすルーチンp&  l edge (ledge、 nledge、  n+rml、  nrr  v gl。
&              pk、  npk、  nwinl、  nw in2)C*孝零傘*傘*零傘ネ零ネ傘傘傘***ネ電傘傘傘**零零ネ傘** 零傘ネネ孝ネ傘零零孝本*傘傘零傘零*ネ*傘ネネ**ネ**本c calcu late edge 1ocations and+elatd pa+ame je+s from paek 1nto。
c pkQ −peak 1ocation+ in IX  vx□c np k −no、of peaksc oinl −no、  of pixels  to the fell of Ii…peak (included )c  nwin2− no、  of pixels to the 「ight  of 1ast peak (included )c  −−−−on  r eturn  −−−−c ledgeQ −1ocation of edg ei in observation vecto!(+x  vs obQ) c nledge −no、  of edgesc n+1ml −no、   of pixel in the observation windowc  ntx  vil −the li+st pixel 1ocation  of the obseryalion window 奄氏@II vsQ。
inleget Iedgefll、 nledge、 nstml、 ntx   vslinlcgc+ pk(1)、 npk、nwinl、 nwin2 nle+Ige =nph do i = 1.  nledge ledge(i> = nwinl + pk(i) −pk (1)en+I d。
nrr  vsl = pk(1) −nvinl÷ln51m1 = nwi nl + pk(npk) −pk(II + nwin2−1nd C傘傘I傘章傘t*H4$幸傘傘傘市傘ネ*零珍重本零零零零市零キ**ネネ零 零ネ傘零**零零拳ネ傘ネ零孝*ネネ零零ネ*零零零C零傘零章零章章章零傘* 傘本本#傘*$傘傘傘宰I章**傘市本市傘*本ネ*傘本ネ孝**本傘零市ネ* 孝傘零*電*傘ネ零掌零零零c lh+esholding the abs  of xQ 、  put delecled 1ocation in 11 c it abs(x(i)) ) lh、 1x(k) = iCI −1n put veclo+ CIIX −dimension oflc  Ib −threshold c lx −detected I 1ocationsc nlx −no  of 1ocation defecledcomplex x(1) integernx、  II(1)、  +111百(cabs(x(i)l 、 ge、 jh) thennlx = nlx + 1 1x(nlx)  コ1 nd1j endd。
+etu+n nd c  $0本*傘##傘!車零季ネを本本零本本傘傘傘***傘本傘*孝傘ネ* 傘市傘*ネ零傘孝零傘零ネ**傘零*孝孝孝****傘零サブルーチンpaus e C章傘傘*零傘零掌*零*tnum零傘傘ネ傘傘傘ネ零傘市本傘キ京*ネ市*傘 を市孝傘ネ*零参ロ傘零市傘傘ネネ傘ネ零零c use plollo cap torqai+7 routine to hold lhc executi on ol progr≠■ c  sc「een displa7 inlegc+ ix、  i! intege「1cha+ call scu+s+(icha+、 ix、 17)Ntun+ nd C孝本傘*零傘零本傘*傘市**傘本ネ傘ネ傘本*傘*ネ市傘*傘本*章ネ傘* 本零孝本本零*ネ傘車ネ零零傘傘ネ本*ネ電零ネ孝傘零$本{本本本 サブルーチンde  ltI (d、a  lt、 slm、ldtfm、n5 1m1.n51m2.  +x、n+xl)Cネ孝傘ネネ傘ネ本*傘市市零傘傘 I零参傘傘#章*ネ傘傘ネ拳傘傘ネ孝本*傘電ネ傘傘ネネ傘傘ネ傘傘本本本零零 零ネ本傘ネネネネ*零傘l零孝ネ c  find  dc  ollselzlQ  =+xQ−tfm  孝  a   lsQ      cl ° DSRCnote  Opage  6 g−75c      T cdlO傘 dO/ n51m1 c a  IsQ −1east 5quare 5olulionc na   Is −dimension of 5olution in a  Is  vect+uc +ImQ −reference Signal matri xc 1dslIIl−leading dimension ol slm  storagec nsfml、  nt1m2−1sl and 2nd d imension of data in slmc +xQ −observ ed datac n+xl −slanting position ol  data in II for the 1east squgre li■ complex d complex a  lt(1)、  slm(ldslm、l)、  r[ I)integer na  Is、Idslm、 nsfml、 n51m2 . n+xldo j= 1. n51m1 d = d + rx(i+n+xl−1)do k = I、  n51m2 d=d−s師(i、k) s a  lt(klendd。
end+In d = d / n5bnl IellHIl nd C*****傘傘傘傘傘nm傘孝傘中零傘*傘中傘傘傘傘傘傘ネ零孝傘傘零ネ孝 ネ孝傘傘零孝傘傘孝傘傘傘傘ネ孝零傘傘傘牟傘零ネ本零傘傘■ サブルーチンs+defl (sl、 nsf、 n5ll、 d2. d3) c *#孝章傘傘###I傘章零傘ネ傘章傘*孝ネ孝孝傘電孝*tネ京傘傘*宰 零*零傘ネ*ネ本零孝車傘零傘傘章傘*ネ$傘本ネ*ネネ本■麻l*ネ c  +1del  (s   +cl    dtlined):  def ine  +elerence  vectorc  II  −s  ref e「ence  vectorc  nsl  −no、  of  ’non e  xeto’  enl+ies  in  slc  nsl  −ce ntu  ol  edge  in  +IQc d2. d3− rang e of composite video to be mapped on t。
c      teler  to note Opage 68−75 fo rnotation and deriyafioncomplex  +1m integer  nsl、   n5ll「eal  d2.  d3 real d6. d7. d8. dO,gc some 1nle+m1d iale signal reIerredc (not 7el imple mcnledlslm = cmplx (−1,8649517,−3,60 1286り5lf2) =cmplx (刊、2861736. −5.659 1640)sl(3) = cmplx (−0,77170418,−If、  897106)II(4) =Cl1lll11 (−6,4951768, −20,4823151s+(5) −cmplx (−19,099679, −25,241158)sl(6) :cmplx (−32,411575, −20,739549)+1(7) =cmplx (−39,163987, −Il、 189711)sl(8) =cmplx (−39,742765 ,−5,1125401)sf(9) = cmplx (−38,90675 3,−3,0546624)d6  =  (40,零d3+100.  傘d 2)/+ 40゜d? = reall(II) i! (real (xl(il)、 gl、 d7) d7 □ real  (sl(i))if (real (sl(i))、 It、 d8) d8  = 「eal (st(i))endd。
d9・(+12*d7− d6td8)/(d7−d8)g・1゜ if  (d?、ne、0.)  theng・(d6−d91/d7 elseil (d8.ne、0.)  theng・(t12−d9)/d8 ndit sl(i) = g$5l(i) + d9tcmpl[、、1,)en+Id 。
C1ton en+1 C本#零ネI零#傘本本本本本本零m本I**零零t1本本*零零傘章零傘本零 零零傘掌章***零本**零本本零***本本**本孝零零Fネ零零本本本ネ零 サブルーチンsg  del (+g、  nhg、  Iaug、 n55g 、  Ifd、 n5ld、 n5ld1. n5sld)c *章章傘*零零 *m傘*傘**零ネ零*中e<車重um*st**n**<*mネネネ孝ネ零拳 *tmum*nutuu*nc generate a ghoiling v ertical 5ync reference loom 5tdOc  w ith edge 1ocated at laugc +gQ −ghost  vecto「Cn1g −n   Ig c laug −edge 1ocation o[sgQ in unit  ol I/Icc n55g −sampling rate in unit  of Icc 5ldQ −++  delailQc n5ld −n   5ld c n5ldl −edge 1ocation in Ifd in uni t of samplec n5ild−n  +ample  tfd in  unit of fccomplex sg(1)、  5ld(1)「ea l taug integernag、  n55g、  n1ld、  n5ldl、  n 5sldn  ++ = n5ldl −n1nl  (laug$1loal  (nss)d))j  =  n   ++  +  (i−1)  本 n 5ifd / n51g百(i、 It、 ]) then 5g(i) −gld(1) elseil (i、gl、n5ld) thensg(il−Ifd(n+I dl lse !g(i)  = Ifd()) ndil nddO Telu+n nd ct市*零零傘*掌傘を零傘傘傘本傘孝本本傘*零零零傘零章傘傘**ネ市*傘 傘傘傘本**零零零傘傘傘*傘*本零零傘零零零本本孝*傘l傘本孝*本ネ傘本 本孝零 サブルーチンgldde((Ifd、 n5ld、 n5fdl、 n5sld lC#零零HtHtHt$*HH$Ht孝ネネm零傘を孝*零傘零京傘*参** 京nanoネ零本零零零零零零零ネ孝傘傘*ネ*零参市京零本c 5fddel  (s  tefdetail  deHned): define refe rence rectorc sld   −s  「eference ve cto「c n5ld −n  Ifd c n+dll −edge 1ocation in ald。
c nold −Sampling +ate in unit of tc+ omplex ildfl) intege「n5ld、  n5ldl、  n5sldn complex vs(256) nssld = 64 vs(64)・cmplr (−]、 4791.  −3.1833 )vs  (65) ・cmplx (−1,5434,−3,2476)vs (66 ) ・cmplx (−1,6720,−3,2797)vs(67) = c mplx  (〜1.6720.  −3.3762 )vi (68) □  cmplx (−1,7363,−3,40841vs (69) □ cmp lx (−1,7363,−3,44051vs (70)  = cmplx   (−1,8006,−3,5048)vi (711−cmplx (−1 ,8006,−3,5370)vs(72) = cmplx (−1,116 50,−3,60+3 )vi (73)  :cmplx  (−1,865 0,−3,6334)vs(74) = cmplx  (−1,865[1, −3,6977)vs (75) −cmplx (−1,8650,−3,7 299)vs (761:cmplx (−1,865(1,−3,82641 vs (77) = cmplx (−1,8650,−4,0193)c    slm ・cmplx(−1,8649517,−3,6012861)vs  (78) =cmplx (−1,8650,−4,1479)vs (79 ) =cmplx (−1,8006,−4,2122)vs (80) =c mplx (−1,8006,−4,30871vs (81) ・cmplx  (−1,8006,−4,3730)vi (82) ・cmplx (−1 ,8006,−4,5016)νh (83) :cmplx (−1,736 3,−4,598] )vl(84) = cmplx  (−1゜6077、   −4.8232 )vs(85)・cmplx (−1,4791,−4, 9518)vs(86)□cmplx(−1,4148,−5,1447)vs  (87) ・cmplx (−1,4148,−5,4341)vs (88 ) ・cmplx (−1,2862,−5,65921vs (89) :c mplx (−1,2862,−5,9486)vi (90) ・cmplx  (−1,1576、−6,1415)vs (91) :cmplx (−1 ,0289,−6,5273)νK (92) −cmplx (−0,900 3,−6,7524)vs (931= cmplr  (−(1,8360, −6,97751vs (94) =CIIIII11 F  −tl、 64 31.  −7.3312 )vh (951=cmplx (−0,6431 ,−7,6527)vs(96) ・cmplx (−0,5788,−8,0 386)vs [971□ cmplx (−0,4502,−8,5531) vs (98) =cmplx (−0,4502,−8,9711)vs(9 9] = cmplx (−0,4502,−9,3891)vs (100)  ・cmplx (−0,4502−9,9035)vs (101) ;cm plx (−0,5145,−10,2251)vs(102):cmplx( −0,5788,−10,6431)vs(1G3):cmplx(−0,64 31,−11,3826)vs(104)=cmplx(−0,7717,−1 1,8971)vs(105’l=cmplx、[−0,9003,−12J7 941v+ (106) :cmplx (−1,0289,−12,7974 )vh (107) = cmplx (−1,0932,−13,2154) vi (1081・cmplx (−1,2219,−13,72991vs  (109) :cmplx (−1,3505,−14,5016)vs (1 10) =cmplx (−1,4791,−15,0482)vtflll) □cmplx(−1,8650,−15,7878)vs(112)・cmpl x(−2,1865,−16,1736)vs (113) =cmplx ( −2,6367、−16,6238)vs(114)・cmplx(−3,08 68,−17,299Q)vs(115]=cmplx(−3,6[113,− 17,6527)vh(116) =cmplx (−3,9871,−18, 2315)vs(117)=cmplx(−4,6945,−19,[1997 )vs(I181=cmplx(−5,0804,−19,5177)vxl1 9)□cmplx(−5,8521,−20,0000)vs (1201・c mplx (−6,4952,−20,4823)vs(121) □ cmp lx (−6,9453,−20,9646)vs (122) = cmpl x (−7,3955,−21,5113)vs(123)・cmplx (− 8,5531,−22,0579)vs (124) =cmplx (−9, 2605,−22,4759)vs (125) =cmplx (−9,96 78,−22,79741vs(126) =cmplx (−10,7395 ,−23,2797)v@(127) :cII+plx (−11,3183 ,−23,6013)vs(12g) :cmplx (−12,2186,− 23,9871)vs(129) =cmplz (−12,9260,−24 ,2765)vs(1301・Cl1191! (−13,8264,−24, 5[116)vh(131) =cmplx (−14,6624,−24,7 588)vs(1321□ cmplx (−15,3698,−24,919 6)vi(182) =cmplx (−39,8714−5,5627)of 183):cmplx(−39,7428−5,2733)vi(184)・c +nplx(−39,7428−5,1125)vi(185)・cmplx  (−39,6785−4,8553)vi (1861= cmplx (−3 9,6785,−4,5981)Ws(+871 ・cmplx (−39,5 498,−4,4051)vi(188) =cmplx (−39,4212 ,−4,3087)vt(189)・cmplx(−39,2926,−4,1 158)vs(190)□cmplx(−39,2926,−3,9550)v sf191)□cmplx(−39,2283,−3,7942)vs(192 Lcmplx(−39,l64Q、   −3,6656)9K (193)  ・cmplx (−39,0354,−3,5370)vi (194) =c mplx (−39,0354,−3,4727)vs (195) :cmp lx (−39,0354,−3,4084)vs(196) = cmplx   (−39,0354,−3,3441)vi(19了1  ・ cmplx   (−39,0354,−3,3119)vi(198) =cmplx   (−38,9711,−3,2154)vt(1991;cmplx  (−3 8,9711,−3,1511)vS(200)・cmplx(−38,911 68,−3,0547)v(f21)I)  = Cmplx  (−311, 842□1.  −3.0517 )vs (202) =cmplx (−3 8,8424,−2,9904)vs(203)  ・ cmplx(−38, 7781,−2,9260)Va (204+ ・cmplx (−38,77 81,−2,8617)vs (205) =cmplx (−38,7781 ,−2,8296)vs(2061= cmplx (−38,7781,−2 ,7331)vi(207)=cmplx(−38,8424,−2,7010 )vi (208) =cmplx (−38,8424,−2,6045)d o i = l、  n+ld +fd(i) ・vs(i+631 cndd。
「eturn nd C掌ネ零傘ネ傘傘傘本****市本***傘*章本傘傘孝**m傘本本本零傘* 零傘傘ネ傘*零章章ネ傘ネ零本傘傘*掌傘零章傘市*ネネ*サブルーチンled ge  laug (taug、 nlaug、n55g、 ledge、 n ledge)c  t$$$を章**電H?H4本*市掌零市本零*本本零市電 1*零傘傘零傘*市傘零傘**零零ネ章**零傘*本*零*ネ零零掌零章P c find IaugQ loom ledgeQ with sampli ng rate n55g in unit of l/l■ 「eal laug(11 integer ntaug、 n55g、 ledge(1)、 nledg entaug = nledge do i = 1.  nlaug laug(i) = float(ledge(i)−1)/float(ns sg)endd。
eturn nd C車零零傘#*零傘ネ傘零ネ****本零It*ネ零tll傘傘*零ネネネ*傘 ネ零車章章零*ネ零*零零零傘掌ネ**市零孝*零掌零零本サブルーチン+gg   del (sgg、n5gg1.n5gg2.ldigg、ns+g、   laug、 nlaug。
&sad、  n5ld、  n5ld1.  n5sld)C傘傘傘*章傘傘 掌傘I傘*零車掌傘*零章章市傘ネ傘牢ネ傘零傘零車傘ネ傘本ネ本零零ネ零車傘 *孝傘零零零零傘傘孝零本市車零零零傘市■■ c generate a ghost reference matrix  f「om +1dQc  with edge 1ocated at jau gQc sggQ −ghost refeIence maHixc n5g g1   [0Wdim、  01 sgg(1(O[# of +ample + in each column)c nssg2− column dim 、 of sggo (o「# of ghosts)c ldigg −le ading dimension of sggQ slo+agec n55 g −8ampling rate in unit of fcc laug O−edge 1ocation of each correspondin g sgQ  in sggQ  in un奄煤@ol I/Ic c ntaug −# ol ghost+c +ldQ −sl delai lQc n5ld −n  5ld c n5ldl −edge 1ocation in sld in uni t of samplec n+sfd −n  sample  +Id i n unit of lccomplex +gg(ldsgg、 l)、山( 1)real laug(1) integer n5gg1. nssg2. ldigg、 n55g、山u g、 ns+d、 n5ldl、 nssldnsgg2− nlaug do ig2 = 1.  n1gg2do 1g1−1.  n5gg1 j= n  +l + (igl−1)傘n5sfd / n55g百〇、 l t、 l) then sgg(igl、 1g2) = sfd[1)el+eii (j、gl、n +Id) thensgg(igl、 1g2) =山(nsld)lse sgg(igl、 1g21 = 5ld(i)ndil endd。
endd。
+elu+n nd C傘零章傘電傘傘傘ネネ傘傘傘を傘tt*ネ***傘I傘を零零掌京傘*本傘零 傘傘傘傘傘*傘**零傘傘傘孝ネ傘傘孝電傘傘傘傘傘*ネ孝磨■ サブルーチン5lddell  (sld、 n5fd、  n5tdl、 n ss[d、 d2.  d3)CネU車章車##章オ重n****傘*零零零零 章拳章零tl傘ネ章掌孝ネ*傘零零章ネ傘孝本零零ネ孝傘ネ傘孝*零零電本ネ本 零零ネネc 5lddel (s  +el  detail delined l: deiine +ele+ence veclo+c  5ldO−s   +ele+ence  veclo+c  n+ld  −n   5ld c n+dll −edge 1ocation in 5ldQc  n5h ld  −+ampliB  「ate  in  unit  of  Ic c d2. d3− 「ange lo be mapped onto lo t the Ieal part in composit■@vide。
C「山[to note Q page 68−75 tot notatio n and de+valioncomplex sld(1) inlege+ n5ld、 n5ld1. n5sld「eal d2.   d3 complex vs(256) real d6.  d7.  d8.  d9.  gn+ild =64 vs(64)□cmplx(−1,4791,−3,1833)vs (65)  :cmplx (−1,5434,−3,2476)vs I 66) =c +nplx (−1,6720,−3,2797)vs(67) = cmpl x (−1,6720,−3,3762)vs (68) ・cmplx (− 1,7363,−3,4084)W! (691・cmplx (−1,736 3,−3,4405)vi (701= cmplx (−1,8006,−3 ,5048)vs(71)=cmplx(−1,8006,−3,5370)v s (12) =cmplx (−1,8650,−3,60131vs(73 ) ・cmplx (−l、865G、   −3,6334)va (741 □ cmplx (−1,8650,−3,6977)vs (75) =cm plx (−1,8650,−3,7299)vs (76) = cmplx  (−1,8650,−3,8264)vs (771:cliplx (−1 ,8650,−4,0193)c    5lcl)□cmplx(−1,86 49517,−3,601286161)vx(7B)=cmplx(−1,8 650,−4,1429]vs (79) =cmplx (−1,8006, −4,2122)vs (80) ・cmplx (−1,8006,−4,3 0871v+ (811・cmplx (−1,8006,−4,3730)v s (82) = cmplx (−!、 8006.  −4.5016 ) vs (831=cmplx (−1,7363,−4,5981)vs(84 )・cmplx (−1,6077、−4,8232)7! (85) ・cm plx (−1,4791,−4,95181vs (86) = cmplx  (−1,4148,−5,1447)va (871=cmplx (−1, 4148,−5,4341)vt(88)・cmplx (−1,2862,− 5,6592)マs (89) :cmplx (−1,21162,−5,9 486)vs(9fll ;cmplx (刊1576、 −6.14+5 ) vs (91) =cmplx  (−1,0289,−6,5273)ws  (92) =cmplx (−0,9003,−6,7524)vs (931 =cmplx  (−0,8360,−6,9775)vs (94) =cm plx (−0,6431,−7,3312)vs (95) :cmplx  (−0,6431,−7,55271v+ f Q61 = cmplx I   −11,57118,−8,0386)vi (97) ・cmpD (−0 ,4502,−8,5531)vs I 981 :cmplx (−0,45 02,−11,9711)vs (99) ・cmplx (−0,4502, −9,3891)va(1001・cmplx(−0,4502−9,H2S) vi (1011・cmplx (−0,5145,−10,2251)vs( 102)=cmplx(−0,5788,−10,6431)vi(103)  =c町1に(−0,6431,−11,3826)vs (104) =cmp lx (−0,7717,−Il、 8971 )vs (105)  □   cmplx  (−0,9003,−12,37941v+ (106) ・c mplx (−1,0289,−12,7974)vs (107) =cmp lx (−1,0932,−13,21541vs (10g) ・cmplx  (−1,2219,−13,7299)vi (109) □ cmplx  (−1,3505,−14,5016)vs(110)=cmplx(−1,4 791,−15,0482)vs(III) ・cmplx (−1,8650 ,−15,7878)vs(1121:cmplx (−2,1865,−16 ,1736)vs (113)  = c+nplx  (−2,6367、− 16,6238)vs (114) = cmplx (−3,0868,−1 7,2990)vs(115)・cmplx(−3,6013,−17,652 7)vs(116):cmplx(−3,9871,−18,2315)vi  (1171=emptx(−4,6945,−19,0997)vs (118 ) =cmplx (−5,0804,−19,5177)vs(119)・c mplx(−5,852]、   −20,0000)vs (120) =c mplx (−6,4952,−20,4823)vs (121) ・cmp lx (−6,9453,−20,9646)vi(1221=cmplx(− 7,3955,−21゜5113 )vs (123) ・cmplx (−8 ,5531,−22,0579)vs (1241−cmplx (−9,2f il15.  −22.4759 )vs (1251= cmplx (−9 ,9678,−22,7974)vs(126) =cmplx (−1(1, 7395,−23,27971y+(127) ・cmplx (−11J18 3.  −23.6013 )vs(128)・cmplx (−12,218 6,−23,9871)vt(129) = cmplx  (−12,926 0,−24,2765)vs(130) = cmplx (i3.8264.   −24.5016 )vt(131) :cmplx (−14,6624 ,−24,7588)vs(132) :cmplx (−15,3698,− 24,9196)vs (133) □ cmplx (−16,0772−2 5,08041vs(134) ・cmplx (−16,9132,−25, 1768)o(135) ・cmplx (−18,0707,−25,209 0)c   5l(5) ・cmplx(−19,H9679,−25,241 158)vs(136) =cmplx (−19,0997,−25,241 2)vs(137) = cmplx  (−19,8714、−25,241 2)vt(138) = cmplx (−20,7074,−25,1768 )vs(139) ・cmplx (−21,8650,−25,1768)v s(140) =cmplx (−23,2797,−25,0482)vs( 141) :cmplx (−24,1801,−24,9196)vs(14 2) ・cmplx  (−24,7588,−24,72671vs(143 ) ・cmplx (−25,4019,−24,5016)vs(144’l  :cmplx (−26,2379,−24,2444)vs(145) : cmplx (−27,0740,−23,92211)vS(146) =c mplx (−27,9743,−23,53701vt(147) ;cmp lx (−28,8746、−23,2476)vs(148) ・cmplr  (−29,4534,−22,7010)vs(149) □ cmplx  (−29,9035,−22,3473)vs(150) = cmplr ( −30,5465,−21,96141vs(151) = cmplx (− 31,4469,−21,5434)vs(152) ・cmply (−32 ,4116,−20,73951vs(153) −cmplx (−33,1 190,−20,1286)vs(154) =cmplx (−33,697 8,−19,511201vs(155) ・cmplx (−34J408.  −19.1640)vsf1561 =cmplx (−34,7910,− 18,6817)vs(157) □ cmplx (−35,2412,−1 7,8457)vs(158)= cmplx (−35,7556,−17, 2669)vs([59) ・cmplx (−36,2701,−16,46 30)vs(160) ・cmplx (−36,7846,−15,9]64  )vs(161) =cmplx (−37,l7Q4. −15.4662  )vs(1621= cmplx (−37,1704,−14,6624) vs(163) =cmplx (−37,9421,−14,1158)vs (164) = cmplx (−38,1994,−13,6977)vs( 165) = cmplx (−38,5209,−13,1190)vs(1 661=cmplx (−38,7781,−12,2830)vs(1671 −cm山(−39,0354,−I+、 7363 )vx(168) ;cm plx (−39,1610,−11,1897)vs(169) □ cmp lx (−39,4212,−10,6752)vs(170) = cmpl x (−39,4212,−10,1929)vs(171)−cmplx(− 39,6785,−9,7106)vs (172) ・cmplx (−39 ,8071,−9,1640)vs(173) = cmplx (−39,9 357,−8,5209)vs (1741= cmplx (−39,871 4,−8,[1707)vs(175) ・cm山(−39,9357,−7, 7814)vs(1761= cmplx  (−40,00Q口、    − 7,4277)vs(177) = cm山(−40,0000,−7,106 1)Ielull end C傘I零*傘中$111章*傘$Ht$$零章市市傘傘傘傘*ネ*孝章本傘孝* **傘零*傘ネ傘**傘本*傘章傘*傘ネ孝零本*零サブルーチンcopyl   + (x、y、nx)C孝章傘零本本#零#零章傘零Uネnun傘傘ネ*孝* 零零零章本傘孝ネ本章を傘零ネ*零章ネ零傘本*零ネ零車牢零ネネネネc co p71−dim real vecto+ 10  to yQ「eal !( 1)、  y(1) integernx y(i) = x(i) endd。
Ielull nd C傘車*I零零傘零****傘零*掌傘ネ傘傘*中零*零零ネ**ネ掌零零*本 零傘ネ京孝ネ傘傘ネ傘孝本孝傘傘*電傘傘傘零*本ネ*サブルーチンlpd ( b、 nb、  er「、  tau、  Is、  to  x、  Ib) ctt幸I零*零傘零I零孝傘n章傘章*市を零ネネ*傘**零*零t1孝** 零零傘本傘本本車零傘車傘本本本零本本孝孝本本本c 1inea+ p+ed ilion coellicienls and euo+ calculat ionc  s    1pd(nus−1au)  = S  b(il   零 l(note−i+1)c cl: DSRCnote R2for de rivation and notationsc tau −lime de la7 f「om cuuenl samplec Is −sampling  1ntervalc nb −n b (no、 of coetlicie nls)、 nb =< 10c ro  I −exteInal auto cor+elation Iunclionc Ib −band width c −−−−on  return −−−−−c bQ −1inea+ p +edicailon coefficientsc eu −predict ion note −−−−1p −−−−− c cov  x  tpQ −sloragc for cova+1anc e mat「ixc ro  x  lpQ −conelalion vec lo+c c  lp −complex !to。
c  l −nb complex b(1)、  ro  Ireal  tau、  Is、   lb、  minteger nb complexcov  I  tp(Ion)、  to  1  tp(1 0)、  C1pin+ege+ 1.  l  Ip(10)、 m  1p (I[l)I:nb cov    I    1p(i+(i−1)零II  = ro   x( float(i−i)ネIs、 Ib)endd。
(o  X  1p(i) = +o−x(float(i−1)零ts−1a u、 Ib)endd。
callminyc (cow  X  tp、  1. c  lp、 I   lp、 m  lp)百((cabs (c  lp)桂)、 eQ、 1. ) prin口、’ Nl conditioned matrix’b(i)   =b(i)  + coy   I   1p(i+(i−1) 傘l)  傘 re  x   1p(i)endd。
endd。
c  tp = c  lp −eonig(ro  x  tp(i))零b  (i)nddO er「= cabs(c  lpl +ela+n nd C#傘H$$*孝孝I市章ネネ*ネ零ネ零*ネネ**ネネネ孝***掌市傘傘章 *本零零本牟本本*零ネネ傘孝傘**本零孝ネネ章市ネ傘傘complex f unction to  cot (tau、  Ib)C**本*掌*零本傘 *傘**ネネネネネ**零傘I車車ネ傘本ネ*零京*孝傘ネネネネ傘孝*孝車零 *傘ネ*孝零本孝*零ネネ零本本本零■ c autocorrelalion function In a unit  powe「signal withc  +a+ied cosinepow e「spectrum density 1unclion in (0,Ib )c cl page 123 of DSRCnote Oc to  co s −output Iunclion valuec tau −time (11) −1+qoenc7 teal  tau、  Ib complex 「o  fp pi = aco+(−1,1 c dl = lh  (ho「1zontal +can frequenc y in unit of Ic)dl = 4./910゜ nl = n1nl (Ib/dl) 1i = float(i)  章d1ro  lp = to  tp +  (1,+ ([atio−1,)傘(1,−cos (piiti))/2.1 endd。
to  cot = Q。
1i = float(it  掌d+to  cos  =  To  co  S +  (1,+  1alio−1,)  ネ (1,−cos (pi tli))/2. @) &           $  cexp(cmplx(0゜2.零pi$Ii ネtau))endd。
141111n nd c  ttt$ネ!幸本傘本##傘**傘傘*傘*ネ零章*章傘傘傘孝傘*傘車 本*ネ本*傘零傘傘ネ*傘*傘**傘零*傘*本*零ネ傘零■■■P complex function to  cos −「at  (tau、   Ib、  ratio)C本傘零*零傘零傘傘傘本零ネ*ネ*傘孝**** *傘零本零傘零**傘京傘傘***傘傘*孝*零零零傘傘*孝孝傘零零零*傘本 本**傘■■ c autocorrelalion function to a unit  powe「Iignal withc  +asied cosine po wer +pecl+um density function in (0, Ib)c cl:page 123 of DSRCnote Oc to   cos  rat −output function valuec  ta u  −time c Ib −1+quency c 「atio −peak to valle7 ratio in the  powe「5pecl+u+n assumed「eal tau、  Ib 、+aliocreal pi、  Ii complex 「01p pi = acos(−1,) c dl = Ih (horizontal 5can frequency  in unit ol Ic)d(・4./910 nl = n1nl(Ib/di) doi=G、nl +i=flogt(il傘d1 [e   Ip  =  to    lp  +  (1,+  (rati o−1,)  ネ (L −cos (pi 掌1i))/Q゜ en+Id。
1i = 1loal(i) $ dita  cos   rat = to   cot   rat + (1,+ (ratio刊、)傘(1,−cos  (pitli))^2、 &       * cexp (cmplx (0,、2,傘pitlitt au))endd。
to  coi  ral=+a  cos−rat/re  1pIetu+ n nd C零零*零ネ**傘ネ*傘傘傘*本*ネ***傘零零零傘本零本章*零零傘傘零 本#I*本零孝ネ本零零本傘ネネ車ネ零孝本零車*車ネ*本complex f unction to  flat  (tau、  1b)C傘傘傘孝#孝* 傘傘*章零零孝*ネ本章傘傘車ネネ**本零零本傘*傘掌本ネネ*零零零零孝零 ネネ#零零零掌零零本本参傘零零零零孝零c autocorrelalion  function to a un日power signal withc   flat power spectrum density functi on in  (0,Ib)c ro  flat −output func tion valuec  tau  −time C篩−1rquency teal  tau、   Ib 1a−a〔n5f−1,)$lau$lb百(abs (Ia)、 l 1.1 0. e−611hen+o  flat =l。
1ie 10−山t  =  cexp(cmplx(0,、Ia))幸sin  (f a)/Isndif eturn nd C傘傘*ネ*I**零零ネネ本ネt#ネ*車傘ネ章ネを置傘ネ傘本本ネ零市傘中 $本*孝ネ傘傘幸零零車孝傘傘$ネ*零参傘傘傘ネ傘傘傘牟l傘ネ傘拳ネ零孝傘 ネ サブルーチンlpd  rat (b、 nb、  e【「、  tau、   Is、  to  cos  +al、  Ib、  +x狽奄潤j C*#*零本**零章*ネ*傘本傘I*傘傘宰本傘傘零零車ネ傘零傘本本本零ネ 章市孝傘章*ネ市ネ章掌傘傘市*本*零掌零孝*本傘本**{本電零市零傘本本 傘 c 1inea+ p+ed山on coellicient+ and eu or calculationc  s    1pd(nets−1au)   −S  b(il  本 5(nt  1+−i+l)c  cl:  DS RCnote  R2lo【 derivation  and  nolat ion+C−−−−−Hent「y −−−−−c tau −time de lay loom cuuenl +amplec Is −sampling  1ntervalc nb −n b +no、 of coellicie nl+)、 nb =< IQc 「e  x −external aulo co++elalion functionc fb −band width C−−−−On  rellllll  −−−−−c bQ −1inea+  p+edicailon coeflicient+c etr −pred iction errorc cov  X  tpQ −sto+Rge l ot cova+1ance matrixc to  K  lpQ −co ++elelion veclo+c c  lp −complex no。
c  l  −nb complex b(1)、 ta  cos  rat「eal  tau、   Is、Ib、  +al  e++intege「nb complexcov  I  1p(100)、  IQ  X  III( 101,C1pinle+ge+ l、  I  1p(10)、 m  tp (10)I=nb cow  I  tp(i+(i−l)+1)=+o cos  rat(fl oat(i−i)$ls、Ib、raliolendd。
「o  !  1p(ij =+a  cos  rat(float(i−1 )傘1a−tau、 lb、 ratio)endd。
callminyc (cov  x  lp、  1. c  lp、 l   lp、 m  1p)tl ((cab+(c  lpQ珪)、 eq、 1 .1 print t、’ ill conditioned maf+ix’ b(i)  =bfil  +  COV  I   1p(i+(i−1)  ネ1) 本 to  x   1p(i)nddn endd。
c    Ip  =  c  ’  jp  −conjg(「o   r    1p(il)  孝 b(1)endd。
e++ = cabs(c  1p) Ielu+n nd c  tuft車孝##本ネ*零ネ#ネII零t1傘ネ*本ネ零*傘本ネ本*ネ 孝傘***零零ネ零傘*零ネ零本*京*零章孝*ネネネ牟孝帥l零ネネ零ネ零零 零 サブルすチ:/ hilbe+I (hh、 nhhlC*零零本*傘傘傘ネ傘 本*ネ傘零零ネ章ネネ電傘傘章傘*本零本市孝傘*車本*零零零本傘**ネネ零 掌傘車ネ零孝ネ零ネ孝孝ネ傘ネ*孝lネ傘零零傘本ネネ傘 c calculate truncated impulse respon se of the digital Hilbert t窒≠氏翌■盾窒■ c hhQ −impulse +eponsec nnn −posiliv ee time no、 of coefficientsreal hh(1 ) do i = I、 nhh hh(i) = O。
endd。
do i = I、  nhh、  2hh(i) = 2  章 (sin( piilloat(i) /2.)章峠)/(pi傘float(i))end d。
「eturn nd C5UB−ROI;TINE  MINVcCPURPO5E CINVERT A complex MATRIXCll5AGE CCALL MINVc (^、 II、 D、 L、 M+CDESCRIP TION OF PARAI、tETER3CA  −INPIIT  MAT RIX、   DESTROYED  IN  COMPuTATION  A ND  REPLACEc  BY CRESL’LTANT  INVER3ECN −0RDEROF MATR IX ACD −RESIILTANT DETERl、1lNANTCL − WORK VECTOROF LENGTHNCM −WORK VECTOR OF LENGTHNC CREMARKS CMATRIX A MtlST BE A GEIIERAL MATRIX C5UB−ROUTINES AND FIINCTION SUBPROGR AMS REQUIREDCN0NE CMETIIOD CTHE 5TANDARD GAIISS−jORDAN METIIOD  Is tlsED、 TIIE DETERMINANTCIs ALSOCA LCULATED、 A DETERMINANT OF zERo INDI CATES THATCTHE MATRIX Is 5INGULARC、、 、、、、、、、、、,0,、、、、、、、、、、、、、、−、、、、、、、、、 、、、、、、、、、、、、、、、、、、、、、、、、A、。
5tlBROUTINE MINVc(^、 N、 D、 L、 M)D1λI ENsION^m、 L il)、 M (1)C、、、、、、、、、、、、、 、、、、、、、、、、、、、、、、、、、、、、、、,0,1,、、、、、、、 、、、、、、、、、、、、、A、、、、、、。
CIF A DOUBLE PRECISION VER5ION OF TH IS ROUTINE Is DESIRED、  THECCIN C0LU lllN I 5HOULD BE REMOVED FROM THE DO IIBLE PRECISIONCSTATEMENT Wl−11CHFOL LOWS。
complex^、D、BIG^、 ll0LDCTHE  CLI[:ST   ALSOBE  RElitO〜Eft  FROM  DOUIILE   PRECISION  ST`TEMENTS CAPPEARING IN 0T)IERROUTINES USED Il I C0NIUNCTION WITHT)ItsCROliTINE CTIIE 1101:IILE PRECISION VER5ION OF  Tlll55IIB−ROUTINE MIIST ALrO CC0NT^IN DOUBLE PRECISION FORTRAN FI INKTIONS、  ^BS IN STATEMENTC10MUST B E CHANGED To DABS。
C、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、 、、、、、、、、、、、、、、、、、、、、、、、、、、、A、、、、、、、、 、、。
C5EARCHFORLARGF、ST ELEMENTD=1.0 NK=−N Do 80 K・1.N NK=NK+N L (K)・K M (K)・K KK=NK1.K BIG人二^(KK) no 201=に、 N +2=++*(1−1) DO201=に、 N II  =12+1 15  IF(cABS([1lGA)−cABsf^(団)) 15.20. 2018 8IGA=Af団 L (K) = 1 M (K) =1 20  C0NTINUE CINTERcIIANGE ROWSJ=L (KI IFO−K)  35.35.25 )10LD=−A (K 11 11=KI−に+1 A (K l) −A (J I) 30  ^(月)  =lIOLD INTERcHANGE COLUMNS351・M (KI 1F・(1−K)  45.45.3838 1P=N零(1−1) )10LD=−^(IKI A fJK) =^(Ill 40  A(J+) =HOLD CDIVIDE C0Ltllll BY MINliS PIVOT (VA LUE OF PIVOT ELEMENT IsCC0NTAINED IN  BIGA)45百(biga、 ne、 cmplx (0,0,11gol e 4846  D・0.O ETURN 4800551=l、N 1F(1−K)  50.55.50 501に=NK+1 ^(IK)=A(IK)/(−BIGA)55   C0NTINじE CREDUCE MATRIX DO651ゴ1.11 1に:Nに+1 )101.D二^(IKl 11=l−N Do 651=1. N 11=lI+N 1F(1・K)  60.65.60 60 1F(1−Kl  62.65.6262  KI=11−1十K A(団=lIOLD孝^(Kl) +A (l I)65  C0NTINli E CDI〜IDE ROW BY PIVOTKI=に−N DO751=l、 N KI=KI+N 1F(1−1) 7Q、?5,1Q ?OA(Kl)・^(Kll/BIG^75  C0NTINUE CPRODIICT OF PIVOTSD=D傘BIGA CREPLACE PIVOT BY RECIPROCAL^(KK) =1 . Q/BIGA 80  C0NTINtlE CFINAL ROW AND COLUMN INTERcII^NGEK= N 100K・(K−1) IF(K)  150.150.105特表千3−503347 (44) と 国際調査報告 国際調査報告 US 8900410

Claims (1)

  1. 【特許請求の範囲】 1.トレーニング信号成分を有する直接信号、および望ましくない多重通路歪み 成分を形成する遅延信号を有する信号を含んでいる入力信号の源と、出力ポート 、前記入力信号の源に結合される入力ポート、およびプログラム可能なフィルタ 係数値を有する無限インパルス応答濾波手段とを含み、 実質的に多重通路歪みのない前記トレーニング信号を表わす基準信号の源(28 5)と、 前記入力信号の源(220)および前記基準信号の源(285)に結合され、前 記基準信号に応答して前記入力信号にほぼ等しい信号を発化する信号濾波システ ムの数学的モデルを作成し、且つ前記遅延信号を実質的に除去した前記直接信号 を出力信号として前記出力ポートに発生するように前記無限インパルス応答濾波 手段を条件づけるために前記数学的モデルで決まる係数値で前記無限インパルス 応答濾波手段(280)の係数値をプログラミングするための制御手段(282 )とを含んでいることを特徴とする多重通路歪みを補正するシステム。 2.前記無限インパルス応答手段がプログラム可能なタップ付き遅延要素を有し 、さらに前記数学的モデルで決まる遅延値で前記無限インパルス応答濾波手段の タップ付き遅延要素をプログラミングするための手段(654)を含むことを特 徴とする請求項1に記載のシステム。 3.サンプリング用クロック信号源(222)と、前記入力信号源(220)に 結合され、前記サンプリング用クロック信号で決まる時点において前記入力信号 を表わすディジタル・サンプルを発生するアナログ・ディジタル変換手段(23 2)であって、前記無限インパルス応答濾波手段(280)および前記制御手段 (282)に前記ディジタル・サンプルを供給するように結合されている前記ア ナログ・ディジタル変換手段とを含み、前記基準信号源が前記基準を表わすサン プル値でプログラミングされる読出し専用メモリ手段(285)を含んでおり、 前記制御手段が、 読出し・書込みメモリ手段(284)と、前記読出し専用メモリ手段および前記 読出し・書込みメモリ手段に結合され、プログラムにより制御されて前記入力信 号のトレーニング信号成分を表わす複数個の前記ディジタル・サンプルを貯え、 前記信号濾波システムの前記数学的モデルを発生し、前記無限インパルス応答濾 波手段についての前記係数値および前記タップ付き遅延要素の値を発生するデー タ処理手段(282)を含むことを特徴とする請求項2に記載のシステム。 4.前記無限インパルス応答濾波手段(280)が、前記入力信号を表わすディ ジタル・サンプルを供給するための入力ポート(Y1)と、 前記入力ポートに結合され、前記入力信号を表わすディジタル・サンプルと補正 信号を表わすディジタル・サンプルを合成し、前記出力信号を発生する第1の信 号合成手段(610)と、 前記第1の合成手段に結合され、それから供給される信号を遅延させるプログラ ム可能な信号遅延手段(622)であって、複数個の出力ポートを有し、且つ前 記合成手段により供給される信号について第1のプログラミング制御信号により 決まる遅延時間だけ遅延させた信号を前記出力ポートの各々に供給するための回 路を有する前記プログラム可能な信号遅延手段と、 各々が前記プログラム可能な信号遅延手段の出力ポートのそれぞれ異なる1つに 結合される複数個のプログラム可能な係数乗算手段(656,658,660) であって、該係数乗算手段の各々に供給される第2のプログラミング制御信号に より決まる前記係数値の各々を前記プログラム可能な信号遅延手段から供給され る信号に掛ける前記複数個のプログラム可能な係数乗算手段と、前記複数個のプ ログラム可能な係数乗算手段の各々に結合され、それから供給される各信号を合 成して前記補正信号を発生する第2の信号合成手段(606,662)とを含み 、 前記データ処理手段が前記プログラムにより制御され、前記第1のプログラミン グ制御信号と前記第2のプログラミング制御信号を発生することを特徴とする請 求項3に記載のシステム。 5.前記プログラム可能な信号遅延手段が、前記サンプリング・クロック信号源 に結合され、前記第1の合成手段(610)の出力信号を受け取るように結合さ れる入力ポートを有し、前記第2のプログラミング制御信号に応答し、前記第1 の合成手段の出力信号について、これとは前記サンプリング・クロック信号の周 期のそれぞれ異なる整数倍だけ前記第1の合成手段の出力信号に対して遅延させ た複数個の選択された遅延信号を各出力ポートに発生するサンプル・データ信号 遅延手段(622)と、 前記サンプル・データ信号遅延手段の出力ポートの中の少なくとも1つに結合さ れ、前記第2の制御信号に応答し、前記サンプル・データ信号遅延手段の1つの 出力ポートに供給される信号を前記サンプリング・クロック信号の周期のほぼN /M(但し、Mは1より大きい整数であり、Nは前記第2の制御信号により決ま る値をとる正の整数である。)倍に等しい時間量だけ事実上遅延させる信号補間 手段(920)とを含むことを特徴とする請求項4に記載のシステム。 6.トレーニング信号成分を有し搬送波信号を変調する直接信号、および望まし くない多重通路歪み成分を形成する遅延信号を含んでいるビデオ信号源(210 )と、前記ビデオ信号源に結合され、それぞれ同相であり前記搬送波と直角位相 関係にあるベースバンド信号を表わす第1および第2の信号成分に前記ビデオ信 号を復調する復調手段(220,230)と、 前記第1および第2の信号をそれぞれ受け取るように結合される第1および第2 の入力ポートを有し、プログラム可能な複素フィルタ係数値を有する複素無限イ ンパルス応答濾波手段(280)とを含み、多重通路歪みを実質的に含んでいな い前記第1および第2の信号のトレーニング信号成分をそれぞれ表わす第1およ び第2の基準信号源(285)と、前記復調手段と前記第1および第2の基準信 号源(285)に結合され、前記第1および第2の基準信号に応答して前記第1 および第2の信号のトレーニング信号成分にそれぞれほぼ等しい第1および第2 の合成信号を発生する信号濾波システムの数学的モデルを形成し、且つ前記遅延 信号を実質的に除去した前記直接信号の第1の信号成分を出力信号として前記出 力ポートに発生するように前記無限インパルス応答濾波手段を条件づけるために 前記数学的モデルで決まる複素係数値で前記無限インパルス応答濾波手段の複素 係数値をプログラミングするための制御手段(282)とを含んでいることを特 徴とするビデオ信号中の多重通路歪みを補正するシステム。 7.前記無限インパルス応答濾波手段がプログラム可能なタップ付き遅延要素を 有し、 前記制御手段が前記数学的モデルで決まる遅延値を使って前記無限インパルス応 答濾波手段のタップ付き遅延要素をプログラミングするための手段(654)を 含むことを特徴とする請求項6に記載のシステム。 8.サンプリング用クロック信号源(222)と、前記復調手段(220,23 0)に結合され、前記サンプリング用クロック信号で決まる時点において前記第 1および第2の信号を表わす第1および第2のディジタル・サンプルのシーケン スをそれぞれ発生する第1および第2のアナログ・ディジタル変換手段(232 ,234)であって、前記無限インパルス応答濾波手段(280)および前記制 御手段(282)に前記第1および第2のディジタル・サンプルのシーケンスを 供給するように結合されている前記第1および第2のアナログ・ディジタル変換 手段とを含み、前記第1および第2の基準信号源が前記第1および第2の基準信 号を表わすサンプル値を使ってプログラミングされる読出し専用メモリ手段(2 85)を含んでおり、前記制御手段が、 読出し・書込みメモリ手段(284)と、前記読出し専用メモリ手段および前記 読出し・書込みメモリ手段に結合され、プログラムにより制御されて前記ビデオ 信号のトレーニング信号成分を表わす前記第1および第2のディジタル・サンプ ルのシーケンスの各々から複数個のディジタル・サンプルを貯え、前記信号濾波 システムの前記数学的モデルを発生し、前記無限インパルス応答濾波手段(28 0)について前記複素係数値および前記タップ付き遅延要素の値を発生するデー タ処理手段(282)を含むことを特徴とする請求項7に記載のシステム。 9.前記無限インパルス応答濾波手段が、前記第1および第2のディジタル・サ ンプルのシーケンスを供給するための第1および第2の入力ポートと、前記第1 および第2の入力ポートにそれぞれ結合され、前記第1および第2のディジタル ・サンプルのシーケンスとそれぞれ第1および第2の補正信号を表わす別個の第 1および第2のディジタル・サンプルのシーケンスとを合成し、補正サンプルの 第1のシーケンスが前記無限インパルス応答濾波手段の出力信号である第1およ び第2の補正サンプルのシーケンスを発生する第1(610)および第2(62 0)の信号合成手段と、前記第1および第2の合成手段に結合され、それらから 供給される第1および第2の補正サンプルのシーケンスをそれぞれ遅延させるプ ログラム可能な信号遅延手段(622)であって、複数個の出力ポートを有し、 且つ前記第1および第2の補正サンプルのシーケンスについて第1のプログラミ ング制御信号により決まる遅延時間量だけ遅延させたものを前記出力ポートの各 々に供給するための回路を有する前記プログラム可能な信号遅延手段と、各々が 前記プログラム可能な信号遅延手段の出力ポートのそれぞれ異なる1つに結合さ れる複数個のプログラム可能な複素係数乗算手段(656,658,660)で あって、該複素係数乗算手段の各々に供給される第2のプログラミング制御信号 により決まる前記複素係数値の各々を前記プログラム可能な信号遅延手段から供 給される第1および第2の補正サンプルのシーケンスの遅延されたものに掛ける 前記複数個のプログラム可能な複素係数乗算手段と、 前記複素係数乗算手段の各々に結合され、それらから供給されるスケール化され たサンプルの第1および第2のシーケンスを合成し、前記第1および第2の補正 サンプルのシーケンスを発生する第3および第4のサンプル合成手段(662− 668)とを含み、前記データ処理手段が前記プログラムにより制御され、前記 第1のプログラミング制御信号と前記第2のプログラミング制御信号を発生する ことを特徴とする請求項8に記載のシステム。 10.前記プログラム可能な信号遅延手段が、前記サンプリング用クロック信号 源に結合され、補正サンプルから成る第1および第2のシーケンスをそれぞれ受 け取るように結合される第1と第2の入力ポートを有し、前記第2のプログラミ ング制御信号に応答し、補正サンプルから成る前記第1および第2のシーケンス に対して前記サンプリング用クロック信号の周期のそれぞれ異なる整数倍だけ遅 延させた補正サンプルから成る前記第1および第2シーケンスについての複数個 の選択された遅延シーケンスを各出力ポートに供給するサンプル・データ信号遅 延手段(622)と、 前記サンプル・データ信号遅延手段の出力ポートの中の少なくとも1つに結合さ れ、前記第2の制御信号に応答し、該出力ポートから供給され補正サンプルから 成る前記第1および第2のシーケンスについての選択された遅延シーケンスを前 記サンプリング用クロック信号の周期のほぼN/M(但し、Mは1より大きい整 数であり、Nは前記第2の制御信号により決まる値をとる正の整数である。)倍 に等しい時間量だけ事実上遅延させるサンプル値補間手段(920)とを含むこ とを特徴とする請求項9に記載のシステム。
JP1502655A 1988-02-08 1989-02-06 多重通路歪みを補正するシステム Expired - Fee Related JP3051984B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US153,748 1988-02-08
US07/153,748 US4864403A (en) 1988-02-08 1988-02-08 Adaptive television ghost cancellation system including filter circuitry with non-integer sample delay

Publications (2)

Publication Number Publication Date
JPH03503347A true JPH03503347A (ja) 1991-07-25
JP3051984B2 JP3051984B2 (ja) 2000-06-12

Family

ID=22548576

Family Applications (1)

Application Number Title Priority Date Filing Date
JP1502655A Expired - Fee Related JP3051984B2 (ja) 1988-02-08 1989-02-06 多重通路歪みを補正するシステム

Country Status (7)

Country Link
US (1) US4864403A (ja)
EP (1) EP0441778B1 (ja)
JP (1) JP3051984B2 (ja)
KR (2) KR900701117A (ja)
AT (1) ATE107103T1 (ja)
DE (1) DE68916034T2 (ja)
WO (1) WO1989007376A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05167886A (ja) * 1991-12-18 1993-07-02 Nec Corp 映像信号補償回路
JP2005039838A (ja) * 2003-07-16 2005-02-10 Samsung Electronics Co Ltd 並列構造のntsc除去フィルタ及びフィルタリング方法

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4980767A (en) * 1989-11-30 1990-12-25 At&T Bell Laboratories Technique for determining signal dispersion characteristics in communications systems
US5260972A (en) * 1990-03-13 1993-11-09 At&T Bell Laboratories Technique for determining signal dispersion characteristics in communications systems
US5121211A (en) * 1990-10-09 1992-06-09 North American Philips Corporation System for echo cancellation comprising an improved ghost cancellation reference signal
US5119196A (en) * 1990-06-25 1992-06-02 At&T Bell Laboratories Ghost cancellation of analog tv signals
US5065242A (en) * 1990-06-29 1991-11-12 General Electric Company Deghosting apparatus using pseudorandom sequences
US5177611A (en) * 1990-07-31 1993-01-05 Rca Licensing Corporation Method and apparatus for canceling quadrature distortion as for video signals having in-phase and quadrature phase components
FR2665599A1 (fr) * 1990-08-03 1992-02-07 Europ Rech Electr Lab Procede de correction a l'emission de signaux de television et dispositif de mise en óoeuvre.
FR2665603A1 (fr) * 1990-08-03 1992-02-07 Europ Rech Electr Lab Procede de correction de defauts de signaux de television lors de leur transmission et dispositif de mise en óoeuvre.
US5172232A (en) * 1990-10-09 1992-12-15 North American Philips Corporation Method and apparatus for communication channel identification and signal restoration
US5111298A (en) * 1990-10-09 1992-05-05 North American Philips Corporation Method and apparatus for communication channel identification and signal restoration
US5047859A (en) * 1990-10-09 1991-09-10 North American Philips Corporation Method and apparatus for communication channel identification and signal restoration
US5179444A (en) * 1990-10-09 1993-01-12 North American Philips Corporation System for echo cancellation comprising an improved ghost cancellation reference signal
KR960015390B1 (ko) * 1991-12-06 1996-11-11 삼성전자 주식회사 고스트제거방법 및 장치
KR100269004B1 (ko) * 1992-01-08 2000-10-16 카츠넬슨 론 디. 다중채널 텔레비젼 신호 스크램블링과 디스크램블링 시스템 및 방법
US5208596A (en) * 1992-04-10 1993-05-04 Rca Thomson Licensing Corporation DAC distortion compensation
US6937292B1 (en) 1992-04-22 2005-08-30 Samsung Electronics Co., Ltd. Ghost cancellation reference signal with bessel chirps and PN sequences, and TV receiver using such signal
US5384865A (en) * 1992-06-01 1995-01-24 Eastman Kodak Company Adaptive, hybrid median filter for temporal noise suppression
US5278637A (en) * 1992-09-25 1994-01-11 Matsushita Electric Corporation Of America Apparatus and a method for quadrature-phase encoding and transmitting a digital signal in a video signal
US5309235A (en) * 1992-09-25 1994-05-03 Matsushita Electric Corporation Of America System and method for transmitting digital data in the overscan portion of a video signal
US5430487A (en) * 1992-09-30 1995-07-04 Matsushita Electric Corporation Of America Method and apparatus for improving temporal video signal processing using motion vectors transmitted with the video signal
TW231393B (ja) * 1992-12-02 1994-10-01 Samsung Electronics Co Ltd
US5335020A (en) * 1993-03-03 1994-08-02 Rca Thomson Licensing Corporation IIR ghost cancelling system with reduction of filter instability
US5532748A (en) * 1995-03-31 1996-07-02 Matsushita Electric Corporation Of America Hybrid analog/digital television transmission system
KR100269130B1 (ko) * 1997-11-21 2000-10-16 윤종용 단일고스트제거기를갖는디지털/아날로그tv방송공용수신기와고스트제거방법
US7038732B1 (en) * 1998-05-12 2006-05-02 Samsung Electronics Company, Ltd. DTV signal with GCR components in plural-data-segment frame headers and receiver apparatus for such signal
US6426972B1 (en) 1998-06-19 2002-07-30 Nxtwave Communications Reduced complexity equalizer for multi mode signaling
US7915962B2 (en) * 2009-07-06 2011-03-29 Nortel Networks Limited System and method for built in self test for timing module holdover
TW201240445A (en) * 2011-03-17 2012-10-01 Novatek Microelectronics Corp Video signal processing circuit and method applicable thereto

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59211315A (ja) * 1983-05-17 1984-11-30 Toshiba Corp 自動等化器におけるタツプ利得修正方法
JPS62159985A (ja) * 1985-12-24 1987-07-15 アールシーエー トムソン ライセンシング コーポレイシヨン 多ゴースト打消し装置

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4127874A (en) * 1976-05-27 1978-11-28 Tokyo Shibaura Electric Co., Ltd. Apparatus for removing ghost signals from received video signals
GB1595636A (en) * 1977-07-22 1981-08-12 Mitsubishi Electric Corp Apparatus for reducing ghosting of tv pictures
JPS54109720A (en) * 1978-02-16 1979-08-28 Toshiba Corp Ghost signal erasing device
US4328772A (en) * 1980-01-14 1982-05-11 Cummins Engine Company, Inc. Combustion chamber for an internal combustion engine
US4458269A (en) * 1980-05-07 1984-07-03 Rca Corporation Sampling arrangement for a television ghost signal cancellation system
US4359778A (en) * 1981-02-05 1982-11-16 Zenith Radio Corporation Channel equalizer and method for cancelling ghosts
US4413282A (en) * 1981-04-06 1983-11-01 Rca Corporation Television ghost signal detection during the video information interval
US4393397A (en) * 1981-10-05 1983-07-12 Rca Corporation Television ghost signal detector with color burst phase delay control
US4564862A (en) * 1982-08-09 1986-01-14 Edwin Cohen Ghost signal elimination circuit
JPS59211388A (ja) * 1983-05-17 1984-11-30 Toshiba Corp ゴ−スト除去装置
JPS59211313A (ja) * 1983-05-17 1984-11-30 Toshiba Corp 自動等化器
US4542408A (en) * 1983-08-31 1985-09-17 Rca Corporation Digital deghosting system
JPS6094580A (ja) * 1983-10-28 1985-05-27 Hitachi Ltd ゴ−スト除去装置
JPS60199275A (ja) * 1984-03-24 1985-10-08 Toshiba Corp ゴ−スト消去装置
US4694414A (en) * 1984-12-19 1987-09-15 Rca Corporation Digital delay interpolation filter with amplitude and phase compensation
US4698680A (en) * 1985-12-24 1987-10-06 Rca Corporation Digital correlation apparatus as for a television deghosting system
US4727424A (en) * 1986-10-16 1988-02-23 Rca Corporation Sampled data filtering system, including a crossbar switch matrix, as for a ghost cancellation system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59211315A (ja) * 1983-05-17 1984-11-30 Toshiba Corp 自動等化器におけるタツプ利得修正方法
JPS62159985A (ja) * 1985-12-24 1987-07-15 アールシーエー トムソン ライセンシング コーポレイシヨン 多ゴースト打消し装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05167886A (ja) * 1991-12-18 1993-07-02 Nec Corp 映像信号補償回路
JP2005039838A (ja) * 2003-07-16 2005-02-10 Samsung Electronics Co Ltd 並列構造のntsc除去フィルタ及びフィルタリング方法
JP4550502B2 (ja) * 2003-07-16 2010-09-22 三星電子株式会社 並列構造のntsc除去フィルタ及びフィルタリング方法

Also Published As

Publication number Publication date
KR900701117A (ko) 1990-08-17
KR930010427B1 (ko) 1993-10-23
DE68916034T2 (de) 1994-12-15
EP0441778B1 (en) 1994-06-08
JP3051984B2 (ja) 2000-06-12
ATE107103T1 (de) 1994-06-15
US4864403A (en) 1989-09-05
DE68916034D1 (de) 1994-07-14
WO1989007376A1 (en) 1989-08-10
EP0441778A1 (en) 1991-08-21

Similar Documents

Publication Publication Date Title
JPH03503347A (ja) 多重通路歪みを補正するシステム
US5341177A (en) System to cancel ghosts generated by multipath transmission of television signals
US5481316A (en) System, apparatus and method for canceling televison ghost signals
US5461426A (en) Apparatus for processing modified NTSC television signals, with digital signals buried therewithin
US6175389B1 (en) Comb filtered signal separation
JP4559734B2 (ja) テレビ受信機およびライン/フレームパターンと同期して信号サンプルストリームを処理するための方法
JPH06233156A (ja) Tv受信機あるいはビデオレコーダに使用されるゴースト消去基準信号捕捉回路
US20070058759A1 (en) Receiver having decisional feedback equalizer with remodulation and related methods
KR0153610B1 (ko) 영상 캐리어와 직각 위상을 이루는 캐리어를 변조한 신호에서의 고스트 억압 장치
EP0558017B1 (en) Contour restoration apparatus
US4698680A (en) Digital correlation apparatus as for a television deghosting system
JPS628070B2 (ja)
KR950014188B1 (ko) 다중 경로 왜곡 보정 시스템
JP2514854B2 (ja) 映像信号のノイズ除去装置
KR0123189B1 (ko) 파형등화기
Chao Multi-path equalization for NTSC video by using digital IIR filter
US20100208135A1 (en) Method for processing an analog television signal
JPH06225183A (ja) ゴースト除去装置
JP3019876B2 (ja) 波形等化システム
JP2570116B2 (ja) 送信装置
JPH09135402A (ja) 波形等化装置
JPH0239774A (ja) ゴースト除去装置
JPH02174377A (ja) ゴースト除去装置
JPH01188116A (ja) ゴースト除去装置
JPS6310633B2 (ja)

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees