- •Глава 4 оптимизация схем
- •4.9. Сравнение оптимизации в частотной области с оптимизацией методом подбора коэффициентов
- •4.10. Минимизация чувствительности [19]
- •4.11. Метод непрерывного преобразования схем с сохранением эквивалентности [18,31]
- •4.12. Заключение
- •Литература
- •Глава 5 моделирование активных элементов
- •5.1. Введение
- •5.2. Классификация моделей
- •5.3. Разработка моделей приборов по описывающим их дифференциальным уравнениям в частных производных
- •5.3.3. Разработка схемной модели
- •5.4. Моделирование полупроводникового диода
4.11. Метод непрерывного преобразования схем с сохранением эквивалентности [18,31]
В качестве заключения в данной главе рассмотрим один метод машинного проектирования схем, который, хотя и выглядит несколько академически по сравнению с мощными итерационными методами, изложенными ранее, все же может служить одним из примеров нетрадиционного подхода к решению классических задач с помощью современных вычислительных машин.
Пусть заданы: схемная функция T(s), которая должна быть синтезирована, и начальная реализация схемы N1. Требуется найти другую схему, эквивалентную N1, но с другими значениями параметров элементов или, может быть, имеющую структуру, отличающуюся от N1. Например, если схемная функция имеет вид
(4.150)
тогда любая схема N2 реализующая функцию
, (4,151)
эквивалентна N1 с точностью до коэффициента усиления.
Один способ получения класса схем, содержащего N2, состоит в том, что коэффициенты а1 берутся из множества решений
дифференциального уравнения
(4.152)
Так как
(4.153)
то очевидно, что b0= exp kx. Здесь без потери общности можно положить k=1.
Значение уравнения (4.152) состоит в том, что любое аi можно записать как функцию адмиттансов элементов yj. Раскрывая уравнение (4.152), имеем
(4.154)
Или можно написать матричное дифференциальное уравнение
При
m
= n,
если уравнения независимы, с помощью
обращения матрицы можно найти единственные
решений для
.
П р и м е р 4.12
Фильтр, изображенный на фиг. 4.11, имеет функцию передачи по напряжению (с учетом сопротивления источника R1)
Необходимо при G2=const и увеличении Ri сохранить функцию передачи неизменной.
Уравнения (4.152) записываются в данном случае в таком виде:
Или
в матричной записи
Переходя к обычной форме записи, получаем
(4.155)
Эти уравнения можно проинтегрировать методами численного анализа, взяв исходные значения параметров элементов схемы фильтра в качестве начальных.
Указанную методику вряд ли можно рекомендовать для больших схем, поскольку при численном интегрировании нужно брать малый шаг, чтобы поддерживать точность элементов в разумных границах. Поэтому почти всегда предпочитают использовать метод подбора коэффициентов. Тем не менее идея введения параметров сама по себе является полезной.
4.12. Заключение
Укажем некоторые направления современных исследований в области машинного проектирования радиоэлектронных схем:
1. Методы чебышевской аппроксимации частотных и временных характеристик [30, 32]. Предложенные до настоящего времени методы подобного типа являются либо узко специализированными, либо имеют сравнительно малые области сходимости или склонны приводить к локальному минимуму. Однако достоинства методов чебышевской аппроксимации вполне заслуживают дальнейших исследований в этом направлении.
2. Оптимизация (нелинейных) переключательных и импульсных схем [30].
3. Автоматическая трансформация и модификация топологии схемы: исключение, замена и введение новых элементов [13, 30, 31]. Автоматический синтез топологии схемы из схем холостого хода и схем короткого замыкания.
ЗАДАЧИ
4.1. Показать, что в методе наискорейшего спуска направление, определенное на i-й итерации, ортогонально направлению, определенному на
(i— 1) -й итерации, если âi-1 вычислено точно.
4.2. Флегчер и Поуэлл
показывают, что их метод обеспечивает
сходимость к А-1
в том случае, когда минимизируемая
функция квадратичная, доказывая,
что вектор
xj
является линейно независимым от
предыдущих (j
— 1) векторов. Почему из этого следует
сходимость к А-1?
4.3. Разложение в ряд Тейлора с остаточным членом для функции n переменных имеет следующий вид:
,
где
х
= x
—
,
— точка глобального минимума, £ лежит
на отрезке х +/Ах (О^/s^l).
Доказать, что если
является локальным минимумом, тогда
гессиан не может быть положительно
определенной матрицей на всем отрезке.
4.4. Плохую обусловленность в задачах, где функция F имеет эллипсоподобные линии постоянного уровня (как на фиг. 4.4а), можно устранить, реализуя итерации не по хi, а по log хi . Показать на примере, как это делается.
4.5. Дополнить математическое условие унимодальности так, чтобы оно годилось для n-мерного пространства. Изобразите чашеобразный участок поверхности, который удовлетворяет условию унимодальности, но не является выпуклым.
4.6. Можно показать, что в методе поиска с помощью чисел Фибоначчи
начальные
интервалы стремятся к пределу
Если выбрать в качестве начальных эти
интервалы, то какими выражениями
описываются длины интервалов при
последующих итерациях? Почему этот
метод менее эффективен, чем изложенный
в подразд. 4.4.3?
4.7. Показать, что в методе подбора коэффициентов должно выполняться следующее обязательное условие: число варьируемых элементов должно быть Не меньще числа полюсов, которые хотят сохранить на том же месте.
4,8, Высокочастотная
малосигнальная модель транзистора,
изображенная на схеме 34.1, содержит
элементы
у,
которые зависят от тока смещения
эмиттера. В частности,
,
(ток
Iе
в миллиамперах). Мы хотим компенсировать
изменения, вызываемые током эмиттера,
путем варьирования параметра одного
или нескольких элементов. Если требуется
сохранить два главных полюса, то
сколько элементов нужно сделать
варьируемыми, чтобы матрица частных
производных была квадратной? Найти эту
матрицу.
4.9. Пусть ограничение
на элемент
G
выражается неравенством
.
Подобрать
а1
и
в выражении
так, чтобы неравенство
для
G
выполнялось в случае
х,
изменяющейся в интервале [0,
.
Найти схему, реализующую
G,
в которой каждая проводимость
Фиг. 34.1. Полосовой усилитель.
либо
постоянна, либо имеет вид Gix.
Обратите внимание на то, что для
итеративная схема, приведенная в
подразд. 4.8.4.2, гарантирует, что
G
будет ограниченной величиной.
4.10 Минимизировать
функцию F(x)
при ограничениях
с помощью метода,
изложенного в подразд. 4.7.2, положив
где
rk<rk-1
. Пусть
обозначает
ту величину
,
при которой
P(x,rk)
достигает минимума. Показать, что
F(
k+1)
<F(
k)
и
.
(Не обязательно требовать, чтобы функции
F
и
Gj
были выпуклой и вогнутой соответственно.)
Если итерации начинаются вблизи от
границы Gj
=0, а минимум достигается далеко от
границы, то каким образом сумма
может убывать, как это требуется?
4.11.В этой главе рассмотрен
метод вычисления производных от
коэффициентов полинома
который в общем случае является
знаменателем некоторой схемной
функции. Попытайтесь развить аналогичный
метод вычислений производных от
коэффициентов полинома в числителе,
исходя из методики измерений
параметров n-полюсника,
аналогичной методике, отраженной на
фиг. 4.15 (Указание: используйте методы
вычисления частных производных в
частотной и временной областях,
изложенные в гл. 3.) Разработайте способ
рекурсивного вычисления производных
от коэффициентов при помощи методов
пространства состояний. (Указание:
примените метод разбиений или методы
теории n-полюсников,
рассмотренные в гл. 1 применительно к
анализу нелинейных цепей.) Следует
отдельно рассматривать случаи резистивных
и реактивных элементов.
4.12. Метод подбора
коэффициентов имеет тесную связь с
методом подбора корней [28]. Пусть
—
заданный (комплексный) корень; положим
Fi=ri - (i=1, 2, …, m), тогда fi = 0 является искомым решением. Пусть n-степень рассматриваемого полинома. Найти связь между двумя методами: а) при m =n; б) при т<n (случай главных и неглавных корней). Коэффициенты, естественно, подбираются с точностью до постоянного множителя. Используйте материал данной главы по чувствительности корней.
