JPS5845671B2 - 周期的な走査方法によつて得られる時間依存の走査信号のｓ／ｎ比を向上させる方法 - Google Patents

Description

【発明の詳細な説明】 本発明は周期的な走査方法によって得られる時間依存の
走査信号におけるＳ／Ｎ比を向上させるための方法に係
り、境界を有するとともに中心が定められたフィールド
を走査装置によって線状に走査し、瞬間的走査位置から
フィールド中心までの距離に対応しかつ距離の瞬間値ｒ
および最大値Ｒを含む時間依存の追加信号を発生し、順
次所定の角度だけずらしながら線走査を行ない、この線
走査の隣接する線がフィールド境界においては互いに接
し、フィールド内部においては部分的に重畳しまたフィ
ールド中心においては全面的に重畳するような前記線走
査によって得られる時間依存の走査信号におけるＳ／Ｎ
比を向上させるための方法に関する。

本発明の課題ならびに解決法をより良く理解するために
は、前記走査方法を例に基すいて説明する必要がある。

しかしこれは本発明の上位概念を記載例に限定するもの
ではない。

赤外線領域で機能する観察装置においては、観察される
フィールドは光学的に撮像されそしてこの映像面は赤外
線に感じる検出器により走査される。

この検出器は、所定の面積を有し、映像フィールドの直
径に沿って移動することにより映像面において円形の映
像フィールドを走査し、この際時間的に連続する軌道は
、映像フィールドの周辺において丁度全ての映像点の隙
間のない検出が行われるように成る角度だけ異る。

全ての軌道は映像フィールドの中心を通り、軌道の重畳
が生じる。

その結果として、映像フィールド全体の走査を行ってい
るうちに、映像フィールド内の各位置は映像中心までの
距離に応じて複数回走査されることになり、しかもその
走査回数は中心までの距離の逆数にほぼ比例する。

この特性は特定の用途、例えば目標追跡装置においては
極めて有意義であり、この場合、映像上の各位置が複数
回走査されている事実を情報理論に基づいて最適評価す
る方法が存在することを前提としている。

なぜならば、このような評価が行われなければ、検出器
ノイズや増幅器ノイズによって検出器信号は誤解釈され
、また観察されるフィールド内の光学的背景によっても
検出器信号は少なからぬ悪影響を受けるからである。

すなわち従来の方法では情報処理理論に基づく高度な評
価技術を用いなければ、充分にＳ／Ｎ比を高めることが
できなかった。

フィールドの純粋な撮像に対してフィールドの走査に相
応する信号を解釈してさらに使用することに関し改良を
行い、そして例えば信号を閾値検出器によってより良く
評価することを保障するために、映像位置の前記の何回
もの走査を有利に行う方法を提供することが本発明の課
題である。

この課題は本発明によれば次のようにして解決される。

すなわち、走査信号に基づいて一定の序列で発生された
Ｎ個の主信号列のうち１番目の主信号Ａ１は走査信号そ
のものとし、２ないしＮ番目の主信号列の順位ｉ番目の
主信号Ａｉは唯一の時間遅延装置およびこの時間遅延装
置と接続された出力信号の帰還路によって循環法により
発生されること、追加信号を２つの関数発生器にそれぞ
れ印加し、１番目の主信号に対応して配設された一方の
関数発生器から、関数（１−ｅｘｐ（−百））に比例す
る補助信号を発生させ、そしてｉ番目の主信号に対応し
て配設された他方の関数発生器から、関数〔ｅｘｐ
（−ｉ））に比例する補助信号を発生させること、主信
号のそれぞれに対応して配設された２個の乗算装置によ
って前記補助信号と前記主信号とをそれぞれ乗算して各
中間信号を発生させること、および加算装置によって前
記各中間信号を加算し出力信号を発生させること、の各
過程から成る方法によって解決される。

このようにして時間依存走査信号は相関分析処理され、
Ｓ／Ｎが著しく改良され、実質的にわずかな費用で他の
近似分析が行わへ実用上十分に良好な結果が得られる。

次に図面に従い本発明にかかる方法およびその好ましい
実施例を詳細に説明する。

第１図は、赤外線観察装置の円形状映像面の一部を示す
。

Ｚは中心を、Ｕは映像面の周辺を示し、その円の半径は
値Ｒを有する。

映像面には、この場合中心Ｐおよび辺長ｄの四角形とし
て形成されている走査装置である放射検出器Ｓがある。

