tfestimate

전달 함수 추정값

구문

txy = tfestimate(x,y)

txy = tfestimate(x,y,window)

txy = tfestimate(x,y,window,nOverlap)

txy = tfestimate(x,y,window,nOverlap,nfft)

txy = tfestimate(___,"mimo")

[txy,w] = tfestimate(___)

[txy,f] = tfestimate(___,Fs)

[txy,w] = tfestimate(x,y,window,nOverlap,w)

[txy,f] = tfestimate(x,y,window,nOverlap,f,Fs)

[___] = tfestimate(x,y,___,freqRange)

[___] = tfestimate(___,Name=Value)

tfestimate(___)

설명

txy = tfestimate(x,y)는 일련의 주파수에서 평가된 입력 신호 x와 출력 신호 y 간의 전달 함수 추정값을 구합니다.

• x와 y가 모두 벡터인 경우 길이가 같아야 합니다.

• 두 신호 중 하나는 행렬이고 다른 하나는 벡터인 경우, 벡터의 길이는 행렬에 포함된 행 개수와 동일해야 합니다. 이 함수는 벡터를 확장하고 열별 전달 함수 추정값으로 구성된 행렬을 반환합니다.

• x와 y가 행의 개수가 같지만 열의 개수가 다른 행렬인 경우 txy는 모든 입력 신호와 출력 신호를 결합하는 MIMO(다중 입력/다중 출력) 전달 함수입니다. txy는 3차원 배열입니다. x는 m개의 열을 가지고 y는 n개의 열을 가지는 경우, txy는 n개의 열과 m개의 페이지를 가집니다. 자세한 내용은 전달 함수 항목을 참조하십시오.

• x와 y가 같은 크기의 행렬인 경우, tfestimate는 각 열에 대해 연산을 수행합니다. 즉, txy(:,n) = tfestimate(x(:,n),y(:,n))입니다. MIMO 추정값을 구하려면 인수 목록에 "mimo"를 추가하십시오.

txy = tfestimate(x,y,window)는 window를 사용하여 x와 y를 세그먼트로 나누고 윈도우를 적용합니다.

txy = tfestimate(x,y,window,nOverlap)은 인접 세그먼트 간에 nOverlap개의 샘플이 중첩되도록 적용합니다.

txy = tfestimate(x,y,window,nOverlap,nfft)는 nfft개 샘플링 점을 사용하여 이산 푸리에 변환을 계산합니다.

txy = tfestimate(___,"mimo")는 행렬 입력값에 대해 MIMO 전달 함수를 계산합니다. 이 구문은 위에 열거된 구문에 사용할 수 있습니다.

[txy,w] = tfestimate(___)는 전달 함수가 추정된 정규화 주파수로 구성된 벡터 w를 반환합니다.

[txy,f] = tfestimate(___,Fs)는 전달 함수가 추정된 주파수(샘플 레이트 f의 값으로 표현됨)로 구성된 벡터 Fs를 반환합니다. Fs는 tfestimate에 대한 6번째 숫자형 입력값이어야 합니다. 샘플 레이트를 입력하고 위에 열거된 옵션 인수의 디폴트 값을 그대로 사용하려면 이 인수를 빈 값 []로 지정하십시오.

[txy,w] = tfestimate(x,y,window,nOverlap,w)는 w에 지정된 정규화 주파수에서 전달 함수 추정값을 반환합니다.

[txy,f] = tfestimate(x,y,window,nOverlap,f,Fs)는 f에 지정된 주파수에서 전달 함수 추정값을 반환합니다.

[___] = tfestimate(x,y,___,freqRange)는 freqRange로 지정된 주파수 범위에 대한 전달 함수 추정값을 반환합니다. freqRange에 유효한 옵션은 "onesided", "twosided", "centered"입니다.

[___] = tfestimate(___,Name=Value)는 이름-값 인수를 사용하여 추가 옵션을 지정합니다. 옵션으로는 전달 함수 추정량과 전달 함수 추정값을 플로팅할 대상 부모 컨테이너가 포함됩니다. 이러한 인수를 위에 열거된 입력 구문에 추가할 수 있습니다.

