Оценивание параметров эквивалентной схемы многополюсников с помощью многооткликовых моделей

УДК 513.73

С.А.Попов, А.Ф.Корчагин

ОЦЕНИВАНИЕ ПАРАМЕТРОВ ЭКВИВАЛЕНТНОЙ СХЕМЫ МНОГОПОЛЮСНИКОВ С ПОМОЩЬЮ МНОГООТКЛИКОВЫХ МОДЕЛЕЙ

A method is presented to evaluate the parameters of the multiterminal networks and their statistical characteristics on the basis of multiresponce regression models. The task of the multiterminal network parameters calculation is formulated in this case as a problem of coefficients estimates determination of the multiresponce model on the experimental data.

Введение

Большинство параметров эквивалентных схем определяются в статическом режиме расчетным путем по результатам измерений вольтамперных характеристик (ВАХ).Однако используемые до настоящего времени методики определения этих параметров не позволяют получать все параметры с достаточной точностью, а для некоторых параметров вообще не разработаны методики их оценивания по экспериментальным данным. Основные недостатки известных методов оценивания параметров состоят в том, что они не являются универсальными, не обеспечивают достаточную точность, не позволяют оптимизировать число измерений, не дают статистических характеристик полученных оценок параметров. Кроме того, некоторые методики оценивания параметров для упрощения расчетов используют различного рода допущения, например деление ВАХ на участки, в пределах которых некоторыми параметрами пренебрегают. В ряде ситуаций такой подход приводит к большим погрешностям оценивания параметров. Отсутствие статистических характеристик получаемых оценок фактически не дает возможности сравнения величин параметров индивидуальных приборов и не позволяет проводить сравнение параметров по партиям приборов.

Описываемый метод оценивания всех статических параметров эквивалентных схем и их статистических характеристик с помощью многооткликовых моделей позволяет получать оценки всех параметров и их статистические характеристики [1,2]. Задача расчета параметров в этом случае формулируется как задача определения статистических оценок параметров многооткликовой регрессионной модели по экспериментальным данным.

Многооткликовые модели для идентификации параметров многополюсников

Использование многооткликовых моделей для расчета параметров многополюсников требует решения нескольких проблем, связанных с конкретными свойствами исследуемых приборов. Эти проблемы включают анализ вида вольтамперной характеристики данного прибора и ее расчет, анализ ошибок наблюдений напряжений или токов, построение плана эксперимента для расчета параметров математической модели и выбор начальных приближений этих параметров.

Структурная математическая модель ВАХ многополюсника, получаемая на основании его эквивалентной схемы, представляется системой уравнений Кирхгофа в виде

F

где U = {u1,и2,...,ик}т — вектор независимых напряжений; I = {ix,i2,к,im}т — вектор наблюдаемых токов; B = {b1, ..., Ь}т — вектор параметров модели; F(B,U,I) — вектор неявных функций ВАХ, известных с точностью до параметров.

Система нелинейных уравнений (1) решается методом Ньютона [3] путем линейного разложения нелинейной функции F (B, X, Y) в виде итерационной процедуры

Is+1 = Is - F

„

где 0\В,и,I) =--------------------------матрица производных уравнении Кирхгофа по зависимым

д1

переменным.

Для получения максимально правдоподобных оценок параметров модели (1) использовался метод Гаусса [4] в виде:

В*+1 - в* + V* £р

]=1

*

где * — номер итерации; п — количество наблюдений; I — вектор расчетных значений тока в соответствии с выражением (2); Ув — ковариационная матрица оценок параметров; УЕ — ковариационная матрица ошибок наблюдений.

д1

Матрица производных — от всех зависимых переменных по параметрам эквива-

дв

лентной схемы определяется как производная неявной функции:

дв - Р(В, и, I ) = др (ви, 1) [<3(В, и, I )]-1. дВ дВ

*

Оценка ковариационной матрицы ошибок наблюдений в точке I рассчитывается по г повторяющимся наблюдениям как

1 г

-1 £ (I, -1*) -1').

г --  1

г-1

В качестве оценки ковариационной матрицы ошибок наблюдений УЕ в формуле (3) использовалась средневзвешенная ковариационная матрица по к наблюдениям в виде

^ 1 к

Е~ к^ Е к ] -1

Для проверки адекватности модели использовалась Т-статистика в виде

Т2 = —1— К],

тп -1"

] -1

где п — число наблюдений, а Я — вектор остатков модели.

При выполнении гипотезы о нулевых математических ожиданиях остатков модели эта статистика приблизительно имеет распределение Фишера со степенями свободы / - п -1 / т и /2 - а -1.

