Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Микроминиатюризация элементов радиоэлектронной аппаратуры

..pdf
Скачиваний:
19
Добавлен:
24.10.2023
Размер:
11.88 Mб
Скачать

T, -

Im

( Д и )

=

{

А.2

-

В,

В2)

(А,

В,

-f Л 4 B3)

f

я

А4 -

-

В3

В,

+

ZI)

(А, В2 •{- А2

В,)

+ 1\

1

В, + А,

В,)

4-

 

 

 

 

 

+

Ц

3 В,

\-

 

А,

В,)

;

 

 

 

 

 

 

7,

=

Re (Ä)

=

Т9

Л,, -

Т, 5-

4- Ц4

 

Z] 4- Л2 2* 4-

 

 

+

Л 2

3

Л4

-

£ 3

Я4 )

-

Вг

3

5 4

4

Л 4

Я,)] ;

 

 

Г( і

-

Im (Д)

=

7 4

Л,

f

Г, В.,

г Z\

4

ZI - f В2

Ц

+

 

 

+

Л 2 а

В,

f

Л4 Д8 ) +

# 2

3

Л 4

-

В,

В,)}

;

 

 

 

 

В,

-

со (С,

4

С,

4- С5 )

 

/ г ? 4 - ( ^ , ) 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ 2

= ш (С2

4- Св

4- С6)

 

 

+

(<к2 )3 ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 з = - и , ( С 3

4 С5 + С7 )

 

 

folg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л , =

щ 4- (»z;

)2

'

 

 

 

 

 

 

 

 

 

 

 

< -

 

 

 

 

 

 

 

Ä> =

coCy;

i

=

1,

2,

 

3;

y

-= 4,

5,

6, 7 .

 

Математическая

задача

оптимизации

параметров

может

быть сведена к экстремальной нелинейной задаче программи­ рования или к задаче решения несовместной системы нелиней­ ных уравнений. Эти две задачи эквивалентны и двойственны,

т. е. одна задача может

быть переведена в другую [4].

 

Для сведения данной задачи к системе уравнений необхо­

димо заданный интервал рабочих частот

Л <

ш < В разбить

на m дискретных значений, определяемых

соотношением

 

А =» ш, < ш, <

. . . < с); <

. . . <

ЯІ =

В .

(8)

В этом случае требуемая частотная характеристика в со­

ответствии с формулой

(2)

на

каждой частоте

будет пред­

ставлять число

 

 

 

 

 

 

 

 

 

ср

(ü>;)

-

bt ,

 

 

 

(9)

а функция (I) становится функцией только

параметров

эле­

ментов

 

 

 

 

 

 

 

 

Уі=/і(Хих2

 

 

XJ

 

*п) •

 

(10)

100

Приравнивая значения этих двух характеристик на часто­ тах, определяемых соотношением ( 8 ) , получаем систему

/ І

(xït

х,_ , . . . ,

Xj , .

. . , х„) -

bt = 0 ,

^

і

1, m, j —

\, п,

т у п

,

 

которая содержит m уравнений с п неизвестными.

Следует отметить, что данная система является в общем случае нелинейной и многоэкстремалыюй, т. е. можег иметь несколько решений.

Необходимо остановиться на соотношении между числом уравнений и числом переменных. Число переменных п четко определяется выбранной структурой фильтра. В решаемой за­ даче число n--13. Число т, как видно из формулы ( 8 ) , может быть выбрано произвольно. Очевидно, что чем больше //;, тем в больших точках будут совпадать требуемая и получаемая час­ тотные характеристики и тем точнее будет решена задача со­ впадения получаемой характеристики с требуемой. Следова­ тельно, с технической точки зрения целесообразно выбрать число m достаточно большим. Однако с его увеличением ус­ ложняется решение математической задачи, что и ограничива­ ет увеличение количества уравнений.