検出器Ｓの中心Ｐおよび映像面中心２間の間隔は値ｒを
有する。

映像面には図示されていない光学手段から光線が放射さ
れ、放射検出器Ｓはこの光線を受光してその光量に比例
した走査信号ｆ（１）を発生する。

そして放射検出器Ｓは中心Ｚからの距離ｒに対応した追
加信号ｒ（、）を発生する。

ここにこの実施例では説明を簡単にするためにｒ（ｔ）
を距離ｒに等しく設定している。

放射検出器Ｓは、定速で直線上の軌道上の映像を走査し
、中心Ｐは映像面の直径を例えばＤのように移動しそし
て周辺から周辺へ即ち直径りの一端から他端へ映像面を
横断する。

それぞれの線走査が終った後、放射検出器Ｓはその軌道
の終りから軌道の始めに戻りそして軌道方向は、新しい
位置もしくは方向の軌道が前の位置もしくは方向の軌道
へ周辺Ｕにおいて丁度液するようにＺを中心として所定
の角度だけ変化する。

第１図において近似されているように、軌道方向を変化
する角度は、値（２ｔａｎ−１（−） ）に相Ｒ 当する。

軌道に沿った放射検出器の帰還の全過程ならびに軌道方
向の回転の全過程は、線走査の時間に対して非常に短く
、かつそれ故重要でない時間内に行われる。

このことから全映像面の走査は、π 線走査の時間の倍数即ち 倍ある２ｔａ
ｎ”（±） Ｒ いは近似的に（πＲ／ａ ）倍の時間期間に相当する時
間内に行われるようになる。

前記から、軌道における放射検出器の移転は、時間経過
について周期的なのこぎり波関数に相当し、これにより
線走査の周期およびフィールド走査の周期が規定され得
ることが分る。

映像走査の前記説明は、本発明の詳細な説明のためのみ
のものであり本発明を限定するものではない。

観測フィールドが前記例において同じ方法で走査される
たびに等しい方式が生れる。

例えば、放射検出器は静止しそして光線が可動鏡システ
ムにより偏向され、あるいはフィールドがレーザもしく
はマイクロ波ビームにより同じ走査方式で照射されそし
て反射されたビームが静止した放射検出器により検出さ
れる。

しかし課題およびその解決法を理解するためには、第１
図に示す場合の説明で十分である。

第１図における放射検出器Ｓの互に隣り同志重畳した軌
道において、中心Ｚから距離ｒを有する映像点がフィー
ルド走査の半周期の間即ちπだけの軌道方向の回転の間
、この距離ｒに応じて１つから多数検出される。

距離Ｒを有する映像点は正確に一度検出される。

何故なら映像円の周囲に在すモして放射検出器Ｓの軌道
はそこで互に丁度接触するからである。

これは例えばｑ＋１で表わす軌道によってだけ検出され
る映像点Ａに対する場合である。

中心までの距離が値旦である映像点は、例えば（ｑ＋２
）および（ｑ＋３ ）で表わす軌道の映像点Ｂのように
正確に２回検出される。

例えば（ｑ＋１ ） 、 （ｑ＋２ ）および（ｑ＋３
）で表わす軌道の映像点Ｃのように中心から値具に相
当する映像点は、正確に３回検出される。

一般的に、中心Ｚから値旦の距離に相当する映像点は、
フィ−ルド走査の半周期中、Ｎが整数だとするとＮ回検
出される。

Ｎが整数でなければ、観測される映像点は映像フィール
ドにおける位置に応じて次の高いか低い整数に応じて検
出される。

何故なら分数の検出は物理的に観測されないからである
。

最終的に、中心Ｚまでの距離が放射検出器の幅ｄよりも
小さくそしてその軌道がそれぞれの走査ごとに検出され
る現像点は、例えば中心Ｚそれ自身になる。

フィールド走査の半周期ごとの検出の最大ト 的にＮｍａｘ＝πＲ／ｄである。

本発明の課題は、走査信号のＳ／Ｎ比を改善するために
、適切な方法により過剰な走査、即ちフィールド走査の
半周期中映像フィールドにおけるある位置の多数回の走
査を利用することである。

この課題の解決につながる技術的思想は、相関方法にお
いて種々の線側期中に発生する走査信号を適当な係数を
乗算しそして適当な数の線周期数だけ遅延させた後、加
算することに在る。

