УДК 513.73
РАСЧЕТ ПАРАМЕТРОВ ЭКВИВАЛЕНТНОЙ СХЕМЫ IGB ТРАНЗИСТОРА
С.А.Попов, И.С.Васильев ESTIMATION OF EQUIVALENT-CIRCUIT PARAMETERS OF IGB TRANSISTOR
S.A.Popov, I.S.Vasil'ev
Политехнический институт НовГУ, [email protected]
Предлагается метод оценивания параметров эквивалентной схемы IGB транзистора. Рассмотрен метод расчета статических параметров эквивалентной схемы с использованием многооткликовой модели. Приводится процедура последовательного планирования эксперимента, обеспечивающая заданную точность оценивания.
Ключевые слова: IGB транзистор, оценивание параметров эквивалентной схемы, многооткликовая модель, последовательное планирование эксперимента
The estimation technique for equivalent-circuit parameters of IGB transistor is offered. The calculation method for static parameters of equivalent circuit using multiresponse model is presented. The procedure of consistent design of experiments ensuring required accuracy of estimation is provided.
Keywords: IGB transistor, estimation of equivalent circuit parameters, multiresponse model, consistent design of experiments
IGBT (Insulated-Gate Bipolar Transistor) — это мощный транзистор с изолированным затвором, специально разработанный для импульсного режима работы. Такие транзисторы применяются в мощной электронике, особенно в импульсных источниках питания, серводвигателях и трехфазных моторах, поскольку они обеспечивают широкий динамиче-
ский диапазон, низкий уровень акустического шума и повышенную эффективность. Однако ЮВ транзисторы имеют довольно сложную эквивалентную схему, определение параметров которой представляет определенные трудности. В статье [1] описан метод оценивания параметров эквивалентных схем многополюсников по результатам измерений вольт-
амперных характеристик (ВАХ). Этот метод позволяет оценивать статические параметры эквивалентных схем и их статистические характеристики на основе построения многооткликовых моделей. Расчет параметров эквивалентной схемы ЮВ транзистора и их статистических характеристик осложняется тем, что в этом случае нет возможности измерить все требуемые для расчета электрические режимы.
Статические режимы и параметры эквивалентной схемы
Эквивалентная схема представлена на рис. 1.
ЮВ транзистора [2]
Рис.1. Эквивалентная схема ЮВ транзистора
Данная эквивалентная схема позволяет получить систему уравнений Кирхгофа в следующем виде:
А -KpVl(l1-0,5Vl)
-1 R ^
"с
ехР
V2-V3 1-1
ф.
т
+а р1р
ехР
V -V
-1
V -у
R
ч
- -а1р
ехр
Фт
ехр
V -V
4 'з |-1
V4
-А_3+а I
ехР
V2_-l |-1
-1 -1
VCC ^
ехр
и 1-1 фт
=0(1)
Здесь а=(1-ар), Vl,V2,Vъ,V4 — напряжения в соответствующих узлах, причем напряжения V и V4 задаются: Vl= Е1 и V = Е2, а и 1усс — токи соответствующих источников напряжений.
Статический режим IGB транзистора зависит от параметров эквивалентной схемы ар,а„, 1с„,К
Rс,R1,R2,1Е8 (всего 8 параметров) и задаваемых источников напряжений: ЕрЕ2. Для расчета этих параметров и их статистических характеристик используется модель ЮВ транзистора, получаемая на основании системы уравнений (1) в виде:
F(X,Y,B)=0, (2)
где F(X,Y,B)={f1(X,Y,B),f2(X,Y,B),...,fя(Х,У,Б)}Т — вектор неявных функций, определяемых правой частью системы уравнений (1), X={хр х2,..., х^}Т — вектор задаваемых (независимых) электрических режи-
мов,
Y= {у
1У 2
Т
ут) — вектор зависимых режимов,
Б= {Ь1,Ь2,...,Ь1}Т — вектор параметров элементов эквивалентной схемы.
Задача расчета параметров с использованием модели (2) формулируется как задача определения статистических оценок параметров многооткликовой регрессионной модели по экспериментальным данным. Параметры эквивалентной схемы изменяются в технологическом процессе с помощью тестовых структур. Этот метод не позволяет определять параметры отдельных транзисторов и, соответственно, статистические характеристики этих параметров по партии приборов. Предлагаемый в данной работе метод расчета параметров эквивалентной схемы основан на результатах измерений доступных зависимых электрических режимов при заданных независимых источниках напряжения. Задача расчета параметров в этом случае формулируется как задача определения статистических оценок параметров многооткликовой модели по экспериментальным данным. Для оценивания параметров эквивалентной схемы проводится эксперимент по измерению зависимых электрических режимов при различных значениях задаваемых режимов. По результатам этого эксперимента путем построения многооткликовой регрессионной модели можно получить максимально правдоподобные оценки параметров эквивалентной схемы и их ковариационную матрицу с помощью следующей процедуры [3]:
Б = Б;
(3)
1=1
где — расчетные значения зависимых переменных, УБ — ковариационная матрица оценок параметров, УЕ — ковариационная матрица ошибок измерений, : — номер итерации, п — количество экспериментов.
Расчетные значения зависимых режимов определяются на основании системы уравнений (1) в виде [4]:
= -|д(в, )]~^(в, ), (4)
c'F(B X Y)
где Q(B,X,Y)=— ду' — матрица производных
системы уравнений (1) по зависимым электрическим режимам.
Матрица производных G(B,X,Y) от всех зависимых режимов по параметрам эквивалентной схемы
т
т
определяется как производная неявной функции следующим образом:
G(B, ВДГ1. (5)
Ковариационная матрица оценок коэффициентов (3) равна [2]:
Vb=|£G(B,X.,Y;И .
(6)
Оценка ковариационной матрицы ошибок наблюдений УЕ рассчитывается в соответствии с [3] по
остаткам модели следующим образом:
V =
m
E mn -1
■]Г [y - Y (B, X)]Y - Y (B, X)[,
(7)
i=i
где Y(B, X) — оценка зависимых переменных. получаемая по формуле (4).
В случае ЮВ транзистора могут быть заданы два независимых источника напряжения х1 = Е1 и
х2 = Е2 (независимые переменные) и могут наблюдаться два тока у1 = 1М05 и у2 = 1уСЕ (зависимые переменные); таким образом. k = 2. т = 2. I = 8.
Планирование и реализация эксперимента
При использовании равномерного расположения точек контроля зависимых режимов величина определителя ковариационной матрицы оценок параметров оказывается очень большой. что не позволяет оценивать параметры с заданной точностью. Критерий ^-оптимальности требует такого расположения точек в области планирования Ех. при котором определитель ковариационной матрицы оценок коэффициентов detУБ минимален [5]. D-
оптимальный план определяется следующим образом:
det УБ (п° )= УБ (П)=
mindetj^ р(в, X. К1[р(В , X. f|
(8)
Л=1
где П° = {х°.х2.....Х°} — оптимальный план в смысле критерия D-оптимальности.
Для повышения точности оценивания параметров использовалась процедура последовательного планирования [5]. В этом случае задается некоторый заранее определенный план эксперимента с небольшим количеством точек контроля. В качестве такого плана можно использовать. например. равномерный план эксперимента. При использовании этого плана выполняются измерения и рассчитываются оценки коэффициентов модели и определитель ковариационной матрицы этих оценок detУБ. На основании полученных оценок отыскивается следующая точка контроля. обеспечивающая наилучшую точность оценивания. Количество измерений увели-
чивается до достижения необходимой фактической точности оценивания параметров данного конкретного прибора. Алгоритм построения последовательного плана может быть представлен следующим выражением:
det VB (п+Xo,B)= mindet VB(n+X,B), (9)
XeiX
где Xo — добавляемая в план П точка.
Для получения экспериментальных данных использовался программно-аппаратный комплекс, позволяющий проводить измерения в автоматизированном режиме [6]. Комплекс состоит из ЭВМ, на которой установлена управляющая программа, управляемых нагрузок, управляемых источников питания, блока преобразований, блока коммутации, измерительной аппаратуры. В комплексе используется универсальная промышленная измерительная аппаратура с интерфейсом для подключения к ЭВМ. Блок коммутации управляется специализированным контроллером, в свою очередь получающим команды от управляющей программы посредством Ethernet LAN-интерфейса.
Программное обеспечение, установленное на
ЭВМ:
1) задает режимы источников питания E1 и E2 в соответствии с планом эксперимента;
2) проводит измерение токов I S и I
для
МО5 УСЕ
заданных режимов;
3) рассчитывает величины параметров эквивалентной схемы (3);
4) рассчитывает величину определителя ковариационной матрицы det УБ (6);
5) если величина det УБ больше заданной. то
рассчитывает следующую точку последовательного плана эксперимента (9).
В результате обеспечивается заданная точность определения параметров эквивалентной схемы.
Результаты эксперимента
Описанный метод был реализован для оценивания параметров эквивалентной схемы ЮВ транзисторов типа IRGPC40F. Были проведены измерения ВАХ при изменении напряжения на затворе в диапазоне 0.1^15 V и на коллекторе в диапазоне 1^600 V. Для каждого прибора выполнялось 100 и 250 измерений тока соответственно относительной погрешности при равномерном распределении точек контроля для получения значимых величин оцениваемых параметров.
На рис.2. показаны зависимости величины логарифма определителя ковариационной матрицы оценок параметров эквивалентной схемы detУБ (6)
от количества наблюдений для равномерного и последовательного планов эксперимента.
Как видно на рис. 2. использование последовательного планирования повышает точность оцени-
log(detVB) -68
-72 -76 -80
-84
\\ \ \ \ 41
\ \ ч ч \ 4
\ \ ----- 4 4 2 4 4 4 4 4 4 N
s N N V N 4 V V к ._______
V. N, •s, •ч «4 ««
100
200
500
1000
2000
Рис.2. Зависимость логарифма определителя ковариационной матрицы detVB оценок параметров модели ВАХ транзистора от
количества наблюдений для равномерного плана (сплошная линия) и для последовательного плана (пунктирная линия) для относительной ошибки наблюдения равной 0,05% (1) и 0,01% (2)
вания в десятки раз по сравнению с равномерным планом при одинаковой погрешности оценивания. Кроме того, как показывают расчеты, некоторые параметры эквивалентной схемы вообще невозможно оценивать с достаточной точностью без использования последовательных планов эксперимента.
Попов С.А., Корчагин А.Ф. Оценивание параметров эквивалентной схемы многополюсников с помощью мно-гооткликовых моделей // Вестник НовГУ. Сер.: Естеств. и техн. науки. 2001. № 19. С.175-179.
Попов С.А., Корчагин А.Ф. Использование многооткли-ковых моделей для расчета параметров электронных приборов // Измерительная техника. 2003. № 4. С.47-51. Бард Й. Нелинейное оценивание параметров. М.: Статистика, 1979. 349 с.
Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1968. 720 с. Федоров В.В. Теория оптимального эксперимента. М.: Наука, 1971. 312 с.
Васильев И.С. Автоматизация тестирования радиоэлектронных компонентов, выпускаемых малыми сериями // Вестник Новгородского государственного университета. Сер.: Техн. науки. 2013. № 5. Т.1. С.98-100.
References
Popov S.A., Korchagin A.F. Otsenivanie parametrov ekviva-lentnoi skhemy mnogopoliusnikov s pomoshch'iu mnogootklikovykh modelei [Estimation of equivalent-circuit parameters of multipoles using multiresponse models]. Vest-nik Novgorodskogo gosudarstvennogo universiteta. Seriia: Estestvennye i tekhnicheskie nauki - Vestnik of Yaroslav the Wise Novgorod State University. Issue: Natural Sciences and Engineering, 2004, no. 19, pp. 175 -179. Popov S.A., Korchagin A.F. Ispol'zovanie mnogootklikovykh modelei dlia rascheta parametrov elektronnykh priborov [Using multiresponse models for estimation of electronic devices' parameters]. Izmeritel'naia tekhnika - Measurement Technique, 2003, no. 4, pp. 47-51.
Bard I. Nelineinoe otsenivanie parametrov [Nonlinear estimation of parameters]. Moscow, "Statistika" Publ., 1979. 349 p. Korn G., Korn T. Spravochnik po matematike dlia nauchnykh rabotnikov i inzhenerov [Mathematics reference book for researchers and engineers]. Moscow, "Nauka" Publ., 1968. 720 p. Fedorov V.V. Teoriia optimal'nogo eksperimenta [Optimal experiment theory]. Moscow, "Nauka" Publ., 1971. 312 p. Vasil'ev I.S. Avtomatizatsiia testirovaniia radioelektronnykh komponentov, vypuskaemykh malymi seriiami [Automating testing of radioelectronic components produced by limited editions]. Vestnik Novgorodskogo gosudarstvennogo universiteta. Seriia: Tekhnicheskie nauki - Vestnik of Yaroslav the Wise Novgorod State University. Issue: Engineering Sciences, 2013, no. 75, vol. 1, pp. 98 -100.
1.
2.
1.
3
2
4
3
5.
4
6
5
6.