Для проверки однородности ковариационных матриц (4) рассчитывается статистика

Т22 - тахТ^Е^Е1}

которая при однородных ковариационных матрицах имеет приблизительно распределение Фишера с числом степеней свободы /1 - г -1 и /2 - к(г -1).

Реализация метода построения многооткликовых моделей и планирование последовательного эксперимента

Исследования, проведенные методом статистического моделирования, показали, что итерационная процедура (2) для исследуемого случая имеет высокую скорость сходимости и обеспечивает устойчивое и точное решение системы (1) в диапазоне реальных значений параметров В. Особенности реализации итерационной процедуры (3) для расчета 8Р1СЕ-параметров состоят в следующем.

На первых шагах итерационной процедуры оценки параметров могут изменяться в значительной степени, что может приводить к переполнению при расчете экспонент и к останову процедуры при вычислении логарифмов.

Параметры имеют ограничения, обусловленные их физическим смыслом. Если эти ограничения нарушаются, то процедура расчета останавливается.

В связи с этим вводятся ограничения на возможные значения параметров моделей в виде системы неравенств. Эти ограничения должны быть максимально широкими, чтобы они включали возможные значения параметров. Наличие таких широких ограничений не позволяет задавать начальные приближения оценок параметров достаточно точно, а в этом случае итерационный процесс для определения оценок часто либо расходится, либо заканчивается в локальном минимуме. Поэтому для расчета оценок параметров в этом случае предложена гибридная схема расчета, совмещающая метод статистического моделирования и метод Гаусса следующим образом.

1. На основе априорной информации задаются интервалы допустимых значений параметров.

2. Методом статистического моделирования генерируются начальные приближения параметров в заданных пределах в соответствии с равномерным распределением.

3. Полученные значения используются в качестве начальных приближений для расчета оценок коэффициентов методом Гаусса.

4. Для полученных оценок рассчитывается значение статистики (5), которое используется в качестве критерия получения реальных оценок. Для адекватной модели величина

Т2 (5) не должна превышать 1 + 3,36/4п (при доверительной вероятности 0,99).

5. Операции 2-4 повторяются, пока величина Т12 не окажется менее 1 + 3,36/>/и .

Оценивание значений отклика является основной целью построения математической модели. Ковариационная матрица ошибок оценивания тока в точке и 0 задается выражением

VI - РТ

Для обеспечения необходимой точности оценок параметров после получения оценок для заданного количества наблюдений выполняется процедура последовательного планирования эксперимента. В качестве критерия оптимальности плана использовалась величина определителя ковариационной матрицы оценок параметров [5]. Порядок измерений и расчетов для определения оценок параметров с заданной точностью с использованием процедуры последовательного планирования может быть представлен следующим образом.

1. При заданной дисперсии относительной ошибки оценивания тока в центре плана эксперимента и при заданном распределении параметров определяется требуемая величина определителя ковариационной матрицы оценок параметров.

2. Для заданной дисперсии ошибки наблюдения тока методом статистического моделирования определяется минимально необходимое количество наблюдений п, обеспечивающее устойчивость итерационного процесса расчета оценок параметров при равномерном плане эксперимента.

3. Задается начальный равномерный план эксперимента О* (* - 0) с количеством точек плана п.

4. В соответствии с этим планом выполняются измерения токов, рассчитываются оценки параметров модели, соответствующие ковариационные матрицы и критерий оптимальности плана.

5. К плану добавляется еще одна точка в середине диапазона изменения измеряемых величин.

6. Отыскивается (п + 1) -я точка эксперимента (л = л +1) ип+* такая, что добавление

этой точки в имеющийся план эксперимента Ол+1 обеспечивает максимальное повышение точности оценивания, т.е.

