
книги из ГПНТБ / Смолов, В. Б. Аналого-цифровые и цифро-аналоговые нелинейные вычислительные устройства
.pdfэто обстоятельство учитывается подстановкой суммарных коэф-
/ |
/ |
i |
i |
фициентов |
2 Ж , |
^ \d 2s и |
^ d 3s.' |
S=1 |
S=1 |
S=1 |
S=1 |
Степень аппроксимирующего полинома может быть повышена за счет использования в цепях обратных связей операционных усилителей линей ных и квадратичных цифровых управляемых сопротивлений (0 — ЦУС).
На рис. 6-7 изображена схема полиномиального аппроксиматора, позво-
4
ляющая получить полином 4-й степени и тлх— ^ А ф 11 на одном операцион- fc=o
ном усилителе.
Токовый квадратичный элемент ТКЭ, включенный во входную цепь
операционного усилителя ОУ, обеспечивает изменение тока, выходящего из
узла а по зависимости вида / а = £/0 (а0 + афЗ + а202)- Сопротивление Rao
цепи обратной связи ОУ изменяется в зависимости от кода. 0 по закону квад ратичного полинома:
axfa + 6) (g2+ |
1~ ° ) |
“ bx6 4 |
“ 6202), |
(6-26) |
Ваб : |
— ^max (2>o 4 |
|||
■Rmax (Ci + 1 + |
сг) |
|
|
|
где 60, bly b2 — постоянные коэффициенты, определяемые величинами от носительных сопротивлений сх и с2:
ь0= Cl(1 + |
Са) -, ьх= |
, |
I |
1 |
(6-27) |
||
1+ С1 4" с2 |
|||||||
1 4 - |
ci |
4 " са |
1 4~ ci 4* с2 |
|
|
||
Таким образом, выходное напряжение схемы £/ВЬ1Х изменяется по зави- |
|||||||
симости вида |
|
|
|
|
|
|
|
t /в ы х = |
- |
1аЯаб = и о {А0 4 - А 10 4- А 202 4 - А а03 4 - А 40 * ), |
(6-28) |
где постоянные коэффициенты Aj (j = 0,4) определяются параметрами ре
зисторов входной цепи и цепи обратной связи ОУ.
6-2. Дробно-рациональные цифро-аналоговые аппроксиматоры
В тех случаях когда полиномиальная аппроксимация не обеспе чивает заданных характеристик цифро-аналогового функциональ ного преобразования (точность, габаритно-весовые показатели, стоимость, надежность и т. д.), следует проанализировать возмож
ности применения дробно-рациональной аппроксимации |
(ДРА), |
при которой функция F (X ) заменяется рациональной |
дробью |
Пп
Q (Х )= 2 AkX k! 2 BkX k, так чтобы коэффициенты Ak, Bk послед- fe=i k=\
ней удовлетворяли заданной ошибке на всем диапазоне изменения аргумента [Хнач, Х кон]:
2 |
Ak* k |
(6-29) |
em (X) = min шах * ( Х ) ~ |
-------- |
|
S |
BkXk |
|
fc=l |
|
150
Кривая |
ошибки |
em (X) |
для |
рассматриваемого |
наилучшего |
|||||
в смысле Чебышева приближения функций F (X) и Q (X) согласно |
||||||||||
теореме |
Ахиозера |
[2] |
имеет |
на |
интервале |
[Хнач, Хкон] не менее |
||||
т + п + |
2 |
альтернирующих |
экстремумов |
в |
критических |
точках |
||||
|
|
■Хнач |
^0 |
|
|
Хт+п+1 |
Х кою |
(6-30) |
||
где т и п — степени полиномов числителя |
|
т |
|
|||||||
Р х (Х )= 2 |
АtXl и зна- |
|||||||||
|
|
|
|
П |
|
|
|
г=о |
|
|
|
|
|
|
|
|
|
|
|
|
|
менателя |
Р а (Х) = |
1-f- |
|
|
|
|
|
|
||
|
|
|
|
/==i |
|
|
|
|
и В/ |
состав |
Для нахождения неизвестных коэффициентов At |
||||||||||
ляется нелинейная система уравнений |
|
|
|
|
||||||
|
|
ет (X/) = extr, |
|
|
| |
|
|
|||
|
|
F (X ,) - Q (X /) = ( - l ) ' E , |
|
|
|
(6-31) |
/ = 0, 1, 2, . . . , tn-j- я -f-1, j
где Е — амплитуда ошибки ет (X,).
Эта система не имеет единственного решения, но проектиров щиков интересует решение, не имеющее полюсов на интервале [Хнач, Хкон]. Такому решению наиболее часто соответствует наи меньшее из всех возможных по абсолютной величине значение Е, поэтому ноль является хорошим начальным приближением. Иссле дованию сходимости при различных вычислительных схемах ал горитмов решения системы (6-31) посвящены многие математиче ские работы [2, 36], в которых алгоритм сводится к двухступенча тому итерационному процессу. На первом этапе («внутренняя ите рация») методом итерации по известным приближениям для крити ческих точек решается система (6-31). На втором этапе находятся экстремумы г1к кривой ошибки em(X) = F (X) — Q (X), прини маемые за новые приближения для критических точек.
Вработе [36 ] доказано, что сходимость итерационного процесса определяется удачным выбором начальных приближений в крити ческих точках.
Вприложении к работе [80] приводится программа на языке АЛГОЛ-60 дляЦВМ«Минск-22»,составленная поалгоритму
Фрэзера |
и Харта. Итерационная схема |
в этом случае |
имеет вид |
|||||
т |
п ' |
■* |
|
|
|
|
|
|
ЦА1+1Х]к + 2J [(— l / |
Euw (Xjk)— F (X//()] X%Brp" + |
|
|
|||||
t'=0 |
i=1 |
|
|
|
|
|
|
|
+ (— 1 ),t£» (Xlk) Erkv1= F (X/k),/ = 0 , 1 , . |
|
. |
. ,m + n + l - |
(6-32) |
||||
В уравнении (6-32) |
k — номер внешней |
итерации, |
г — номер |
|||||
внутренней итерации, w (X) — весовая |
функция. Если |
амплитуда |
||||||
ошибки Ek известна, то на шаге г + |
1 |
(6-32) |
становится |
линейной |
||||
системой |
относительно неизвестных А;+1, |
ВрИ, Егк+\ Внутренняя |
||||||
итерация |
продолжается до тех пор, |
пока |
значения Ek |
и |
Ek+l не |
151
Таблица 6-4
Дробно-рациональные аппроксимации тригонометрических функций
Вид аппроксимации
|
. |
я |
|
п |
|
1,7970 |
|
|
tg — |
|
9 я ----------------- |
1,147 — 02 |
|||
|
|
2 |
|
|
|
||
. |
п |
|
|
|
|
|
0,639 0 |
tg я0 к 0,59Л -\----------------- |
О,25 — 02 |
||||||
|
|
|
|
|
|
|
|
|
. |
2 я |
|
„ |
а |
|
1,28 0 |
|
t g ----- |
3 |
|
0 |
0 , § 2 — 02 |
||
|
|
|
|
|
|||
|
|
, |
па |
« |
1,93 0 |
||
|
arc tg |
2 0 |
0,98 + 02 |
||||
|
|
|
|
|
|
||
|
. |
я |
„ |
|
4 0 ( 2 - 0 ) |
||
|
|
2 |
|
|
|
5 — 29 + 02 |
|
. |
„ |
ж — 6 |
,а0 0 н |
17,180 |
|||
sin Л0 |
1,82 + 02 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
„ |
ж |
1 - |
402 |
|
|
COS JT0 |
1 + 0 2 |
||||
|
|
|
|
|
|
||
|
. |
лГ— А |
|
3,28 0 |
|||
arc sin у |
|
2 |
0 ж --------------- |
2,34 — 02 |
|||
|
|
|
|
|
|
|
|
arc sin 20 « |
1,97 0 + |
0,284 — 02 |
|||||
|
|
|
|
|
|
|
|
|
|
|
я |
|
. |
|
2,18 |
|
sec — |
0 ж |
---------------2,184 — 02 |
||||
|
|
|
3 |
|
|
Пределы
П • Д min’ °max
о сл |
о сл |
I |
|
—049; 0,49
1 О сл |
о сл |
о1 сл |
о сл |
0; |
1 |
—0,5; |
0,5 |
1 О СЛ |
о сл |
о сл |
о сл |
—0,5; |
0,5 |
— 1 ,0 ; |
1 , 0 |
Относительная
приведенная ошибка б, %
0,08
0 , 2
0 ,2 2
0,33
0,17
0,06
0,17
0, 2
1 , 1
1, 0
t g J L 0 s |
>’ 792 0 |
—0,9; 0,9 |
0,25 |
3 |
3,792— 02 |
|
|
152
будут совпадать с заданной точностью. Как было отмечено ранее,
удобно принять Е0 = 0.
В табл. 6-4 приведены аппроксимации ряда элементарных функ ций F (X), полученные по алгоритму Фрезера и Харта.
Рациональные приближающие дроби Q (X) могут быть обра зованы прежде всего при помощи пассивных схем (рис. 6-8, а), содержащих по аналогии с известными потенциометрическими [55] или широтно-импульсными [61] схемами цифровые линейные
проводимости Yq = Kmax0 |
(сопротивления |
RQ= i?max0) и п0_ |
стоянные проводимости Yj |
(сопротивления |
Rj). Выходное напря- |
Рис. 6-8 . Пассивные схемы ДРА
жение и вых пассивных схем ДРА изменяется |
в общем случае по |
|||
зависимости |
|
|
|
|
Uвых = и 0 |
Ао -f- A -f- А 262 |
(6-33) |
||
B0+ B XQ+ BJP |
||||
|
|
’ |
||
где |
|
|
|
|
Ak — Gk{d1. |
d2, d3, blt b2, ba, bt), |
|
||
Bk = Lk(bi, |
b2t b3, fr4), |
(6-34) |
||
|
|
|
||
Uо |
|
tfmax= r 0(2 " - l) - |
||
^max |
|
Приближающая рациональная дробь, воспроизводимая пассив ными схемами, имеет взаимосвязанные коэффициенты Ак, Вк, и поэ'тому возможности ее применения весьма ограничены. Так, на-
153
пример, при |
Ьг = |
оо, Ь2 = 0, b3 = 1, |
d3 = — 1, d4 = 0 и Ь4 + 0 |
(схема рис. |
6-8, б) |
рациональная дробь |
имеет вид |
k'm,!.ч —0,56V 2&,° |
|
(6-35) |
|||
|
1+ 2 6 4— 02 |
|
|||
при й4 = 0,896 и [ 0 |< ; 1 |
с приведенной |
ошибкой не более 0,25% |
|||
приближающийся к функции тангенса: |
|
|
|
||
Uвых= 0,5U0 |
1,7920 |
0,6061 - ^ - tg 0 . |
(6-36) |
||
3,792 — 62 |
|||||
|
|
2 |
|
Схема 6-8, в пригодна для воспроизведения секансной функции:
^вых |
2,18 |
(6-37) |
sec 0 |
||
|
2.184 — 02 |
’ |
|
Рис. 6-9. Схема активного ДРА |
|
|
|
|
||||
так как при Ьг = 0,592, |
Ъг = |
0, Ь3 = 1, 64 = о о , |
= |
0, |
d a = |
1, |
|||
d3 = 0 рациональная дробь (6-37) |
обеспечивает |
ошибку |
прибли |
||||||
жения не более 0,7% для |0 |< ;.1 . |
d3 = 0, |
Ьг — с о , |
b2 = 0,844, |
||||||
Наконец, при dx = |
0, d2 — 0, |
||||||||
Ь3 = со, |
bi — 0,067 схема рис. 6-8, г моделирует зависимость |
|
|||||||
|
UВЬ1Х= и й л [ |
1,791 |
« и о1,1067 - |
9 |
|
(6-38) |
|||
|
вых |
V |
и 1 , 9 1 1 - 0 |
|
V |
' |
|||
с ошибкой не более 1% |
для 0,04 < 0 < 1 . |
|
|
|
|
|
|||
Во всех рассмотренных схемах 0 |
= RJRm&-^ т. е- используются |
||||||||
цифровые |
последовательные |
управляемые |
резисторы |
А+ЦУС), |
однако подобные схемы легко выполняются с применением цифро вых параллельных управляемых резисторов ( 5Д-ЦУП). Более широкие возможности с точки зрения класса воспроизводимых за висимостей имеют активные схемы ДРА.
Одна из возможных активных схем, изображенная на рис. 6-9,
содержит во входной цепи и |
цепи обратной связи |
ОУ сложные |
|
цифро-управляемые' цепи из |
У^-ЦУП (А+ЦУС) |
и |
постоянных |
резисторов (У*, Rk). Если |
dl k = Y 1A/Y mах, |
d2k ^ . Y J Yтах, |
154
О |
= |
Y J Y mах, |
то сложные проводимости |
Ух (х) и У2 (х) равны со |
|||||||||||
ответственно: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
d u d i 2 |
(1 |
4 ~ d 1 з 4 ~ d 14) -f- 0 |
[с/ц (1 -f- d 1 3 |
-)- d 14) — |
d u d n ] |
^ ц 02 __ |
|||||||
|
|
du (dn + di2) (1 + |
d13) -f- 0 [rfM(1 + |
di3) — du (du 4- ^12)] dH02 |
|||||||||||
|
|
|
|
|
|
|
a io 4 ~ Д ц 9 4~ fli 202 |
|
( 6 - 3 9 ) |
||||||
|
|
|
|
|
|
|
1>ю 4 ~ & ц 0 4~ ^ 2i 0a |
|
|||||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
__ d 2 1 d 2 2 |
(1 |
4~ d ~23 -f- d.M) -f- Qdg| (1 -f- d 23 4 ~ d 2i ■— dm) — |
^216" |
|
|
||||||||
|
|
d%i (d2 1 + |
d 22) (1 + d 23) -(- Qd2i (l -|- d23— |
d2L— |
d22) |
d2$ |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
% |
|
|
|
|
|
|
|
|
|
|
|
^20 4 ~ fl210 "I” ^220^ |
|
(6-40) |
||||||
|
|
|
|
|
|
|
^20 4- b2lQ4" b2S i |
|
|
||||||
|
Так как для идеализированного ОУ (Ки — |
|
|
I' = |
0» U’ = 0) |
||||||||||
U = |
/ а. то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U. |
|
- и пYi (*) |
У |
(«10 4 ~ а ц 0 4 ~ а 1202) |
(^20 4~ ^210 4~ ^ 220г) |
_ |
|||||||||
|
|
Y*(*) |
|
(^10 4- Ьц& 4” |
^ 1202) |
( а 2 0 + |
а 210 4 “ а 2202) |
|
|||||||
|
|
|
|
|
|
|
= |
- |
t / 0- i^ ------. |
|
(6-41) |
||||
|
|
|
|
|
|
|
|
|
|
|
2 |
B/0/ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
;=о |
|
|
|
|
|
Аналогичная |
|
зависимость |
имеет |
место |
для |
случая |
dlk = |
|||||||
= |
RlJRlmax, |
d2k = R2k/R2mx, |
6 = R J R max И |
СОВМеСТНОГО |
В К Л Ю - |
||||||||||
чения |
цифровых |
управляемых |
проводимостей |
и |
сопротивлений |
в цепь обратной связи. Полученная рациональная дробь четвертой степени (т = 4, п = '4) имеет более широкие возможности для аппроксимации элементарных функций, однако, как и ранее, ко эффициенты А/, B j являются взаимосвязанными. Для обеспечения взаимонезависимых коэффициентов А/, В/ любого знака исполь зуется универсальная схема ДРА [60] на рис. 6-10, содержащая три операционных усилителя ОУ1, ОУ2, ОУЗ и два токовых квад ратичных элемента ТКЭ1, ТКЭ2\ напряжения f/x—Ue задаются при помощи линейных потенциометров ЛП1—ЛП6. Сопротивления Д х—Re служат для задания расчетных значений коэффициентов рациональной дроби. Очевидно, что токи / х и / 2, проходящие по квадратичным токовым цепям, равны:
/ 1 = |
^ A fi> = UsYa |
U\Yx -\- U2 У2 |
(6-42) |
(^ ш а х -■Yx), |
|||
|
/=о |
У 2 4 " У ш ах |
|
/2 = |
l , B f l l = UeY e |
{/4У4+ £ /5К5 |
(6-43) |
(^шах — ■Yx). |
|||
|
/=о |
У ь 4 " У ш ах |
|
155
Решая систему уравнений |
|
. /х = / 2, |
|
d/ = 4 r > |
• / = 1> 2> 3’ |
сг = - ^ - , |
г = 4, 5, 6. |
С^ВЫХ |
|
относительно и вих с учетом |
соотношения |
2 ^ -0 / ^вых — Uо /=0
2 * ^
/= 0
где
— dscs |
^2^2 |
Во — d6ce - |
&ьсь |
||
1 + с2 |
1 + съ |
||||
д^ _ rfi + d2c2 |
_ |
dj —dbcb |
|||
1 + сг |
1 |
1+ |
cs |
||
л 2= |
|
B,= |
|
d4 |
|
|
|
|
|
||
( + C2 |
|
1 |
+ |
c 5 |
|
Уi |
/ = 2, |
3, 4, |
5, |
6. |
|
|
p,
(6-44)
Утах0, получают
(6-45)
(6-46)
156
Расширение функциональных возможностей схемы рис. 6-10 осуществляется путем введения в схему пассивных цифро-анало говых делителей напряжения и суммирующих ОУ.
Схеме рис. 6-10 присуще свойство многофункциональности, заключаю щееся в возможности получения различных аппроксимирующих зависимо стей РВых = UeF (0) за счет автоматического включения соответствующих
номиналов постоянных проводимостей без изменения структурных связей блоков.
В качестве примеров ниже приводятся схемы некоторых тригонометри ческих дробно-рациональных преобразователей, использующих аппрокси мирующие зависимости, приведенные в табл. 6-4.
Пример 6-1. Тангенсный множительный цифро-аналоговый преобразо ватель, реализующий зависимость UBHX — £/Bxtg (л0/2) путем ее аппрокси
мации рациональной дробью
1,7949
t/вых — Uв 1,147 — 02 (6-47)
с ошибкой
|
|
|
|
1,7940 |
|
|
|
|
|
|
^шах — |
1.147 — 02 |
|
|
|
||
|
|
tg (я0/2) |
< 0,0008 |
|
|
|||
|
|
|
|
|
|
|
|
|
при |
| 0 | <0,5, изображен на рис. 6-11. |
|
FBX = |
1,794 Y max и |
||||
|
Схема |
имеет постоянную входную проводимость |
||||||
выходную |
7 Вых ~ |
аппроксимирующая |
рациональная дробь |
(6-47) яв |
||||
ляется корнем неявной функции |
|
|
|
|
||||
|
|
U ' = 1,7941/вх0 + |
О,1471/Вых + |
( 1 - 0 2) ^ в ы х « О , |
|
(6-48) |
||
описывающей равновесное состояние схемы |
рис. 6-11 |
при U’ = |
0. |
|||||
|
Пример 6-2. Тангенсный множительный цифро-аналоговый преобразо |
|||||||
ватель для реализации функции UBых= Ubx tg (2л0/3) при | 0 | |
< |
0,5 с ошиб |
||||||
кой |
бгаах < 0,214% |
использует |
аппроксимирующую |
рациональную дробь |
||||
вида |
|
|
1,289 |
|
|
|
||
|
|
|
|
|
|
(6-49) |
||
|
|
|
^ВЫХ = [/„ |
|
|
|
||
|
|
|
|
0,6 2 — 0 2 ’ |
|
|
|
157
которая является корнем неявной функции |
|
U ' = - 1,28(Увх0 - 0,38t/Bblx + (1 - 0=) UBbIX и 0, |
(6-50) |
описывающей равновесное состояние схемы рис. 6-1 2 .
Схема содержит три операционных усилителя и имеет постоянные вход-, ную и выходную проводимости.
Рис. 6-12. Множительный |
цифро-аналоговый |
преобразователь для |
функции |
2д |
|
и вы х= fAixtg----- |
0 |
|
|
3 |
|
Пример 6-3. Тангенсный множительный цифро-аналоговый преобразо ватель, воспроизводящий функцию
UBm- = U BXtg(nQ/2) |
(6-51) |
Рис. 6-13. Схема для воспроизведения функции тангенса
при помощи аппроксимирующей зависимости
|
^вых — Ug- |
0(2,8 —( |
|
(6-52) |
||
|
|
"(1,8 + 0) ( 1 - 0 ) ’ |
|
|
||
являющейся корнем неявной функции |
|
|
|
|||
f l ' |
_ п ^>8У тах0 |
, |
,, |
1 ,8 ^ т а х ( 1 — 9) |
(6-53) |
|
и |
— о'вх — |
+ |
ь 'в ы х |
-------------- *--------- |
0 |
|
|
1,8 + 0 |
|
|
2,8 — |
|
описывающей равновесное состояние схемы рис. 6-13.
158
При | 0 | < 0,95 ошибка аппроксимации ие превышает 1%, схема имеет
переменную входную проводимость (0 < FBX < 1,8 |
Hmax : 2,8) и постоянную |
|
выходную ( Квь,х ~ °°)- Изменение начала отсчета |
цифрового аргумента на |
|
я/2 означает в рассмотренной схеме подачу кода 0 = |
1—0 во входную цепь |
|
ОУ и кода 0 — в цепь его обратной связи. Поэтому |
рассмотренные схемы |
|
пригодны для реализации котангенсной функции. |
|
|
Пример 6-4. Синусный множительный цифро-аналоговый преобразова |
||
тель, реализующий приближение к зависимости |
|
|
Uвых— — 0.25С/вх sin |
Я0 |
(6-54) |
|
рациональной дробью |
0) |
0 (2 - |
|
^вых — — |
(6-55) |
5 — 20 — 02
которая при 0 < 0 < 1 обеспечивает ошибку 6тах < 0,17%.
Схема, воспроизводящая рациональную дробь (6-55) в виде корня неяв ной функции
U' = U. |
0(2 — 0) |
0 (2 - 0) |
0, |
(6-56) |
+ ^в |
|
изображена на рис. 6-15.
Пример 6-5. Арксинусный множительный цифро-аналоговый преобразо
ватель, имеющий характеристику |
вида |
|
|
|
|
3,28 — 0 |
(6-57) |
|
|
2,34 — 03 |
|
|
|
’ |
|
изображен на рис. 6-15. |
|
|
|
При | 0 | < 0,5 |
зависимость (6-57) приближается к функции |
||
|
Uвых = |
£/вх arc sinyT O |
(6-58) |
С ошибкой 6п1ах < |
0.2%. |
|
|
159