задача скалярной оптимизации будет сформулирована в виде задачи минимаксной аппроксимации, а именно, подставляя (2.4) в (2,5) имеем
min max e( X , ). |
(2.6) |
XE
Вобщем случае аналитическое решение задачи (2.3) неизвестно. Для численных
методов расчета характеристик устройства необходимо перейти к дискретной задаче путем
покрытия области E |
-сетью с дискретными значениями |
1, |
2,…, |
Р. Если обозначить |
f j ( X ) Hi ( X , l ) , Z1 j |
Hнi ( l ) , Z2 j Hвi ( l ) , [1, Pi ], |
i |
[1,P], |
где Pi - количество |
точек в поддиапазоне E i , то соответственно условия работоспособности будут определяться соотношениями
Z1 j f j ( X ) |
Z2 j , j [1, M ], |
(2.7) |
и |
|
|
Qj (X ) |
e j (X , ) , |
(2.8) |
а дискретная минимаксная задача, следующая из исходной непрерывной задачи (2.3), может быть сформулирована в виде задачи скалярной оптимизации Ошибка! Источник ссылки не найден., где
|
|
|
Q( X ) |
max Qj ( X ) |
(2.9) |
||
|
|
|
|
|
j [1,M ] |
|
|
или в виде |
|
|
|
|
|
|
|
|
|
|
min max Qj ( X ) |
(2.10) |
|||
|
|
|
|
X Ew |
j [1,M ] |
|
|
Таким образом, задача векторной оптимизации фактически сведена к задаче |
|||||||
скалярной оптимизации, решение которой связано с минимизацией M функции максимума |
|||||||
(2.9), заданной на EN. |
|
|
|
|
|
|
|
В случае, когда критерий Q(X ) задан на всем пространстве EN, имеем следующие |
|||||||
выражения для непрерывных и дискретных задач скалярной оптимизации |
|
||||||
|
|
|
|
min Q( X ), |
(2.11) |
||
|
|
|
|
X |
D |
|
|
|
|
|
min max e( X , ), |
(2.12) |
|||
|
|
|
X |
D |
E |
|
|
|
|
|
|
min max Q j ( X ). |
(2.13) |
||
|
|
|
|
X D |
j [1,M ] |
|
|
Нередко в качестве обобщенного критерия выбирается среднестепенная норма |
|||||||
|
|
|
|
|
M |
|
|
|
|
|
Q( X ) |
|
Q ( X ) q |
(2.14) |
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
i |
1 |
|
где Q j ( X ) e j ( X ) |
|
, j |
[1, M ]; |
H 0 j - заданная "эталонная" функция; |
|
||
j |
H 0 j f j |
|
|||||
q 1 - целое положительное число.
94
Когда требования ЧТЗ заданы в виде коридоров допусков, то в [2] предложено вместо обобщенного критерия (2.14)использовать критерий вида
|
|
|
|
Qj ( X ) s |
1s |
|
|||
Q( X ) |
|
Qmax |
|
|
|
, |
(2.15) |
||
j J Qmax |
|||||||||
|
|
|
|
|
|
||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
s qsign(Qmax ) , |
|
|
||||
|
|
|
|
|
|||||
J |
{ j |
Q j ( X |
0)}, при Qmax |
0 |
|
||||
1 : M , при Qmax 0 |
|
|
|||||||
|
|
|
|||||||
|
|
Qmax |
max Q j ( X ) |
|
|
||||
|
|
|
|
[1,M ] |
|
|
|||
При положительном максимальном отклонении |
Qmax обобщенный критерий (2.15) |
||||||||
формируется из положительных функций-отклонений и наоборот. |
|
||||||||
Функции f j (X ) , Q j ( X ) , |
|
j [1: M ] могут быть как линейными, так и нелинейными |
|||||||
функциями параметров X. Исходя из |
условий работоспособности, |
решение ЗВО может |
|||||||
считаться удовлетворительным, если все локальные критерии Q j ( X ) принимают отрицательные значения. Таким образом, задача векторной оптимизации УЦОС сведена к задаче скалярной оптимизации. В качестве целевой функции последней задачи выступает обобщенный критерий оптимальности, формируемый на основе совокупности локальных
критериев Ошибка! Источник ссылки не найден. или функций f j (X ) , j |
[1: M ] . Задача |
||||||
скалярной оптимизации имеет в рассматриваемом случае вид |
|
||||||
|
|
|
min Q( X ). |
|
|||
|
|
|
X |
D |
|
|
|
Мультипликативный и аддитивный обобщенные критерии соответственно |
|||||||
определяются выражениями |
|
|
|
|
|
|
|
|
|
|
m1 |
|
|
|
|
|
|
|
|
iQi ( X ) |
|
||
|
Q( X ) |
|
i 1 |
|
|
|
(2.146) |
|
|
m2 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
k Qk ( X ) |
|
||
|
|
|
k m1 |
1 |
|
|
|
|
m1 |
|
|
m2 |
|
|
|
Q( X ) |
|
iQi ( X ) |
k Qk ( X ), |
(2.17) |
|||
|
i 1 |
|
|
k m1 |
1 |
|
|
где первые m1 функций должны уменьшаться, а остальные m2-m1 - увеличиваться.
Недостатками критериев вида (2.16), (2.17) являются [7, 10, 30, 48]: неограниченная возможность компенсации уменьшения качества по одному локальному критерию
95
увеличением качества по другому (другим), что в ряде случаев маскирует внутренние технические противоречия оптимизируемого УЦОС.
Данные недостатки устранимы при использовании максиминного обобщенного критерия. В максиминном критерии, предложенном в [54], вводятся запасы работоспособности, представляющие собой относительные (безразмерные) оценки выполнения каждого из условий работоспособности
Qj |
j min |
f j ( X ) Z1 j |
1, |
Z |
2 j |
f j |
( X ) |
1 , |
(2.18) |
|
|
|
|
|
|||||
j |
|
|
j |
|
|||||
|
|
|
|
|
|
|
|
||
где j - оценка рассеяния значений j-той функции. |
|
|
|
|
|
|
|||
Величина j задается равной допуску или статистическому разбросу |
f j (X ) . Если эти |
||||||||
величины неизвестны, ими надо задаваться исходя либо из предварительных расчетов, либо из опыта. Например, если путем расчетов или опытным путем установлено, что функция
f j (X ) распределена по нормальному закону с дисперсией |
j, |
то полагая, что в процессе |
||
оптимизации дисперсия |
j меняется мало, можно считать j=3 |
j. |
|
|
При известных оценках рассеяния задачу векторной оптимизации можно поставить и |
||||
решить как максиминную |
|
|
|
|
|
max min Qj ( X ) |
|
(2.19) |
|
|
X D j [1:M ] |
|
|
|
или минимаксную |
|
|
|
|
|
min max |
Q j ( X ) |
|
(2.20) |
|
X D j [1:M ] |
|
|
|
При отсутствии |
оценок рассеяния |
j, j [1: M ] задача |
(2.20) сводится к задаче |
|
Ошибка! Источник ссылки не найден..
Критерий (2.18) также, как и критерий Ошибка! Источник ссылки не найден., не является гладким, что существенно усложняет задачу оптимизации и требует применения специальных алгоритмов [54, 55]. В ряде случаев применяется подход, основанный на процедуре сглаживания критерия (2.20) с последующим обращением к методам гладкой оптимизации [3, 5, 7, 11, 28, 48].
Так как
arg min Qj |
arg max e |
Q j |
|
||
j |
j |
|
то задача (2.20) эквивалентна задаче
max e Q j min . |
(2.151) |
j |
|
96
Здесь к (2.21) применима среднестепенная свертка Ошибка! Источник ссылки не найден.. В результате приходим к следующему обобщенному критерию оптимальности:
M |
|
qQ j ( X ) |
|
|
|
Q( X ) |
e |
min (q 1,2, ). |
(2.22) |
||
|
|||||
j |
1 |
|
X |
|
|
|
|
|
На
Расчет частотных (временных)
|
|
|
|
характеристик H i |
X , |
|
и их |
граничных значений |
|||
|
|
|
|
H Hi , H Bi |
,i |
1, P |
|
Формирование |
|
вектора |
функций отклонения (ЛКО) |
||
|
|
|
e1(X, |
, e2 X , |
,..., eP X |
|
|
|
Выбор номера обобщенного критерия оптимальности из базового набора
Нормирование
функций
ei X , , i 1, P
Расчет численного значения обобщенного
критерия
97
Ко
Рис 2.1 Схема алгоритма расчета обобщенных критериев оптимальности УЦОС
В результате, при решении задач векторной оптимизации УЦОС, когда требования ЧТЗ к ЧХ заданы либо в виде односторонних, либо в виде двусторонних неравенств, в
качестве обобщенных критериев оптимальности целесообразно использовать целевые функции вида Ошибка! Источник ссылки не найден., (2.14) – (2.20) и (2.22). Алгоритм формирования базового набора обобщенных критериев оптимальности УЦОС приведен на рис.2.1.
Рассмотренные выше целевые функции вошли в библиотеку целевых функций подсистемы оптимального проектирования УЦОС [62]. Базовый набор целевых функций,
введенный в состав библиотеки алгоритмов подсистемы, позволяет сформировать различные требования к проектируемому устройству, и обеспечивает решение широкого круга задач оптимального проектирования УЦОС в различной постановке.
2.2 Алгоритмы экономичных вычислений частотных характеристик
цифровых фильтров
Оптимальное проектирование УЦОС крайне затруднено даже с применением
высокопроизводительных ПЭВМ и инженерных рабочих станциях (ИРС). Так как:
-в обобщенные критерии входит большое количество характеристик УЦОС, расчет которых требует значительных машинных и временных ресурсов;
-на каждом оптимизационном шаге используется многократный расчет всего устройства в целом, при этом основное время и ресурсы занимает расчет передаточных функций цифровых фильтров, входящих в состав УЦОС.
Рассмотрим основные приемы сокращения вычислительных затрат при решении задач
оптимизации УЦОС, когда в их структуре имеется либо нерекурсивный , либо рекурсивный цифровые фильтры (НЦФ, РЦФ) [15, 33, 72].
Прежде всего рассмотрим алгоритм расчета частотных характеристик (ЧХ) НЦФ.
Практически реализуемую ЧХ НЦФ можно записать в виде [72]
H (k) A(k)e j (k ) R(k) Y (k ) , |
(2.23) |
где k [1: M ] , M - число точек дискретизации на оси частот. Тогда процедуру расчета ЧХ НЦФ условно можно представить следующей схемой:
98