tf_estimate（MathScript RTモジュール関数）

所有クラス: spectral analysis

要件：MathScript RTモジュール

構文

tf_estimate(x, y)

tf_estimate(x, y, fftsize, fs, window)

tf_estimate(x, y, fftsize, fs, window, noverlap)

tf_estimate(x, y, fftsize, fs, window, noverlap, flag)

H = tf_estimate(x, y)

H = tf_estimate(x, y, fftsize, fs, window)

H = tf_estimate(x, y, fftsize, fs, window, noverlap)

H = tf_estimate(x, y, fftsize, fs, window, noverlap, flag)

[H, f] = tf_estimate(x, y)

[H, f] = tf_estimate(x, y, fftsize, fs, window)

[H, f] = tf_estimate(x, y, fftsize, fs, window, noverlap)

[H, f] = tf_estimate(x, y, fftsize, fs, window, noverlap, flag)

従来名: tfe

説明

刺激または応答から伝達関数を予測します。出力を指定しない場合、この関数は現在のプロットウィンドウで予測された伝達関数をプロットします。

詳細

サンプル

入力

名前	説明
x	刺激信号を指定します。
y	xによって刺激された応答信号を指定します。length(y)はlength(x)と等しくなければなりません。
fftsize	FFTの長さを指定して、データの各セグメントに実行します。デフォルトは256の最小値、xの長さになります。
fs	サンプリング周波数を指定します。デフォルトは2です。
window	窓を指定して、xおよびyに適用します。窓の長さはfftsize以下でなければなりません。デフォルトはfftsizeの長さのハニング窓です。
noverlap	データの隣接するセグメント間のオーバーラップするデータ要素の数を指定します。noverlapはwindowの長さより短くなければなりません。デフォルトは0です。
flag	windowを適用する前にxおよびyを前処理するのに使用するメソッドを指定します。flagは以下の値を受け入れる文字列です。 'linear' LabVIEWはデータの各セグメントの最良線形フィットを削除します。 'mean' LabVIEWはデータの各セグメントの平均を削除します。 'none'（デフォルト） LabVIEWはxおよびyを前処理しません。

出力

名前	説明
H	伝達関数を返すxとyが実数の場合、Hの長さは偶数fftsizeに対してfftsize / 2 + 1、奇数fftsizeに対して(fftsize + 1) / 2になります。XとYが複素数の場合、Hの長さはfftsizeになります。
f	Hに対応する周波数を返します。

詳細

この関数はウェルチメソッドを使用して、xのパワースペクトルとxとyのクロスパワーを求めます。LabVIEWは以下の公式を使用して、H(f) = P_xy(f)/P_xx(f)の伝達関数を計算します。

次の表には、この関数のサポート特性 が記載されています。

LabVIEWランタイムエンジンでサポートされる	はい（出力を要求した場合）
RTターゲットでサポートされる	はい（出力を要求した場合）
RTで制限付きの実行時間に適する	特性なし

サンプル

fs = 1000;
T = 2;
t = 0:1 / fs:T-1 / fs;
fftsize = 512;
window = win_hann2(fftsize);
noverlap = fftsize / 2;
dflag = 'none';
xn = randnormal(1, length(t));
h = fir_win(5, 0.2, 'low', win_hann2(6));
yn = filter(h, 1, xn);
[H, f] = tf_estimate(xn, yn, fftsize, fs, window, noverlap, dflag);

関連項目

coherence
crosssd
filter
fir_win
randnormal
win_hann2