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

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

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

ния хп через 'сумму более низких

степеней

оба автора

используют смещенные полиномы

Чебышева Т*п(х) =

= Т п(2х— 1), приведенные к виду т*я (х) =

Т*п(х)!2гп~х

в интервале [0, 1]. Как известно, в единичном интервале из всех полиномов п-то порядка со старшим коэффи­

циентом единица наименьшее возможное уклонение от нуля, равное 2~2n+1, имеют полиномы т*п(*). Пониже­ ние степени полинома при экономизации производит­ ся на основе приближен­ ного равенства

хпя*хп—т*п(я)

с наибольшей абсолютной

погрешностью

± 2 _2п+1.

Для

аппроксимации АЧХ

такая

замена

представ­

ляет определенное неудоб­ ство, так как т*п(0) фО и

после нескольких последовательных снижений степени л; накапливается некоторая постоянная составляющая и тогда М (0) Ф 1. Во избежание этого перенесем начало

координат и изменим масштаб так, чтобы новый поли­ ном Т°п(х) удовлетворял равенствам

Рп(0)=0; Р я(1)=7,*п(1).

Определим правила перехода от полинома Т*п(х*) к Т°п(хР). Нули несмещенного полинома Чебышева Тп(х) располагаются в точках

х , = cos ( 2 / - 1)'^-.

 

 

- шЛ

Поэтому нули смещенного полинома Т*п(х*) имеют

абсциссы

 

 

г*

* /+ ! _ *

+ C0S(2/- 1)^

Х } — ~О

о

С увеличением интервала растет уклонение полинома тп(х) от нуля и точность аппроксимации снижается.

Чтобы потеря точности была минимальной, перенесем Начало координат в ближайший нул^

у,

• а ТС

Л==53Ш w

71

в растянем масштаб оси абсцисс в

/

1

--------—

раз.

 

1 - **п

7Z

Г

 

 

cos2w

 

 

Тогда получим полином

 

 

р „(* 0) = т * п[х<>(\—х*п) + Х * Л

(2-23)

график которого показан на рис. 2-3. Нули полинома расположены в точках

п

COS (2/

1)

4" COS 2п

c o s /o T

cos ( / —

 

1) 2л"

x*i =

 

%

 

п

 

 

l,+C0S5F

 

COS2 Ап

 

(2-24)

 

 

 

 

 

 

а коэффициент при х п [’равен 2,л"1 cosartj

Наибольшее

уклонение от

нуля полинома

т:0л(л0) =

-----Т°п

- бу-

 

 

 

 

22Я-1 COS2rt Ап

дет

 

 

 

 

 

 

 

Дмакс = niz ---------------

 

(2-25)

По сравнению с полиномом т*п(х*) уклонение воз­ растает в cos-2n(я/4/г) раз. Однако ухудшение точности

.аппроксимации незначительно, так как полиномы при­ меняются при п ^ 4 . Например, для п= 4 уклонение со­

ставляет 0,00912 вместо 0,00781, т. е. возрастает на 17%. С увеличением п это различие сглаживается.

В дальнейшем в процессе экономизации будем исполь­ зовать полиномы Т°п(х°) и считать, что

Хп« Хп—Х°п(х) = Qn ) .

Коэффициенты q™ полиномов Qn (x) до п—9 и ма­

ксимальные уклонения Лпмакс приведены в табл. 2-3. Впервые метод экономизации для исследования широко­ полосных усилителей применен в ,[29].

Установим зависимости (2-20) между коэффициен­ тами числителя и знаменателя АЧХ, используя полино­ мы Qn(x)- В качестве примера, иллюстрирующего ме-

72

Р 4 4 ? 4 в ) А

Таблица 2-3

п

3

4

5

6

7

8

9

Я\п

—0,4308

0.1835

-0,0698

0,0247

-0.00828

0,00266

-0,00084

Я2п

1,3923

—0,1135

0,6781

—0,3496

0,16106

-0,06833

0,02724

ЯЗп

 

1,9209

—0,0477

1,6094

—1,04262

0,58990

—0,30184

Я4п

 

 

2,4373

—3,2330

3,10301

—2,41419

1,61706

Явп

 

 

 

2,9480

—4,66885

5,28380

-4,79266

ЯВп

 

 

 

 

3,45556

-6,35508

8.27703

Я7п

 

 

 

 

 

3,96119

—8,29160

Я8п

0,03848

0,00912

 

0.00054

0,00013

0,00003

4,46555