Система (10), в которой т>п, несовместна и в алгебраиче­ ском смысле не имеет решения. Решение такой системы рас­ сматривают как задачу приближения, т. е. задачу определения таких значений х ь х2 , . . . , * „ , при которых величина невязки, определяемая левой частью уравнений системы (11), будет ми­ нимально возможной. Величину невязки можно оценивать по чебышевскому или среднестепенному критериям.

В основе предлагаемого метода решения несовместной си­ стемы нелинейных уравнении лежит идея работы [5], однако используют не чебышевский, а среднеквадратичный критерий приближения. Обусловлено это тем, что применение последне­ го коитерия позволяет более просто реализовать алгоритм на ЦВМ.

Следуя идее работы [5], заменяем исходную задачу после­

довательностью трех более простых задач.

х*,...,

1. Задача аппроксимации нелинейных функций /,'(*'ь

хп) исходной системы (11) в виде суперпозиций функций,

каж­

дая из которых зависит только от одной переменной, т. е..каж­

дая функция fi (х)

заменяется функцией

 

Ф,

= £ ацѴ/(х})-

• (12)

101

Решение этой задачи выполняется таким образом, чтобы от­

личие между множеством

аппроксимируемых

функций / / (х)

и множеством аппроксимирующих функций Ф,(*)

было мини­

мально возможным по среднеквадратичному критерию [8].

2. Задача среднеквадратичного приближения

для несов­

местной системы линейных уравнений. Эта система

получается

из исходной системы (11),

если функции fi(x)

заменить ап­

проксимирующими функциями и выполнить замену перемен­ ных

' O U ' , ) О .

Новая система

уравнений

 

п

аи Г, ' bt = 0, і^-ТГт

 

V

( 1 3 )

является несовместной, но линейной. Методы ее решения из­ вестны как для среднеквадратичного [6], так и для чебышевскоіо критерия [7].

3. Задача определения действительных корней п независи­ мых алгебраических или трансцендентных уравнений

 

 

Wj(xj)-

r;

=

o ,

 

 

( H )

где Г* — решение

несовместной

системы

линейных

уравне­

ний

(13). •

 

 

 

 

 

 

Следует отметить, что характер уравнений

(14) определяет­

ся функциями

Wj(xj),

выбранными на первом этапе.

 

Если в качестве

I ' , - (Xj)

принять обыкновенные полиномы,

определяемые

соотношением

 

 

 

 

 

 

Wj

(xj)

= (i0 f dx

x j \

 

. . . i- f/v

X) ,

(15)

то определение их коэффициентов d0,

dlt...,

dw

сводится к ре­

шению несовместной системы линейных уравнений. Определив функции Wy- (xj) и задав определенное множество дискретных значений параметров хи х^... ,х„, можно свести задачу опре­ деления коэффициентов а0- также к решению несовместной си­ стемы линейных уравнений. Таким образом, определение мно­ жества Ч/Xj), определение матрицы [А] и решение второй за­ дачи, определяемой соотношением ( 1 3 ) , сводится к решению одной и той же задачи — решению несовместной системы ли­ нейных уравнений.

Применение среднеквадратичного критерия позволяет со­ здать для решения этих задач достаточно простую и легко реа­ лизуемую на ЦВМ стандартную процедуру. Эта процедура определяется соотношением

102

II Г*

i| = Il A+ II X II * Ii

(16)

и состоит в нахождении

псевдообратной матрицы || А 4

|| от ко­

эффициентов уравнений и умножении ее на столбец свободных членов II b II [6].

Поскольку уравнения (14) независимы, то искомые значе­ ния х2 состоят из набора возможных корней уравнений. Сле­ довательно, возможные решения исходной системы определяют множеством

 

X. =

{*!*, ,---,Xjk.

xnkn)

;

kj

= 1,

$j ,

 

(17)

іде

 

 

 

 

 

 

 

 

 

 

XjK

kj-ü

корень /-го уравнения

(14);

 

 

 

 

ßy- — число действительных корней /-го уравнения.

 