Б = тт йег[ув

В Е^и

7. В полученной точке плана выполняются измерения.

8. Критерием останова служит заданная предельная величина целевой функции — величина определителя ковариационной матрицы оценок параметров Бтах. Если

Б > Бтах, то в план добавляется следующая точка.

Результаты экстракции параметров 8Р1СЕ-модели биполярных транзисторов

Описанный метод определения параметров эквивалентной схемы использовался для экстракции параметров 8Р1СБ-модели биполярных транзисторов. Всего определялись 18 параметров. Методом статистических испытаний строилась зависимость максимально допустимой дисперсии относительной ошибки наблюдения тока от количества наблюдений для получения реальных оценок параметров (рис.1). Использовался равномерный план эксперимента.

100 200 500 1000 2000

Рис.1. Зависимость дисперсии относительной ошибки наблюдений от количества наблюдений п для различных моделей ВАХ транзистора: 1 — полная модель, включающая все 18 параметров; 2 — модель без параметра пк; 3 — модель без параметров пк и УЛР; 4 — модель без параметров nk, Уар и Уае

Рис.2. Зависимость определителя ковариационной матрицы О оценок коэффициентов модели ВАХ транзистора от количества наблюдений для равномерного плана

На рис.1 приведены зависимости для различных видов модели ВАХ: для полной модели и для моделей, предполагающих, что некоторые, наименее точно определяемые параметры, известны априорно. Эти зависимости используются для определения минимально необходимого количества наблюдений для начала процедуры последовательного планирования.

На рис.2. показаны зависимости величины Б (7) от количества наблюдений для равномерного и последовательного планов эксперимента.

На рис.3 показана зависимость дисперсии относительной погрешности оценивания тока ^ в центре плана от величины определителя ковариационной матрицы Б, построенные на основании выражения (6). При этом в качестве ковариационной матрицы в выражении (6) использовалась диагональная матрица, диагональные элементы которой принимались равными априорным дисперсиям оценок параметров модели ВАХ. Рис.3 позволяет по заданной максимально допустимой дисперсии прогнозирования тока определить требуемую величину определителя ковариационной матрицы оценок параметров, а по рис.2 можно оценить необходимое количество наблюдений.

Как видно, подобный подход снижает необходимое количество измерений в десятки раз по сравнению с равномерным планом. Для каждого прибора после расчета коэффициентов модели выполнялась проверка адекватности модели в соответствии с критерием (5). Для этого на основании расчет-

и _ гр 2

ной величины Tl определялась величина Р-уровня по формуле

^ = Ъ/

где FР, 1 / — квантиль распределения Фишера, соответствующий вероятности Р со степенями свободы / = п -1 / т и

/2 = к(а -1).

На рис.4 приведено распределение Р-уровня по партии приборов объемом 38 приборов КТ315А. Для получения нормального распределения строилось распределение величины 1п(Р).

Рис.3. Зависимость величины определителя ковариационной матрицы оценок параметров й в центре плана эксперимента от дисперсии относительной ошибки д, оценивания тока транзистора

Рис.4. Распределение величины !п(Р) по партии приборов

-15 -НО -0$ 0Я 0.5 ио г

Рис.5. Распределение величины г по партии приборов

Полученное распределение хорошо описывается нормальным законом, что подтверждается критерием Колмогорова — Смирнова. На основании нормального распределения рассчитывался нижний предел 95%-го доверительного интервала величины Р, который составил 0,906, что означает адекватность всех моделей на уровне значимости 0,05. Кроме проверки адекватности проводился анализ остатков модели и выполнялась проверка нормальности их распределения для всех параметров, строились диаграммы рассеяния и графики распределения остатков для визуального контроля нормальности остатков и выполнялась проверка однородности ковариационных матриц ошибок наблюдений.

Для сравнения результатов экстракции расчетные значения параметров хс сравнивались с экспериментальными значениями хт, которые можно было измерить непосредствен-

но. По этим данным рассчитывалась статистика г = -

где .* — стандартное отклоне-

ние. На рис.5 приведено распределение величины г, полученное по партии приборов для параметра 13.

хт ~ с

*

Заключение

Описываемый метод определения параметров эквивалентных схем с использованием многооткликовых регрессионных моделей позволяет рассчитать все параметры в едином процессе с заданной точностью. Для этого разработана двухступенчатая процедура расчета оценок параметров. На первом этапе этой процедуры формируется равномерный план эксперимента, а начальные приближения параметров задаются методом статистического моделирования при заданных ограничениях, что дает возможность получить реальные значения параметров. На втором этапе используется последовательное планирование эксперимента до достижения требуемой точности оценок параметров. Использование процедуры последовательного планирования позволяет повысить точность оценивания в десятки раз.

1. Попов С.А., Корчагин А.Ф. // Измерительная техника. № 4. 2003. С. 47-51.

2. Попов С.А., Корчагин А.Ф. // Сб. тр. VII Междунар. семинара «Современные проблемы прочности» им. В.А.Лихачева. В.Новгород, 2003. С.184-190

3. Г.Корн, Т.Корн. Справочник по математике для научных работников и инженеров. М.: Наука, 1968. 720 с.

4. Бард Й. Нелинейное оценивание параметров. М.: Статистика, 1979. 349 с.

5. Федоров В.В. Теория оптимального эксперимента. М.: Наука, 1971. 312 с.