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

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

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

тьет точность решения исходной задачи. Однако с увеличени­ ем m не только возрастает объем вычислений, но может изме­ ниться характер задачи. Действительно, при яг>н-система (4) становится несовместной и в алгебраическом смысле не имеет решения. В этом случае решение задачи рассматривают как задачу приближения по чебышевскому пли среднестепенному критерию.

Один из методов решения несовместной системы нелиней­

ных уравнений

предложен в работе

[ 7 ] . Этот метод

положен в

основу алгоритма решения системы ( 4 ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЛИТЕРАТУРА

1. Ф е л ь д ш т с й н

А. Л., Я в и ч

 

Л. Р.

Синтез

четырехполюсников и

восьмиполюсников на СВЧ. М., «Связь», 1971.

 

 

 

 

 

2.

Л и т в и и е н к о О. Ы., С о ш н и к о в

В. И.

Колебательные

системы

из отрезков неоднородных линии. М., «Сов. радио»,

1971.

 

 

 

3.

Ку р и л и н

Б. I I . Колебательные

системы из отрезков фидерных ли­

ний. К., «Техщка»,

1969.

 

 

 

 

 

 

 

 

 

 

4.

Р е м е з

Е. Я. Общие

вычислительные

методы

чебышевского

прибли­

жения. Изд-во АН УССР, 1957.

 

 

 

 

 

 

 

 

5.

Б е р с з н n

И. С.,

Ж и д к о в

Н. П.

Методы

вычислений,

т.

1 и 2.

Фпзматгиз,

1959 и

1960.

 

 

 

Y.

Ci. Il а с

 

 

 

 

6.

B r a

у t o r i

К. К.,

G u s t a v s o n

h te I G.

1).

 

A New

F.ssecient Algorithm for Solving Differential Algebraic

Systems

Lsin^

Im ре -

lieit

Backward

Differention

Formulas

Pros. IEEE. v. GO, № I .

 

 

 

7.

К у р и л и "

Б. И.

К

решению чебышевской

задачи приближения для

несовместной системы нелинейных уравнений. «Журнал вычислительной ма­ тематики и математической физики», МЬ 1, 1970, с. 3.

УДК 621.3.089.52

Б. И. Курилин, В. П. Пуганов

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

В

статье предложены

метод и алгоритм

минимиза­

ции

параметрической

чувствительности

устройств. В

ос­

нове

метода лежат

свойства полиномов

Чебышева

и

принципы, чебышевского

приближения

для

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

системы нелинейных

уравнений.

 

 

 

Вопросы минимизации параметрической чувствительности различных устройств в последнее время приобретают все боль-

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

Одним из основных элементов ряда функциональных узлов РЭА являются резонансные системы. Они определяют рабо­ чий диапазон, электромаіннтную совместимость и некоторые другие характеристики передающей, приемной, измерительной аппаратуры. Поэтому минимизация отклонения характеристик резонаторов при изменении параметров элементов, особенно активных, является актуальной задачей.

Разработаны [3, 4 ] методы решения дайной задачи, позво­ ляющие получить заданное уменьшение чувствительности тою пли иного показателя функционирования РЭА к отклонению одного или нескольких параметров.

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

При решении задачи минимизации чувствительности исхо­ дим из того, что определена аналитическая зависимость харак­

теристик устройства Г,-

от п внутренних

параметров

х

 

УІ - ( * j

-г, ,

**

Xj

, *J

( О

и заданы пределы изменения N = n—А' внутренних переменных параметров Xj(j — k+\,n)

dj <xj<:П}.

(2)

Необходимо определить такие значения к постоянных внут­ ренних параметров л\ , • . . , д", при которых, величина макси­ мального отклонения выходных характеристик была бы мини­ мально возможной

АК^ = min А К; , л — IA'J , . . . , xw , . . . , хк) .

X

При решении задачи минимизации чувствительности в ин­ тервале (2) используем свойство полиномов Чебышева мини­ мально отклоняться от нуля [5] в интервале z — \1, 1]. При этом