Кроме решений системы (11), в множестве

(17)

появляются

лишние корни, которые должны быть исключены.

Для

этого

используют свойство выпуклости суммарной функции

 

 

 

 

 

m

 

 

 

 

 

 

 

 

М А - )

=

е

/ , ( * ) .

 

 

(is)

в соответствии с которым для

х+,

являющимся решением си­

стемы

(11), выполняется условие

 

 

 

 

 

 

 

h \>У

•+ (1 - ).) х"\ <

\h

(х')

+

(1 -

/,) Д

(х") ,

(19)

где х' и х" — значения внутренних

параметров в

достаточно

малой окрестности исследуемого решения

хл,

0 < Х < 1 .

 

Именно это свойство и используют в алгоритме для опреде­ ления всех возможных решений исходной системы (11). Прове­ ряют на выпуклость суммарную функцию (18) и набор пара­ метров, удовлетворяющих выпуклости функции, записывают в

множество х'~. Следовательно,

множество

х+ содержит

все

те наборы х\, Х2,...,хп,

которые

являются

решением системы

(11). Среди множества

х' необходимо найти оптимальное

ре­

шение или глобальный экстремум исходной задачи. Оптималь­

ным решением будет такое значение

х* , при котором обеспе­

чивается минимальная

невязка

[5]. Для

определения глобаль­

ного

экстремума используют

метод

моделирования,

т. е.

сравниваются величины суммарной

невязки

для

различных

,ѵ+,

принадлежащих

множеству х ; ,

и искомым

решением

х'а

будет то значение,

которое

удовлетворяет

условию

мини­

мума суммарной невязки. Для этого производится случайный

выбор

A"j из множества х 1 и сравнивается

невязка для зна­

чения

Х'1' с невязкой для значения xu : _,. Если

Z.(.v+)<Z.(x+,),

103

т. е. если погрешность

для исследуемого значения

ха

мень­

ше, чем для предыдущего значения параметра

" то соот­

ветствующее значение

л; запоминается. В итоге будет опреде­

ли j иіхѵЗних

 

 

!7Г

 

 

 

.

f

 

 

 

 

i

z

 

*"* Определение

ry 7-T~~

z

I

n

v

t

u

 

Г

 

 

 

trÛ

ïpêtsIf/OtrrtU

 

Решение

it

 

 

 

Рис. 2

лено значение дгк'еЛГ+, которое обеспечивает минимальное значение суммарной невязки, т. е.

L (xl) => min L (xa ) .

xj « .Y

Алгоритм оптимизации ФСС представлен на рис. 2.

104

Таким образом, предложенный алгоритм реализует в себе два метода и его можно назвать комбинированным. Причем первый метод, основанный на идеях статьи [5], позволяет ло­ кализовать экстремумы, т. е. определить достаточно грубо ок­ рестности решения исходной системы уравнений, а вторая часть метода, использующая свойства выпуклых функций, по­ зволяет находить достаточно быстро решение с требуемой точ­ ностью. Именно поэтому сочетание этих двух методов и позво­ ляет найти глобальный экстремум задачи с заданной точ­ ностью. Применение методов выпуклого программирования или методов случайного поиска для решения сформулирован­ ной задачи, а также других известных методов гарантирует нахождение только локального экстремума. В этом принципи­ альное отличие предлагаемого метода от известных методов решения нелинейных экстремальных задач.

Алгоритм реализован на ЦВМ «М-222».

 

 

 

 

 

 

ЛИТЕРАТУРА

1.

С о б с І І и и Я. М-

Расчет полиномиальных фильтров. Связьиздаг, 19G3.

2.

M а т т е и

Д. Л.

и

др.

Фильтры СВЧ, согласующие цепи связи.

М„ «Связь», 1971.

 

 

 

 

3. С и г о р с к я й Рі. П.

Анализ электронных схем. К., «Техаіка», 1963.

4.

В о л к о в

В. М.,

К у р и

л и и Б. И. Методы оптимального проекти­

рования сложных

систем

(на

украинском языке). К., «Техщка», 1971.

5. К у р и л и и

Б. I I . К решению чебышевской задачи приближения для

несовместной системы линейных уравнений. «Журнал вычислительной ма­

тематики

и математической физики», т. 10, № 1,

1970.

 

 

6.

Г a

H

т м а X с р

Ф. Р. Теория матриц. М.,

Гостсхпздат, 1953.

7

Р е м е з fi.

Я.

Общие вычислительные методы чебышевского прибли­

жения. Изд. АН УССР. 1957.

ил и H

 

 

 

 

 

8.

В о л к о в

В. M , К у р

Б. И.,

M и р о н ы ч е в Я

Л., О р е ­

х о в

Е. Ф. Оптимизация параметров и допусков элементов радиоприемных

устройств.

Сб. «Некоторые

вопросы

оптимальной

обработки

сигналов и

улучшения характеристик РЛС». Изд. КВИРТУ, 1971.

 

 

9.

К а л а X а н Д.

Современный синтез

цепей. М.,

«Энергия», 1966.

УДК 621.375.019.3

Е. К. Вервейко, В. М. Волков

АЛГОРИТМ ОПТИМАЛЬНОГО РАСЧЕТА ДОПУСКОВ ПАРАМЕТРОВ РАДИОТЕХНИЧЕСКИХ УСТРОЙСТВ

В статье методом статистических испытаний решена задача оптимального распределения допусков параметров элементов радиотехнического устройства по заданному допуску на его выходной параметр.

Аналитические методы решения подобных экстремальных задач удается проводить далеко не всегда. Для преодоления математических трудностей приходится вводить ограничения и допущения, и чем их больше, тем меньше доверия конечному результату. А главное, довольно часто отсутствует явная ана­ литическая зависимость между выходными и входными пара­ метрами устройства.

В этом отношении обладают большей универсальностью алгоритмические методы, не дающие конечную формулу, а лишь указывающие путь отыскания экстремума.

Довольно широко внедряется в практику расчета электрон­ ных устройств метод статистических испытаний. Он особенно эффектен для нелинейных цепей, для цепей с неявной аналити­ ческой зависимостью, в условиях произвольных законов рас­ пределения и неограниченного разброса входных параметров. Большие перспективы имеет этот метод для расчета ИС и БИС.

В статье предлагается алгоритм оптимального распределе­ ния допусков на входные параметры устройства при заданном допуске на выходной параметр. Алгоритм ограничен для уст­ ройств с 99% вероятностью безотказной работы и со статисти­ чески независимыми входными параметрами.

Предполагается, что в нашем распоряжении имеется дискрет­ ный ряд допусков на параметры элементов и, что вполне есте­ ственно, чем меньше допуск, тем выше стоимость элемента.

Постановка задачи.

Много факторов

влияет на точность

выходных параметров

радиотехнических

устройств, и многие

из них носят случайный характер. Усложнение задач

практики

постоянно повышает требования к точности выходных

парамет­

ров радиотехнических устройств, а это достигается все больши­ ми материальными затратами. Вот почему актуальной задачей

является достижение заданной точности на выходной

параметр

с помощью минимальных затрат.

 

Нам известны: топология схемы, номинальные

значения

элементов, дискретные значения допусков на элементы устрой-

106

ства, критерий для характеристики устройства, допуск на вы­ ходной параметр, данные стоимости элемента в единицу до­ пуска. Задана максимально возможная сумма расходов Сн.

Решение. Самый простой, но и самый неоптимальный путь для решения задачи — перебор всех возможных сочетаний до­ пусков на параметры элементов, статистическая проверка на надежность и выбор затем оптимального сочетания. Для боль­ шого количества элементов, даже с использованием самых со­ временных ЦВМ, решение этой задачи потребует много време­ ни. Необходимо искать пути совершенствования решения.

1. Ограничение по критерию «наихудшего случая».

Радиотехнические элементы промышленность выпускает с различными допусками на параметры. Например, резисторы имеют допуски на сопротивление ±20; 10; 5; 2; 1%. Все ли до­ пуски на параметры элементов надо рассматривать при реше­ нии нашей задачи? Допуск на выходной параметр А)> будем определять через допуски на входные параметры { по фор­ муле

О)

где

КІ — коэффициент, характеризующий вес /-го элемента в определении выходного параметра;

а — количество элементов.

Из этой формулы можно найти наихудший возможный до­ пуск для параметра /-го элемента &х\. Определим его в пред­ положении, что все остальные элементы имеют номинальные значения, а лишь /-й элемент обеспечивает допуск на выходной параметр. Очевидно, что в дальнейшем нужно рассматривать только допуски, меньшие Ах',.

Составим таблицу возможных значений допусков на пара­ метры элементов, располагая их в порядке убывания абсолют­ ных значений (табл. 1). Здесь Рі— параметр /-го элемента, s — возможное количество допусков на параметр /-го эле­ мента.

Рассмотрим структуру многокорневого дерева (рис. 1). Точ­ ки первого столбика представляют собой элементы S\ — полу­ ченные допуски для параметра р\.

От каждой точки первого столбика идут ответвления, веду­

щие ко всем элементам s2

полученным допускам для рг

и

т. д. И, наконец, последний

столбик

имеет

все' элементы

sn,

повторяющиеся при каждом

выборе

(и—1)

предшествующих

допусков элементов.

107

Таблица 1

 

....

Pi

P

 

 

 

 

 

дат,,

 

л:Х.л

 

л . х а

-

 

 

f

 

 

 

- -

 

 

Л*,. S

 

Л:>',.:

tarn .»„

Рис. i

108

 

Существует С1} (і— 1, 2, ... , п) точек в последнем

