книги / Оптимизация параметров импульсных и широкополосных усилителей
..pdfотносительной погрешности в заданном интервале или же максимального интервала изменений аргумента, в котором погрешность не превышает допустимую. Мож но существенно расширить интервал аппроксимации без увеличения погрешности, если воспользоваться инте гральным методом приближения функций [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,316х2 |
(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
A°SP\L-
А“ 1 = — агзМ 2и— аззРз/3и, |
Рис. 1-14. |
а абсолютная погрешность
А- —A“i = (СС23—Ct22)'P2^2H+Ct33p3^3u.
С учетом условия Pi=0 это дает:
А- —А- 1=0,1404 (а2—b2) 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