JPH0783977A - Response characteristics measuring apparatus - Google Patents
Response characteristics measuring apparatusInfo
- Publication number
- JPH0783977A JPH0783977A JP23022793A JP23022793A JPH0783977A JP H0783977 A JPH0783977 A JP H0783977A JP 23022793 A JP23022793 A JP 23022793A JP 23022793 A JP23022793 A JP 23022793A JP H0783977 A JPH0783977 A JP H0783977A
- Authority
- JP
- Japan
- Prior art keywords
- input
- output signal
- input signal
- time
- signal
- 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
Abstract
Description
【０００１】[0001]
【産業上の利用分野】本発明は、計測対象へ入力された
入力信号とその計測対象から出力された出力信号とに基
づいてその計測対象の伝達関数もしくはインパルス応答
関数を求める応答特性計測装置に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a response characteristic measuring apparatus for obtaining a transfer function or an impulse response function of an object to be measured based on an input signal input to the object to be measured and an output signal output from the object to be measured. .
【０００２】[0002]
【従来の技術】従来より、計測対象の入力信号と出力信
号とのクロススペクトルを求めることにより、計測対象
の伝達関数を求め、あるいはその伝達関数を逆フーリエ
変換することによりその計測対象のインパルス応答特性
を求める手法が知られており、また広く用いられてい
る。2. Description of the Related Art Conventionally, a transfer function of a measurement object is obtained by obtaining a cross spectrum of an input signal and an output signal of the measurement object, or an impulse response of the measurement object is obtained by performing an inverse Fourier transform of the transfer function. Techniques for obtaining characteristics are known and widely used.
【０００３】図７は、伝達関数ないしインパルス応答関
数の求め方の説明図である。計測対象１０への入力信号
をｘ（ｔ）とする。ここで、ｔは時刻を表わす。その入
力信号ｘ（ｔ）のフーリエ変換をＸ（ω）とする。ωは
角周波数である。またこの計測対象１０のインパルス応
答関数をｈ（ｔ），そのフーリエ変換すなわち伝達関数
をＨ（ｋ）とする。FIG. 7 is an explanatory diagram of how to obtain a transfer function or an impulse response function. The input signal to the measurement target 10 is x (t). Here, t represents time. Let the Fourier transform of the input signal x (t) be X (ω). ω is the angular frequency. Further, the impulse response function of the measurement object 10 is h (t), and its Fourier transform, that is, the transfer function is H (k).
【０００４】この計測対象１０からは雑音ｎ（ｔ）が混
入した出力信号ｙ（ｔ）が出力される。この出力信号ｙ
（ｎ）の離散フーリエ変換をＹ（ｋ）とする。このとき Ｙ（ｋ）＝Ｘ（ｋ）・Ｈ（ｋ）＋Ｎ（ｋ） ……（１） が成立する。両辺にＸ* （ｋ）（* は複素共役を表わ
す）を掛けると、 Ｘ* （ｋ）・Ｙ（ｋ） ＝Ｘ* （ｋ）・Ｘ（ｋ）・Ｈ（ｋ）＋Ｘ* （ｋ）・Ｎ（ｋ） ……（２） 入力信号ｘ（ｔ）と雑音ｎ（ｔ）は無相関であるため、
多数回に亘って（２）式の演算を行い、それらの平均値
＜ ＞を求めると、＜Ｘ* （ｋ）Ｎ（ｋ）＞は零に近づ
き、 ＜Ｘ* （ｋ）・Ｙ（ｋ）＞＝＜Ｘ* （ｋ）・Ｘ（ｋ）＞・Ｈ（ｋ） ……（３） が成立する。したがって、式 Ｈ（ｋ）＝＜Ｘ* （ｋ）・Ｙ（ｋ）＞／＜Ｘ* （ｋ）・Ｘ（ｋ）＞ ……（４） により計測対象１０の伝達関数Ｈ（ｋ）が求められ、ま
たこの求められた伝達関数Ｈ（ｋ）を逆フーリエ変換す
ることにより計測対象１０のインパルス応答関数ｈ
（ｔ）が求められる。An output signal y (t) mixed with noise n (t) is output from the measuring object 10. This output signal y
Let the discrete Fourier transform of (n) be Y (k). At this time, Y (k) = X (k) · H (k) + N (k) (1) holds. Multiplying both sides by X * (k) ( * represents a complex conjugate), X * (k) * Y (k) = X * (k) * X (k) * H (k) + X * (k) N (k) (2) Since the input signal x (t) and the noise n (t) are uncorrelated,
When the equation (2) is calculated over a number of times and the average value <> thereof is obtained, <X * (k) N (k)> approaches zero and <X * (k) · Y (k )> = <X * (k) .X (k)>. H (k) (3) holds. Therefore, the transfer function H (k) of the measurement object 10 is calculated by the formula H (k) = <X * (k) · Y (k)> / <X * (k) · X (k)> (4). The impulse response function h of the measurement object 10 is obtained by inverse Fourier transforming the obtained transfer function H (k).
(T) is required.
【０００５】図８は、フーリエ変換のための時間窓を表
わした図である。（Ａ），（Ｂ）はそれぞれ入力信号，
出力信号の時間窓であるが、ここでは同一の時間窓を表
わしている。実際の演算にあたっては、入力信号ｘ
（ｔ），出力信号ｙ（ｔ）を所定の微小時間幅Δｔ毎に
サンプリングすることにより入力信号ｘ（ｔ），出力信
号ｙ（ｔ）の各サンプリング値ｘ（ｎ），ｙ（ｎ）を求
め、これらのサンプリング値ｘ（ｎ），ｙ（ｎ）に基づ
いてフーリエ変換等の演算が行われる。ここで、ｎは、
微小時間間隔Δｔ毎にサンプリングした際の各サンプリ
ング点の時間的な順序を表わしており、ここでは、図７
に示すようにｎ＝０，１，…，Ｎ−１のＮ個のサンプリ
ング点のサンプリング値ｘ（ｎ），ｙ（ｎ）、ｎ＝０，
１，…，Ｎ−１に基づいた演算が行われるものとする。FIG. 8 is a diagram showing a time window for Fourier transform. (A) and (B) are input signals,
Although it is a time window of the output signal, the same time window is shown here. In the actual calculation, input signal x
(T), the output signal y (t) is sampled for each predetermined minute time width Δt to obtain each sampling value x (n), y (n) of the input signal x (t) and the output signal y (t). Then, the Fourier transform or the like is performed based on the obtained sampling values x (n) and y (n). Where n is
FIG. 7 shows the temporal order of each sampling point when sampling is performed at every minute time interval Δt.
, N−1, sampling values x (n), y (n), n = 0, of N sampling points of n−1,
It is assumed that the calculation based on 1, ..., N-1 is performed.
【０００６】以下、このサンプリングされた入力信号，
出力信号を、単に入力信号ｘ（ｎ），出力信号ｙ（ｎ）
と称する。それらのフーリエ変換値は、前述と同様にＸ
（ｋ），Ｙ（ｋ）で表わされるが、この場合、ｋは離散
的な角周波数として理解される。また、図８に示すよう
に、時間的に一定な方形窓Ｗ1 を用いることに代えて、
入力信号ｘ（ｎ），出力信号ｙ（ｎ）に時間的に重みが
変化する時間窓Ｗ2 を適用することも広く知られてい
る。この時間的に変化する時間窓Ｗ2 としては、例え
ば、ハニング窓，ハミング窓，リース窓，ブラックマン
ハリス窓等、多数の関数形が知られている。Hereinafter, this sampled input signal,
The output signal is simply the input signal x (n) and the output signal y (n)
Called. Those Fourier transform values are X
Although represented by (k) and Y (k), k is understood as a discrete angular frequency in this case. Further, as shown in FIG. 8, instead of using the rectangular window W 1 which is constant in time,
It is also widely known to apply a time window W 2 whose weight changes with time to the input signal x (n) and the output signal y (n). As the time window W 2 that changes with time, a number of function forms such as a Hanning window, a Hamming window, a lease window, a Blackman Harris window, etc. are known.
【０００７】[0007]
【発明が解決しようとする課題】従来は、上述したよう
に、所定の時間幅の入力信号ｘ（ｔ），出力信号ｙ
（ｔ）を切出し、即ち所定の微小時間幅Δｔ毎にサンプ
リングされたＮ個のサンプリング値ｘ（ｎ），ｙ（ｎ）
を求め、そのように切り出した、即ちそのようにサンプ
リングされたＮ個の入力信号ｘ（ｎ），出力信号ｙ
（ｎ）に基づいて、もしくは切り出した入力信号ｘ
（ｎ），出力信号ｙ（ｎ）に例えばハニング窓等の時間
窓に応じた重み演算を施した信号に基づいて、計測対象
の伝達関数ないしインパルス応答関数が求められてい
る。Conventionally, as described above, the input signal x (t) and the output signal y having a predetermined time width are used.
(T) is cut out, that is, N sampling values x (n), y (n) sampled at each predetermined minute time width Δt.
, And the N input signals x (n) and output signals y thus cut out, that is, so sampled
Input signal x based on (n) or cut out
(N), the transfer function or impulse response function to be measured is obtained based on a signal obtained by performing weighting calculation on the output signal y (n) according to a time window such as a Hanning window.
【０００８】この伝達関数ないしインパルス応答関数を
高精度に求めるためには、インパルス応答の継続時間よ
りも十分に長い時間窓（方形窓を含む）を用いることが
必要であるが、例えば振動の減衰の少ない金属構造物や
広い空間を計測対象としてその金属構造物の振動特性や
その空間の音響特性を求める場合等、インパルス応答の
継続時間が非常に長い場合、十分に長い時間窓を用いる
という条件を満たそうとうするとサンプリング数Ｎが膨
大なものとなり、演算が不可能となり、もしくはその演
算に膨大な時間を要することとなり、十分に長い時間窓
を用いるという条件を満たすことは実際上困難となる。
この場合、求められた伝達関数ないしインパルス応答関
数の誤差が非常に大きくなってしまうという問題があ
る。In order to obtain this transfer function or impulse response function with high accuracy, it is necessary to use a time window (including a square window) that is sufficiently longer than the duration of the impulse response. For example, vibration damping If the duration of the impulse response is very long, such as when the vibration characteristics of the metal structure with a small number of objects or a large space are to be measured and the acoustic characteristics of the space are to be obtained, a condition that a sufficiently long time window is used To satisfy the above condition, the number of samplings N becomes enormous, the operation becomes impossible, or the operation requires an enormous amount of time, which makes it practically difficult to satisfy the condition of using a sufficiently long time window.
In this case, there is a problem that the error of the obtained transfer function or impulse response function becomes very large.
【０００９】本発明は、上記事情に鑑み、特にインパル
ス応答の継続時間に比べ十分に長い時間窓を用いること
ができない場合に、従来と比べ高精度の伝達関数、イン
パルス応答関数を求めることのできる応答特性計測装置
を提供することを目的とする。In view of the above circumstances, the present invention can obtain a transfer function and an impulse response function with higher accuracy than before, especially when a time window that is sufficiently longer than the duration of the impulse response cannot be used. It is an object to provide a response characteristic measuring device.
【００１０】[0010]
【課題を解決するための手段】上記目的を達成する本発
明の応答特性計測装置は、計測対象へ入力された入力信
号とその計測対象から出力された出力信号とに基づいて
その計測対象の伝達関数もしくはインパルス応答関数を
求める応答特性計測装置において、 （１）所定の時刻ｔ0 から所定の第１の時間幅Ｔ1 の間
の入力信号を切り出すとともに上記第１の時間Ｔ1 より
も短い所定の第２の時間幅をＴ2 としたときに時刻ｔ0
＋Ｔ2 から時刻ｔ0 ＋Ｔ1 までの間の入力信号を零に置
換することにより置換入力信号を生成する入力信号切出
し置換手段 （２）上記所定の時刻ｔ0 から第１の時間幅Ｔ1 の間の
出力信号を切り出す出力信号切出し手段 （３）上記入力信号切出し置換手段により生成された置
換入力信号と出力信号切出し手段により切出された出力
信号とに基づいて上記計測対象の伝達関数もしくはイン
パルス応答関数を求める応答特性演算手段 を備えたことを特徴とするものである。A response characteristic measuring apparatus of the present invention that achieves the above object transmits a measurement target based on an input signal input to the measurement target and an output signal output from the measurement target. In a response characteristic measuring apparatus that obtains a function or an impulse response function, (1) cut out an input signal within a predetermined first time width T 1 from a predetermined time t 0 and a predetermined time shorter than the first time T 1. When the second time width of T 2 is T 2 , the time t 0
Input signal cutout replacement means for generating a replacement input signal by replacing the input signal from + T 2 to time t 0 + T 1 with zero. (2) From the predetermined time t 0 to the first time width T 1 Output signal cutout means for cutting out an output signal between (3) The transfer function or impulse of the measurement object based on the replacement input signal generated by the input signal cutout replacement means and the output signal cut out by the output signal cutout means It is characterized in that a response characteristic calculating means for obtaining a response function is provided.
【００１１】[0011]
【作用】インパルス応答の継続時間が長い場合、入力信
号に対する応答が時間窓からはみ出してしまうことによ
り計測誤差が生じる。この場合に、時間窓内の入力信号
であっても、その時間窓の後端側の入力信号成分に対す
る応答ほどその時間窓からはみ出してしまう成分が多く
なる。When the duration of the impulse response is long, the response to the input signal is out of the time window, which causes a measurement error. In this case, even with respect to the input signal within the time window, the more the response to the input signal component at the rear end side of the time window, the more the component is out of the time window.
【００１２】本発明は、上記の点に着目し、入力信号に
対する時間窓を、出力信号に対する時間窓と比べ時間的
に初期の領域だけに限定することに想到し、後述するよ
うにシミュレーションにより従来と比べ誤差が低減化さ
れることを確認し、本発明を成すに至ったのである。本
発明によれば、インパルス応答のｔ0 ＋Ｔ2 〜ｔ0 ＋Ｔ
1 の区間はパイアス誤差（偏り誤差）を生じることな
く、正確に求めることができる。ランダム信号源を用い
て本発明により求めたインパルス応答関数は、ランダム
誤差しか含まないから、クロススペクトル法における平
均回数を大きくするだけで、誤差を小さくすることがで
きる。In consideration of the above points, the present invention contemplates that the time window for the input signal is limited only to the initial region in terms of time as compared with the time window for the output signal. It was confirmed that the error was reduced as compared with, and the present invention was accomplished. According to the present invention, the impulse response from t 0 + T 2 to t 0 + T.
The section 1 can be accurately determined without causing a bias error (bias error). Since the impulse response function obtained by the present invention using a random signal source includes only a random error, the error can be reduced only by increasing the average number in the cross spectrum method.
【００１３】尚、Ｔ2 の値は、インパルス応答の継続時
間が時間窓長に比して長いときには小さな値にすること
が好ましい。それは、信号源側の時間窓内の窓内の最後
の応答が、出力信号を切り取った時間窓の中に入ってい
るとき、バイアス誤差を含まない演算が可能になるため
であり、それよりもインパルス応答の継続時間が長いと
きにはバイアス誤差が生じるためである。The value of T 2 is preferably small when the duration of the impulse response is longer than the time window length. This is because when the last response in the window within the time window on the signal source side is within the time window obtained by clipping the output signal, it is possible to perform calculation without bias error. This is because a bias error occurs when the duration of the impulse response is long.
【００１４】[0014]
【実施例】以下、本発明の実施例について説明する。図
１は、本発明の一実施例の応答特性計測装置の構成ブロ
ック図である。ここでは特に限定されない所定の伝達系
（計測対象）２０の応答特性が計測される。この伝達系
２０に入力信号ｘ（ｔ）が入力されるとともにその伝達
系２０から出力信号ｙ（ｔ）が出力されている。これら
入力信号ｘ（ｔ），出力信号ｙ（ｔ）は、応答特性計測
装置３０に入力される。EXAMPLES Examples of the present invention will be described below. FIG. 1 is a configuration block diagram of a response characteristic measuring apparatus according to an embodiment of the present invention. Here, the response characteristic of a predetermined transmission system (measurement target) 20, which is not particularly limited, is measured. The input signal x (t) is input to the transfer system 20, and the output signal y (t) is output from the transfer system 20. The input signal x (t) and the output signal y (t) are input to the response characteristic measuring device 30.
【００１５】応答特性計測装置３０に入力された、入力
信号ｘ（ｔ），出力信号ｙ（ｔ）は、それぞれＡ／Ｄ変
換器３２，３４に入力され、それぞれ微小時間間隔Δｔ
毎にサンプリングされ、ディジタルの入力信号ｘ
（ｎ），出力信号ｙ（ｎ）が生成される。Ａ／Ｄ変換器
３１から出力された入力信号ｘ（ｎ）は、入力信号切出
し置換手段３３に入力され、後述するようにして置換入
力信号ｚ（ｎ）が求められる。この置換入力信号ｚ
（ｎ）は、応答特性演算手段３５に入力される。The input signal x (t) and the output signal y (t) input to the response characteristic measuring device 30 are input to the A / D converters 32 and 34, respectively, and each of them has a minute time interval Δt.
Digital input signal x sampled for each
(N), the output signal y (n) is generated. The input signal x (n) output from the A / D converter 31 is input to the input signal cutout replacement unit 33, and the replacement input signal z (n) is obtained as described later. This replacement input signal z
(N) is input to the response characteristic calculation means 35.
【００１６】またＡ／Ｄ変換器３２から出力された出力
信号ｙ（ｎ）は、出力信号切出し手段３４に入力され、
この出力信号切出し手段３４において所定の時間幅Ｔ1
の部分、即ち、ｙ（ｎ），ｎ＝０，１，…，Ｎ−１の部
分だけ切出されて応答特性演算手段３５に入力される。
図２は、図１に示す入力信号切出し置換手段３３、出力
信号切出し手段３４の説明のための時間窓を表わした図
である。The output signal y (n) output from the A / D converter 32 is input to the output signal cutting means 34,
The output signal cutting means 34 has a predetermined time width T 1
, That is, y (n), n = 0, 1, ..., N−1 is cut out and input to the response characteristic calculation means 35.
FIG. 2 is a diagram showing time windows for explaining the input signal cutout replacing means 33 and the output signal cutout means 34 shown in FIG.
【００１７】図２（Ａ）は、入力信号ｘ（ｎ）に適用さ
れる時間窓であり、ｘ（０）〜ｘ（Ｍ−１）（但し、Ｍ
＜Ｎ）の間入力信号ｘ（ｎ）が切出され、ｘ（Ｍ）〜ｘ
（Ｎ−１）は零に置換され、全体としてＮ個のサンプリ
ング点からなる置換入力信号ｚ（ｎ），ｎ＝０，１，
…，Ｍ−１，Ｍ，…，Ｎ−１が求められる。図２（Ｂ）
は、出力信号ｙ（ｎ）に適用される時間窓であり、ｙ
（０）〜ｙ（Ｎ−１）の間の出力信号ｙ（ｎ）が切出さ
れる。FIG. 2A is a time window applied to the input signal x (n), where x (0) to x (M-1) (where M is
While <N), the input signal x (n) is cut out and x (M) to x
(N-1) is replaced by zero, and the replacement input signal z (n), n = 0, 1, consisting of N sampling points as a whole.
, M-1, M, ..., N-1 are obtained. Figure 2 (B)
Is the time window applied to the output signal y (n), y
The output signal y (n) between (0) and y (N-1) is cut out.
【００１８】図１に示す応答特性演算手段３５では、上
記のようにして求められた置換入力信号ｚ（ｎ）のフー
リエ変換関数Ｚ（ｋ）、および上記のようにして切出さ
れた出力信号ｙ（ｎ）のフーリエ変換関数Ｙ（ｋ）が求
められ、入出力間のクロススペクトルＺ* （ｋ）・Ｙ
（ｋ）、置換入力信号のパワースペクトルＺ* （ｋ）・
Ｚ（ｋ）が求められる。上記の過程が多数回繰り返され
た後、伝達系２０の伝達関数Ｈ（ｋ），インパルス応答
関数ｈ（ｎ）が、 Ｈ（ｋ）＝＜Ｚ* （ｋ）・Ｙ（ｋ）＞／＜Ｚ* （ｋ）・Ｚ（ｋ）＞ ……（５） ｈ（ｎ）＝Ｆ-1｛Ｈ（ｋ）｝ ……（６） 但しＦ-1は逆フーリエ変換を表わす として求められる。In the response characteristic calculation means 35 shown in FIG. 1, the Fourier transform function Z (k) of the replacement input signal z (n) obtained as described above and the output signal cut out as described above. The Fourier transform function Y (k) of y (n) is obtained, and the cross spectrum Z * (k) · Y between input and output is obtained.
(K), the power spectrum Z * (k) of the replacement input signal
Z (k) is required. After the above process is repeated many times, the transfer function H (k) and the impulse response function h (n) of the transfer system 20 become H (k) = <Z * (k) · Y (k)> / < Z * (k) · Z (k)> (5) h (n) = F −1 {H (k)} (6) where F −1 is obtained as an inverse Fourier transform.
【００１９】次に、本発明の応答特性計測装置における
誤差と従来法においてハニング窓を適用した場合の誤差
との比較を表わすシミュレーション結果について説明す
る。図３〜図６は、このシミュレーション結果を表わし
た図である。図３〜図６は、時間窓長Ｎを、インパルス
応答が６０ｄＢ減衰するまでの時間に相当するサンプル
数Ｔの、それぞれ４倍、２倍、１倍、２／３倍のとき
の、Ｍの値による計測誤差の変化を調べた結果である。
各図中ａは、図８に示すような従来法においてハニング
窓を適用した場合の誤差、ｂは本発明を適用したときの
誤差を示している。Next, a simulation result showing a comparison between an error in the response characteristic measuring apparatus of the present invention and an error when the Hanning window is applied in the conventional method will be described. 3 to 6 are diagrams showing the results of this simulation. 3 to 6 show the time window length N of M when the number of samples T corresponding to the time until the impulse response is attenuated by 60 dB is 4 times, 2 times, 1 time, and 2/3 times, respectively. It is the result of examining the change of the measurement error due to the value.
In each figure, a indicates an error when the Hanning window is applied in the conventional method as shown in FIG. 8, and b indicates an error when the present invention is applied.
【００２０】これらの図により、従来法のように、入出
力信号ともに長さＮのハニング窓で切り取ったときに比
べて、Ｎ／Ｔが４以下のときに、適当なＭの値を選ぶこ
とにより誤差が減少することがわかる。From these figures, select an appropriate value of M when N / T is 4 or less, as compared with the case where both the input and output signals are cut by the Hanning window of length N as in the conventional method. It can be seen that the error decreases.
【００２１】[0021]
【発明の効果】以上説明したように、本発明によれば、
インパルス応答の継続時間に対し、十分長い時間窓長を
確保することができない場合に、従来法よりも高精度
に、伝達関数ないしインパルス応答関数が求められる。As described above, according to the present invention,
When it is not possible to secure a sufficiently long time window length with respect to the duration of the impulse response, the transfer function or impulse response function is obtained with higher accuracy than the conventional method.
【図１】本発明の一実施例の応答特性計測装置の構成ブ
ロック図である。FIG. 1 is a configuration block diagram of a response characteristic measuring apparatus according to an embodiment of the present invention.
【図２】時間窓を表わした図である。FIG. 2 is a diagram showing a time window.
【図３】シミュレーション結果を表わした図である。FIG. 3 is a diagram showing a simulation result.
【図４】シミュレーション結果を表わした図である。FIG. 4 is a diagram showing a simulation result.
【図５】シミュレーション結果を表わした図である。FIG. 5 is a diagram showing a simulation result.
【図６】シミュレーション結果を表わした図である。FIG. 6 is a diagram showing a simulation result.
【図７】伝達関数ないしインパルス応答関数の求め方の
説明図である。FIG. 7 is an explanatory diagram of how to obtain a transfer function or an impulse response function.
【図８】フーリエ変換のための時間窓を表わした図であ
る。FIG. 8 is a diagram showing a time window for Fourier transform.
１０ 計測対象 ２０ 伝達系 ３０ 応答特性計測装置 ３１，３２ Ａ／Ｄ変換器 ３３ 入力信号切出し置換手段 ３４ 出力信号切出し手段 ３５ 応答特性演算手段 10 Measurement Target 20 Transfer System 30 Response Characteristic Measuring Device 31, 32 A / D Converter 33 Input Signal Cutout Replacement Unit 34 Output Signal Cutout Unit 35 Response Characteristic Calculation Unit
───────────────────────────────────────────────────── フロントページの続き (72)発明者 城戸 健一 神奈川県横浜市緑区白山１丁目16番１号 株式会社小野測器テクニカルセンター内 (72)発明者 浦壁 浩司 神奈川県横浜市緑区白山１丁目16番１号 株式会社小野測器テクニカルセンター内 (72)発明者 鈴木 英男 神奈川県横浜市緑区白山１丁目16番１号 株式会社小野測器テクニカルセンター内 ─────────────────────────────────────────────────── ─── Continuation of front page (72) Kenichi Kido 1-16-1 Shirayama, Midori-ku, Yokohama-shi, Kanagawa Ono Sokki Technical Center Co., Ltd. (72) Koji Urabei 1 Shirayama, Midori-ku, Yokohama-shi, Kanagawa 16-1 1-1 Ono Sokki Technical Center Co., Ltd. (72) Inventor Hideo Suzuki 1-1-16 Shirayama Midori-ku, Yokohama-shi Kanagawa Ono Sokki Technical Center Co., Ltd.
Claims (1)
対象から出力された出力信号とに基づいて前記計測対象
の伝達関数もしくはインパルス応答関数を求める応答特
性計測装置において、 所定の時刻ｔ0 から所定の第１の時間幅Ｔ1 の間の前記
入力信号を切り出すとともに前記第１の時間Ｔ1 よりも
短い所定の第２の時間幅をＴ2 としたときに時刻ｔ0 ＋
Ｔ2 から時刻ｔ0 ＋Ｔ1 までの間の前記入力信号を零に
置換することにより置換入力信号を生成する入力信号切
出し置換手段と、 前記所定の時刻ｔ0 から前記第１の時間幅Ｔ1 の間の前
記出力信号を切り出す出力信号切出し手段と、 前記入力信号切出し置換手段により生成された前記置換
入力信号と前記出力信号切出し手段により切出された前
記出力信号とに基づいて前記計測対象の伝達関数もしく
はインパルス応答関数を求める応答特性演算手段とを備
えたことを特徴とする応答特性計測装置。1. A response characteristic measuring apparatus that obtains a transfer function or an impulse response function of the measurement target based on an input signal input to the measurement target and an output signal output from the measurement target, at a predetermined time t 0. When the input signal is cut out for a predetermined first time width T 1 from the time t 2 and a predetermined second time width shorter than the first time T 1 is T 2 , the time t 0 +
An input signal cut replacement means for generating a replacement input signals by substituting zero the input signals between the T 2 to time t 0 + T 1, wherein the predetermined time t 0 a first time width T 1 Between the output signal cutout means for cutting out the output signal between, and the replacement input signal generated by the input signal cutout replacement means and the output signal cut out by the output signal cutout means A response characteristic measuring device, comprising: a response characteristic calculating means for obtaining a transfer function or an impulse response function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP23022793A JP3234365B2 (en) | 1993-09-16 | 1993-09-16 | Response characteristic measuring device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP23022793A JP3234365B2 (en) | 1993-09-16 | 1993-09-16 | Response characteristic measuring device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0783977A true JPH0783977A (en) | 1995-03-31 |
| JP3234365B2 JP3234365B2 (en) | 2001-12-04 |
Family
ID=16904541
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP23022793A Expired - Fee Related JP3234365B2 (en) | 1993-09-16 | 1993-09-16 | Response characteristic measuring device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3234365B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017167591A (en) * | 2016-03-14 | 2017-09-21 | オムロン株式会社 | Simulation apparatus, simulation method, control program, and record medium |
| JP2017167592A (en) * | 2016-03-14 | 2017-09-21 | オムロン株式会社 | Simulation apparatus, simulation method, control program, and record medium |
-
1993
- 1993-09-16 JP JP23022793A patent/JP3234365B2/en not_active Expired - Fee Related
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017167591A (en) * | 2016-03-14 | 2017-09-21 | オムロン株式会社 | Simulation apparatus, simulation method, control program, and record medium |
| JP2017167592A (en) * | 2016-03-14 | 2017-09-21 | オムロン株式会社 | Simulation apparatus, simulation method, control program, and record medium |
| US10909278B2 (en) | 2016-03-14 | 2021-02-02 | Omron Corporation | Simulation device, simulation method, control program and recording medium |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3234365B2 (en) | 2001-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rife et al. | Transfer-function measurement with maximum-length sequences | |
| US5874916A (en) | Frequency selective TDOA/FDOA cross-correlation | |
| US7880672B1 (en) | Generating nonlinear FM chirp radar signals by multiple integrations | |
| JPH06317619A (en) | Group delay estimation system using least-square adaptation for phase-response inclination | |
| US7630432B2 (en) | Method for analysing the channel impulse response of a transmission channel | |
| US5168213A (en) | Swept signal analysis instrument and method | |
| Challis et al. | A wide bandwidth spectrometer for rapid ultrasonic absorption measurements in liquids | |
| JPH0783977A (en) | Response characteristics measuring apparatus | |
| US5917738A (en) | Removing the gibbs phenomenon in fourier transform processing in digital filters or other spectral resolution devices | |
| JPH09510802A (en) | Method and apparatus for determining mask limits | |
| US5732044A (en) | System and method for compensating for doppler shifts in signals by downsampling | |
| JP2867769B2 (en) | Sound measurement method and device | |
| US5610847A (en) | Ratiometric fourier analyzer | |
| JPS6244620B2 (en) | ||
| Van der Auweraer et al. | Accurate modal analysis measurements with programmed sine wave excitation | |
| Nguyen et al. | Digital interpolating peak locator | |
| White et al. | Signal processing techniques | |
| US7561639B2 (en) | Method and device for estimating channel properties of a transmission channel | |
| Wan et al. | Design of periodic signals using FM sweeps and amplitude modulation for ocean acoustic travel-time measurements | |
| JPH07120517A (en) | Response measuring apparatus | |
| JPH07311081A (en) | Method and device for acoustic measurement | |
| Gupta et al. | Time and frequency domain analyses of single-degree-of-freedom systems | |
| Wolf et al. | Amplitude and frequency estimator for aperiodic multi-frequency noisy vibration signals of a tram gearbox | |
| Ramakrishnan et al. | Generation of desired signals from acoustic drivers | |
| JPH07120516A (en) | Response characteristics measuring apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20010911 |
|
| LAPS | Cancellation because of no payment of annual fees |