Сызықты динамикалық жүйелерді идентификаттау. Үйірткі теңдеуінен салмақтық функцияны анықтау Уақыттың шектелген аралығында сызықты стационарлы жүйенің кіріс және шығыс сигналдарын бақылау нәтижелері бойынша оның салмақтық функциясын (импульстік өтпелі функциясын) анықтау керек.
Сурет 13.2 – Сызықты стационарлы жүйе
Кірісі мен нөльдік бастапқы шарттарда жүйенің шығыс сигналы үйірткі интегралымен өрнектеледі:
(13.23)
, болғанда деп жорамалданады.
Енді уақыттың w(t)кіріс функциясын қадаммен N нүктелерде құрама-тұрақты функциямен аппроксимациялауды енгізейік, сонымен қатар, :
при (13.24)
Бөлу нүктелердің арасындағы h(t) –ны тұрақты етіп қабылдаймыз:
nt=n болғанда (13.23) интеграл және сатылы аппроксимациялау терминдерінде жуықтап келесі түрде жазылады:
(13.25)
Шығысты (өлшемі N) бақылау векторын YT(T) = [y() y(2) ... y(N)] және салмақтық функция мәндері: арқылы белгілейік.
(13.25) теңдеуді векторлы-матрицалық түрде қайта жазайық:
Y(T) = Wh(T) (13.26)
W матрица келесі теңдікпен анықталады:
W – сол жақты үшбұрыш матрица, диагоналінде w(0).
Енді есеп (13.26) теңдеуден фиксациялау нүктелерінде салмақтық функция мәндерінің һ векторын анықтауға келтірілді. болғандықтан w – өзгешеленбеген және .
Сондықтан (13.26) теңдеу шешімін формальды келесі түрде жазуға болады:
(13.27)
W матрицаның сол жақты үшбұрышты формада болуының арқасында һ үшін өрнекті рекуррентті түрде қайта жазуға болады:
(13.28)
мұнда:
Қарастырған тәсілдің құндылығы – кез-келген кіріс сигналдарды пайдалану мүмкіндігі. Арнайы тестілік сигналдарды пайдалану қажеттілігі болмағандықтан жүйені қалыпты пайдалану үрдісінде алынған жүзеге асыруларын пайдалануға болады.
Егер кіріс сигнал бірлік секірудің функциясы болатын болса, (13.28) алгоритмі едәуір жеңілдейді. Бұл жағдайда барлық і үшін w(i) = 1, (13.28) келесі түрге ие болады:
.
шаманы анықтап аламыз, бұл кезде Hn = Hn-1 +hn-1.