Верификации прогнозов АРПСС-моделей временных рядов применяемых для прогнозирования долговечности интегральных схем ИС
В предыдущих публикациях автора показано, что модели цифровых фильтров, модели временных рядов и нейронные сети могут быть использованы для прогнозирования процесса деградации контролируемых электрических параметров ТТЛ ИС при испытаниях на долговечность.
Цифровые адаптивные фильтры способны строить лишь одношаговые прогнозы и не пригодны для прогнозирования времени наступления параметрических отказов. Цифровые фильтры эффективны в задачах слежения за процессом деградации электрических параметров ИС. Модели временных рядов (модель авторегрессии проинтегрированного скользящего среднего или модель АРПСС(p,d,q), где p — число параметров авторегрессии, d — порядок дифференцирования, q — число параметров скользящего среднего) позволяют прогнозировать время наступления параметрического отказа по траектории процесса деградации контролируемого параметра [1, 2]. Нейронные сети способны эффективно решать более широкий круг задач, в том числе задачи слежения и прогнозирования процесса деградации [5].
В данном разделе рассматривается графическая верификация прогнозов модели АРПСС(0,1,2), подогнанной к ряду деградации наихудших (экстремальных) значений параметра UOL (выходное напряжение низкого уровня) ТТЛ ИС типа 106ЛБ1 выпуска января 1980 года при испытаниях на долговечность продолжительностью 120 тыс. ч.
Главная особенность такого ряда — отсутствие корреляции между изменением параметра и временем испытаний: UOL = 0,3371–0,0007×t, где t — время испытаний. Коэффициент детерминации R2 = 0,0329 (коэффициент детерминации показывает, что в данной ситуации только 3,29% общей вариабельности объясняется изменением значений параметра UOL под воздействием времени испытаний, в то время как на все остальные факторы приходится 96,71%), коэффициент корреляции r = –0,1812. Таким образом, налицо вероятностная природа ряда деградации (процесс деградации развивается по законам теории вероятности).
Прогнозирование осуществляется с использованием модуля «Анализ временных рядов и прогнозирование» системы Statistica for Windows компании StatSoft [6].
Ряд замеров параметра UOL ТТЛ ИС типа 106ЛБ1 при испытаниях на долговечность в течение 120 тыс. ч имеет вид: 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; 40; 50; 60; 70; 80; 90; 100; 110; 120 тыс. ч.
Ряд деградации параметра UOL содержит данные с пропусками. Число отсчетов ряда — 40 и они получены в неравные интервалы времени. Необходимое число отсчетов для построения АРПСС-модели — не менее 50 [7–10], оптимальное число замеров — 150, а в некоторых случаях — более 1000. Тем не менее, система Statistica for Windows позволяет построить АРПСС-модели при числе замеров не менее 15. Однако в этом случае использование АРПСС-модели носит демонстрационный характер. Для получения недостающих значений необходимо использовать анализ данных с пропусками. Система Statistica for Windows предлагает 5 стандартных методов заполнения пропусков: заполнение пропусков средним значением ряда; интерполяция по соседним точкам; среднее по N соседним точкам; медиана N ближайших значений ряда; заполнение пропусков прогнозами линейной регрессии. Однако могут быть использованы и другие методы получения недостающих значений, например, аппроксимация кубическими сплайнами или аппроксимация недостающих значений двухслойной нейронной сетью с радиальными базисными элементами.
Выбор того или иного метода для аппроксимации недостающих значений ряда деградации контролируемого электрического параметра ИС — достаточно сложная задача. Пользователь должен самостоятельно выбрать любой из подходящих методов.
Например, ряд замеров параметров ТТЛ ИС типа 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 тыс. ч, то есть общее число отсчетов равно 45. На рис. 1а показан ряд деградации параметра UOL ИС типа 133ЛА8 при испытаниях на долговечность в течение 150 тыс. ч с пропусками, заполненными методом интерполяции, на рис. 1б — с пропусками, заполненными прогнозами линейной регрессии.
Видно, что ряд с пропусками, заполненными методом интерполяции, соответствует нашему представлению о природе временного ряда. Ряд с пропусками, заполненными прогнозами линейной регрессии, имеет необычную форму, резко прослеживается линейный тренд в данных между опорными точками. На рис. 1в–д показаны сглаженные ряды с использованием скользящего среднего с интервалом 5 точек и скользящей медианы с интервалом 5 точек.
В данной работе предлагается использовать два метода получения недостающих значений — аппроксимация кубическими сплайнами (система Matlab) и метод интерполяции (система Statistica for Windows). Это позволяет для ряда деградации параметра UOL ТТЛ ИС типа 106ЛБ1 получить 121 отсчет с интервалом 1000 ч. При этом какая-либо дополнительная проверка на адекватность использования выбранного метода получения недостающих значений ни в системе Matlab, ни в системе Statistica for Windows не проводится.
Основная идея графической верификации прогнозов заключается в том, что сначала на основе некоторого рационального критерия или критериев необходимо разделить на подмножества имеющийся массив данных о деградации контролируемых параметров ИС, полученный на основании анализа измерений экстремальных значений на интервале времени [0, tu], где tu — длительность испытаний, а затем использовать одну или несколько частей этих данных для построения «предсказывающей» модели процесса деградации параметров ИС, оставшуюся же часть или части использовать для «проверки» («экзамена») этой модели, то есть для оценки того, насколько хорошо предсказывает эта модель. Прогнозирование долговечности ИС до наступления параметрических отказов произведем с использованием метода Бокса-Дженкинса [7–10]. В данном случае рациональным критерием, позволяющим разбить исходный ряд деградации на части, является периодичность замеров: до 30 тыс. ч замер параметров осуществляется через 1000 ч, а после 30 тыс. ч — через каждые 10 тыс. ч. Для того чтобы природа ряда деградации как можно меньше искажалась, применим методы получения недостающих значений лишь к той части временного ряда, где существуют пропуски.
а — интерполяция;
б — прогнозы линейной регрессии;
в — сглаживание скользящим средним ряда деградации, пропуски которого заполнены методом интерполяции;
г — сглаживание скользящим средним ряда деградации, пропуски которого заполнены прогнозами линейной регрессией;
д — сглаживание медианной ряда деградации, пропуски которого заполнены методом интерполяции
С этой целью ряд деградации параметра UOL ТТЛ ИС типа 106ЛБ1 разбивается на две части: от 0 до 30 тыс. ч— исходные данные, с 30 до 120 — аппроксимация кубическим сплайном. Например, между 30 и 40 тыс. ч для выбранного шага 1000 ч необходимо 9 недостающих значений.
Затем ряд продолжительностью 120 тыс. ч, составленный из исходных и аппроксимированных значений, разбивается на 4 ряда. Для ряда 1 (30 тыс. ч испытаний) строится прогноз с 30 до 120 тыс. ч, для ряда 2 (60 тыс. ч испытаний) строится прогноз с 60 до 120 тыс. ч, для ряда 3 (90 тыс. ч испытаний) — прогноз с 90 до 120 тыс. ч, а для ряда 4 (120 тыс. ч испытаний) — прогноз со 120 тыс. ч до наступления параметрического отказа. К каждому ряду подгоняется своя АРПСС-модель, с помощью которой и осуществляется прогнозирование. На рис. 2 показана графическая верификация прогнозов модели АРПСС(0,1,2) для ряда деградации параметра UOL.
Подогнанные модели к четырем рядам процесса деградации параметра UOL с учетом оценок дисперсии белого шума σ2a сведены в таблице 1. В графе «Идентифицированная модель» в скобках обозначены стандартные ошибки оцениваемых параметров АРПСС-модели. В таблице 1 также приведена сводка совокупного критерия согласия Q, примененного к остаточным ошибкам модели АРПСС(0,1,2). Сравнивая параметры θ1 и θ2 модели АРПСС(0,1,2) и их стандартные ошибки для рядов 1, 2 и 3, видим, что они устойчивы на временном диапазоне продолжительностью 120 тыс. ч.
Эвентуальная прогнозирующая функция модели АРПСС(0,1,2) для упреждений l ≥ 2 показывает, что прогнозы модели АРПСС(0,1,2) рассматриваемых временных рядов 1, 2, 3, 4 в момент времени t будут следовать прямой линии, параллельной временной оси. В соответствии с теорией временных рядов лучшими будут два первых прогнозных значения, поэтому на рис. 2а прогнозы модели АРПСС(0,1,2) представлены прямыми линиями. На рис. 2а также показан ряд деградации параметра UOL с прогнозами моделей АРПСС, подогнанных к его частям (ряды 1–3) и 90%-ными доверительными интервалами. За параметрический отказ принято достижение верхней 90%-ной границей доверительного интервала отказового уровня, установленного в технических условиях (0,35 В).
Сравнение полученных прогнозов с рядом 4 (рис. 2) показало, что прогнозы, построенные с использованием моделей АРПСС, дают удовлетворительные результаты. Таким образом, анализ рис. 2а показывает, что модель АРПСС(0,1,2) адекватно описывает процесс деградации параметра UOL в течение 120 тыс. ч испытаний. Прогнозирование долговечности ТТЛ ИС типа 106ЛБ1 до наступления параметрических отказов следует вести с использованием модели АРПСС(0,1,2), идентифицированной для ряда 4.
Рассмотрим проверку на «избыточность» (перепараметризация модели) модели АРПСС(0,1,2). Испытание модели АРПСС(0,1,2) для ряда 4 путем введения избыточных параметров скользящего среднего θi, i>2 согласно рекомендациям работ [6–9] показало, что воспроизводимость значений ряда деградации параметра UOL не улучшается. Поиск альтернативных (конкурирующих) моделей АРПСС идентифицированной модели АРПСС(0,1,2), связанный с тем, что малые изменения оператора скользящего среднего могут быть почти точно скомпенсированы соответствующими изменениями оператора авторегрессии [6–9], дал положительный результат. Выявлена конкурирующая модель АРПСС(1,0,0): Xt = 0,264+0,669Xt–1+at. Все параметры конкурирующей модели статистически значимы. Данная модель более экономична, чем основная модель АРПСС(0,1,2).
а) 1— прогноз для ряда 1; 2— прогноз для ряда 2; 3— прогноз для ряда 3;
4 — прогноз для ряда 4;
б) сравнение прогнозов модели АРПСС (0, 1, 2) и конкурирующей модели АРПСС (1, 0, 0) с константой;
КТ — контрольная точка (наихудшее значение параметра после 130 тыс. ч испытаний
Для подтверждения или опровержения достоверности прогнозов деградации параметров UOL и UOH ТТЛ ИС типа 106ЛБ1 АРПСС-моделей был проведен контрольный замер параметров этих же ИС по прошествии еще 10 тыс. ч. Таким образом, фактическая наработка этих ИС при испытаниях на долговечность составила 130 тыс. ч. Отказов за это время не зафиксировано.
Сравнение точечного прогноза процесса деградации параметра UOL на момент времени 130 тыс. ч, построенного с использованием модели АРПСС(0,1,2) по наработке 120 тыс. ч, с контрольной точкой (КТ) на момент времени 130 тыс. ч показало, что прогнозное значение составило 0,300 В, и оказалось завышенным на величину 0,018 В (рис. 2б). Однако КТ укладывается в 90%-ный доверительный интервал модели АРПСС(0,1,2). В целом можно считать, что получено хорошее согласие точечных прогнозов с КТ.
На рис. 1б показано сравнение прогнозов модели АРПСС(0,1,2) с конкурирующей моделью АРПСС(1,0,0): Xt = 0,264+0,669Xt–1+at. По точечным и вероятностным прогнозам конкурирующей модели параметрический отказ не наблюдается. Точечный прогноз с использованием конкурирующей модели на момент времени 130 тыс. ч дает величину 0,265 В, которая оказывается заниженной на величину 0,017 В относительно КТ (0,282 В). Однако если строить долгосрочные прогнозы на глубину, в 2–3 раза превышающую глубину исходного ряда деградации, то они могут стать небезопасными и противоречить нашему представлению о процессе деградации, так как точечные прогнозы начинают носить асимптотический характер, то есть прогнозные значения перестают изменяться.
В таблице 3 приведены результаты прогнозирования процесса деградации параметра UOL по данным испытаний на долговечность ТТЛ ИС типа 106ЛБ1 в течение 130 тыс. ч на глубину 30 тыс. ч с использованием АРПСС-моделей. В случае аппроксимации кубическими сплайнами используются только наихудшие значения параметра UOL в выборке, а для метода интерполяции используются наихудшие, средние и наилучшие значения.
а) ряд деградации без пропусков;
б) недостающие значения ряда получены методом интерполяции:
1 — наихудшие значения;
2 — средние значения;
3 — наилучшие значения
На рис. 3а показано прогнозирование процесса деградации параметра UOL ТТЛ ИС типа 106ЛБ1 на глубину 30 тыс. ч для ряда деградации без пропусков, составленного из наихудших значений с использованием основной модели АРПСС(0,1,2): ∇1Xt = at–0,379at–1–0,599at–2 и конкурирующей модели АРПСС(1,0,0): Xt = 0,260+0,403Xt–1+at. Обе модели дают приемлемый результат. На рис. 3б показано прогнозирование процесса деградации параметра UOL по наихудшим, средним и наилучшим значениям в выборке. Недостающие значения ряда заполнены методом интерполяции. Все прогнозы носят асимптотический характер.
Анализируя таблицу 3, можно сделать вывод, что общий вид модели АРПСС(0,1,q) для рядов деградации параметра UOL, составленных из наихудших, наилучших и минимальных значений параметра, в выборке не изменился, меняется лишь число параметров скользящего среднего. Разница между последним точечным прогнозным значением на глубину 30 тыс. ч для ряда деградации без пропусков с использованием модели АРПСС(0,1,2): ∇1Xt = at–0,379at–1–0,599at–2 и для ряда деградации с пропусками, заполненными методом интерполяции с использованием модели АРПСС(0,1,2): ∇1Xt = at–0,291at–1–0,557at–2 составляет 0,025 В.
Представляет интерес рассмотрение достоверности «грубых» прогнозов АРПСС-модели (долгосрочных прогнозов до наступления параметрического отказа), построенные к ряду деградации параметра UOL ТТЛ ИС типа 106ЛБ1 при испытаниях на долговечность. На рис. 4 показано поведение параметра UOL ТТЛ ИС типа 106ЛБ1 при испытаниях на долговечность (наихудшие значения). Ряд деградации с пропусками заполняется методом интерполяции (кривая 1) и аппроксимацией кубическим сплайном (кривая 4). Подбирается адекватная модель АРПСС(0,1,2) к ряду деградации с пропусками, заполненными методом интерполяции, параметры модели оцениваются методом максимального правдоподобия. Далее осуществляется построение прогноза к части ряда с 30 до 130 тыс.ч, прогнозирование назад (кривая 2) с использованием модели АРПСС(0,1,2), подогнанной по целому ряду (кривая 1). Видим, что точечный прогноз (кривая 2) модели АРПСС(0,1,2) достаточно хорошо отражает поведение ряда в будущем.
Для сравнения качества прогнозов модели АРПСС(0,1,2), построенных с 30 до 130 тыс. ч, то есть на 100 тыс. ч вперед (глубина прогнозирования почти в 3 раза превышает исходную глубину ряда), рассмотрим статистически возможные траектории (кривые 5, 6, 7, 8) модели АРПСС(0,1,2) ∇1Xt = at–0,291at–1–0,557at–2, полученные методом Монте-Карло, путем генерации значений белого шума at с нормальным законом распределения (μ = 0, σ2 = 0,005).
Из рис. 4 видно, что вероятностные траектории ряда укладываются в 90%-ный доверительный интервал модели АРПСС(0,1,2). Кроме того, видно, что чем на большую глубину будем прогнозировать, тем сильнее отстоят друг от друга вероятностные траектории, и тем быстрее будут раскрываться 90%-ные доверительные интервалы АРПСС-модели. Из рис. 4 также видно, что в районе 100 тыс. ч по верхней границе 90%-ного доверительного интервала модели АРПСС(0,1,2) фиксируется параметрический отказ, тогда как по экспериментальным данным он не фиксируется (кривая 1).
1 — ряд деградации с пропусками, заполненными методом интерполяции, и прогноз с использованием модели АРПСС(0,1,2), построенный к части ряда (от 30 до 130 тыс. ч) прогнозированием назад (2) с 90%-ным доверительным интервалом (3);
4 — аппроксимация ряда кубическим сплайном;
5, 6, 7, 8 — статистически возможные траектории ряда 1, полученные методом Монте-Карло путем генерации значений белого шума модели АРПСС(0,1,2)
Эвентуальная прогнозирующая функция модели АРПСС(0,1,2) для ряда деградации параметра UOL ТТЛ ИС типа 106ЛБ1 показывает, что модель вычисляет «точно» только первые два прогнозных значения, далее прогнозы будут следовать прямой линии, параллельной временной оси. Графическая верификация показывает, что «долгосрочные» точечные прогнозы модели АРПСС(0,1,2) для ряда деградации параметра UOL ТТЛ ИС типа 106ЛБ1 вполне состоятельны. Сравнение с экспериментальными данными показывает, что пересечение доверительным интервалом границы параметрического отказа, фиксируемое как параметрический отказ, показывает, что степень доверия к точечным прогнозам упала до некого критического уровня, а долгосрочные точечные прогнозы не могут быть признаны состоятельными. Для ТТЛ ИС типа 106ЛБ1 для параметра UOL по рядам деградации без пропусков параметрический отказ не фиксируется; по 120 тыс. ч испытаний в случае заполнения пропусков аппроксимацией кубическими сплайнами прогнозируемое время составляет 150 тыс. ч, а по 130 тыс. ч испытаний в случае использования метода интерполяции для заполнения пропусков прогнозируемое время составляет более 1000 тыс. ч. По уравнению линейной регрессии параметрический отказ так же не фиксируется.
Литература
- Строгонов А. Прогнозирование деградации выходных параметров ТТЛ ИС. Часть I // Компоненты и технологии. 2005. № 8.
- Строгонов А. Прогнозирование деградации выходных параметров ТТЛ ИС. Часть II // Компоненты и технологии, 2005, № 9.
- Строгонов А. Использование цифровых фильтров для моделирования деградации выходных параметров ТТЛ ИС в системе Matlab/Similink // Компоненты и технологии. 2005. № 8.
- Строгонов А. Использование нейронных сетей для прогнозирования деградации выходных параметров ТТЛ ИС в системе Matlab/Simulink // Компоненты и технологии. 2006. № 1.
- Строгонов А. Использование линейной нейронной сети в задачах прогнозирования деградации выходных параметров ИС // Компоненты и технологии. 2006. № 2.
- Боровиков В. STATISTICA. Искусство анализа данных на компьютере: Для профессионалов. 2-е изд. (+CD). СПб.: Питер. 2003.
- Справочник по прикладной статистике. В 2-х т. Т. 2: Пер. с англ. М.: Финансы и статистика. 1990.
- Бокс Дж., Дженкинс Г. Анализ временных рядов. Прогноз и управление. М.: Мир. 1974.
- Боровиков В. П., Ивченко Г. И. Прогнозирование в системе STATISTICA в среде Windows. Основы теории и интенсивная практика на компьютере: Учеб. пособие. М.: Финансы и статистика. 1999.
- Грешилов А. А., Стакун В. А., Стакун А. А. Математические методы построения прогнозов. М.: Радио и связь. 1997.