所定の映像位置に対する放射検出器Ｓにおいて映像フィ
ールド走査の際発生される走査信号は時間に依存した信
号振幅ｆ（ｔ）として表現され、かつＮ個のｆ（ｔ）の
平均値を求める際に線型的に加算されるのに対して、相
互に独立な雑音信号は平方和として加算されるので、定
常的な雑音の場合には平均化に際してＳ／Ｎ比はＮ倍に
改善される。

それ故走査信号ｆ（ｔ）は、値ｒ＝是で表わす位置に対
して、経過中に位置が検出されるそれぞれ線周期にわた
る平均化により処理されて改善された出力信号Ｆ（１）
になる。

これは、Ｔを線周期、ｉを整数指数とすると次の式（１
）のようになる。

式（１）は旦の整数値にのみ使用できる。

しかしこれは、Ｂ−山中間値に対しても計算されねばな
らない。

何故なら映像フィールドにおける走査された位置は連続
になっているからである。

本発明は、次の確認を基にしている。

即ち第１図により整数Ｎについては（Ｎ’ＷＩ）くｒ’
４により規定される間隔において映像フィールドの個々
の位置がＮ回、しかし残りの位置は（Ｎ＋１）回フィー
ルド走査の半周期ごとに検出さへそして前記間隔におい
て増大する値ｒについて検出される位置の数Ｎ回と検出
される位置の数（Ｎ＋１）回との比は１からＯへ向う（
ｒが増大すると位置はより少く検出される）ことである
。

他方、Ｎは近似的にＲ＝＝Ｎとなるので前もって知られ
た整数値ではない。

そのため、式（１）により計算を行う場合加数の数はＮ
でなく、より大きなしかし周知の値ＮｒＴ１ａｘに限定
される。

この場合処理には伺ら貢献しないが追加的な雑音を導入
するような加数が加算される。

この欠点を除くために、個々の加数に、望ましくない加
数をしゃ断する重み係数を付加する。

整数値Ｒ−Ｎを有するそれぞれの位置においては、所属
の指数ｉ ＝ １〜Ｎを有する加数には重み１を与え、
一方所属の指数１＝（Ｎ＋１）〜Ｎｍａｘを有する残り
の加数には重み零を与える。

整数でない値−を有する位置の場合、”＝ （Ｎ＋ ｘ
）に値するとき（Ｎは整数、Ｘは端数（０＜Ｘ＜１）
）、所属の指数ｉ ＝ １〜Ｎを有する加数には重み１
を、そして所属の指数１＝（Ｎ＋２）〜Ｎｍａｘを有す
る加数には重み零が与えらへ一方所属の指数１＝（Ｎ＋
１）を有する加数には、零および１の間の重み、例えば
重みＸが与えられる。

式の次の組合せ（２）により規定されるランプ関数Ｓ（
ロ）について考える。

５６ｖ）＝ｏ ｗ＜ｏについて Ｓ（ｗ）＝Ｗ ＯくＷく１について ・・・・・・
（２）Ｓ （ｗ）＝ Ｉ Ｗ＞ １について ランプ関数を本発明の重み係数として使用するために、
独立変数Ｗとして例えばＷ＝（”−ｉ＋１★挿入する。

式（１）から、重み係数の挿入および限界値として７〜
Ｎｍａｘを有する全ての可能な旦の値に拡大することに
より、次の式（３）が生れる。

Ｆ（ｔ）＝ｉ’Σ：、：、Ｎｒｒｌａ Ｘ ｆ （１−
（ｓ−１）・Ｔ）・（Ｓ（旦−ｉ＋１）） ・・・・
・・（３）この式は、可変係数を有する擬似−ローパス
フィルタの関係である。

整数値旦については式（３）および（１）が等しくなる
ことが容易に確認できる。

走査信号ｆ（ｔ）を式（３）による出力信号Ｆ（１）と
なるように処理するために、第２図のブロック図の方法
を次に提案する。

この方法は同一出願人にかかる特願昭５０−１２６９５
７号に開示されるもので、以下のような構成を有するも
のである。

走査信号ｆ（ｔ）から、それぞれ所属の指数ｉ＝１〜Ｎ
を有する多数のＮ 個の主信号Ａｉが形ａＸ 戒される。

各主信号Ａｉは、主信号（Ａｉ−１）に対して線走査の
周期Ｔだけ遅延され、このとき、主信号Ａ１は当然例外
的に、遅延されない。

これは、第２図でＴ２・・・・・・Ｔｉ・・・・・・Ｔ
Ｎ で示す周知の遅延装置により行われる。

