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
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/14—Picture signal circuitry for video frequency region
- H04N5/21—Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03012—Arrangements for removing intersymbol interference operating in the time domain
- H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
- H04L25/03038—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception with a non-recursive structure
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/14—Picture signal circuitry for video frequency region
- H04N5/21—Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
- H04N5/211—Ghost 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.トレーニング信号成分を有する直接信号、および望ましくない多重通路歪み 成分を形成する遅延信号を有する信号を含んでいる入力信号の源と、出力ポート 、前記入力信号の源に結合される入力ポート、およびプログラム可能なフィルタ 係数値を有する無限インパルス応答濾波手段とを含み、 実質的に多重通路歪みのない前記トレーニング信号を表わす基準信号の源(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に記載のシステム。
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)
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)
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)
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)
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 |
-
1988
- 1988-02-08 US US07/153,748 patent/US4864403A/en not_active Expired - Lifetime
-
1989
- 1989-02-06 JP JP1502655A patent/JP3051984B2/ja not_active Expired - Fee Related
- 1989-02-06 KR KR1019890701840A patent/KR900701117A/ko not_active IP Right Cessation
- 1989-02-06 KR KR1019890701840A patent/KR930010427B1/ko active
- 1989-02-06 DE DE68916034T patent/DE68916034T2/de not_active Expired - Fee Related
- 1989-02-06 WO PCT/US1989/000410 patent/WO1989007376A1/en active IP Right Grant
- 1989-02-06 EP EP89902863A patent/EP0441778B1/en not_active Expired - Lifetime
- 1989-02-06 AT AT89902863T patent/ATE107103T1/de not_active IP Right Cessation
Patent Citations (2)
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)
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 |