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

книги / Оптимизация параметров импульсных и широкополосных усилителей

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
17.93 Mб
Скачать

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

П

р„ (х)= 2 а‘х‘’

/=о

аппроксимирующего функцию у(х) на интервале [а, Ь], начинается с приближения не самой функции, а ее п-й производной. Старший коэффициент ап полинома опре­ деляется из равенства

ь

ь

*/(") (х) dx =

| Р^1) (л:) dx = п\ап — а),

а

а

представляющего собой условие выравнивания площадей, ограниченных кривыми у№ (х) и Р^1) (л:) и осью абсцисс

на интервале аппроксимации. Отсюда

 

у(п-\) (Ь) — у(п-\) (л)

 

а“—

л! (i> - а)

-

Коэффициент ап-1 определяется из аналогичного усло­

вия

ь

ь

 

 

 

 

(х) d x = ^

Р^я-1) (*) d x =

 

а

а

 

 

=

{п— 1)! (Ь — а)

 

+ -J- (b* - а?)ап.

В результате получаем:

 

 

 

a«_i

__ //(«“ 2) (6)_f/(rt-2)

(Л)

~2~ФЛ~Я) а,х'

(Л— i)j

— л)

Так же находятся остальные коэффициенты, за исключением свободного члена, определяемого услови­ ем Р п( а ) = у ( а) • Нетрудно показать, что тогда выпол­ няется и равенство Рп(Ь) —у(Ь).

Таким образом, значения функции и аппроксимирую­ щего полинома совпадают на концах интервала аппрок­ симации. Значения их производных в рассматриваемом

31

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

Воспользуемся изложенным методом для аппрокси­

мации экспоненты

ех

в интервале [0,

а] полиномом

Р2(х) =ao+aiX-ba2X2.

Коэффициент

определим так,

чтобы выполнить условие

 

а

а

 

 

£(**)" dx =

^ ? '\ (x )d x , откуда

=

о

о

 

 

Аналогично найдем коэффициент сц, потребовав вы­ полнения равенства

С(е*)' dx = f Р\ (х) dx.

О

о

 

Отсюда получаем;

 

 

 

еа— 1

<?“ — 1. а0 = е° = 1.

 

 

~ ~ а ~

2 *

 

Следовательно, экспоненту ех в интервале (0, а]

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

 

 

 

х + ^ х > .

Например, при а = — 1 имеем:

 

е* =5=1+

0,9482л:— 0,312

(1-22)

Так как в конце интервала погрешность равна нулю,

то интегральной

аппроксимацией можно

пользоваться

и при | я | > | а | .

Например, полином (1-22)

дает относи­

тельную погрешность, не превышающую 2% в интерва­ ле [— 1,1; 0]. Сводная таблица коэффициентов полино­ мов, используемых в дальнейшем для аппроксимации экспоненты, приведена в табл. 1-6. Там же указана гра­ ница I интервала аппроксимации •[—/, 0] и наибольшие относительные погрешности |е макс1 % .

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

32

 

 

 

 

 

Т аблица 1-6

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

Г

а1

а2

а3

«4

1

1

 

 

 

 

 

Емакс||

 

 

 

 

 

 

 

%

Маклорена

2

1 , 0 0 0 0

0 ,5 0 0 0

 

___

0 ,4 5

2

 

3 1 , 0 0 0 0

0 ,5 0 0 0

0,1667

0 ,7 8

3

 

4 1 , 0 0 0 0

0 ,5 0 0 0

0,1667

0,0 4 1 7

1,07

3

Интегральная

2

0 ,9 4 8 2

0,3161

0 ,0 7 7 3

 

1.Ю

2

 

3

1,0063

0 ,4 4 0 5

2 , 0 0

3

 

 

 

 

 

 

 

 

4

1,0174

0,5 0 8 7

0,1 3 7 7

0,0 1 5 3

2 ,5 0

3

Модифицирован-

2

0 ,8 8 3 3

0 ,2 4 5 9

 

1 .5

2 .3

ная чебышевская

3

0 ,9 5 8 7

0 ,3 8 6 3

0,0617

2 , 0

0,6