それぞれの遅延装置はＴだけ遅延を行う。

全ての遅延装置は直列に接続され、その内の最初のＴ２
は走査信号ｆ（ｏを供給される。

そこで遅延装置の出力ごとに主信号が生じそして主信号
は、ステップ状に増大する遅延がそれぞれ周期Ｔごとに
見出される一連の序列を形成する。

強制的なものではないが、最初の主信号Ａ１として走査
信号自身を使用すると合目的である。

これにより遅延装置Ｔ１を節約できる。したがってこれ
は第２図に示さない。

並列に接続されそして走査信号ｆ（ｔ）を供給され、そ
れぞれが新値の遅延Ｔ 、２Ｔ 、３Ｔ等を行う一群の
（ＮｍａＸ−１）個の全ての遅延装置により等価の方法
で主信号が発生され得る。

導線２０は走査信号ｆ（ｔ）を案内しそして導線２１は
所属の指数１＝（Ｎｒｒｌａｘ−１）を有する主信号を
案内する。

追加信号ｒは、時間に依存しそしてそれ故ｒ（ｔ）とし
て表わされる。

第２図において、追加信号ｒ（ｔ）が、導線２２を通し
てｑ・・・・・・Ｇｉ・・・・・・ＧＮｍａｘで表わす
多数のＮｍａＸ個の関数発生器に導入される。

それぞれの関数発生器には、それ故ｉ＝１およびｉ＝Ｎ
ｍａｘ間の所定の指数ｉが所属しそしてこの指数に対応
して所属の指数ｉを有する関数発生器において関数（ｒ
、 Ｓ （Ｒ−１＋１））が発生される。

値Ｒは、Ｒがｒ（ｔ）の最大値なので予め定められてい
る。

独立変数Ｗのランプ関数Ｓ＜＝＞は、前に定義した。

関数ｒから旦を形成すること、独立変数（Ｂ−−ｉ＋１
）およびランプ関数Ｓの発生ならびにランプ関数Ｓと値
貴との乗算は、現在の演算技術において周知であるため
詳しく説明しない。

関数発生器０１〜ＧＮｍａｘにおける追加信号ｒの処理
の結果として、所定の指数ｉを含む多数のＮｍａＸ個の
補助信号Ｂｉが関数数発生器Ｇｉ −ＧＮ の出力
側に生じる。

ここにＡｉ＝ｆ （ｔ−（ｉ−１）・Ｔ）、ＢｉＲ （豆・５（７−ｉ＋１））である。

導線２３は所属の指数ｉ＝（Ｎｍａｘ−１）を有する補
助信号を案内する。

同じ所属の指数ｉを有するそれぞれ補助信号および主信
号がそれぞれ周知の乗算装置において互に乗算され、し
かして同じ指数ｉが所属するそれぞれ中間信号Ｃ・＝Ｌ
−ｆ（を−（ｉ−１）−Ｔ）−（Ｓ（旦１−ｉＲ＋１）
）ヵ３８よあゎ６゜相応する多くのＮｍａｘ個の乗算装
置は第２図でＭ□・・・・・・Ｍｉ・・・・・・ＭＮｍ
ａｘで表わす。

所属の指数ｉ−（Ｎｍａｘｌ）を有する中間信号を案内
する導線は、第２図において２４を表わす。

全中間信号は、それ自体周知のアナログ−加算装置にお
いて加算さへしかして式（３）による出力信号Ｆ（１）
が形成される。

第２図で加算装置はΣで表わし、出力信号はＦ（、）は
導線２５を通して案内される。

このような、第２図に示す方法では、理想的な近似値が
得られるが、反面それ相応の費用が必要となる。

そこで、所要費用を低減しながら実用上十分である良好
な結果が得られる方法を以下に開示する。

すなわち、時間遅延装置の出力信号をその入力にフィー
ドバックすることにより、相互にステップ状に増加する
遅延時間を有する一連の信号が発生するという確認を基
にして、僅かな費用でありながら、式（３）の近似のみ
を達成する非常に有利な循環法を利用する方法を提案す
る。

これは、前記および第３図に示す方法において、指数ｉ
＝ ２〜Ｎｍａｘが属する全ての部材を指数Ｋを有す
る唯一の循環部材と置換するのと等価であり、関数発生
器ＧＫならびにＧ１で発生される関係は、循環法に適合
する。