max I Тп (?) | = 1;

min Тп (г) = 0 .

г

z

Отклонения Л К,, не будут превышать заданной величины, если положить

121

где Xj — текущее значение переменного параметра, приведен­ ное к стандартному интервалу [—1, 1] на основе соотношения

 

 

 

_

 

 

2xj

- dj

-

Dj

 

 

 

 

 

 

 

X j

=

 

Dj

-

dj

 

'

 

 

 

a величина

удовлетворяет условию

 

 

 

 

 

 

 

 

О <

Ay <

 

2BtJ

 

'

 

 

 

 

 

 

-^jzTjr

 

 

 

 

где

B,j

— допустимое

отклонение

àYy

i-'k

характеристики

при изменении ;'-го параметра.

 

 

 

 

 

 

 

 

С учетом условия (3) задачу можно

свести

к решению си­

стемы уравнений вида:

 

 

 

 

 

 

 

 

 

 

 

 

 

/ > ( * ) - Ѵ = о ,

il =

ТГл*\

 

 

(4)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fv- = fi (xJ)

~

fi

(xj

±

AXJ) .

 

 

 

если

Axj

достаточно велико и

 

 

 

 

 

 

 

 

 

 

 

F

=•-

S •

Vi

 

 

 

 

 

 

 

 

 

dx7

 

 

 

 

 

 

если Ах,- мало, £i a .=const.

 

 

 

 

 

 

 

k неизвестных

и получена

из равен­

Система

(4) содержит

ства

(3)

при условии, что в стандартном

интервале z для каж­

дого

Xj выделено pj значений

х,

 

удовлетворяющих

неравен­

ству

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 < xjy

<

. . . <

xje