Чебышевский метод приближения позволяет наилуч­ шим (в указанном смысле) образом аппроксимировать непрерывную функцию f(x) на замкнутом и ограничен­ ном множестве Е с весом р(х) полиномом Р п(х)•

Эффективным методом решения этой задачи является второй полиномиальный алгоритм Ремеза [7]. При на­ хождении аппроксимирующего полинома за основу взят этот алгоритм с некоторой модификацией [8]. Были по­ строены полиномы Р п(х), у которых одинаковые макси­

мумы относительного отклонения

f (х) Р п(х) fix)

достигаются в п+ 1 точке в отличие от чебышевского подхода, когда требуется выполнить это условие в п+ 2 точках. Требование максимального отклонения в п + 2

точках заменено условием f(0 )= P n(0). Это уменьшает относительную погрешность при расчете коротких импульсов ценой очень небольшого возрастания макси­ мальной погрешности в интервале.

Полином нулевого приближения Р^0) (л:) строится как

наилучшее приближение функции f(x) на множестве п+ 1 точек, выбираемых произвольно внутри интервала.

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

3—195 33

ные относительные отклонения 6o(*t) полинома от аппроксимирующей функции и точки Хи в. которых они

достигаются.

Находится поправка А ^ 0) (*) к полиному нулевого приближения и образуется полином первого приближения

( * ) + ЛРГ (•*)•

Этапыповторяются до выравнивания модулей экстремальных относительных отклонений с заданной

точностью. 'Процесс выравнивания показан на рис. 1-12.

На рис. 1-13 показаны зависимости максимальной отно­ сительной погрешности приближения и коэффициентов полиномов Рг(х) и Рг{х) от длины интервала аппрокси­

мации. В табл. 1-6 приведены данные для двух полино­ мов. По сравнению с интегральной аппроксимацией полиномом второй степени при |еМакс|=2% наблюдает­ ся расширение интервала в 1,29 :раза, а полиномом третьей степени при |еМакс| =3% — в 1,47 раза.

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

34

относительная погрешность определения ПХ равна:

 

 

т

 

 

 

т

 

т

 

 

 

 

 

2

АпеР" ‘

(!+ < /.)- 2

А^ Рп‘

2

АпеРп‘

 

 

^общ

п= 1

 

 

П= 1

 

П=1

 

 

 

 

 

т

 

 

 

т

 

 

 

 

 

 

 

2 'АпёV

 

 

2

л,*'"'

 

 

 

 

 

/ 2 = 1

 

 

 

/1=1

 

 

 

 

Если все

корни

близки

по величине, то

81= 82= ...

 

= еп и е0бщ =еп. Когда это условие не соблюдено, то

а12

а 22

 

 

6,%ajj

 

 

 

 

°73

а23

106

0,51

 

 

12

t

 

 

 

 

 

 

 

 

ПН

 

 

 

 

106 051

tooШ5

 

 

10 0)2)

 

 

 

 

1,00- 0/5

О,Oh

 

 

 

 

1

”^4.

/ J

 

0J9h- 0,33

 

 

 

V *

 

 

 

\ V

 

 

 

ом 0,33

 

 

 

 

 

\

/

 

0,33

 

 

 

 

 

 

 

 

 

 

\* 2 2 \J

 

 