この概念を利用することにより、第２図から次に述べる
第３図が生れる。

関数発生器Ｇ１および乗算装置Ｍ１ならびに加算装置Σ
は第２図から受継がれ、一方全ての関数発生器％〜ＧＮ
ｍａｘにかえて唯一の関数発生器ＧＫが、また全ての乗
算装置Ｍ２〜ＭＮｍａｘにかえて唯一の乗算装置ＭＫが
用いられている。

さらに全ての時間遅延装置Ｔ２〜Ｔ にかえて唯ｍ
ａｘ −の時間遅延装置ＴＫが用いられ、この時間遅延装置に
は走査信号ｆ（ｔ）ではなく出力信号Ｆ（１）が入力さ
れる。

しかして循環法が成立する。走査信号ｆ（１）はそれ故
導線３０を通して乗算装置Ｍ１にのみ案内さへ一方導線
３１〜３５は導線２１〜２５と同じ機能を有し、さらに
導線３６は出力信号Ｆ（、）を時間遅延装置ＴＫへ帰還
する。

関数発生器Ｇ１およびＧＫにおいて発生されるｒ（１）
に依存する関数をそれぞれｖｌおよび１Ｋにて表わすと
、第３図の方法は次の式（４）により次の数式で表わさ
れる。

Ｆ（ｔ）＝１に−Ｆ（ｔ−Ｔ）＋ｒ１・Ｆ（１）・・・
・・・（４）式（４）は周波敷写において−Ｔ ／ 、
ｆｆｎ （ｆ Ｋ ）の擬似一時定数およびｙ、／（ｉ
−ｆＫ）の減衰を有する１次の擬似−ローパスフィルタ
の特性を表わす。

瞬間的な走査信号を走査信号の旦の先の値との相関によ
り処理することが望まれるので、Ｔ一旦の時間が検出さ
れ、この時間および擬似一時定数の等化により値ｆＫ＝
ｅｘｐ（Ｒ）が得られそして値？、＝（１−ｅｘｐ（−
昔））が最大感度に対する条件として得られる（走査信
号における周波敷写において減衰−１）。

したがって、式（４）は次の式（５）として書き換えら
れ得る。

Ｆ（ｔ）＝ｅｘｐ （Ｒ） ・Ｆ （ｔ Ｔ ） ＋
（１−ｅＸｐ（−！−））・ｆ（１） ・
・・・・・（５）この式（５）による循環法において、
出力信号は、循環グループ（３６，ＴＫ、３１２ＭＫ、
３４゜Σ、第３図）を通してのそれぞれの通過の際に更
新されて周期Ｔだけ遅延され、係数ｅｘｐ（Ｒ）だけ減
衰される。

前の線走査の瞬間的な出力信号への貢献は、それ故遅延
時間の増加△ｔと共に素早く即ちｅｘｐ（−Ｎ・了）で
減少する。

つまり循環法が収斂する。

この循環法においてはそれ故一方では走査信号ｆ（１）
と一致する所属の指数ｉ＝１を有する主信号が発生さへ
他方では所属の指数ｉ≧２を有する全ての主信号を集め
そして包含する所属の指数Ｋを有する主信号が発生され
る。

所属の指数ｉを有する主信号に相当する時間遅延は、そ
れ自体周知の時間遅延装置ＴＫにより行われる。

追加信号ｒ（１）から、それぞれ関数発生器Ｇ１もしく
はＧＫにおいてそれぞれ補助信号が発生さへここでこれ
らの補助信号のそれぞれは関数ｙ、−（１ｅｘｐ（−Ｌ
））もしくはＳ’ ＝（ｅｘｐ（−ｉ））に比例する
。

現在の演算技術ではこのような値の代数組合せおよびこ
のような独立変数の指数関数値を形成することは周知で
あり、詳述する必要はないであろう。

同じ所属の指数ｉ＝１もしくはＫを有するそれぞれ補助
信号および主信号は、それぞれそれ自体周知の乗算装置
Ｍ１もしくはＭＫにおいて互に乗算され、しかして同じ
指数ｉ＝１もしくはＫか所属するそれぞれの中間信号が
発生される。

この中間信号は、それ自体周知の加算装置Σにおいて積
算され、しかして式（５）による出力信号Ｆ（１）が形
成される。

第２図における直接法および第３図における循環方法に
よる走査信号ｆ（１）の出力信号Ｆ（１）への処理を比
較すると、循環法は本質的に低費用になることが直ちに
分るはずである。

