JPH06110919A - 近似直線計算方法 - Google Patents

近似直線計算方法

Info

Publication number
JPH06110919A
JPH06110919A JP27951292A JP27951292A JPH06110919A JP H06110919 A JPH06110919 A JP H06110919A JP 27951292 A JP27951292 A JP 27951292A JP 27951292 A JP27951292 A JP 27951292A JP H06110919 A JPH06110919 A JP H06110919A
Authority
JP
Japan
Prior art keywords
straight line
approximate straight
data
coordinate data
calculated
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.)
Pending
Application number
JP27951292A
Other languages
English (en)
Inventor
Osamu Hosoi
修 細井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tektronix Japan Ltd
Original Assignee
Sony Tektronix Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sony Tektronix Corp filed Critical Sony Tektronix Corp
Priority to JP27951292A priority Critical patent/JPH06110919A/ja
Publication of JPH06110919A publication Critical patent/JPH06110919A/ja
Pending legal-status Critical Current

Links

Landscapes

  • Complex Calculations (AREA)

Abstract

(57)【要約】 【目的】最小二乗法を利用して、複数の座標データの任
意の領域に対して高速に近似直線を計算する方法を実現
すること。 【構成】N個のXY座標データ(Xi,Yi)(iは、1
≦i≦Nの整数)から最小二乗法を用いて近似直線(Y
=A1X+A0)を迅速に計算する際に、 SXk = Σ(i=1, k)Xi SYk = Σ(i=1, k)Yi SXYk = Σ(i=1, k)XiYi SX2k = Σ(i=1, k)Xi2 (但し、kは、1≦
k≦Nの整数) で定義されたデータを夫々計算してメモリ26に予め記
憶する。次に、これらの記憶データに基づいて、オペレ
ータが指定した任意の領域(i=j〜k)の座標データ
に関する近似直線のパラメータA1及びA0を求める。 【効果】計算量が大幅に低減し、指定された任意の領域
の座標データの近似直線を迅速に計算出来る。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、最小二乗法を用いて複
数の座標データを近似する直線を高速に計算する方法に
関する。
【0002】
【従来技術】デジタル・オシロスコープや半導体測定装
置のように測定値をデジタルデータに変換し、座標デー
タとしてメモリに記憶し、これら記憶データを読み出し
て陰極線管等のスクリーン上に測定波形を表示する場
合、座標データの所望部分の近似直線を表示することが
行われている。この近似直線を表示する方法としては、
(a)データの2点(Xi,Yi)及び(Xj,Yj)を結
ぶ直線を表示する方法(b)データの1点を通る直線を
表示し、その直線の傾きを表示波形に合うようにオペレ
ータが調整する方法、(c)オペレータが目視で任意の
領域の波形に最も近似した直線を表示するように調整す
る方法、(d)最小二乗法に基づいて座標データを計算
して近似直線のパラメータを計算する方法等がある。
【0003】
【発明が解決しようとする課題】上述の(a)の方法で
は、2点のみのデータに依存するので適正な近似直線を
得ることは困難であり、(b)及び(c)の方法では、
オペレータによる人為的な誤差が生じる上に煩雑であ
る。また、(d)の従来の最小二乗法によれば、点(X
j,Yj)から点(Xk,Yk)までの(k−j+1)個の
データ点に近似した直線の式をY=A1X+A0とする
と、パラメータA1及びA0は、以下の式で表される。な
お、この式のΣ記号は、添字「i」についてjからkま
での和をとるものとする。 A1={Σ(1)Σ(XiYi)-Σ(Xi)Σ(Yi)}/{Σ(1)Σ(X
i2)-[ΣXi]2} A0={Σ(Xi2)Σ(Yi)-Σ(XiYi)Σ(Xi)}/{Σ(1)Σ(X
i2)-[ΣXi]2
【0004】この(d)の最小二乗法を用いる方法で
は、任意の領域(例えば点jから点kまで)が指定され
た場合、各点のデータ値をメモリから逐一読み出して計
算するので、読み出し動作及び計算動作の両方に多大の
時間がかかり、近似直線を迅速に求めることが困難であ
った。従って、本発明の目的は、最小二乗法を利用して
近似直線を求める際に、複数の座標データの任意の領域
に対して高速に正確な結果を求めることが可能な近似直
線計算方法を実現することである。
【0005】
【課題を解決する為の手段】本発明の近似直線計算方法
は、N個のXY座標データ(Xi,Yi)(iは、1≦i
≦Nの整数)から最小二乗法を用いて近似直線(Y=A
1X+A0)を迅速に計算する。まず、 SXk = Σ(i=1, k)Xi SYk = Σ(i=1, k)Yi SXYk = Σ(i=1, k)XiYi SX2k = Σ(i=1, k)Xi2 (但し、kは、1≦k≦Nの整数) で定義されたデータを夫々計算してメモリに予め記憶す
る。次に、これらの記憶データに基づいて、オペレータ
が指定した任意の領域(i=j〜k)の座標データに関
する近似直線のパラメータA1及びA0を求める。
【0006】
【実施例】図2は、本発明を適用するのに好適な素子特
性測定装置の構成を示すブロック図である。CPU10
はバス12を介して、ROM(リード・オンリ・メモ
リ)14内のプログラム・データ又は入力装置16の操
作による入力データを受け取り、これらのデータに従い
バス12を介して他の構成ブロックを制御する。入力装
置16は、押しボタン及び制御ノブ等を有する操作パネ
ル、表示画面及び項目選択用押しボタン等を組み合わせ
た入力装置、又は表示画面を利用したタッチスクリーン
その他の周知の入力装置で良い。測定回路18は、CP
U(中央処理装置)10により制御されて可変電圧源、
可変電流源、電圧測定回路及び電流測定回路として機能
すると共に測定点の接地又は開放を制御する測定ユニッ
トを少なくとも3個含んでいる。CPU10に制御され
て、測定回路18は、第1セットの出力端子からDUT
(被測定素子)20に試験信号を供給する。図2では、
DUT20はバイポーラ・トランジスタであるが、その
他の種々の素子でも良い。例えば、このトランジスタの
コレクタ・エミッタ間電圧VCE−コレクタ電流IC特性
を測定するためには、コレクタに定電圧、ベースに階段
波状バイアス電流を供給し、エミッタを接地する。これ
により、バイアス電流の各ステップ毎にVCE及びICに
夫々相当するX及びYアナログ信号が測定回路18の第
2セットの出力端子から出力される。
【0007】ADC(アナログ・デジタル変換器)22
は、CPU10により制御されるクロック発生器24か
らのサンプリング・クロックfsのタイミングで、X及
びYアナログ信号をX及びYデジタル・データに夫々変
換する。ADC22からのX及びYデジタル・データ
は、CPU10によりアドレス制御されてメモリ26に
記憶される。表示制御回路28は、陰極線管の如き表示
装置30上に、メモリ26から読み出されたX及びYデ
ジタル・データを夫々水平軸用及び垂直軸用データとし
て、バイポーラ・トランジスタであるDUT20のコレ
クタ・エミッタ間電圧VCE−コレクタ電流IC特性を表
示する。この様な素子特性測定装置では、この特性測定
の他に、DUT20の各端子に接続された測定ユニット
の機能を制御することにより、ICBO−VCBO特性、β−
IC特性、IC−VBE特性その他の多様な特性を測定して
表示できる。
【0008】また、表示制御回路28は、周知の文字発
生器を含んでおり、表示に必要な文字、カーソル、直線
等を表示することが出来る。本発明の近似直線表示方法
を利用する場合には、直線Y=A1X+A0のパラメータ
A1及びA0を後述する式に基づいて計算し、この直線の
式に従って文字発生器により従来より周知の方法で直線
を表示装置30のスクリーンに表示する。
【0009】図1は、本発明の近似直線計算方法の一実
施例の手順を示す流れ図である。先ず、ステップ100
において、N個の点i(1≦i≦N)に対応するXY座
標データ(Xi,Yi)をメモリに記憶する。次のステッ
プ102では、 SXk=Σ(i=1, k)Xi ・・・ (1) をk=1〜NまでのN個の整数についての値を計算して
メモリに記憶する。なお、Σ(i=1, k)は、添字i=1
〜Kまでのk個の整数について総ての項の和をとること
を意味する。また、SX1=X1であり、SXm=SXm-1
+Xmであるので、X1から順次加算することにより、容
易にSXkを計算することが出来る。次のステップ10
4では、 SYk=Σ(i=1, k)Yi ・・・ (2) をk=1〜NまでのN個の整数についての値を計算して
メモリに記憶する。この計算も(1)式の場合と同様に、
SYk=Y1から順次加算することにより容易に得ること
が出来る。次のステップ106では、 SXYk=Σ(i=1, k)XiYi ・・・ (3) をk=1〜NまでのN個の整数についての値を計算して
メモリに記憶する。この計算も上述のように、SXY1
=X1Y1から順次加算することにより、容易に計算出来
る。次のステップ108では、 SX2k=Σ(i=1, k)Xi2 ・・・ (4) をK=1〜NまでのN個の整数についての値を計算して
メモリに記憶する。これも上述のように、SX21=X1
2から順次加算することにより容易に計算出来る。
【0010】以上の(1)〜(4)までの計算結果をメモリ
に予め記憶しておくことにより、オペレータが指定する
任意の領域の座標データに近似した直線の式を極めて高
速且つ正確に計算することが出来る。すなわち、Σ(i=
j, k)Xi=SXk−SXj-1と表され、この式の右辺の
項の値は既にメモリに記憶されているので、計算が極め
て容易となる。これら予め記憶されたデータを用いれ
ば、ステップ110において、オペレータが指定した任
意の領域(例えば、データの添字j〜kまでの領域)の
座標データに近似した直線の式をY=A1X+A0とする
と、この直線を決定する2つのパラメータA1及びA0
は、以下の式で求めることが出来る。
【数1】
【数2】
【0011】上述のように、(1)式から(4)式までの計
算を予め実行してデータをメモリに記憶しておけば、オ
ペレータが任意に指定した領域の座標データに近似する
直線を即座に求めることが出来る。従って、オペレータ
が指定領域を変更した場合でもその変化する領域に迅速
に追従して近似直線を表示する、所謂リアルタイム応答
が可能になる。
【発明の効果】本発明の近似直線計算方法によれば、最
小二乗法の高速計算に必要なデータを予め計算し、メモ
リに準備しておくことにより、従来より遥かに高速に任
意の領域の座標データに関する近似直線を計算すること
が可能になる。
【図面の簡単な説明】
【図1】本発明の一実施例の手順を示す流れ図である。
【図2】本発明を利用するのに好適な素子特性測定装置
の構成を示すブロック図である。
【符号の説明】
10 CPU 16 入力装置 22 ADC(アナログ・デジタル変換器) 26 メモリ 28 表示制御回路 30 表示装置

Claims (1)

    【特許請求の範囲】
  1. 【請求項1】 N個のXY座標データ(Xi,Yi)(i
    は、1≦i≦Nの整数)から最小二乗法を用いて近似直
    線(Y=A1X+A0)を計算する近似直線計算方法にお
    いて、 SXk = Σ(i=1, k)Xi SYk = Σ(i=1, k)Yi SXYk = Σ(i=1, k)XiYi SX2k = Σ(i=1, k)Xi2 (但し、kは、1≦k≦Nの整数) で定義されたデータを夫々計算してメモリに予め記憶
    し、これらの記憶データに基づいて上記近似直線のパラ
    メータA1及びA0を求めることを特徴とする近似直線計
    算方法。
JP27951292A 1992-09-24 1992-09-24 近似直線計算方法 Pending JPH06110919A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP27951292A JPH06110919A (ja) 1992-09-24 1992-09-24 近似直線計算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP27951292A JPH06110919A (ja) 1992-09-24 1992-09-24 近似直線計算方法

Publications (1)

Publication Number Publication Date
JPH06110919A true JPH06110919A (ja) 1994-04-22

Family

ID=17612074

Family Applications (1)

Application Number Title Priority Date Filing Date
JP27951292A Pending JPH06110919A (ja) 1992-09-24 1992-09-24 近似直線計算方法

Country Status (1)

Country Link
JP (1) JPH06110919A (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58105372A (ja) * 1981-12-16 1983-06-23 Fujitsu Ltd 最小二乗法による線図作図方式
JPS623428B2 (ja) * 1979-08-09 1987-01-24 Tokyo Shibaura Electric Co
JPS6336387A (ja) * 1986-07-30 1988-02-17 Nec Corp グラフ表示方式
JPS63259693A (ja) * 1987-04-17 1988-10-26 ソニー株式会社 音源正規化方法
JPH01189765A (ja) * 1988-01-25 1989-07-28 Hisaya Kotajima 関数演算装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS623428B2 (ja) * 1979-08-09 1987-01-24 Tokyo Shibaura Electric Co
JPS58105372A (ja) * 1981-12-16 1983-06-23 Fujitsu Ltd 最小二乗法による線図作図方式
JPS6336387A (ja) * 1986-07-30 1988-02-17 Nec Corp グラフ表示方式
JPS63259693A (ja) * 1987-04-17 1988-10-26 ソニー株式会社 音源正規化方法
JPH01189765A (ja) * 1988-01-25 1989-07-28 Hisaya Kotajima 関数演算装置

Similar Documents

Publication Publication Date Title
US5397981A (en) Digital storage oscilloscope with automatic time base
US6681191B1 (en) Frequency domain analysis system for a time domain measurement instrument
EP0701138B1 (en) Method of signal analysis employing histograms
JPH10333798A (ja) 測定信号波形スケーリング・システムおよびその方法
JPH0543993B2 (ja)
US20020103651A1 (en) Voice-responsive command and control system and methodology for use in a signal measurement system
US4801873A (en) Waveform measuring apparatus with marker zone displaying function
US4761640A (en) Positioning cursors at specific points on a waveform display
JP2533813B2 (ja) スペクトラム・アナライザ
EP1111396B1 (en) Frequency domain analysis system for a time domain measurement instrument
JPH09510783A (ja) 波形分析のための方法および装置
JPH06110919A (ja) 近似直線計算方法
JPH06258356A (ja) 信号波形の最大,最小値測定装置
GB2132058A (en) Waveform display apparatus
JP2807709B2 (ja) 期待値表示機能を有する掃引同調試験装置と期待値表示方法
JPH0360439B2 (ja)
US6646428B2 (en) Sweep synchronization testing
JPH05240906A (ja) 素子特性測定装置のカーソル測定方法
JP2928867B2 (ja) 測定装置
JPH06281706A (ja) 素子特性測定装置の測定方法
JP3071875B2 (ja) Icテスト装置
Guido et al. Automatic pulse parameter determination with the computer augmented oscilloscope system
JPS6281573A (ja) カ−ソル測定機能を有するオシロスコ−プ
JP3410035B2 (ja) 傾斜カーソル表示方法
JP2592744B2 (ja) 素子特性測定装置