Радіотехнічні кола і сигнали
УДК621.372.061
ЕФЕКТ ПРИСКОРЕННЯ ПРИ ОПТИМІЗАЦІЇ ЕЛЕКТРОННИХ КІЛ
Проектування аналогових електронних кіл з вибраною топологією зводиться до процесу параметричної оптимізації при одночасному аналізі кола на кожному кроці оптимізаційної процедури. Час проектування в цьому випадку швидко зростає при збільшенні розмірів і складності кола. Відомі традиційні підходи скорочення часу аналізу кола, пов’язані з використанням властивостей розріджених матриць і техніки декомпозиції. Нові ідеї, пов’язані з використанням методів теорії управління, були запропоновані авторами в ряді робіт [1-2]. Ідея відмови від точного дотримання законів Кірхгофа, яка запропонована в роботі [3], була узагальнена на основі застосування теорії оптимального управління. Це призводить, як показано в [1-2], до появи множини різних стратегій проектування і можливості побудування на основі цієї множини стратегії, що наближується до оптимальної за швидкодією. Цей підхід узагальнює процес проектування і дозволяє визначити задачу пошуку оптимального за часом алгоритму як задачу теорії управління. Основна концепція нового підходу полягає у введенні спеціальних керуючих функцій, які узагальнюють процес проектування і дозволяють керувати цим процесом для досягнення оптимуму цільової функції проектування за мінімально можливий машинний час. Ця можливість з’являється внаслідок практично необмеженого числа різних стратегій проектування, які існують в межах запропонованого підходу. В межах цієї концепції традиційна стратегія проектування є лише одним з можливих представників множини різних стратегій проектування.
У даній роботі аналізується новий ефект, що виникає в процесі проектування кіл методами теорії управління. Дослідження проводиться для одно- , дво- і трьохкаскадного транзисторних підсилювачів. Всі розглянуті приклади демонструють наявність ефекту прискорення процесу оптиміза-ції кола, який виникає внаслідок різниці у поведінці траєкторій проектування, що мають різні керуючі функції. Ми вважаємо тут і далі, що вектор параметрів математичної моделі кола X може бути розділений на дві частини: X=(Х, X"), де X 'є RK є вектором незалежних змінних і K - число незалежних змінних, X "є RM є вектором залежних змінних і M - число залежних змінних, при цьому N = K+M.
Процедура оптимізації і модель кола відповідно до запропонованої методології [1] задаються рівняннями:
Земляк А.М., Маркіна Т.М.
Постановка задачі
х.
(1)
12
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2010.-№40
Радіотехнічні кола та сигнали
І1 - uj )gj (x ) = 0, j = 1,2,..., M (2)
де s - номер ітерації, t s - ітераційний параметр, тобто крок оптимізаційної процедури і ts є R1, U є керуючий вектор, що складається з M компонент, U = (щ,U2,...,um), де Uj є Q; Q = {0;1}. Функції f (Х,U), що представляють собою компоненти вектора спрямованого руху, визначаються вибраним методом оптимізації. Ці функції, наприклад для градієнтного метода, визначаються формулами [1]:
f (XU) = -^TF(X•U), i = 1’2— K (3)
f (X,U) = -u-K -^-F(X,U) +
(1- Ui-K )r
dt
-xi +^i(X)
, i = K + 1, K + 2,..., N
де F (X,U) нова узагальнена цільова функція, що задається формулою F (X ) = C (X )+ф( X), де C (X) - цільова функція процесу оптимізації кола, а структура додаткової штрафної функції ф( X) може бути задана як
ф
/ \ 1 / ч
(XS) =1 £UJ ■ gj2 (XS)
(4)
Оператор
47
визначається формулою і^ф(X ) =
Эф^) K+M 5ф(X)cxf
8 xt т Г J' 8xi dxi
+ £
p=K+1 dxp dxi
Сенс керуючих змінних uj наступний: рівняння номер j присутнє в системі
(2), а член g j(X) вилучений з правої частини формули (4) у випадку, коли uj = 0, і навпаки, рівняння j вилучено з системи (2), а відповідний член
8
з’являється в правій частині формули (4), коли uj = 1. Всі керуючі змінні uj є функціями поточної точки процесу оптимізації.
Чисельні результати
Розглянемо процес проектування підсилювача на транзисторах (Рис. 1), що був проаналізований раніше в роботі [1], але без врахування цього можливого ефекту прискорення.
В якості моделі транзистора використовувалася статична модель Еберса-Молла, реалізована в системі SPICE [4]. Цільова функція визначена як сума квадратів різниць між заздалегідь заданими і поточними значеннями напруг на переходах транзисторів.
У роботі аналізується процес оптимі-зації кола для трьох різних випадків: одно-, дво- і трьохкаскадного підсилювачів. У першому випадку можна визначити три незалежних змінних, провідності y1,y2, y3, (К=3), и три залеж-
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2010.-№40
13
Радіотехнічні кола і сигнали
них змінних, вузлові напруги Vx,V2,V3, (М=3), в другому - є п’ять незалежних змінних y1, y 2, y3, y 4, У5, (К=5), и п’ять залежних змінних
Vj, V2, V3, V4, V5 (М=5), та в третьому - сім незалежних змінних
Уі,У2,Уз,У4,У5,Уб,У7, (К=7) и сім залежних змінних Vj,V2,V3,V4, V5,V6,V7, (М=7). Алгоритм оптимізації включає систему шести рівнянь виду (1) в першому випадку, десяти рівнянь - в другому випадку и чотирнадцяти - в третьому. Шаг процедури оптимізації ts вибирався постійним. Модель схеми (2) визначається трьома нелінійними рівняннями в першому випадку, п’ятьма рівняннями в другому, і сімома - в третьому. Базис стратегій проектування містить 2м різних стратегій, і в першому випадку - це 8 стратегій, в другому - 32 і в третьому - 128. Час оптимізації кожного з розглянутих підсилювачів на основі традиційної стратегії проектування (ТСП), для якої U = 0, V7=U..M, становив відповідно в першому випадку
13,79 с для градієнтного метода і 5,55 с для метода Давідона-Флетчера-Пауелла (ДФП), в другому - 257,37 с для градієнтного метода і 62,64 с для метода ДФП і в третьому - 14052,58 с для градієнтного метода і 1421,05 с для метода ДФП.
Для отримання ефекту прискорення початкова точка вектора змінних X в процедурі оптимізації вибиралася з умови негативної напруги на виході підсилювача.
Є безліч різних квазіоптімальних рішень при різному виборі стратегій з базису і різних початкових значеннях вектора X. Вибір найбільш перспективних стратегій проектування здійснювався на основі аналізу проекцій траєкторій на площині У3 - V3 для першого випадку, У5 - V5 для другого і У7 - V7 для третього. Згідно із загальною ідеєю побудови квазіопті-мальной за часом траєкторії [5], для першої частини траєкторії обиралася одна із стратегій, яка відноситься до групи стратегій, що мають траєкторії, подібні модифікованої традиційної стратегії проектування (МТСП), для якої Uj =1,V/=1,2,...M. Для здійснення перемикання вибиралися стратегії, що
мають траєкторії, подібні ТСП. В усіх трьох випадках квазіоптимальна стратегія включала два перемикання в керуючому векторі. Дані, відповідні квазіоптимальним стратегіям, при оптимізації всіх трьох підсилювачів наведені в таблиці для двох різних методів оптимізації - градієнтного метода и метода ДФП. При оптимізації однокаскадного підсилювача і використанні квазіоптимальної стратегії отримано виграш у часі в порівнянні з ТСП більше ста разів для обох методів оптимізації. Для двокаскадного підсилювача виграш становить майже тисяча для градієнтного методу і більше 600 для методу ДФП. Для трьохкаскадного підсилювача виграш становить більше 7000 для градієнтного методу і більше 4000 для методу ДФП.
14
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2010.-№40
Радіотехнічні кола та сигнали
Таблиця
Число каскадов Метод Оптимальный вектор управляющих функций U Число итераций Точки переключ. Время (сек) Выигрыш во времени
1 Градиентный (111); (000); (111) 1341 5 6 0,083 166
ДФП (111); (000); (111) 230 2 3 0,036 154
2 Градиентный (11111); (01000); (11111) 499 1 2 0,261 986
ДФП (11111); (00000); (11111) 222 1 2 0,098 639
3 Градиентный (1111111); (0000000); (1111111) 5090 10; 11 1,82 7721
ДФП (1111111); (0000000); (1111111) 462 2; 3 0,33 4306
Результати оптимізації цих активних кіл представлені на рис. 2 у вигляді залежності виграшу в часі G від числа каскадів підсилювача.
Видно, що отриманий виграш значно зростає при збільшенні розміру та складності кола.
Висновки
Проведені дослідження показують, що структура квазіоптімальной траєкторії проектування повинна базуватися на додатковому ефекті прискорення процесу. Однією з умов появи цього ефекту є правильний вибір початкового наближення, тобто наявність у початковому значенні вектора змінних X принаймні однієї негативної компоненти. Другою умовою є правильний вибір послідовності стратегій проектування в процесі побудови оптимальної траєкторії. Тут мається на увазі вибір на початковому етапі однієї зі стратегій, що має траєкторію подібну МТСП, і потім зміну її на стратегію, що має траєкторію подібну ТСП. Додаткові зміни керуючого вектора дозволяють ще більш наблизитися до оптимальної за часом траєкторії. Третя і найістотніша умова - це визначення точного моменту перемикання вектора керуючих функцій.
Сформульовані умови можна визначити як достатні для отримання додаткового ефекту прискорення процесу проектування. Результати, отримані в роботі, є важливим кроком до виявлення структури найкращої стратегії оптимізації кола і служать основою для пошуку оптимального за часом алгоритму проектування.
Література
1. Zemliak A.M. Analog System Design Problem Formulation by Optimum Control Theory, IEICE Transactions on Fundamentals of Electronics, Communication and Computer Sciences, 2001, Vol. E84-A, № 8, pp. 2029-2041.
2. Земляк А.М., Маркина Т.М., Оптимизация электронных цепей как управляемый динамический процесс, Автоматизированные системы управления и приборы автоматики. Вып. 146, 2009, стр. 62-69.
3. Каширский И.С., Трохименко Я.К. Обобщенная оптимизация электронных схем. - Киев: Техника. - 1979.- 192 с.
4. Massobrio G., Antognetti P. Semiconductor Device Modeling with SPICE.- N.Y.: Mc. Graw-Hill, Inc.- 1993.
5. Zemliak A. Comparative Analysis of the Different Design Strategies for the Analog Electronic System Design//WSEAS Trans. on CS, V. 3, No. 2, 2004, pp. 294-299.
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2010.-№40
15
Радіотехнічні кола і сигнали
Земляк О.М., Маркіна Т.М. Ефект прискорення при оптимізації електронних кіл.
Проаналізовано ефект додаткового прискорення процесу оптимізації електронного кола. Цей ефект виникає внаслідок різної поведінки траєкторій проектування при зміні вектора керуючих функцій. Ефект проявляється для всіх кіл, що аналізуються, и дозволяє додатково скоротить витрати машинного часу на кілька порядків.
Ключові слова: оптимізація електронних кіл, електронні кола.____________________
Земляк А.М., Маркина Т.М. Эффект ускорения при оптимизации электронных цепей. Проанализирован эффект дополнительного ускорения процесса оптимизации электронной цепи. Этот эффект возникает вследствие различного поведения траекторий проектирования при изменении вектора управляющих функций. Эффект проявляется для всех анализируемых цепей и позволяет дополнительно сократить затраты машинного времени на несколько порядков.
Ключевые слова: оптимизация электронных цепей, электронные цепи_________________
Zemliak A.М., Markina T.M. Effect of acceleration in optimization of the electronics circuits. An additional acceleration effect has been analyzed for the electronic circuit optimization process. This effect has been appeared due to the different behavior of the design trajectories whenever the vector of the control functions was changed. Effect was obtainedfor all analyzed circuits and it can reduce the computer time to some powers.
Key words: optimization electronics circuits, electronics circuits, design , electronics circuits
УДК 621.372.061:391.266
ОЦІНКА НАДІЙНОСТІ АУТЕНТИФІКАЦІЇ ОСОБИ ЗА ДИНАМІЧНО ВВЕДЕНИМ ПІДПИСОМ ПРИ НАЯВНОСТІ
АДИТИВНОГО ШУМУ
Рибін О.І., Кузьменко О.М.
Задача аутентифікації особи за динамічно введеним підписом є частковим випадком задачі класифікації сигналів. Таку задачу розв’язують за допомогою різноманітних методів, до яких належить класична лінійна (погоджена) фільтрація [1], нелінійні косинусна [2] та нормальна [3-7] фільтрація, а також погоджена фільтрація на базі перетворення Карунена - Лоє-ва [8-10]. В роботах [9,10] було запропоновано модифікований алгоритм з використанням перетворення Карунена - Лоєва, який показав надзвичайну надійність [10] саме при розв’язанні задачі аутентифікації за динамічно введеним підписом. Достатня для оптимістичних висновків, але недостатня для статистичних досліджень вибірка підписів тієї самої особи (для багатьох осіб) не дозволяє встановити границі надійної аутентифікації при зростанні випадкових відхилень реалізацій підписів від математичного очікування (що, наприклад, має місце при варіації підпису через тривалий період часу для тієї самої особи).
Теоретичні положення
У відповідності до алгоритму[10] аутентифікації особи для кожного класу підписів формується кореляційна матриця
Cor
1 M --- --т
—V Ах-Ах M “1
(1)
де Ах = х - m - стовпець відхилень значення кожного відліку дискретизо-ваного підпису з певного класу(тобто певної особи) від математичного очі-
16
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2010.-№40