tfestimate(___)에 출력 인수를 지정하지 않으면 현재 Figure 창에 전달 함수 추정값을 플로팅합니다.

예제

두 시퀀스 간 전달 함수

두 시퀀스 x와 y 간의 전달 함수 추정값을 계산하고 플로팅합니다. 시퀀스 x는 백색 가우스 잡음으로 구성되어 있습니다. y는 $$0.2\pi$$ rad/sample의 정규화된 차단 주파수를 갖는 30차 저역통과 필터를 사용하여 x를 필터링한 결과입니다. 사각 윈도우를 사용하여 필터를 설계합니다. 전달 함수 추정값에 대해 샘플 레이트가 500Hz이고 길이가 1024인 해밍 윈도우를 지정합니다.

h = fir1(30,0.2,rectwin(31));
x = randn(16384,1);
y = filter(h,1,x);

fs = 500;
tfestimate(x,y,1024,[],[],fs)

전달 함수가 필터의 주파수 응답에 근사하는지 확인합니다.

freqz(h,1,[],fs)

변수로 전달 함수 추정값을 반환하고 추정값의 절댓값을 데시벨 단위로 플로팅하여 동일한 결과를 구합니다.

[Txy,f] = tfestimate(x,y,1024,[],[],fs);

plot(f,mag2db(abs(Txy)))

SISO 전달 함수

간단한 SISO(단일 입력 단일 출력) 시스템에 대한 전달 함수를 추정하고 이를 정의와 비교합니다.

1차원의 이산시간 진동 시스템은 단위 탄성 상수를 갖는 스프링으로 벽에 연결된 단위 질량 $m$(단위: kg)으로 구성됩니다. 센서가 질량의 가속도 $a$를 $F_s=1$Hz로 샘플링합니다. 댐퍼는 감쇠 상수 $b=0.01$kg/s를 사용하여 속도에 비례하는 힘을 가해 질량의 운동을 저해합니다.

2000개의 시간 샘플을 생성합니다. 샘플링 간격을 $\Delta t=1/F_s$로 정의합니다.

Fs = 1;
dt = 1/Fs;
N = 2000;
t = dt*(0:N-1);
b = 0.01;

시스템은 다음 상태공간 모델로 설명할 수 있습니다.

$$\begin{array}{c}x(k+1)=Ax(k)+Bu(k),\\ y(k)=Cx(k)+Du(k),\end{array}$$

여기서 $x={\left[\begin{array}{cc}r& v\end{array}\right]}^T$는 상태 벡터로, $r$과 $v$는 질량의 위치와 속도를 각각 나타내고, $u$는 구동력이며, $y=a$는 측정 출력입니다. 상태공간 행렬은 다음과 같습니다.

$$A=\exp (A_c\Delta t),B=A_c^{-1}(A-I)B_c,C=\left[\begin{array}{cc}-1& -b\end{array}\right],D=1,$$

$I$는 $2\times 2$ 단위 행렬이고, 연속시간 상태공간 행렬은 다음과 같습니다.

$$A_c=\left[\begin{array}{cc}0& 1\\ -1& -b\end{array}\right],B_c=\left[\begin{array}{c}0\\ 1\end{array}\right].$$

Ac = [0 1;-1 -b];
A = expm(Ac*dt);

Bc = [0;1];
B = Ac\(A-eye(size(A)))*Bc;

C = [-1 -b];
D = 1;

질량은 측정 구간의 절반에 대한 임의의 입력값으로 구동됩니다. 상태공간 모델을 사용하여 모든 초기 상태값이 0인 시공간 시스템의 변화를 계산합니다. 질량의 가속도를 시간의 함수로 플로팅합니다.

rng("default")

u = zeros(1,N);
u(1:N/2) = randn(1,N/2);

y = 0;
x = [0;0];
for k = 1:N
    y(k) = C*x + D*u(k);
    x = A*x + B*u(k);
end

plot(t,y)

시스템의 전달 함수를 주파수의 함수로 추정합니다. 2048개의 DFT 점을 사용하고 형태 인자가 15인 카이저 윈도우를 지정합니다. 인접 세그먼트 간에 디폴트 중첩값을 적용합니다.