/*33

 

 

10.82-0,27

ОМ•ом

 

 

h

0,06

 

/ \

\

 

 

 

 

0,76- 0,21

 

А

2 от

 

 

1 -L"

1

0.21

0,7010,15

 

4

 

 

 

 

^ 1

0,15

7

О

0,02

 

 

 

 

 

2t

7

 

 

 

hi

 

 

 

 

 

 

Рис.

1-13.

 

 

 

 

 

Еобщ может

отличаться

от еМакс

как в большую,

так и

в меньшую сторону. Это вызвано тем, что коэффициен­ ты Ап зависят от расположения всех полюсов и нулей и

могут значительно отличаться по величине, а также иметь разные знаки. Поэтому необходимым условием хорошей точности при использовании аппроксимирую­ щих формул является использование полиномов, обеспе­ чивающих близкое приближение к ПХ в возможно боль­ шем интервале времени. Однако этого недостаточно. Для оценки погрешности расчета можно после нахожде­ ния нормированной длительности импульса определить искажения вершины с помощью используемого полинома и полинома более высокого порядка. Если их разность превышает допустимую погрешность расчета 8ДОп, то

3*

3

длительность импульса следует уменьшить. Если раз­ ность меньше еДОш то полученный результат достаточно близок к точному, даже если для наибольшего по абсо­ лютной величине корня произведение \pn\tvi превышает

величину допустимого интервала аппроксимации.

В этом случае большая относительная погрешность

определения еРп и мало влияет

на общую погрешность

Л

.

из-за малости произведения Апе

б) Коррекция вершин одиночных прям оугольны х им пульсов

1. Коррекция одним элементом. Для синтеза параме­ тра коррекции аппроксимируем ПХ усилительного ка­ скада полиномом второй степени

т

от

 

 

Н (t) — 2 Ап6

^

Ап (1“|” а 12/?л^“|~ &ггр2п12)

П—\

Л=1

 

 

или, иначе

 

 

(1-23)

Н (/)« 1 +

“Ь a&fct2.

Точность, обеспечиваемая

формулой (-l-SS), вполне

достаточна

во м’нотих практических случаях. Так, например, в схеме коррекции вершин импульсов с помощью цепи Сф#ф в ламповом усилителе

 

Pi — о>\Ь\=

1

—т

 

 

 

т

 

 

 

 

 

 

.

, , ,

,

,

=

— b— m

рг =

Дг — Ьг b\\fl\ Ь\)

~ --------t

 

 

 

 

 

т2

где

Сф/?а

 

 

/?а

 

 

 

 

 

 

” = а д Г : ь = Лф

:

Я (/) ^

1 + 0.8833

1 +

0,2459

V. (1-24)

Результаты расчетов по формуле (il-24) достаточно точно совпа­ дают с графиками, приведенными в приложениях к монографиям [9, 10] и рассчитанными по точным зависимостям. Подобные же гра­ фики используются для расчета транзисторных усилителей (11].

Рассмотрим такой вариант коррекции, в котором па­ раметры каскада выбираются из условия компенсации линейной составляющей ПХ, и определим искажения, создаваемые квадратичной составляющей. В этом слу­ чае ПХ каскада имеет вид, изображенный на рис. 1-14.

36

 

Спад

вершины

импульса А~ определяется

выраже­

нием

 

т

 

 

 

 

 

А

1

Н (^н) ^

tl—- ^22^2^0»

(1“25)

 

 

 

Л=1

 

при этом необходимо выполнение условия

 

т

(1-26)

 

2 АпРп — pi — О*

 

п—1

 

На

основании формул (1-20) выражения

(1-25) и

(1-26)

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

 

A - = a 22t21i(b2—a2)\

a i—bi.

Равенство (1-28) дает уравнение величины корректирующего пара­ метра. Из формулы (1-27) по за­ данному спаду определяется нор­ мированная длительность импуль­ са /л= тп/т, где т — постоянная времени, по которой производи­ лось нормирование при нахожде­ нии изображения (1-17).

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

(1-27)

(1-28)

для определения

W )

7

Зг

tu t

SP\L-

А“ 1 = — агзМ 2и— аззРз/3и,

Рис. 1-14.

а абсолютная погрешность

А- —A“i = (СС23—Ct22)'P2^2H+Ct33p3^3u.

С учетом условия Pi=0 это дает:

А- —А- 1=0,1404 2b2) t2n +

+ 0,0617 [ 3—Ь3) bi (a2—b2) ] t\u

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

37

трению одного каскада. Такая задача значительно про­

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

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

Н (р)Н\ (р) Нг (р) ... Hs (р) =

___ {pm + a u p nt-

1 +

fl2ljPm ~ 2 +

••») { р т +

CL\^Pm ~ X +

 

 

{рт - \ - Ь п р т ~ х _j_ Ьглр™ - 2 +

...)

\ р т +

Ьхг р т ~ х +

 

 

- f а 2г р м - г +

... )

(р т +

g isp m - i -{ -a asjpM- 2 4 . ... )

 

 

Ьйърт ~ 2

... )

(р т -{-

b \s p tn~ x -|- bzsprn~ z -J- . . . )

*

 

Это изображение

может

быть

представлено

в виде

И ( п \ ~ — />sm +

g '* /?sm~ l +

a / 2l

s w - 2 +

... + a rsm

n

nqx

pS'n +

bripS"x- x+

b'zpS'n-2 +

t" + b'sm’

K

'

где

e'«= 2 a><-; 6' - = S 6ir;

r = 1

 

r ~ 1

s

s

s

Q 2 ==2 ^2r

S

2

r= 1

 

r = \ q=l

s

s

s

 

b's= 2 b*r+

2

2

r Ф я*

r ~ \

r= 1 <7= 1

 

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

ство

a'i=&'i

или, иначе,

S

S

 

2 a , r = 2 V

(1-30)

Г—1

г=1

 

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

а1г=&1Г. (1-31)

38

Тогда выражение (1-27), определяющее квадратич­ ную составляющую спада вершины, с учетом (1-31) зна­ чительно упрощается и принимает вид:

А

----- 0&22^2Ц 2 (р2Г -----

dzf)»

Г=I

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

А = Aj -J-Д2 —|—... —|—As .

(1-32)-

Если усилитель состоит из одинаковых каскадов, то

A“ = sAy\

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

Я (0)= #(г„).

(1-33)

 

 

Тогда ПХ приобретает форму,

 

 

показанную да рис. 1-15. При

 

 

этом такой способ коррекции по

 

 

сравнению

с рассмотренным вы­

 

 

ше позволяет при тех же допусти­

 

 

мых искажениях Д увеличить дли­

 

 

тельность

импульса или

умень­

 

 

шить низкочастотную постоянную

. Рис.

1-15..

времени усилителя.

 

 

 

Определим условия коррекции для этого случая. При

применении

квадратичной

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

вершина

импульса считается параболической с максимумом в мо­ мент времени t = t ul2 . Тогда условие (1-33) принимает

РИД:

тт

2

==2 ^п^ —1—

i а82Рал^*и)»

/1=1

/1=1

 

рткуда определяется длительность импульса

^

<*12р|_______ Oi2_____________til bj___________

(1-34)

**

a 22p8

&22 tti bz~—Ь\ {tl\ bi'j

 

39

Величина подъема вершины Д+ определяется следую­

щим образом:

т т

Д+ = Я Л у^ — Я (0) =

а12 -у*

Апрп 4 ”

Чг

AnPin

 

 

 

 

 

л=1

л=1

 

с учетом (1-34)

 

 

{ах — М2

 

 

 

 

 

ttgl2

 

 

 

(1-35)

 

Д+ = — 4а22

CLZ bz b\ (flj'— b i)

 

 

Формула (1-35) позволяет по заданным искажениям

определить

численное

значение параметра коррекции,

 

 

 

 

 

входящего

в коэффициен­

 

 

 

 

 

ты изображения

ПХ. То­

 

 

 

 

 

гда с помощью равенства

 

 

 

 

 

(1-34) находится норми­

 

 

 

 

 

рованная

длительность

 

 

 

 

 

импульса /и, для которого

 

 

 

 

 

выполняется

условие

 

 

 

 

 

(1-33).

 

 

 

 

 

 

 

 

При больших подъемах

 

 

 

 

 

вершины

или

глубокой

 

 

 

 

 

коррекции больших спадов

ся от параболической

 

 

форма импульса отличает­

(кривая 1 -на рис. 1-Ь6).. В этом слу­

чае величина

подъема

по-прежнему достаточно

точно

определяется

с помощью формул (1-27) или (1-35)

(кри­

вая 2 )< а

длительность

импульса'выходит ,за*'пределы

интервала

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

При этом формула.... (1-34)

дает всегда заниженное значение t u . Для

более-точного

определения t u запишем условие (1-33) с точностью до

членов третьего порядка:

тт

2 ^ п

2

4” ^13Рп^И~Н ®2зР2«^2И

 

й-зф*п№н)'з

п = \

Л=1

 

 

 

откуда

ai3Pi+ (Х23р2^и+ азэРз^2и= 0

 

 

 

 

 

И

 

 

 

 

^

__—агзРг i 1^а22зР22 — 4агзазз51Рг

(1-36)

 

 

 

 

Длительность

импульса определяет

меньший поло­

жительный корень уравнения (1-36).

 

 

40