|длмакс|

0,00221

0,00001

тодику расчета и возможную область применения мето­ да, рассмотрим АЧХ системы, имеющей изображение с двумя нулями и четырьмя полюсами. В единичном; интервале ее можно представить формулой

М (х) =

\ / ______ 1-f- AiX + Ax2_____ #

 

 

1 + В ,х +

Bzx* + В 3х 3 + В4х*

 

На основании

выражения

(2-21) запишем

равен­

ство

 

 

 

(I+ A ix +Агх2) (1 +рцх+рих2+рз1х3+рглх1) =

(2-26)

= l + B lx + B 2Xi + B 3x^+B ixi.

Множитель х6 в левой части заменим полиномом

 

6

 

 

 

<з,

qi'x ‘'

 

/=i

а *5 в полученном выражении — полиномом

5

q i % x i *

/=i

Затем приравняем коэффициенты при одинаковых «степенях в левой и правой частях полученного прибли­ женного равенства. Тогда получим систему уравнений

* + 0 +^44^1*) А + (Р34?15+

Р**Я*ъ9\Ь”}“P44<7le) А

(2-27)

Р*4 ~|“ (р*4 “{-Р**Я*й) А ”1" (Р*4

-J- р » *Я + РиЯмЯм + = В*.

73

Выражения (2-27) связывают коэффициенты Лг* и £* Простой линейной зависимостью. Однако при рассмотре­ нии АЧХ в качестве обобщённых параметров целесооб­ разнее использовать коэффициенты изображения ПХ. Во-первых, при использовании зависимостей между ко­ эффициентами числителя и знаменателя изображения ПХ полностью унифицируется методика расчета пара­ метров конкретных схем при исходных требованиях, заданных во временной или частотной областях. Во-вто­ рых, коэффициенты изображения проще выражаются через физические параметры схемы, чем коэффициенты АЧХ. Поэтому заменим коэффициенты Лг- и B i их вы­

ражениями через коэффициенты дважды нормированно­ го изображения:

A i= (g2i—2gz)q; A2—g\q2\ B i= \ d h —2d2)q\

B2=•'(dS,—2did$H* 1)q2] £ 3— (d2s2d2) p3; В 4='p4.

Тогда последнее уравнение системы (2-27) будет со­ держать только одно неизвестное q\

qkg\ (р24 р34р45+ р44^56^45+ Р44^4б) р2—

— (Рз4+Р44р4б) (gh—2g2)q—Р44-О ,

если заданы коэффициенты gi и g2. Остальные уравне­

ния (2-27) образуют систему

d‘, — 2 d ,= — ;

Ч

d*s — 2 d ,d ,= ~ — 1;

d U - Ы г — Цг

q*

С тремя неизвестными di, d2, d2. Эта система легко ре­

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

ле

низких степеней

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

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

di.

В

зависимости от

величины коэффициентов

qa,

рц>

а также выбранных значений gi и

g2 могут

существо­

вать несколько решений. Из них отбираются те, которые удовлетворяют условиям устойчивости [30]. Полученные результаты приведены в табл. 2-5 в виде функциональ­

74

ных зависимостей d{—fi(gu £2) при АЛ1= ±0,01, а так­

же для оптимальных монотонных АЧХ.

При больших допустимых отклонениях АЧХ метод экономизации вносит существенную погрешность. На рис. 2-4 показана АЧХ системы четвертого порядка для ДМ= ±0,05 . Исходная АЧХ при gi = 0 имеет требуемые отклонения от единичного уровня. С увеличением gi

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

личины и при g i= l,5 достигает MMai<c=U. В случае ДМ>0,05 погрешность экономизации еще больше. По­ этому коэффициенты di для ДМ—±0,05 и ±0,1, приве­

денные в той же таблице, получены численными мето­ дами, изложенными в § 2-3 и свободными от указанного недостатка.

б) Подвижный интервал аппроксимации

В практике исследования и расчета широкополосных усилителей часто встречается необходимость рассмо­ треть небольшое число АЧХ, выраженных полиномами высокого порядка. Определение экстремумов АЧХ и граничной частоты при этом затруднительно, а исполь­ зование ЭВМ нецелесообразно из-за малого количества вариантов. В таком случае весьма эффективно приме­ нение подвижного интервала аппроксимации.

Пусть полином п-го порядка Р п{х)

со

старшим

ко­

эффициентом единица уклоняется от нуля

на |ДПмакс|

в единичном интервале. Переведем его

в