nfs = 2048;
wind = kaiser(N,15);

[txy,ft] = tfestimate(u,y,wind,[],nfs,Fs);

이산시간 시스템의 주파수-응답 함수는 단위원에서 계산된, 시스템의 시간 영역 전달 함수에 대한 Z 변환으로 표현할 수 있습니다. tfestimate로 계산된 추정값이 이 정의와 일치하는지 확인합니다.

[b,a] = ss2tf(A,B,C,D);

fz = 0:1/nfs:1/2-1/nfs;
z = exp(2j*pi*fz);
frf = polyval(b,z)./polyval(a,z);

plot(ft,mag2db(abs(txy)))
hold on
plot(fz,mag2db(abs(frf)))
hold off
grid
ylim([-60 40])

tfestimate의 내장 기능을 사용하여 추정값을 플로팅합니다.

tfestimate(u,y,wind,[],nfs,Fs)

MIMO 시스템의 전달 함수 추정

MIMO(다중 입력/다중 출력) 시스템에 대해 전달 함수를 추정합니다.

두 질량 양쪽에 스프링과 댐퍼를 연결하면 이상적인 1차원 이산시간 진동 시스템이 구성됩니다. 시스템 입력 배열 u는 질량에 가해지는 무작위 구동력으로 구성됩니다. 시스템 출력 배열 y는 초기 기준 위치로부터 관측된 질량의 변위를 포함합니다. 시스템은 40Hz의 Fs 레이트로 샘플링됩니다.

MIMO 시스템 입력, 시스템 출력 및 샘플 레이트가 포함된 데이터 파일을 불러옵니다. 예제 Frequency-Response Analysis of MIMO System에서는 이 예제에 사용된 데이터를 생성한 시스템을 분석합니다.

load MIMOdata Fs u y

시스템 데이터와 함수 tfestimate를 사용하여 시스템의 주파수 영역 전달 함수를 추정하고 플로팅합니다. 인접 세그먼트 간의 중첩이 50%가 되도록 신호를 5000개 샘플 세그먼트로 나눕니다. 각 세그먼트에 해닝 윈도우를 적용하고 8192개 주파수 점을 사용하여 신호 세그먼트의 이산 푸리에 변환을 계산합니다. "mimo" 옵션을 선택하여 네 개의 모든 전달 함수를 생성합니다.

g = hann(5000);
nfft = 8192;
ol = 2500;

[tXY,ft] = tfestimate(u,y,g,ol,nfft,Fs,"mimo");

tiledlayout flow
for jk = 1:2
    for kj = 1:2
        nexttile
        plot(ft,mag2db(abs(tXY(:,jk,kj))))
        grid on
        ylim([-120 0])
        title("Input "+jk+", Output "+kj)
        xlabel("Frequency (Hz)")
        ylabel("Magnitude (dB)")
    end
end

UI Figure 창을 만들고 좌표축과 패널 컨테이너를 추가합니다.

uif = uifigure(Position=[100 100 720 540]);
ax = uiaxes(uif,Position=[5 280 450 240]);
p = uipanel(uif,Position=[180 5 520 270], ...
    Title="Transfer-Function Estimate in Panel Container", ...
    BackgroundColor="white");

시스템에 대한 네 개 전달 함수의 추정값을 좌표축에 플로팅하고 시스템의 각 입력-출력 쌍을 패널 컨테이너에 플로팅합니다. (R2026a 이후)

