Main Content

시계열 모델 분석하기

이 예제에서는 시계열 모델을 분석하는 방법을 보여줍니다.

시계열 모델에는 입력이 없습니다. 하지만 해당 모델에서 여러 가지 응답 계산 명령을 사용할 수 있습니다. 잡음원 e(t)는 (암시적으로) 측정된 입력으로 취급됩니다. 따라서 step(sys)는 계단 입력이 잡음 채널 e(t)에 적용되었다고 가정하여 계단 응답을 플로팅합니다.

시계열 모델을 처리하는 방식에서의 모호성을 방지하기 위해 noise2meas를 사용하여 명시적으로 이를 입력-출력 모델로 변환할 수 있습니다. 이 명령은 잡음 입력 e(t)를 측정된 입력으로 취급하고, Ny 출력을 갖는 선형 시계열 모델을 Ny 출력과 Ny 입력을 갖는 입력-출력 모델로 변환합니다. 결과로 생성된 모델을 bode, nyquist, iopzmap 등의 명령과 함께 사용하여 H 전달 함수의 특성을 조사할 수 있습니다.

시계열 모델을 추정합니다.

load iddata9
sys = ar(z9,4);

시계열 모델을 입력-출력 모델로 변환합니다.

iosys = noise2meas(sys);

H의 계단 응답을 플로팅합니다.

step(iosys);

Figure contains an axes object. The axes object with title From: e@y1 To: y1 contains an object of type line. This object represents iosys.

H의 극점과 영점을 플로팅합니다.

iopzmap(iosys);

Figure contains an axes object. The axes object with title From: e@y1 To: y1 contains 2 objects of type line. This object represents iosys.

입력-출력 모델로 변환하지 않고 시계열 스펙트럼을 직접 계산하여 플로팅합니다.

spectrum(sys);

Figure contains an axes object. The axes object with title From: e@y1 To: y1 contains an object of type line. This object represents sys.

이 명령은 시계열 스펙트럼 진폭 Φ(ω)=H(ω)2를 플로팅합니다.

참고 항목

|

관련 항목