интервал

[а,

75

b] с

помощью

замены

переменных

*i= а + /х , где /==

= 6 —а. Уклонение при этом не

изменяется,

а коэффи­

циент

при xni

станет

равным

1Jln.

Чтобы

получить

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

максимальное уклонение изменяется во столько же раз. Поэтому полином nn(X i)= lnP n (Xi) уклоняется от нуля

на величину

| Ал макс | fа,Ь] 1п | А п макс

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

Сужать интервал аппроксимации, естественно, мож­ но лишь в том случае, когда хотя бы грубо приблизи­ тельно известно положение искомой экстремальной точ­ ки или граничной частоты. Известно, что граничные ча­ стоты по уровню 0,7 систем 3—6-го порядков, имеющих АЧХ вида

1 + 2 Dix‘,

укладываются в интервал [1, 5; 2, 5]. Поэтому если по­ рядок системы определен, то положение экстремумов и граничная частота известны с точностью, достаточной для применения подвижного интервала аппроксимации. Тогда в области исследуемого экстремума АЧХ аппро­ ксимируется выражением

а в области граничной частоты (рис. 2-5) выражением

+ G",*,

о +

76

из которых легко определяются значения xig, х1в и Мд.

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

менения

масштаба. Коэффициенты

G'0, D'Q и величина

М (а) не равны единице. По­

 

этому

требование,

чтобы

 

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

полино­

 

мы вначале

интервала

об­

 

ращались в нуль, не являет­

 

ся обязательным.

 

В даль­

 

нейшем воспользуемся поли­

 

номами

т*П'(я),

обеспечи­

 

вающими наименьшую

по­

 

грешность.

 

использо­

 

Правильность

Рис. 2-5.

вания метода

экономизации

 

проверяется попаданием рассчитанной точки в интервал аппроксимации. Наибольшая погрешность при экономи­

зации возникает, когда

x3i заменяется квадратичным

полиномом, a x2i — линейным. Рассмотрим

применение

такой методики к АЧХ вида

 

 

 

м (х,) =

г

 

1-f- G\X\

(2-28)

D

Z)2#2i~b X3I

 

 

 

для определения х^, Ммй1 С и Хщ. Сначала определим

интервал, в котором возможен подъем АЧХ. На верхней и нижней границах этого интервала соблюдается ра­ венство

_____1 -f- GjXi________i Q 1-}- D\%\-f- Z?2X2i-p x»i

которое приводит к уравнению

(Gi—Di) Xi—D2x2i—x3i= 0.

Отсюда

_

n