< . . . < xJ{

< 1 .

(5)

При этом общее число уравнений в системе

(4)

равно

 

 

 

M

=

тК,

К =

П

Pj •

 

 

 

Число Pj

необходимо выбирать так, чтобы M >

k.

 

Система

уравнений

(4) может быть

как совместной, так и

несовместной. Методы решения совместных систем уравнений известны [6J.

Решение несовместной системы принято [5] рассматривать как задачу приближения по чебышевскому или среднестепенному критерию.

122

Фѵ ѵ )
(8)

Решение по чебышевскому критерию приближения состо­

ит в нахождении таких значений х*=

(х'ѵ...,х*к),

при которых

максимальная невязка

 

 

А = max I Ff, (х) ~

bt, \ ,

(6)

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

А = Д * = m i n А .

X

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

 

А =

M

(/V (х) - Ь,Г .

 

 

(7)

 

М.=І

 

 

 

 

V

 

 

 

 

 

На практике наиболее часто используют частный случай

этого критерия

— среднеквадратичный,

соответствующий

? = 2.

 

 

 

 

 

 

 

Чебышевский

критерий

более предпочтителен

[7],

однако

решение по среднестепенному критерию

проще

и менее гро­

моздко при реализации на ЦВМ.

 

 

 

 

Известные методы решения этой задачи

обладают

рядом

существенных недостатков, в частности,

не

всегда

обеспечива­

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

Можно предложить следующий алгоритм (рис. 1) решения

системы

(4) по среднеквадратичному критерию. В основе алго­

ритма лежит

идея [8] замены исходной задачи последователь­

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

Алгоритм содержит следующие основные операторы.

Оператор

1 формирует

исходную информацию и вводит ее

в машину.

 

 

 

Оператор

2

выбирает,

случайно или детерминированно,

функцию

Fv,

(х)

б / .

 

Оператор 3 выбирает аппроксимирующие функции

/Ѵ(*) = £ Ф*Ы

ѵ = 1

с помощью стандартной подпрограммы [9].

Оператор 4 определяет элементы матрицы А по формуле

û,v = у F,, (АѴ) Ч\ ( j C v T ) ,

(9)

123

полученной из условия минимизации среднеквадратичной ошибки приближения:

V

«... 'Г.

)

dx

(10)

 

 

 

 

I

'

при разбиении интервала t/, < хѵ

<

определяемого усло­

виями технической реализации, на

* (•: -=*•

к

1)

фиксирован­

ных точек.

' >.. :<

au

у-; ,

 

и с te i'/.'t.'.

 

boni!

 

 

) ыосф

 

ici no H Н на

F,uc{>!(z.f

 

nut

 

 

Определение

 

 

Z T Z

 

il

ônfitoejttHtje

9

 

 

вычисление

Л

 

д tek)

5

 

Реал: >>ч/е

 

 

Сися-гемм (а)

 

 

 

_..

f

4 <.д jor?.

* ЬыЬор Ъ€ X *

7

On1>ь t. •' -if ние ii»О* -іа X*

Р.ІГ. 1

Оператор 5 решает несовместную систему

ь

 

 

V г/,,. Г ~ / ѵ

0 ,

(11)

полученную из формулы (4 ) путем замены функций FIL (х) аппроксимирующими функциями ( 8 ) и перехода к новым пе-

121

ременным Y». — Чт^ (.%). Решение выполняется на основе среднеквадратичного критерия.

Оператор 6 решает M уравнений

«V M

- )Ѵ

(12)

и находит множество X — {Хх

, . . . , Л'г }, состоящее из Г =

С*

сочетаний но k параметров

 

 

 

л '

-

І

Л', .

 

где Л', — число корней

ѵ-го уравнения.

 

Оператор 7 определяет

множество А"~ 6 АГ сочетаний,

яв­

ляющихся решением системы уравнений (12). проверяя выпук­ лость всех функций системы,

Оператор 8 выбирает сочетание,

определяющее

значение

вектора X, принадлежащего

множеству

А'",

т . е .

 

А", б

А 4 .

Оператор 9

вычисляет

среднеквадратичную ошибку

(10)

для выбранного значения А'».

 

 

 

 

 

 

Оператор 10 сравнивает

величину ошибки А (Х„),

соответ­

ствующей Х„

с ошибкой A (АѴі)

для предыдущего резуль­

тата Ха-І.

 

 

 

 

 

 

 

 

Если А (Х„)

< A (АѴіК

то сочетание Ха

запоминается

оператором 11,

иначе управление

передается

оператору

12,

который определяет, все ли

сочетания

множества

Xh

прове­

рены. Если вычислены ошибки для всех сочетаний

Л'~;~, то опе­

ратор 13 сравнивает минимальную ошибку, записанную в опе­

раторе

I I , с допустимой ошибкой. Если А (Х„) <

A t r n , то ре­

зультат

(Л"*, А) выводит на печать оператор 15. В

противном

случае управление передается оператору 14, который опреде­ ляет, проанализированы ли все M функции при решении зада­

чи аппроксимации.

Если

р-0<Лі,

то управление

передается

оператору 2

для

выбора

нового

класса

аппроксимирующих

функций, иначе результат выводится на печать.

 

 

 

 

 

 

ЛИТЕРАТУРА

1. А й э и н о в

М.М. Избранные вопросы теории

сигналов

и теории це­

пей. М., «Связь»,

1971.

 

 

 

 

 

2.С и г о р с к и й В. П. Содержание и методы теории электронных це­ пей. «Известия вузов, Радиоэлектроника», № 7, 1967.

3.К а л а X а н Д. А. Методы машинного расчета электронных схем. М., «Мир», 1970.

4. С и г о р с к и й В. П.,

П е т р е н к о А. И. Алгоритмы анализа элек­

тронных схем. К., «Техніка»,

1970.

5. Р е м е з Е. Я. Общие вычислительные методы чебышевского прибли­ жения. Изд-во АН УССР, 1957,

125

6.Б е p е з и н И. С , Ж и д к о в Н. П. Методы вычислений, т. 1. М., Физматгиз, 1959.

7.П е р в о з в а н е к и й А. А. Критерий равномерного приближения в задачах оптимального управления. «Оптимальные системы. Статистические

методы». Труды I I I Всесоюзного совещания по автоматике. М , «Наука», 1967.

8. К у р и л и и Б. И. К решению чебышевской задачи приближения для несовместной системы нелинейных сравнений.-Журнал вычислительной ма­ тематики и математической физики, 197Ö, 10, № 1.

9. Ш у р а - Б у p а М. Р. Аппроксимация функций многих переменных функциями, каждая из которых зависит от одной переменной. В сб. «Вы­ числительная математика». Вып. 2, Изд-во АН СССР, 1957.

УДК 621.396.6

И. Н. Бобров, H. М. Шапкин

СИГНАЛЬНЫЕ ХАРАКТЕРИСТИКИ ТУННЕЛЬНОГО ПЕРЕХОДА

На основании

эквивалентной схемы

и

аппроксимаций

вольт-амперной

и

кулон-вольтной

характеристик

тун­

нельного перехода найдена зависимость его

входной

про­

водимости от напряжения смещения,

частоты и

ампли­

туды гармонического

напряжения,

действующего

на тун­

нельный переход,

и электронных

свойств

собственного

перехода.

 

 

 

 

 

 

 

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

Эквивалентная схема такого'перехода показана на рис. 1, На этой схеме g — дифференциальная проводимость туннель­ ного перехода; С — емкость перехода; г — сопротивление по­ терь перехода.

Входная проводимость туннельного перехода

_

г

( g2

2

+ со*2

 

С

2

)

+

g

,

, _

 

 

 

 

 

 

 

С

2

)

 

в *

г2

( g

+ < u C

2

H - 2 r g + 1

ь J

г- {g- Ч со2

 

- f 2rg - f 1

 

 

 

 

 

 

126

Выражая г через предельную частоту перехода

п р

 

 

 

 

 

 

 

 

 

+ 8'

 

 

 

С

 

 

 

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которая определяется из равенства нулю

активной входной

проводимости, получаем

 

 

 

со

 

 

.со

 

С

 

о

 

со' С -

S- tuf

 

е2

£

2

- f

4

2

2

пр

 

о

 

 

 

 

1 +

пр

" п

р С 2

 

пр

+

 

 

 

 

2

* > 2

 

С

2« С

 

 

 

. te2

2

+

2

<оп2р С 2 )

 

 

 

4 У

g

 

w

 

" п р С

 

 

 

 

 

 

Введем безразмерные параметры

5 —

со сопр

I нормированную проводимость туннельного

перехода и

— нормированную частоту.

. - 1 7

Рис.

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

У

_

(S + ! ) ( ! -

Q)

,

,.

(5 4- 1)2 соС

0)

7

вх

S

С О

I

 

 

'

У

59

4- 1

 

 

 

 

50

4- 1

 

 

 

 

 

Для реальных туннельных переходов при

w <

сопр

 

 

 

 

 

5 « 1 ,

 

а <

1 .

 

 

 

Поэтому выражение

(1)

можно представить в виде:

 

 

Гвх =

О в х f

у Я м

^

g

(1

-

Q) +

jmC .

 

(2)

127

Дифференциальную проводимость и емкость туннельного перехода определим из его вольг-амперной и кулон-вольтной характеристик [1].

С достаточной для практики точностью кривые емкости перехода на падающем участке вольт-амперной характерис­ тики можно аппроксимировать выражением

 

 

 

 

Сп

— С,

 

I

 

 

 

 

С — С a

yr~ZZ'\J~

 

 

 

где

 

 

 

 

 

 

 

 

 

 

Сп — емкость перехода,

соответствующая

максимуму

(пи­

 

ку) вольт-амперной характеристики;

 

 

 

С Б

емкость

 

перехода,

соответствующая

минимуму

(впа­

 

дине) вольт-амперной характеристики;

 

 

и — напряжение на переходе.

 

 

 

 

Если и = ѵ + Е представляет собой сумму

высокочастотного

сш нала ѵ и напряжение смещения П, то

 

 

 

 

С=*С„-~

 

-[ѵ~(Ѵп-Е)},

 

 

(3)

где

 

 

 

 

 

 

 

 

 

 

АС = С п

С„ ,

 

 

 

 

 

 

АѴ = Св

-

 

Ѵв .

 

 

 

 

 

 

Если

о =

Ѵт

cos (at+

?).

то (3)

принимает

вид

 

 

 

 

 

С = С0

f

cos

4- ?)

,

 

(4)

где

 

 

 

 

 

 

 

 

 

 

C0 = C„

AC

 

(£" — V'„)

— постоянная составляющая и

 

 

Д 1/

 

 

 

 

 

 

 

 

д с

 

 

 

 

 

 

 

 

 

С! =

~

 

 

— переменная составляющая

емкости пере-

хода.

Вольт-амперную характеристику туннельного перехода с достаточной для практики точностью можно представить поли­ номом 5-й степени [2] вида

і - - ^ г | 5

(и -

Ѵп) (и

- Ѵ,У - (и -

Ѵ.у\ + / в ,

(о)

где

= /„ -

/

 

;

 

 

 

Д /

в

 

 

 

/„

— ток

максимума вольт-амперной характеристики;

/ в

— ток

минимума

вольт-амперной

характеристики,

 

128

При и-^ѵ + Е выражение (5) принимает вид

 

А /

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

і

=

 

( 4 ^

-

5 [ 3 ( Ѵ „ - £ ) - ( £ - С / п )

г;* 4

2 0 ( V B - £ )

[V.

-

-Е-(Е-Ѵп)}ѵ*-\0

 

 

 

 

ш ~ £)'-[ Ѵв

- 3 ( £ - Ѵа)\ гі?

-

 

-

20

(V, -

ЕУ ( £

- ]/„) г. - f

(К. -

£ ) 4 [V. -

£

+

 

 

 

 

 

 

 

 

 

+ 5 ( £ - Ѵ „ ) ] } -1- / 8 .

 

 

 

 

 

 

Дифференциальная

проводимость туннельного

перехода

 

 

g -

- ^ - - 2 0 - ^ 5 - { г ^ - [ 3 ( Ѵ в - £ ) - ( £ » Ѵ в ) 1 г . Ч -

4

3 (Ѵв

-

£ )

 

[Va -

Е -

(Е -

Ѵп)}

V* -

в -

Еу-

[V, -

£

-

 

Если

 

_ 3

( Е - У п ) ] г > - ( Ѵ в - £ ) » ( £ - Ѵ п ) } .

 

 

 

 

 

 

 

 

 

 

 

 

 

ѵ—Ѵт

 

cos (со£ - f ?),

то дифференциальную

проводи­

мость можно представить как

 

 

 

 

 

 

 

 

 

 

g = ßa

 

4- £i

cos

(w* +

?)

+

g a cos 2 (со* + <p)

+

 

 

 

 

 

4

g 3

cos 3 (ш* -f- ?)

4

g4

cos

4 (wf

4 <p) ,

 

 

(6)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

go - 20 - ^ { J - Ѵ*я + }

F

e (V. -

E)

- (V„ -

£ ) •

( £ -

Vn )} ;

- £ ~ 3 ( £ - V'n )l Vm) ;

£ 3 = 10

{ ^ . + 3 (V. - £ ) [V. - £ - ( £ -

{ѴШ-Е-{Е-ѴЯ)\-

Ѵв ) V%) ;

Л e - 5 д ^ Г ^ [3 (V, ~ E ) - ( E - V„)\ ;

Подставляя

(4)

и (6)

в

(2),

получаем

 

Ув*

- Go

-f

G) cos

(ш* +

<p) +

G2 cos 2

4 ?) 4

4 Os cos

3 (иИт)

4 0A

cos

4 (<of+<p)

4 j [ ß 0 + ß i

cos (w*+<p)],

9. Зак. 205.

129

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