matlab에 : 전달 함수

Question

내가 S2 = S1 ' 은 사실 내가 S2 * S1

s = tf('s') 
Suu = -1.6/((s-4)*(s+4)) 
Sux = -0.8/((s+4)*(s-4)*(s^2 + 0.1*s + 1)) 
Sxx = 0.3*(s - 4.163)*(s + 4.163)/((s+4)*(s-4)*(s^2 - 0.1*s + 1)*(s^2 + 0.1*s + 1)) 
Sxu = Sux' 
SxdSdx = Sxx - (Sxu*Sux)/Suu 

SXD을 확인하는 방법 그리고 제품을 같은 방법으로 2 개 전송 기능으로 내 전송 기능을 분할해야

를 인수 분해하는 방법 Sdx = Sdx '이면? 아무도 도와 줄 수 없나요?

Answer 1

당신이하려고하는 것은 스펙트럼 인수 분해입니다. 스펙트럼 인수 분해 알고리즘에 대한 좋은 설문지는 H. Sayed와 T. Kailath, "Spectral Factorization Methods의 조사", Numer. 선형 대수학 Appl., vol. 8, PP. 467-496,

당신의 전달 함수는 합리적인 함수로 표현 될 경우 2001 년 (즉, S에서 다항식의 비율) 당신이 S1을 선택할 수 있습니다 다음 인수 분해 할 수 있습니다 누구의 제로 합리적인 기능 할 수 (분자 다항식의 근원)과 극 (분모 다항식의 근원)은 모두 왼쪽 평면에 있습니다 (즉, 실수 부분이 0보다 작음). 실제 루트가있는 경우 쌍으로 표시되어야하며 각 쌍의 멤버 중 하나를 S1에 할당 할 수 있습니다. 이렇게하면 S1 * conj (S1)이 원래 전송 함수가됩니다.

전달 함수를 요소로하는 수치 적 방법은 전달 함수의 극점과 영점 (몇 개나 그 근처에 있는지 등), 전달 함수를 아는 방법 그것의 다항식 계수를 통해, 당신은 공간의 선택된 점에서의 값을 통해 그것을 안다?), 얼마나 정확하게 그것을 알고 있는가 (전달 함수 계수를 얼마나 정확하게 나타낼 것인가), 그리고 여러분의 인수 분해가 얼마나 정확한지를 알 수 있습니다.