- O ,± /Д « ,- 4 (О,- О , Г .

Лц

\J9 *^12,11 —■

 

2

Точки

с ординатами

A f(*i)>l существуют, если

 

D h >4(D i— Gi).

 

Амплитудно-частотная характеристика

имеет мини­

мум и максимум, когда Xi2> 0

и лг1з> 0, что выполняется

при условиях D i> G i и D2< 0.

Тогда за

приближенное

77

значение Хъ, при котором Af=MMaKC, можно принять ве­

личину

^

#12 + Хи

^2

 

Хгэ

2

1Г*

 

Пусть задано Gi = 2 ,5 ,

Di = '2 ,9 8 4 ,

Dz==—2,046. Тогда Xia—’1,008.

Для уточнения применим

метод экономизации

выражения (2-26)

в интервале [0,5; 1,5]. Полином т*з’(х)

переведем

в этот интервал и

воспользуемся выражением

 

 

 

х3~х3—

0,5) = 'Зх2—2 ,8 1'25х+0,81(25.

Тогда АЧХ принимает вид:

М {х г) Y

1+ GiXi

1,8125 + [Di*— 2,8125} Xj + (D2 + 3) х»,

12,5xi

-Г - 1,8125-f 0,1715x,+ 0,984x2,

Отсюда из условия Af,*(xio) = 0 получаем Xi3=0,99 и M23=ll,'179. Расчет без применения экономизации требует решения кубического уравнения и дает Xi3=0,994, Л423=И,4794. Для определения XiD пе­ реведем полиномы т*з(х) и x*2'(^) в интервал ['2, 3]. Тогда в этом интервале получаем:

лг3^7,5л:2—18,5625а:+ 15,1563;

х2~5х —6,125;

,,, ч т /

1+ 2,5xi

M {x i)sssy

— i7 ,433 -f- П ,842xi

Отсюда при М(х\п) = 0 ,7 0 7 получаем X m = 2 ,8 4 . Точный расчет требует решения не линейного, а кубического уравнения и дает хщ = = '2 ,8 4 5 . Погрешность при определении xi3 и Хщ очень мала, поэтому подвижный интервал аппроксимации удобно использовать для на­ хождения х\ъ при изменениях параметров схемы, .нахождения коэф­ фициентов АЧХ, при которых М0 имеет заданную величину, и т. д.

2-3. СИНТЕЗ АМПЛИТУДНО-ЧАСТОТНЫХ ХАРАКТЕРИСТИК С ЗАДАННЫМИ ЭКСТРЕМУМАМИ

Как уже отмечалось, при допустимых отклонениях АЧХ от единичного уровня на величину ДМ^О,05 не удается получить равноволновую характеристику с по­ мощью метода экономизации. Нахождение равноволно­ вой АЧХ является задачей, родственной задаче чебышевского приближения, хотя не тождественной ей. При нахождении чебышевского приближения минимизирует­ ся максимальное отклонение, в данном же случае это

78

бТкЛонёние фиксируется. Отыскиваются коэффициенты дробно-рациональной функции, позволяющие получить заданные максимальные отклонения внутри выбранного интервала и удовлетворяющие условиям физической реализуемости. При этом число точек, образующих альтернанс1, не превышает п+ 1 независимо от степени

полинома числителя, так как отклонение

пт

2 B/xi - 2 Aixi

1 - Af* { х ) = ^ ---------------------

1 + 2 B,xi

/=1

не может иметь более п нулей внутри конечного интер­

вала. Кроме того, АЧХ .принято нормировать так, чтобы выполнялось условие М (0) = 1. В конце интервала исхо­

дя из допустимых искажений также задается значение ординаты АЧХ. Для упрощения расчета в дальнейшем принято условие М(1) = 1. Это не отражается на конеч­ ных результатах, так как после получения АЧХ с задан­ ными отклонениями независимо от первоначальной нор­ мировки она представляется в виде

м (х.)=

у

1 -f- G 1X 1 -{- G2Xsi - j - . . . -j- GmXmi

1 -J-D\%\ -j- DzX^i

... -f* xni

 

 

Указанные условия накладывают два дополнитель­ ных требования на форму АЧХ, поэтому альтернанс можно получить только в п—1 точке.

Итак, необходимо решить следующую задачу: по за­ данным значениям коэффициентов Агнайдем коэффи­ циенты £ 3-, обеспечивающие в единичном интервале п—1

экстремум функции

т

1 + 2 Aix‘

М (X) = — ^ ------.

(2-29)

1 + 2 B/xi

/=1

1 Альтернансом называется совокупность точек, в которых раз­ ность между аппроксимируемой /*(*) и аппроксимирующей f(x) функцией A= f ( x ) —*P{x) принимает равные по абсолютной величине Амако и чередующиеся по знаку экстремальные значения.

79

равный l+AM/t. Для этого воспользуемся следующим алгоритмом.

1. Задается начальное приближение коэффициентов обеспечивающих существование п— 1 экстре­

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

2. Вычисляется п—1 точка экстремума. Проверяет­

ся, равны ли экстремальные отклонения заданным величинам AMh (k=\, ..., п— 1). Если равенство выпол­

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

3. Вычисляются -новые значения коэффициентов BWj, при которых аппроксимирующая кривая проходит

через точки с абсциссами, равными абсциссам найден­ ных экстремумов, и ординатами, равными 1 +A M h. Ко­ эффициенты находятся из системы, содержащей п ли­

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

 

т

 

 

1 + 2

Л/А

 

Ч ------- =

(1+ДЛЪ)»;

1 + 2

в/*1*

 

 

/=1

 

til

 

п

 

2

^

==2 ^

•••>п ~~ i)«

/=1

 

/=i

 

4. Повторяется второй этап.

Итерационный процесс продолжается до тех пор, по­ ка не будет получена АЧХ с заданными отклонениями. В качестве начального приближения используются ко­ эффициенты Bj, найденные аналитически в § 2-2,в при

А{—0.

Изложенный алгоритм представляет собой модифи­ кацию второго алгоритма Е. Я. Ремеза [7] и относится к той же группе модификаций, что и алгоритм А. А. Ланнэ [31], отличаясь от него классом аппроксимирующих функций и требованиями к форме характеристик. Алго­ ритм очень эффективен: подтягивание происходит за 2—4 итерации в зависимости от начального приближе­ ния и порядка системы.

80