столбике

и полное количество путей

от левой

части дерева

к правой,

Каждый путь связан с выбором вектора допуска ДЛ\

следова­

тельно, имеет соответствующую стоимость

 

 

С =

V

С, ( \Х*)

,

 

г де

Ci (àX,-) — функция стоимости допуска /-го элемента;

 

k — номер столбика табл. 1.

 

 

 

Из всего многообразия

возможных

путей дерева

необходи­

мо найти искомый. Будем исходить из предположения, что если допуски на параметры элементов удовлетворяют условию (1), то в этом случае будет обеспечена высокая вероятность безот­ казной работы устройства. Бесспорно, что это утверждение име­ ет высокую степень достоверности. Обозначим область воз­ можных путей, удовлетворяющих условию (1), через х*.

2. Метод определения самого дешевого пути.

Задача определения самого дешевого пути сводится к опре­ делению пути через дерево слева направо так, чтобы удовлет­ ворялись два условия:

— путь входил в область допустимых путей

х*;

— стоимость

 

пути

 

 

 

 

 

 

 

 

п

 

 

 

должна

быть

меньше

Сѵ .

 

 

 

И вот тогда уже из оставшихся путей будет выбран опти­

мальный, т. е. путь, имеющий минимальную стоимость.

Алгоритм определения такого пути представлен на рис. 2.

Начнем слева

по табл. 1 и возьмем для р\ самый большой

допуск

АЛ"ы не больше

аХ'и

Будем двигаться через дерево,

выбирая допуски

k-x элементов

при

 

 

где функция

Т/,А (bXt)

характеризует

путь до k-ro элемента

по условию (1). После каждого

такого выбора

подсчитываем

частичную стоимость Ск

для пройденной части пути.

Нас может остановить только два обстоятельства:

— допуски на параметры элементов пройденного пути уже

превысили заданный допуск на

выходной

параметр по усло­

вию (1);

 

 

 

 

 

 

109

Соседние файлы в папке книги из ГПНТБ