この利点は、価格のために信号処理が劣化する犠牲を払
っている。

直接法においてはそれぞれの線走査について専ら旦の線
走査が考慮され、しかも一定の重みで考慮され、これに
対し循環法においては基本的に全ての先の線走査が考慮
され、しかも減少する重みで考慮される。

直接法では、重みは関数ｒ（１）、即ち走査された位置
からフィールドの中心迄の距離に依存しない（式（３）
において因子畳は平均値形成に役立ち、重み付けには役
立たない）。

これに対し循環法においては、重みは関数ｒ（１）に依
存する。

循環法による走査信号の処理によってのＳ／Ｎ比の向上
は、本質的に費用を多く要する直接法による処理のよう
には良好でない結果をもたらす。

本発明による方法は、アナログ−関数装置を使用してア
ナログ−信号の処理に特に使用される。

しかしながら単独あるいは全信号をディジタル形式で処
理することならびに単独あるいは全装置をディジタルあ
るいは混合した関数方法で挿入することを除外するもの
でなく本発明に含まれるものである。

以上の本発明の詳細な説明中、簡単のために多数の主信
号に対する符号Ｎｍａｘの代りに符号Ｎを使用した。

【図面の簡単な説明】

第１図は本発明を説明するための映像面における走査方
法の図、第２図は直接法の実施例を説明するためブロッ
ク図および第３図は本発明にかかる循環法の実施例を説
明するブロック図、なお、図中の主な参照符号の対応は
次の通りである。 Ｚ・・映像面中心、Ｕ・・・映像面周辺、Ｒ・・・映像
面半径、ｒ・・・瞬間値、Ｓ・・・放射検出器、Ｐ・・
・検出器中心、ｆ（ｔ）・・・走査信号、ｒ（１）・・
・追加信号、Ｆ（１）・・・出力信号、Ｃ・・・・中間
信号、Ｔｉ・・・時間遅延装置、Ｇｉ・・・関数発生器
、Ｍｉ・・・乗算装置、Σ・・・加算器。

Claims (1)

1. 【特許請求の範囲】 １ 周期的な走査方法によって得られる時間依存の走査
   信号におけるＳ／Ｎ比を向上させるための方法であって
   、境界を有するとともに中心が定められたフィールドを
   走査装置によって線状に走査し、瞬間的走査位置からフ
   ィールド中心までの距離に対応しかつ距離の瞬間値ｒお
   よび最大値Ｒを含む時間依存の追加信号を発生し、順次
   所定の角度だけずらしながら線走査を行ない、この線走
   査の隣接する線がフィールド境界では互いに接しフィー
   ルド内部では部分的に重畳しまたフィールド中心では全
   面的に重畳するような前記線走査によって得られる時間
   依存の走査信号におけるＳ／Ｎ比を向上させるための前
   記方法において、走査信号ｆ（１）に基づいて一定の序
   列で発生されるＮ個の主信号列Ａ１．Ａ２・・・・・・
   Ａｉ・・・・・・ＡＮのうち１番日の主信号Ａ１は走査
   信号そのものとし、２ないしＮ番目の主信号列Ａ２・・
   曲Ａｉ・・・・・・ＡＮの順位ｉ番目の主信号Ａｉは、
   唯一の時間遅延装置ＴＫおよびこの時間遅延装置ＴＫと
   接続された出力信号Ｆ （ｔ）の帰還路によって循環法
   により発生されること、 前記追加信号ｒ （ｔ）を２つの関数発生器Ｇｌ、ＧＫ
   にそれぞれ印加し、前記１番目の主信号に対応して配設
   された一方の前記関数発生器Ｇ１から、関数 （１−ｅｘｐ（−百）〕 に比例する補助信号？□を発生させ、そして前記ｉ番目
   の主信号に対応して配設された他方の前記関数発生器Ｇ
   Ｋから、関数 （ｅｘｐ（−五）〕 に比例する補助信号ｆｆＫを発生させること、前記主信
   号のそれぞれに対応して配設された２個の乗算装置Ｍｌ
   、ＭＫによって前記補助信号？１．ｆＩＫと前記主信号
   とをそれぞれ乗算して各中間信号Ｃ１，ＣＫを発生させ
   ること、および、加算装置Σによって前記中間信号Ｃ，
   ，ＣＫを加算し出力信号Ｆ（ｔ）を発生させること、の
   各過程から成ることを特徴とする前記方法。