CC BY
62
УДК 81.255.2
ОЦЕНИВАНИЕ SPICE-ПАРАМЕТРОВ С ЗАДАННОЙ ТОЧНОСТЬЮ
С.А.Попов, И.С.Васильев
ESTIMATION OF SPICE-PARAMETERS WITH REQUIRED ACCURACY
S.A.Popov, LS.Vasilyiev
Политехнический институт НовГУ, stanislav.popov@novsu.ru
Рассматривается метод оценивания Spice-параметров и их ковариационной матрицы на основании многооткликовой модели эквивалентной схемы. Полученные параметры используются для расчета и коррекции электрического режима электронных устройств. Предлагается метод построения последовательного плана эксперимента, позволяющего получить необходимую точность оценивания.
Ключевые слова: оценивание Spice-параметров, многооткликовая модель, точность оценивания, план эксперимента
Spice-parameters and their covariance matrix are estimated on the basis of multiresponse model of equivalent circuit. The obtained parameters are used for calculation and correction of electrical mode of electron devices. The method is offered for construction of sequential design of experiment that allows to provide the required accuracy of estimation. Keywords: Spice-parameters estimation, multiresponse model, accuracy of estimation, design of experiment
Для прогнозирования выходных электрических режимов электронных устройств, состоящих из транзисторов и других компонентов, величины параметров которых можно точно измерить, необходимо знание SPICE-параметров транзисторов. Задача оценивания SPICE-параметров транзисторов и их статистических характеристик с учетом мно-гооткликового характера моделей этих приборов формулируется как задача определения статистических оценок параметров многооткликовой модели по экспериментальным данным. Такие задачи включают расчет вольт-амперных характеристик (ВАХ), анализ ошибок наблюдений и построение плана эксперимента для расчета параметров математической модели [1,2]. Затем эти параметры используются для расчета электрических режимов электронных устройств и сравнения их с допустимыми значениями.
Модель ВАХ транзистора с однородной ошибкой и расчет SPICE-параметров
Многооткликовая модель ВАХ транзистора с абсолютной ошибкой представляется в виде [1]
I = F (В, и, I)+Я, (1)
где F(B,U, I) — система уравнений Кирхгофа, описывающая эквивалентную схему транзистора, I={il,...,ikf — вектор токов, U={up...,um}T — вектор напряжений, B={pl,...,b} — вектор коэффициентов модели, которые и являются определяемыми SPICE параметрами, R= {rp...,rm}T — абсолютная
ошибка наблюдения.
При измерении тока в очень широких пределах абсолютная ошибка наблюдения также изменяется. На практике принимается, что абсолютная ошибка измерения пропорциональна величине измеряемой переменной. Таким образом, зависимая от величины тока абсолютная ошибка наблюдения тока R выражается через постоянную относительную ошибку наблюдения Е в виде R=IET. Тогда
ln( I+R)=ln(l+IET )=ln I+ln(J+E),
где J={1,...,1}T .
Разлагая функцию ln( J+E) в ряд Тейлора возле точки Е = 0 и ограничиваясь одним членом для малых Е, получим
ln(J+E)« E =Д. « E.
dE
1+E
Таким образом, для приведения ошибки наблюдений к однородной можно принять, что
1п( I+Я)=1п(/+1Ет )= 1п I+1п(/+Е)«1п I+Е.
Тогда модель ВАХ с однородной ошибкой наблюдений выражается в виде
1п I = 1п[^ (В,и, 1)]+Е, (2)
где Е={в1,.,вт}т — однородная ошибка наблюдений с ковариационной матрицей VE .
Для расчета оценок параметров используется итерационная процедура, основанная на модифицированном методе Гаусса, который является достаточно быстродействующим и относительно устойчивым в вычислительном отношении. Оценки параметров, обеспечивающие максимум функции правдоподобия, в этом методе можно выразить в следующем виде [3]:
В8+1 = В8 +
Iр(В 8 ,и., I, НВ 8 ,и., I, Г
1_ 1=1
IФ8,и1,I1 ^I/ -1пЪI
(3)
1=1
где 8 — номер итерации, п — количество наблюдений ВАХ, I* — наблюденное значение тока, II —
расчетное значение тока по модели при заданном напряжении в соответствии с выражением (2).
Векторная функция Р(В,и, I) определяется как производная неявной функции (2) следующим образом:
Р(В,и, I)=-1 ^ (ВВи, ) [и(В,и, I)]1, (4)
где L — единичная матрица, а и (В,и, I) = 'дF (В,и, I)
дI
—L
Ковариационная матрица VB оценок параметров рассчитывается по формуле
УВ НТР(В,и,,I/ УЕР(В,и/,I/Т
(5)
Расчет модели ВАХ и предсказание по модели
Для реализации процедуры расчета (3) коэффициентов модели (2) на каждой итерации необходимо рассчитывать значение тока I. Расчет выполняется путем решения системы уравнений ВАХ в виде
F (В,и, I)-1=Ф(В,и,!)=0. (6)
Эта система нелинейных уравнений решается методом Ньютона [4] линейным разложением нелинейной функции Ф(В,и, I) в окрестности точки I в виде
-г
Ф(и, В, I)« ф(В,и, I*)
дФ(В,и, I)
_ дI I=т8 -
(I -1 *).
Тогда итерационная процедура расчета вектора токов имеет вид
18+1 = 18-ф(В,и, 18 рад!)
- ^
I=18
(7)
Исследования, проведенные методом статистического моделирования, показали, что итерационная процедура (7) имеет высокую скорость сходимости и обеспечивает устойчивое и точное решение системы (6) в диапазоне реальных значений параметров В внутри цикла (3).
Точность оценивания и планирование эксперимента
Основной целью построения математической модели транзистора является предсказание значений параметров-критериев годности транзисторных устройств (например, транзисторных ключей). Оценка
токов I при заданных напряжениях определяется решением системы уравнений в виде
н(B,Z,U,I)-I=0, (8)
где Н(В,Z,и,I) — система уравнений Кирхгофа,
описывающая электронное устройство, В — оценки SPICE-параметров транзисторов, получаемые по формуле (3), Z — известные компоненты электрической схемы, параметры которых можно измерить.
Система уравнений (8) решается с помощью итерационной процедуры, аналогичной (7). Если признать, что модель транзистора сама по себе правильна, то существуют два возможных источника ошибок предсказания выходных параметров устройства: ошибки в оценках коэффициентов и ошибки наблюдения откликов. Тогда ковариационная матрица ошибок предсказания выходных параметров рассчитывается следующим образом:
vY=Q(B, ад гУв (Q(B, Z,и, I)j+VE, (9)
где Q(B, Z ,и, /) — матрица производных, получаемая
по аналогии с (4).
Доверительные интервалы оцениваемых по модели (8) токов i. выражаются в виде
г .4 - V. <г <г .+t ..IV..
1 ау ] ] ] аУ ] -
(10)
где г. — оценка .-го тока по формуле (7), V. — дисперсия оценки тока (/-й диагональный элемент ковариационной матрицы У) tа — квантиль распределения Стьюдента на уровне значимости а и с числом степеней свободы /=п -1 /k, /=1Д
В качестве критерия оптимальности плана в этом случае принимается определитель ковариационной матрицы det VY [5]. Для расчета этого критерия
оптимальности необходимо для заданного плана эксперимента рассчитать ковариационную матрицу оценок коэффициентов УВ, а затем по формуле (9) рассчитать ковариационную матрицу оценок откликов VY и ее определитель. Оптимальные планы, полученные для этого критерия оптимальности, называются G-оптимальными. Известно [5], что при непрерывном планировании для построения G-оптимальных планов вместо критерия оптимальности в виде det VY можно
использовать критерий оптимальности в виде det УВ
-1
X
(^-оптимальные планы). Это в значительной степени ускоряет процедуру планирования, поскольку упрощает вычисления критерия оптимальности.
В случае последовательного планирования задается некоторый начальный заранее заданный план эксперимента с небольшим количеством точек наблюдений. Используя этот план, выполняются измерения и рассчитываются оценки коэффициентов модели. На основании полученных оценок отыскивается следующая точка наблюдения, обеспечивающая наилучшую точность оценивания. Количество измерений увеличивается до достижения необходимой точности оценивания параметров данного конкретного прибора. Алгоритм построения последовательного плана может быть представлен следующим выражением:
detVB (П+X0, в)= (П+X, В), (11)
Хе1х
где П — построенный план эксперимента, X0 — следующая добавляемая в план П точка.
При разработке электронных схем, например транзисторных ключей, в качестве параметров-критериев годности задаются точки на ВАХ и соответствующие допустимые значения, например, для токов /'. в виде
i ^ i ,
mj
(12)
нижнее
где гт. — измеренное значение тока, г.
допустимое значение тока годных изделий.
Чтобы сделать обоснованный (на уровне значимости а) вывод о годности прибора, необходимо, чтобы нижняя граница доверительного интервала (10) была больше допустимого значения г'., т.е. требуется
выполнение следующего неравенства:
-и^ч-. (13)
Если выполняется неравенство
i +1 Jv. < i ,
J J lJ'
(14)
то прибор не соответствует ТУ по данному параметру. Если оба неравенства (13) и (14) не выполняются, то эксперимент нужно продолжать.
При выполнении неравенства (14) требуется коррекция технологического процесса с целью подгонки SPICE-параметров для получения требуемых значений параметров-критериев годности. Выражения (8) и (9) позволяют получать точечную и интервальную оценки параметров-критериев годности электронного устройства при заданных SPICE-параметрах транзистора, что дает возможность корректировать технологических процесс для обеспечения выполнения ТУ.
1. Попов С.А., Корчагин А.Ф. Использование многоот-кликовых моделей для расчета параметров электронных приборов // Измерительная техника. 2003. №4. С.47-51.
2. Попов С.А., Васильев И.С. Расчет параметров эквивалентной схемы IGB транзистора // Вестник НовГУ. Сер.: Физ.-мат. науки. 2014. №80. С.33-36.
3. Бард Й. Нелинейное оценивание параметров. М.: Статистика, 1979. 349 с.
4. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1968. 720 с.
5. Федоров В.В. Теория оптимального эксперимента. М.: Наука, 1971. 312 с.
References
1. Popov S.A., Korchagin A.F. Ispol'zovanie mnogootklikovykh modelei dlia rascheta parametrov elektronnykh priborov [Using multiresponse models for calculating the parameters of electron devices]. Izmeritel'naia tekhnika - Measurement Techniques, 2003, no. 4, pp. 47-51.
2. Popov S.A., Vasil'ev I.S. Raschet parametrov ekvivalentnoi skhemy IGB tranzistora [Estimation of equivalent-circuit parameters of IGB transistor]. Vestnik NovGU. Ser. Fiziko-matematicheskie nauki - Vestnik NovSU. Issue: Physico-Mathematical Sciences, 2014, no. 80, pp. 33-36.
3. Bard Y. Nonlinear Parameter Estimation. Academic Press, New York and London, 1974. (Russ. ed.: Bard I. Nelineinoe otsenivanie parametrov. Moscow, "Statistika" Publ., 1979. 349 p.).
4. Korn G.A., Korn T.M. Mathematical Handbook for Scientists and Engineers. McGraw-Hill, New York, 1968. (Russ. ed.: G. Korn, T. Korn. Spravochnik po matematike dlia nauchnykh rabotnikov i inzhenerov. Moscow, "Nauka" Publ., 1968. 720 p.).
5. Fedorov V.V. Teoriia optimal'nogo eksperimenta [Optimal experiment theory]. Moscow, "Nauka" Publ., 1971. 312 p.
