CC BY
185168 www.finestreet.ru САПР
Андрей СТРОГОНОВ, к. т. н.
Использование цифровых фильтров
для моделирования деградации выходных параметров ТТЛ ИС в системе MatLab/Simulink
Моделирование процесса деградации параметров ТТЛ ИС с использованием цифровых фильтров
Применение цифровых фильтров для изучения процесса деградации параметров ТТЛ ИС есть результат подхода к исследуемым объектам как к «черному ящику». Метод «черного ящика» — кибернетический. Объект исследования представляется как некоторая кибернетическая система, и она может быть описана своим функциональным оператором. Такой подход преследует цель посредством построения некоторой модели установить изоморфизм не с внутренней структурой и ее функционированием, а с внешними проявлениями ее информативных параметров. Недостатком данного подхода к моделированию процесса деградации параметров ТТЛ ИС является то, что построенная модель не основана на физических представлениях.
Из теории цифровых фильтров известна их связь с линейными моделями временных рядов. Если модели временных рядов успешно используются для построения прогнозов деградации параметров ИС уже много лет, то представляет определенный интерес построить прогнозы с применением цифровых фильтров и посмотреть, насколько они полезны в данном случае.
Воспользуемся авторегрессионными моделями [1-6] и системой MatLab с пакетом моделирования динамических систем Simulink — библиотекой моделирования цифровых фильтров DSP Blockset [7] для моделирования процесса деградации параметров ТТЛ ИС. На практике AR-модель применяется в основном для построения цифровых фильтров, связанных с обработкой сигнала, и для идентификации динамических систем. Так как цифровые фильтры не строят прогнозы за пределы ряда деградации, в случае прекращения подачи сигнала на его вход, будем называть прогнозирование в пределах ряда — моделированием, зачастую отождествляя два термина «моделирование» и «прогнозирование».
В качестве базового описания линейной системы с аддитивными случайными возмущениями используется соотношение: y(t) = G(q)u(t) + H(q)e(t), где G(q), H(q) —передаточные функции системы, q — оператор сдвига, ut —входной и yt — выходной сигналы (временные ряды) в моменты времени t = 1,2,...N, e(t) — белый шум, воздействующий на систему. Произведение G(q)u(t) записывается в виде:
G(q)u(t)=fäg(k)u(t-k),
к=1
где к — задержка. На практике часто работают со свободной от шумов или детерминированной моделью: y(t) = G(q)u(t) [1].
ARX-модель (сочетание AR относится к авторегрессионной (AvtoRegressive) части A(q)y(t), X показывает наличие внешнего (eXternal) входа, на который поступает сигнал (ut) и записывается в виде [1]:
A(q)y(t) = B(q)u(t) + e(t), y(t) = G(q,Q)u(t) + H(q,Q)e(t) =
и (t) + -^—e(t)
Aq)
где Б(с¡) и А(с) — полиномы: А(с) = 1 + + а^-1 + ... + атс[ш и Б(с) = ^ + Ь^1 + + ... + Ьпас[па+1, па и пЬ — порядки полиномов. Когда внешний вход системы отсутствует, получим AR-модель, то есть и( = 0, У = UOL(t), пЬ = 0:
А(сМ*) = еМ.
Предсказатель для ARX-модели:
у(Щ = Б(с)и(г) + [1- А(с)]уМ,
или У^/0) = фт(^)0,
где 0 = [а1,...,апа, Ь1,...Ьпь]т вектор параметров, подлежащих отысканию, ф(£) = [-у^ - 1),..., -у^-па ),и^ - 1),...,и(^ - пь)]т. Предсказатель представляет скалярное произведение известного вектора данных ф($ и вектора параметров 0.
Для исследования процесса деградации выберем наихудшие значения параметра UOL (выходное напряжение низкого уровня) по результатам испытаний на долговечность в течение 150 тыс. ч выборки из 20 ТТЛ ИС типа 133ЛА8. Временной ряд замеров параметров ИС типа 133ЛА8 при испытаниях на долговечность в течение 150 тыс. ч имеет вид: 0; 1; 2; 3; 4; 5; 6; 7; 8; 9; 10; 11; 12; 13; 14; 15; 16 17; 18; 19; 20; 21; 22; 23; 24; 25; 26; 27; 28; 29; 30 35; 40; 45; 50; 60; 70; 80; 90; 100; 110; 120; 130; 140; 150 тыс. ч. Общее число замеров N = 45, средняя продолжительность между замерами 150 тыс. ч/45 составляет 3333 ч. Отказов за время испытаний (достигнута величина гарантийной наработки) зафиксировано не было. Временной ряд деградации параметра UOL имеет вид:
0.283 0.286 0.399 0.287 0.399 0.280 0.290 0.295 0.321 0.313 0.283 0.392 0.379 0.399 0.339 0.286 0.282 0.296 0.400 0.400 0.348 0.284 0.285 0.281 0.390 0.348 0.287 0.391 0.400 0.342 0.351 0.287 0.280 0.308 0.314 0.287 0.298 0.361 0.271 0.387 0.270 0.264 0.258 0.325 0.248
Предположим, что процесс деградации параметра UOL может быть описан моделью авторегрессии второго порядка (АР(2)-моде-лью или ARX-моделью с отсутствующим внешним входом системы и na =2). Используя редактор рабочей области сессии (переменные в MatLab хранятся в особой области памяти, именуемой рабочей областью, то есть в виде бинарных файлов, с расширением (*.mat) или Mat-files), необходимо отредактировать любой из имеющихся в системе Mat-файлов. Сохраним файл с именем DATA_IN.mat. С помощью команды load загрузим его в рабочую область. Подав команды в командном окне y = iddata(in); m = arx(y,[2]) вычислим коэффициенты ARX-модели: A(q) = 1 - 0.526 qA-1 - 0.4625 qA-2 (рис. 1). Запишем полученное уравнение авторегрессии в виде: y(t) - 0.526y(t- 1) - 0.4625y(t- 2) = e(t), где yt — выходной сигнал в моменты времени t = 1,2,...N, e(t) — белый шум, воздействующий на систему. Параметры модели оцени-
ваются линейным МНК посредством метода QR-декомпрессии для решения системы линейных уравнений.
В работах [2, 3] показано, что коэффициенты, минимизирующие дисперсию ошибки предсказания, совпадают с коэффициентами AR-модели формирования сигнала. Поэтому для линейного предсказания воспользуемся вычисленными коэффициентами AR-модели.
На рис. 2 показано линейное предсказание сигнала КИХ-фильтром (фильтр с конечной импульсной характеристикой, или finite impulse response, FIR) в системе MatLab/Simulink. Предварительно необходимо загрузить ряд деградации параметра в рабочую область (рис. 1) и задать в блоке параметров КИХ-фильтра коэффициенты: 0.526 и 0.463. Показано два примера использования одного и того же КИХ-фильтра, но с различным способом построения. Первый способ — используется встроенная функция (прямая реализация фильтра, Direct Form FIR) из пакета DSP Blockset. Второй — с использованием набора элементарных функций: задержка на такт, умножитель, сумматор. Эллипсом на рис. 2 обозначен КИХ-фильтр из элементарного набора функций. В качестве коэффициентов КИХ-фильтра возьмем коэффициенты ARX-модели. В первом способе коэффициенты задаются внутри блока, во втором — реализуются на умножителях. На рис. 3а показаны исходный сигнал (ряд деградации параметра UOL), предсказанный сигнал КИХ-фильтром типа DF из библиотеки DSP BlockSet и ошибка предсказания сигнала. На рис. 3б представлен предсказанный сигнал КИХ-фильтром из элементарного набора функций. В обоих случаях, независимо от способа построения фильтров, предсказанный сигнал одинаков. Более подробно с описанием блоков библиотеки DSP BlockSet можно ознакомиться на математическом образовательном сайте www.exponenta.ru или с использованием справочной информации Signal Processing BlockSet for use with Simulink. User's Guide. Version 6.
Существует соотношение между параметрами АР-модели и автокорреляционной функцией. Это соотношение известно под названием уравнений Юла-Уокера (формула внизу): где RXX(0),...,RXX(p) оценки автокорреляционной функции. Алгоритм Левинсона-Дербина обеспечивает эффективное решение выше приведенной системы уравнений. Алгоритм позволяет рекурсивно вычислять набор
параметров (я11, с12}, {a21, a22 af
ap2,.
pv
,app, ap}, где p — порядок авторегрессии.
Рис. 1. Загрузка ряда деградации параметра и0[_ в рабочую область сессии
Предсказанное значение КИХ-фильтром
Рис. 2. Линейное предсказание сигнала КИХ-фильтром в системе MatLab/Simulink
Параметры {а21,а22,...,а^^} называют коэффициентами отражения и обозначают через
{ВД,...,.^}.
На практике решетчатые структуры фильтров используются в системах линейного предсказания. На рис. 4а показано применение алгоритма Левинсона-Дербина для оценки коэффициентов отражения решетчатого КИХ-фильтра. На рис. 4б показано прогнозирование КИХ-фильтром типа DF (прямая фор-
R
Rxxiі) R
jsar(O) Rxx( 1)
ляг(О)
• Rxx(-P) Rxx(-(P-1))
Ахх (p) ^XX
(p-1) ...Яж(0)
1 а2
«1 = 0
.ар. 0
ма), коэффициенты которого предварительно вычислены с использованием ARX-модели.
Иногда в задачах прогнозирования требуется вычислять не прогнозные значения временного ряда, а ошибки прогнозирования. На рис. 4в показано предсказание ошибки решетчатым нерекурсивным КИХ-фильтром со структурой скользящего среднего (Moving Average). Навхо-ды фильтра подаются ряд деградации и коэффициенты отражения, на выходе имеем ошибку предсказания. Такие фильтры называют «обеляющими» или «обратными». На рис. 4г показано восстановление сигнала решетчатым рекурсивным БИХ-фильтром, со структурой
170
САПР
Рис. 3. Временные диаграммы работы КИХ-фильтра: а — исходный сигнал (ряд деградации параметра U0¡), предсказанный сигнал КИХ-фильтром из DSP BlockSet, ошибка предсказания; б — предсказанный сигнал КИХ-фильтром из элементарного набора функций
Рис. 4. Обработка сигнала цифровыми фильтрами: а — применение алгоритма Левинсона-Дербина для оценки коэффициентов решетчатых фильтров; б — прогнозирование поведения сигнала с использованием КИХ-фильтра; в — предсказание ошибки решетчатым КИХ-фильтром со структурой скользящего среднего; г — восстановление сигнала решетчатым БИХ-фильтром с авторегрессионной структурой
АР-модели. На входы фильтра подаются ошибки предсказания и коэффициенты отражения. Выходом фильтра является восстановленный сигнал (выходной сигнал). Согласно АР-модели выходной сигнал формируется путем пропускания дискретного белого шума (в нашем случае ошибки предсказания «обеляющего» фильтра) через «чисто рекурсивный» или «чисто полюсный» фильтр.
Вывод результатов моделирования осуществляется в векторы-столбцы. Следует заметить, что коэффициенты КИХ-фильтров, вычисленные с использованием ARX-модели и рекурсивного алгоритма Левинсона-Дербина, несколько отличаются, что объясняется различными методами нахождения параметров авторегрессионной модели АР(2) формирования сигнала. На рис. 5 представлены параметры блоков, используемых на рис. 4.
Использование Калмановской фильтрации для моделирования процесса деградации параметров ТТЛ ИС
Для большинства реальных систем формирование моделей на основе физических соображений проще осуществлять в непрерывном, а не в дискретном времени. Калмановская фильтрация — метод предсказания поведения сигнала по имеющимся наблюдениям у, содержащим ошибку , так что у + V. Термин «фильтрация» указывает на стремление удалить шум V(из наблюдений так, чтобы получить наилучшую оценку сигнала Самым известным применением фильтра Калмана было использование его в исторической миссии «Аполлона-11» в качестве системы слежения за космическим кораблем и лунным модулем. Типичное применение — обработка данных метеорологических наблюдений.
В пространстве состояний связь между входными сигналами, шумами и выходными сигналами записывается в виде системы дифференциальных или разностных уравнений первого порядка посредством введения вспомогательного вектора состояний х.
В основе метода Калмановской фильтрации лежит общая линейная модель (ОЛМ) с дискретным временем, описывающая состояние системы в фазовом пространстве. Для системы с входом и имеем: х+ = А(0)х( + + Б(0)и( + где х(-(пх1) — вектор фазовых
переменных состояния системы, w (-(пх1) — вектор шума системы (объекта).
С уравнением системы связано уравнение измерения: у( = С(0)х( + V, где у( - (тх1) —вектор измерений, С(0) - (тхп) матрица измерений, V(- (пх1) — вектор ошибок измерений, А(0), Б(0), С(0) — параметры, описывающие систему (матрицы), 0 — вектор параметров, и>р V (— интерпретируются как помехи, могут являться сигналами известной физической природы. Задача состоит в том, чтобы оценить фазовую переменную х из зашумленных наблюдений у1,у2,...,у(.
Block Parameters: Overlap Analysis Windows
Buffer (mask) (Ink)
Convert scalar samples to a frame output at a low« sample rate You can also convert a frame to a smaler or larger see with optional overlap For calculation of sample delay, see the rebuff erdelay function
Parameters
Output tnifer size (per channel):
Buffer overlap
Cancel I Help I _____________________________I
Block Parameters: All-pole Filter
Digital Filer (mask) (ink)
Independently filter each channel of the input over time using a specified digital Her implementation Fie# coefficients can be specified using either tunable mask dialog parameters or separate input ports (useful for itfne-varyng coefficients)
T me-varymg coefficients are supported at one of I wo possible rates of coefficient update n frame-based processng. The coefficients can update at a rale of one filer per frame (te. they stay constant for the duration of an nput frame and change for the next frame! Oi, the coefficients can update at a rate of one filer per sample (le. they change with every sample n the nput frame)
Parameters
T ransfer function type; [j|
Fit! structure: (Died lorn
¡3
■3
Coefficient source: | Specify via dialog Numéralo« coefficients:
13
I[0526 0 453)
|7 First denominator coefficient « 1 remove 1 /a0 term n the structure
Coefficient update rater |0ne fier pet frame Initial condtons:
Help
Block Parameters: Autocorrelation
Autocorrelation (mask) (ink)
Computes the autocorrelation of vectors or frame based matrices For frame-based inputs, it computes along the inpul's columns
When set to 'Compute al non-negative lags', it computes using lags in the range (0, lengthfnput)-1) Otherwise, it computes using lags in the range (0, maxLag). where you specify the value of maxLag in "Maximum non-negative lag'
Parameters
Г” Compute al non-negative lags
Maxmum non-negative lag (less than inpul length)
Scaling: | Biased Computation domain: prime
Caned I Help
Apply
б)
Block Parameters: FIR
Digital Fie* (mask) (Ink)
Independently lier each channel of the nput over I roe using a speofied cigital fier implementation Fier coefficients can be speeded usng «(her tunable mask ciaJog parameters or separate input ports (usefvi for time-varying coefficients)
Time-varying coefficienls are supported al one of two possible rates of coefficient update in frame-based processing The coefficients can update at a rate of one filter per frame (Le. they stay constant for the duration of an input frame and change for the next frame) Or, the coefficients can update at a rate of one fier per sample (i e. they change with every sample n the nput frame)
Parameters
Transfer function type- [
Fier structure. I Lattice MA Cixttrxrt taiar | Input patts)
"3
"3
{7 Гпи denominator coefficient _ ; / remove 1/a0 term in the structure Coetfaent update rate: | One (*et pel Irame 3
Initial condtoons:
P
Cancel I Help I App1,' |
Block Parameters: Levinson-Durbin
Levinson-Oubm (mask) (ink)
Solve Hermitian Toepitj system of equations using the Levrison-Durbn recursion Input is lypicaly a vector of autocorrelation coefficients with lag 0 as the first element Outputs polynomial coefficients A, reflection coefficients K. and/or the jaediclion error power P
Parameters Outpi*(s) [i
13
W Output predction error power (P)
W li the value of lag 0 is zero, A»(1 zeros], К «{zeros], P-0
Help
Block Parameters: IIR
Digital Fier (mask) (Ink)
Independently fie* each channel of the input over time using a specified rigital fier implementation Fiei coefficients can be specified using either tunable mask cfcalog parameters or separate input ports (useful for rime varying coefficients).
T*ne varying coefficients are supported at one of two possible rates of coefficient update n frame-based processing T he coefficients can update al a rale of one fier pet frame (ie. they slay constant for the duration of an inptf frame and change for the next frame) Or, the coefficients can update al a rate of one fier per sample fie. they change with every sample n the input frame)
Parameters
T ramier function type: I Fier structure I Lattice АЯ
13
■З
Coefficient lomee I Input poct(s) 3
!✓ Fічі denominator c-oefficren* - I, remove 1 /«9 term m Ihe structure
Coefficient update rate: | One fier p Initial condtions.
"З
Help
Рис. 5. Параметры блоков, используемых на рис.4:
а) — формирование буфера глубиной 45 отсчетов на канал (фрейм данных); б — блок параметров автокорреляционной функции, лаг равен 2;
в) — блок параметров алгоритма Левинсона-Дербина, сконфигурирован на вычисление коэффициентов A и K и мощности ошибки предсказания P;
г) — блок параметров КИХ-фильтра (чисто нулевая модель), структура фильтра — прямая форма, коэффициенты фильтра 0.526 и 0.463 (числитель функции передачи);
д) — решетчатый нерекурсивный КИХ-фильтр со структурой скользящего среднего (Moving Average); е) — решетчатый рекурсивный БИХ-фильтр со структурой авторегрессии
Калмановская фильтрация используется для построения адаптивных фильтров. Адаптивные устройства обработки сигнала действуют по принципу замкнутого контура (обратной связи). Входной сигнал фильтруется или взвешивается в программируемом фильтре для получения выходного сигнала ¥р который затем сравнивается с полезным, стандартным или обучающим сигналом У( для нахождения сигнала ошибки е. Затем этот сигнал ошибки используется для корректировки весовых параметров обрабатывающего устройства (обычно итеративным методом) с целью постепенной минимизации ошибки.
Рассмотрим применение скалярного фильтра Калмана. Предположим, что процесс деградации параметра UOL описывается моделью АР(1). Устройство Калмановской оценки реализует процесс параметрического оценивания, основанный на авторегрессивной (АР) модели процесса генерации сигнала. АР-модель процесса первого порядка данно-
го типа показана на рис. 6а, соответствующая модель измерений представлена на рис. 6б. Модель измерений (уравнение измерения) представляет усилительное звено и источник аддитивного белого шума
Y
Модель генерации сигнала (уравнение системы) описывается АР(1)-моделью:
X
aXt-l + w-l,
где — обусловлен несовершенством представления Х( через Xt_1, константа а считается известной. Уравнение системы можно применять для прогноза X: X = aXt_1.
В общем виде наилучшая линейная несмещенная оценка X, формируется взвешиванием линейной комбинации двух имеющихся оценок:
X = L X-1 + ад,
где L (и К — зависящие от времени весовые матрицы, которые выбираются так, чтобы оценка в каждый момент времени была несмещенной и имела максимальную дисперсию. В нашем случае рекурсивная формула оценки первого порядка имеет вид:
X = ь X _!+ ад,
где коэффициенты Ь, к( изменяющиеся во времени, отыскиваются линейным методом наименьших квадратов, путем дифференцирования среднеквадратичной ошибки Е[X - X]2 по Ьи к. Проведя преобразования получим:
X = а^ _! + к( [У( - асХ_1].
Первый член aXf-1 предсказывает Xt, а второй член acXj-1 корректирует X(на основании оценки ошибки с учетом Калмановского коэффициента к . Структурная схема этого устройства оценки показана на рис. 6г. Обоб-
а)
в)
г)
д)
е)
172
САПР
Рис. 6. Калмановская фильтрация:
а — рекуррентная модель генерации сигнала первого порядка; б — модель схемы измерения данных; в — обобщенная структурная схема рекурсивного устройства оценки первого порядка; г — блок-схема фильтра Калмана первого порядка
щенная структурная схема рекурсивного устройства оценки первого порядка представлена на рис.бв. В общем виде задача прогнозирования сигнала Х( записывается так [1]:
Х(* + 1,0) = Л(0)Х(*,0) + В(0)и(£) + +К(0)№) - С(0)Д<,0)],
У(*,0) = С(0)Х(*,0),
где К(0) определяется как:
К(0) = [Л(0)Р(0)]Ст(0) +
+ ^12(0)][С(0)Р(0)СТ(0) + Я2(0)]-1,
Р — матрица ковариации ошибки в оценке состояния
-(0) = Ё[х(£) - х(*,0)][х(*) - 5с(£,0)]т,
где Я1, Я2, Я12 — ковариации шума объекта и шума измерений: Ем'(ґ)'мт(ґ) = Л1(0), Еу(ґ)ут(ґ) = Л2(0), Ем/(ґ)ут(ґ) = Л12(0).
В справочной системе Ма^аЬ фильтр Кал-мана в матричной форме записывается в следующем виде:
К{п) = К{п-1)—§(п)ин(п)К(п -1) + <2Р, К(п-\)и(п)
эффициентов (отводов) КИХ-фильтра на шаге п; у(п) — профильтрованные выходные значения фильтра Калмана на шаге п; е(п) — ошибки на шаге п; й(п) — эталонные значения; (м — матрица корреляций шума измерений и (>р — матрица корреляций шума объекта. Матрицы (м и (>р формируются путем заполнения диагонали скалярными значени-
ями дисперсии белого шума vt уравнения измерения и дисперсии белого шума wt процесса генерации сигнала. Общая идея построения фильтра Калмана показана на рис. 7а, также показан адаптивный фильтр в общем виде (б) и использование адаптивного фильтра для построения прогнозирующего фильтра (в). Кроме адаптивного алгоритма Калмана, в качестве адаптивных алгоритмов используются простейшие алгоритмы по критерию наименьшего среднеквадратичного отклонения (СКО) (Least Mean Square, LMS), нормализованные LMS (NLMS)-алгоритмы, рекурсивные алгоритмы по критерию наименьших квадратов (Recursive Least Squares, RLS) и быстрые RLS-алгоритмы [8]. Данные алгоритмы различаются качеством и вычислительной сложностью. Более высокое качество требует применения более сложных алгоритмов. Адаптивные фильтры также могут быть использованы для идентификации систем и для удаления шума из сигнала.
Из теории цифровых фильтров известно, что фильтр Калмана эффективен в случае, если полезный сигнал (в нашем случае ряд деградации параметра UOL) — медленно меняющаяся функция времени (синусоида), а помеха — некоррелированный шум. Рассмотрим применение фильтра Калмана для прогнозирования процесса деградации параметра UOL ИС типа 133ЛА8. На рис. 8 показан блок параметров фильтра Калмана. Назначения портов блока фильтра следующие: In — вход
g{n) =
и н {п)К (и - 1)м(и) + Qh
у(п) = и11 (n)w(n), e(n) = d(n)-y{n), w(n +1) = w{n) + e(n)g(n),
где n — текущая итерация; u(n) — значение буферизованного входного сигнала In на шаге n; K(n) — матрица корреляций ошибки в оценке состояния; g(n) — вектор усиления Калмана на шаге n; w(n) — вектор оценок ко-
Шум
объекта
1/Z
Буффер
w{n)
Произведение
Матрица переходов
О
Out
ж
Шум
измерений
Эталонный сигнал d(n)
Рис. 7. Адаптивные фильтры: а — общая идея построения фильтра Калмана; б — адаптивный фильтр в общем виде; в — использование адаптивного фильтра для построения прогнозирующего фильтра
фильтра, скалярная величина, которая буферизуется в вектор u(n); Out — выход фильтра, скалярная величина y(n); Err - e(n), скалярная величина или оценки ошибок на каждой итерации; Taps - w(n), вектор оценок коэффициентов КИХ-фильтра, которые меняются на каждой итерации n. Задаются следующие параметры: порядок КИХ-фильт-ра — 2, дисперсия белого шума vt уравнения измерения — 0.03, дисперсия белого шума wt процесса генерации сигнала — 0.05, первоначальные значения коэффициентов КИХ-фильтра (вектор отводов фильтра) — 0, первоначальное значение корреляционной матрицы ошибок в оценке состояния — 0.01 (первоначальное значение ).
На рис. 9 показано прогнозирование процесса деградации параметра UOL с использованием Калмановского адаптивного фильтра. Также показаны временные диаграммы работы фильтра Калмана, динамически обновляемые частотная характеристика фильтра и коэффициенты КИХ-фильтра второго порядка (импульсная характеристика).
Поскольку заранее известно, что процесс деградации параметра UOL описывается АР(2)-моделью, то в Калмановском фильтре используется КИХ-фильтр второго порядка. Значения дисперсий шумов заданы без физического истолкования. Следует заметить, что погрешность измерений параметров ИС типа 133ЛА8 при испытаниях на долговечность согласно ТУ на ИС данной серии была не хуже 3% для параметра UOL. На рис. 10 показано сравнение результатов моделирования процесса деградации с применением адаптивного фильтра Калмана и КИХ-фильтра второго порядка. В первом случае (рис. 10а) для моделирования процесса деградации используется библиотека DSP Blockset/Adaptive Filters, во втором
Q
Kalman Adaptive Fier (mask) (Ink)
One-step Kalman predctor algorithm ten adaptive FIR fAerng of nput «gnat
If 'Adapt input' is checked and the Adapt signal is zero. the algorithm teenporarity cSsables adaptation of the fief coefficients. A non-zero signal enables fie« adaptation
If the Tiese* input' is enabled and the Reset signal transitions appropriately, the lie« coefficients are reset to the initial conditions.
Parameter;
OR fie« length
IT-
Measurement noise variance:
|q<¡3
Process noise vaiiance
|ôœ
Initial value of fie« taps
Initial error correlation matrix (ÔÔÎ
Г Adapt input
Reset input: I None 3
OK ~| Cancel I Help I _______________________________________I
Рис. 8. Блок параметров фильтра Калмана
Block Parameters: Kalman Adaptive Filter
Рис. 9. Прогнозирование с использованием Калмановского адаптивного фильтра в MatLab/Simulink: а — блок-схема метода; б — временные диаграммы; динамически обновляемые частотная характеристика фильтра (в) и коэффициенты КИХ-фильтра второго порядка (импульсная характеристика) (г)
Рис. 10. Сравнение результатов моделирования процесса деградации с использованием адаптивного фильтра Калмана и КИХ-фильтра второго порядка, построенных: а — с использованием библиотеки DSP Blockset/Adaptive Filters; б — с использованием функции adaptkalman и ARX-модели второго порядка (АР(2))
174
(рис. 10б) — функция adaptkalman. Судя по поведению прогнозных значений и по ошибкам работы фильтра Калмана, его работу можно признать удовлетворительной, однако прогнозные значения фильтра пересекают верхнюю границу параметрического отказа 0,4 В. Налицо явное усиление сигнала фильтром Калмана. Сравнивая результаты, видим, что КИХ-фильтр лучше моделирует процесс деградации. Отчасти это может быть объяснено тем, что ковариационная матрица шумов процесса генерации сигнала заранее неизвестна и число отсчетов явно недостаточно. По мере того как фильтр Калмана обучается прогнозированию, выходной сигнал у приближается к входному сигналу х, ив районе 40-45 отсчетов прогнозы начинают совпадать, а сигнал ошибки е начинает уменьшаться. Пример использования функции adaptkalman показан ниже.
d=x;
delay = zeros(1,1);
z= delay; w0 = zeros(1,2);
K0 = 0.01*eye(2); Qm = 0.03;
%Ряд деградации %о(эталонные значения) %Задержка на один шаг %(нулевой элемент)
%оПервоначальные значения %коэффициентов КИХ-фильтра %2-го порядка %Корреляционная матрица %ошибок в оценке состояния
%Дисперсия белого шума %уравнения измерения %(скалярная величина)
Ор = 0.05*еуе(2); %Ковариационная матрица
%дисперсии белого шума %процесса генерации сигнала т=^ d(1:44)]; %Объединение задержки
%(один отсчет) и сигнала х %(44 отсчета) для создания %сигнала т
Б = mitkalman(w0,K0,Qm,Qp); %Инициализация параметров %функции фильтра Калмана [у,е,Б] = adaptkalman(in,d,S); %Фильтр Калмана
Выводы
Использование цифровых фильтров для моделирования процесса деградации параметров ТТЛ ИС эффективно только в случае, если происходит непосредственное наблюдение за процессом деградации и сам процесс адекватно описывается авторегрессионными моделями. Цифровые фильтры, несмотря на их связь с линейными моделями временных рядов, не могут строить прогнозы за пределы временного ряда, поэтому они не могут быть использованы для прогнозирования с упреждением, а только для моделирования.
Работа выполнена по программе гранта РФФИ 05-08-01225-а.
Литература
1. Льюнг Л. Идентификация систем. Теория для пользователя: Пер. с англ. / Под ред. Я. З. Цып-кина. М.: Наука, 1991. 432 с.
2. Кей С. М., Марпл С. Л. Современные методы спектрального анализа: Обзор. ТИИЭР, 1981, т. 69, № 11, с. 5-43.
3. Цифровая обработка сигналов / А.Б. Сергиен-ко. СПб.: Питер, 2003. 608 с.
4. Бокс Дж., Дженкинс Г. Анализ временных рядов. Прогноз и управление. М.: Мир, 1974. 406 с.
5. Справочник по прикладной статистике. В 2-х т. Т. 2: Пер. с англ. / Под ред. Э. Лойда, У. Ледер-мана, С. А. Айвазяна, Ю. Н. Тюрина. М.: Финансы и статистика, 1990. 526 с.
6. Уидроу Б., Стирнз С. Адаптивная обработка сигналов: Пер. с англ. М.: Радио и связь, 1989. 440 с.
7. Дьяконов В. П., Абраменкова И. В., Круглов В. В. MatLab 5.3.1 с пакетами расширений. Под ред. проф. В. П. Дьяконова. М.: Нолидж, 2001. 880 с.
8. Filter Design Toolbox. For use with MatLab.