tfestimate(u,y,g,ol,nfft,Fs,"mimo",Parent=ax)
legend(ax,"Input "+[1 2]+", Output "+[1 2]',Location="best")

tfestimate(u,y,g,ol,nfft,Fs,Parent=p)

입력 인수

x — 입력 신호

벡터 | 행렬

입력 신호로, 벡터나 행렬로 지정됩니다.

예: cos(pi/4*(0:159))+randn(1,160)은 백색 가우스 잡음을 포함하는 정현파를 지정합니다.

데이터형: single | double
복소수 지원 여부: 예

y — 출력 신호

벡터 | 행렬

출력 신호로, 벡터나 행렬로 지정됩니다.

데이터형: single | double
복소수 지원 여부: 예

window — 윈도우

정수 | 벡터 | []

윈도우로, 정수로 지정되거나 행 벡터 또는 열 벡터로 지정됩니다. window를 사용하여 신호를 여러 세그먼트로 나눕니다.

• window가 정수이면 tfestimate는 x와 y를 길이가 window인 세그먼트로 나누고 각 세그먼트에 해당 길이의 해밍 윈도우를 적용합니다.

• window가 벡터이면 tfestimate는 x와 y를 벡터와 길이가 같은 세그먼트로 나누고 window를 사용하여 각 세그먼트에 윈도우를 적용합니다.

x와 y의 길이를 nOverlap개의 샘플이 중첩된 정수 개수의 세그먼트로 정확히 나눌 수 없는 경우, 신호가 이에 맞게 적절하게 잘립니다.

window를 빈 값으로 지정하면, tfestimate는 해밍 윈도우를 사용하여 x와 y를 nOverlap개의 중첩 샘플을 가지는 8개의 세그먼트로 나눕니다.

사용 가능한 윈도우 목록은 윈도우 항목을 참조하십시오.

예: hann(N+1) 및 (1-cos(2*pi*(0:N)'/N))/2는 모두 길이가 N + 1인 핸 윈도우를 지정합니다.

데이터형: single | double

nOverlap — 중첩 샘플의 개수

양의 정수 | []

중첩(Overlap) 샘플의 개수로, 양의 정수로 지정됩니다.

• window가 스칼라이면 nOverlap은 window보다 작아야 합니다.

• window가 벡터이면 nOverlap은 window의 길이보다 작아야 합니다.

nOverlap을 빈 값으로 지정하면 tfestimate는 세그먼트 간에 50% 중첩을 만드는 값을 사용합니다.

데이터형: double | single

nfft — DFT를 적용할 지점의 개수

양의 정수 | []

DFT를 적용할 지점의 개수로, 양의 정수로 지정됩니다. nfft를 빈 값으로 지정하면 tfestimate는 이 인수를 max(256,2^p)으로 설정합니다. 여기서, 길이가 N인 입력 신호의 경우 p = ⌈log₂ N⌉이고 ⌈ ⌉ 기호는 천장 함수(Ceiling function)를 나타냅니다.

데이터형: single | double

Fs — 샘플 레이트

양의 스칼라

샘플 레이트로, 양의 스칼라로 지정됩니다. 샘플 레이트는 단위 시간당 샘플 개수입니다. 시간 단위가 초이면 샘플 레이트의 단위는 Hz입니다.

• Fs를 빈 값([])으로 지정하면, tfestimate 함수는 입력 신호 x의 샘플 레이트를 1Hz로 가정합니다.

• Fs를 지정하지 않으면, tfestimate 함수는 입력 신호 x의 샘플 레이트를 2π rad/sample로 가정합니다.

w — 정규화 주파수

벡터

정규화 주파수로, 적어도 2개의 요소를 가진 행 벡터나 열 벡터로 지정됩니다. 정규화 주파수는 rad/sample을 단위로 사용합니다.

예: w = [pi/4 pi/2]

데이터형: double | single

f — 주파수

벡터

주파수로, 적어도 2개의 요소를 가진 행 벡터나 열 벡터로 지정됩니다. 주파수는 단위 시간당 주기를 단위로 사용합니다. 단위 시간은 샘플 레이트 fs로 지정됩니다. fs의 단위가 샘플/초이면 f는 Hz 단위입니다.

예: fs = 1000; f = [100 200]

데이터형: double | single

freqRange — 전달 함수 추정값에 사용할 주파수 범위

"onesided" | "twosided" | "centered"

전달 함수 추정값의 주파수 범위로, "onesided", "twosided" 또는 "centered"로 지정됩니다. 디폴트 값은 실수 값을 갖는 신호의 경우 "onesided"이고, 복소수 값을 갖는 신호의 경우 "twosided"입니다.

• "onesided" — 실수 값을 갖는 두 개의 입력 신호 x와 y 간 전달 함수의 단측 추정값을 반환합니다. nfft가 짝수이면 txy가 nfft/2 + 1개의 행을 가지고 구간 [0,π] rad/sample에 대해 계산됩니다. nfft가 홀수이면 txy가 (nfft + 1)/2개의 행을 가지고 구간은 [0,π) rad/sample이 됩니다. Fs를 지정하는 경우, 이에 대응하는 구간은 nfft가 짝수이면 [0,Fs/2] cycles/unit time이고, nfft가 홀수이면 [0,Fs/2) cycles/unit time입니다.

• "twosided" — 실수 값 또는 복소수 값을 갖는 두 개의 입력 신호 x와 y 간 전달 함수의 양측 추정값을 반환합니다. 이 경우, txy는 nfft개의 행을 가지고 구간 [0,2π) rad/sample에 대해 계산됩니다. Fs를 지정하는 경우 구간은 [0,Fs) cycles/unit time이 됩니다.

• "centered" — 실수 값 또는 복소수 값을 갖는 두 개의 입력 신호 x와 y 간 전달 함수의, 중심이 맞춰진 양측 추정값을 반환합니다. 이 경우, txy는 nfft개의 행을 가지고, nfft가 짝수이면 구간 (–π,π] rad/sample에 대해 계산되며, nfft가 홀수이면 구간 (–π,π) rad/sample에 대해 계산됩니다. Fs를 지정하는 경우, 이에 대응하는 구간은 nfft가 짝수이면 (–Fs/2, Fs/2] cycles/unit time이고, nfft가 홀수이면 (–Fs/2, Fs/2) cycles/unit time입니다.

주파수로 구성된 벡터 w 또는 f와 freqRange를 입력값으로 지정하면 tfestimate 함수는 freqRange를 무시합니다.

이름-값 인수

선택적 인수 쌍을 Name1=Value1,...,NameN=ValueN으로 지정합니다. 여기서 Name은 인수 이름이고 Value는 대응값입니다. 이름-값 인수는 다른 인수 뒤에 와야 하지만, 인수 쌍의 순서는 상관없습니다.

예: txy = tfestimate(x,y,Estimator="H1",Parent=h)는 "H1" 추정량을 사용하여 입력 x와 출력 신호 y를 가진 시스템의 전달 함수 추정값을 계산하고 대상 부모 컨테이너 h에 플로팅합니다.

Estimator — 전달 함수 추정량

"H1" (디폴트 값) | "H2"

전달 함수 추정량으로, "H1" 또는 "H2"로 지정됩니다.

• 잡음이 입력 신호와 상관 관계가 없는 경우 "H1"을 사용합니다.

• 잡음이 출력 신호와 상관 관계가 없는 경우 "H2"를 사용합니다. 이 경우, 입력 신호의 개수는 출력 신호의 개수와 같아야 합니다.

자세한 내용은 전달 함수 항목을 참조하십시오.

Parent — 대상 부모 컨테이너

[] (디폴트 값) | Axes 객체 | UIAxes 객체 | Panel 객체

R2026a 이후

대상 부모 컨테이너로, Axes 객체, UIAxes 객체 또는 Panel 객체로 지정됩니다.

Parent를 지정하면 tfestimate 함수는 출력 인수를 사용하여 함수를 호출하는지 여부와 관계없이 전달 함수 추정값을 지정된 대상 부모 컨테이너에 플로팅합니다.

MATLAB^® 그래픽스의 대상 컨테이너와 부모-자식 관계에 대한 자세한 내용은 그래픽스 객체 계층 구조 항목을 참조하십시오. 앱 설계 시 UIAxes 객체와 Panel 객체에서 Parent를 사용하는 방법에 대한 자세한 내용은 Plot Spectral Representations of Signal in App Designer 항목을 참조하십시오.

출력 인수

txy — 전달 함수 추정값

벡터 | 행렬 | 3차원 배열

전달 함수 추정값으로, 벡터, 행렬 또는 3차원 배열로 반환됩니다.

w — 정규화 주파수

벡터

정규화 주파수로, 실수 값 열 벡터로 반환됩니다.

f — 순환 주파수

벡터

순환 주파수로, 실수 값 열 벡터로 반환됩니다.

세부 정보

전달 함수

입력값 x와 출력값 y 간의 관계는 선형 시불변 전달 함수 txy로 모델링됩니다. 주파수 영역에서 Y(f) = H(f)X(f)입니다.

• SISO(단일 입력 단일 출력) 시스템의 경우, 전달 함수의 H₁ 추정값은 다음과 같이 지정됩니다.

  $$H_1(f)=\frac{P_{yx}\left(f\right)}{P_{xx}\left(f\right)},$$

  여기서 P_yx는 x와 y의 상호 파워 스펙트럼 밀도이고, P_xx는 x의 파워 스펙트럼 밀도입니다. 이 추정값은 잡음이 시스템 입력값과 상관관계가 없다고 가정합니다.

  MIMO(다중 입력/다중 출력) 시스템의 경우, H₁ 추정량은 m개 입력값과 n개 출력값에 대해 다음과 같이 됩니다.

  $$H_1(f)=P_{YX}(f)P_{XX}^{-1}(f)=\left[\begin{array}{cccc}{P}_{y_1{x}_1}(f)& P_{y_1{x}_2}(f)& \cdots & P_{y_1{x}_m}(f)\\ P_{y_2{x}_1}(f)& P_{y_2{x}_2}(f)& \cdots & P_{y_2{x}_m}(f)\\ ⋮& ⋮& \ddots & ⋮\\ P_{y_n{x}_1}(f)& P_{y_n{x}_2}(f)& \cdots & P_{y_n{x}_m}(f)\end{array}\right]{\left[\begin{array}{cccc}{P}_{x_1{x}_1}(f)& P_{x_1{x}_2}(f)& \cdots & P_{x_1{x}_m}(f)\\ P_{x_2{x}_1}(f)& P_{x_2{x}_2}(f)& \cdots & P_{x_2{x}_m}(f)\\ ⋮& ⋮& \ddots & ⋮\\ P_{x_m{x}_1}(f)& P_{x_m{x}_2}(f)& \cdots & P_{x_m{x}_m}(f)\end{array}\right]}^{-1}$$

  여기서,

  • P_yixk는 k번째 입력값과 i번째 출력값의 상호 파워 스펙트럼 밀도입니다.

  • P_xixk는 k번째 입력값과 i번째 입력값의 상호 파워 스펙트럼 밀도입니다.

  입력값과 출력값이 각각 두 개인 경우, 추정량은 다음과 같은 행렬입니다.

  $$H_1(f)=\frac{\left[\begin{array}{cc}{P}_{y_1{x}_1}(f)P_{x_2{x}_2}(f)-P_{y_1{x}_2}(f)P_{x_2{x}_1}(f)& P_{y_1{x}_2}(f)P_{x_1{x}_1}(f)-P_{y_1{x}_1}(f)P_{x_1{x}_2}(f)\\ P_{y_2{x}_1}(f)P_{x_2{x}_2}(f)-P_{y_2{x}_2}(f)P_{x_2{x}_1}(f)& P_{y_2{x}_2}(f)P_{x_1{x}_1}(f)-P_{y_2{x}_1}(f)P_{x_1{x}_2}(f)\end{array}\right]}{P_{x_1{x}_1}(f)P_{x_2{x}_2}(f)-P_{x_1{x}_2}(f)P_{x_2{x}_1}(f)}.$$

• SISO(단일 입력 단일 출력) 시스템의 경우, 전달 함수의 H₂ 추정값은 다음과 같이 지정됩니다.

  $$H_2\left(f\right)=\frac{P_{yy}\left(f\right)}{P_{xy}\left(f\right)},$$

  여기서 P_yy는 y의 파워 스펙트럼 밀도이고, P_xy = P^*_yx는 x와 y에 대한 상호 파워 스펙트럼 밀도의 켤레 복소수입니다. 이 추정값은 잡음이 시스템 출력값과 상관관계가 없다고 가정합니다.

  MIMO 시스템의 경우, H₂ 추정량은 입력값과 출력값의 개수가 같은 경우, 즉 n = m인 경우에만 제대로 정의됩니다. 이 추정량은 다음과 같이 됩니다.

  $$H_2(f)=P_{YY}(f)P_{XY}^{-1}(f)=\left[\begin{array}{cccc}{P}_{y_1{y}_1}(f)& P_{y_1{y}_2}(f)& \cdots & P_{y_1{y}_n}(f)\\ P_{y_2{y}_1}(f)& P_{y_2{y}_2}(f)& \cdots & P_{y_2{y}_n}(f)\\ ⋮& ⋮& \ddots & ⋮\\ P_{y_n{y}_1}(f)& P_{y_n{y}_2}(f)& \cdots & P_{y_n{y}_n}(f)\end{array}\right]{\left[\begin{array}{cccc}{P}_{x_1{y}_1}(f)& P_{x_1{y}_2}(f)& \cdots & P_{x_1{y}_n}(f)\\ P_{x_2{y}_1}(f)& P_{x_2{y}_2}(f)& \cdots & P_{x_2{y}_n}(f)\\ ⋮& ⋮& \ddots & ⋮\\ P_{x_n{y}_1}(f)& P_{x_n{y}_2}(f)& \cdots & P_{x_n{y}_n}(f)\end{array}\right]}^{-1},$$

  여기서,

  • P_yiyk는 k번째 출력값과 i번째 출력값의 상호 파워 스펙트럼 밀도입니다.

  • P_xiyk는 i번째 입력값과 k번째 출력값에 대한 상호 파워 스펙트럼 밀도의 켤레 복소수입니다.

알고리즘

tfestimate는 Welch의 평균화된 주기도 방법을 사용합니다. 자세한 내용은 pwelch 항목을 참조하십시오.

참고 문헌

[1] Vold, Håvard, John Crowley, and G. Thomas Rocklin. “New Ways of Estimating Frequency Response Functions.” Sound and Vibration. Vol. 18, November 1984, pp. 34–38.

확장 기능

C/C++ 코드 생성

MATLAB® Coder™를 사용하여 C 코드나 C++ 코드를 생성할 수 있습니다.

사용법 관련 참고 및 제한 사항:

• 이름 값 쌍을 사용하여 지정되는 인수는 컴파일타임 상수여야 합니다.

GPU 코드 생성

GPU Coder™를 사용하여 NVIDIA® GPU용 CUDA® 코드를 생성할 수 있습니다.

C/C++ 코드 생성 섹션의 사용법 관련 참고 및 제한 사항을 참조하십시오. GPU 코드 생성에도 동일한 사용법 관련 참고 및 제한 사항이 적용됩니다.

GPU 배열

Parallel Computing Toolbox™를 사용해 GPU(그래픽스 처리 장치)에서 실행하여 코드 실행 속도를 높일 수 있습니다.

이 함수는 GPU 배열을 완전히 지원합니다. 자세한 내용은 GPU에서 MATLAB 함수 실행하기 (Parallel Computing Toolbox) 항목을 참조하십시오.

버전 내역

R2006a 이전에 개발됨

R2026a: 지정된 부모 대상에서 함수 출력값 시각화

tfestimate 함수는 Parent 입력 인수를 사용하여 지정한 부모 좌표축 또는 패널 컨테이너에 출력값을 플로팅할 수 있습니다.

R2024a: 단정밀도 가변 크기 창 입력값에 대한 코드 생성 지원

tfestimate 함수는 코드 생성 시 단정밀도 가변 크기 창 입력값을 지원합니다.

R2023b: 단정밀도 데이터와 gpuArray 객체를 사용함

tfestimate 함수는 단정밀도 입력값과 gpuArray 객체를 지원합니다. gpuArray 객체를 사용하려면 Parallel Computing Toolbox™가 있어야 합니다.

참고 항목

함수

• cpsd | mscohere | periodogram | pwelch

도움말 항목

• Plot Spectral Representations of Signal in App Designer