книги / Методы электрических измерений
..pdfДискретные и цифровые фильтры принято делить на два класса: нерекурсивные (НФ) и рекурсивные (РФ). Если в выражении (5.38) все коэффициенты at — 0, то получим нерекурсивный фильтр, реализующий алгоритм
|
|
N - 1 |
|
|
|
у (kT) = S blX[(k-l)T]. |
(5.39) |
|
|
/=о |
|
Если |
в |
выражении (5.38) хотя бы один из |
коэффициентов |
at Ф 0, |
то |
фильтр, реализующий этот алгоритм, |
называют ре |
курсивным. Очевидно, что НФ представляет собой устройство без обратной связи, а РФ — устройство с обратной связью.
Строго говоря, цифровые фильтры представляют собой нелиней ные устройства, к которым не применимы методы анализа и син теза линейных систем. Однако число разрядов в кодовых словах, циркулирующих в этих фильтрах, как правило, достаточно ве лико, чтобы сигналы считать приблизительно дискретными, а фильтры — линейными дискретными. Это позволяет использо вать известные методы анализа и синтеза подобных устройств. Вводимые ниже характеристики (передаточная функция и др.) относятся к линейным дискретным фильтрам, точно реализую щим алгоритм (5.38), а эффекты квантования величии х (kT), у (kT), аь bt и их произведений будут учтены путем анализа про хождения шумов квантования через линейный дискретный фильтр.
Передаточной функцией (или системной функцией) цифровых фильтров называют отношение z-образов выходного Y (z) и вход
ного X (z) сигналов фильтра |
при нулевых начальных |
условиях |
£18]: |
|
|
Н (z) — |
Y (z)!X (z). |
(5.40) |
Здесь z-образ получают в результате г-преобразования, т. е. применения преобразования Лапласа к дискретному сигналу, записанному в виде последовательности 6-функций [32]:
/Ю = |
а НЬТ)6(1-кТ); |
|
|
ОО |
&=0 |
ОО |
(5.41) |
|
|||
F (s) = f |
dt = 2 / (kT) е-*т\ |
|
|
0 |
|
к=0 |
|
если обозначить |
|
|
|
|
e*r = |
z, |
(5.42) |
то из выражений (5.41) получим |
z-образ |
|
|
|
ОО |
|
|
F(z) = Е |
f(kT)z~k. |
(5.43) |
|
|
A=Q |
|
|
ш
Для рекурсивных и нерекурсивных фильтров из выражений
(5.38) и (5.39) с помощью |
(5.43) получаем: |
|
Я р <*) = ( Д |
Ь р -')К 1 + % аГ ‘) ; |
<5-44) |
|
N-1 |
|
Я н (г) = £ biz4. |
(5.45) |
|
|
1—0 |
|
Комплексные частотные характеристики цифровых фильтров представляют собой функции, полученные в результате подста новки z = exp (/©!Г) в передаточные функции (5.44) или (5.45):
Я |
р |
= ( д |
+ Д |
(5.46) |
|
|
|
N-1 |
|
|
|
Я н (е/“Г) = |
S bte - i^ Tl. |
(5.47) |
|
|
|
/=о |
|
Модуль |
комплексной частотной характеристики |
Л (о) = |
||
*= | Я lexp (/©71)] |, |
называемый |
амплитудно-частотной |
характе |
ристикой (АЧХ) фильтра, определяет амплитуду выходного сигнала устойчивого фильтра в установившемся режиме при входном сигнале х (kT) = exp(jkaiT), а аргумент <р (©) = = arg [Я (exp [/©71])], называемый фазо-частотной характерис
тикой |
(ФЧХ) фильтра, — фазу |
выходного |
сигнала |
при |
том же |
||||||
входном сигнале: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Г/N -l N-l |
|
|
\ |
|
|
|||
|
Лр (ш) = |
|
2 |
,Е_ &пА cos (m — k)(x)Tj х |
|
|
|||||
|
|
|
. \ m = Q |
fc=Q |
|
|
|
|
|
||
|
/Л1-Г* |
A f - I |
|
|
|
|
\ - I |
|
|
|
|
|
X ( 2 |
S |
aparcos (p — г) © Г) |
, a, = |
1; |
(5.48) |
|||||
|
s p = Q |
r = Q |
Г N-l |
|
|
I N-l |
|
|
|
||
|
|
|
|
|
|
b, cos laT |
|
|
|||
|
<pp(to) = —arctg I Д |
bi sin l<s>T / Д |
I + |
|
|||||||
|
|
|
M-1 |
|
|
I M-l |
|
1 |
] |
|
|
|
|
|
[ |
2 Qat sin mT I 2 |
ai cos i©T I ; |
|
(5.49) |
||||
|
|
|
/ |
N-l |
N-l |
bmbhcos (m — k) ©Г; |
|
(5.50) |
|||
|
|
|
|
2 |
S |
|
|||||
|
|
|
|
m = 0 |
fc=Q |
|
|
|
|
|
|
|
|
|
|
AM |
|
|
/ AM |
|
-I |
|
|
|
|
|
|
[ |
2 |
^sln/w T / 2 |
bjC os/© rJ. |
(5.51) |
|||
Пример 5.1. Пусть |
Я (z) = |
1/(1 — 0,5z-1). |
Из |
выражений |
(5.44) и (5.46), |
||||||
(5.48) и |
(5.49) получим: |
|
|
|
|
|
|
|
|
|
|
|
Я (е/шГ) = |
( 1 - 0,5е-?й)Г) - 1; |
Л И |
= l / K I ^ - c o s < o f ; |
|
Ф (со) = arctg [0,5 sin ч>ТЦ\ — 0,5 cos wT1)].
Из формул (5.44)...(5.51) следуют основные свойства частот ных характеристик цифровых фильтров с вещественными коэф фициентами:
1) все частотные характеристики представляют собой перио
дические функции частоты то с |
периодом сод = 2я/Т; |
||
2) АЧХ |
фильтра является |
четной функцией, а ФЧХ — не |
|
четной функцией |
частоты со. |
|
|
Из указанных свойств следует, что требования к частотным |
|||
характеристикам |
фильтров при |
постоянном Т задают лишь на |
|
интервале |
10; сод/2], а входные аналоговые сигналы должны |
быть ограничены по частоте сверху величиной <од/2. Импульсная характеристика (весовая функция) h (kT) цифро
вых фильтров представляет собой реакцию фильтра при нулевых начальных условиях на входное воздействие:
к = 0;
к^= 0.
Из этого определения и определения передаточной функции фильтров следует, что
|
h{kT) = Z-1 \Н (2)|; |
Н (z) = Z\h{kT)), |
(5.52) |
|
где Z {•) |
и Z-1 {.} — операторы прямого и обратного |
г-преобра- |
||
зований. |
Из формул (5.52) следует, что h {kT) и Н [ехр (/о>Т) ] |
|||
связаны |
парой преобразований |
Фурье: |
|
|
|
|
+л/т |
|
|
|
h {kT) — |
Г |
H{e^T) e ^ Tda>; |
|
|
|
- п/Т |
|
(5.53) |
|
|
|
|
|
|
Я(е/“Г) = |
S |
h(kT)erl*»T. |
|
|
|
k~Q |
|
|
5.3.2. Некоторые формы реализации цифровых фильтров. На точность филь* тров влияет форма реализации. Из всего многообразия вариантов реализации выделим два: прямую форму и каноническую. Прямая форма соответствует реали зации фильтра согласно выражениям (5.38) или (5.39). Структурная схема фильтра
.в прямой форме (рис. 5.1) содержит элементы задержки (г-1), т. е. регистры па мяти, блоки умножения (У.н) на коэффициенты и сумматор.
Канонической форме (рис. 5.2) соответствует последовательная реализация сначала рекурсивного фильтра, соот ветствующего знаменателю передаточной функции (5.44), а затем реализация не рекурсивного фильтра, соответствующего числителю в формуле (5.44). При этом удается уменьшить число элементов за
держки |
по сравнению с их числом при |
|
прямой |
форме от |
величины Ln — N + |
+ М — 2 до LK~ |
max (N — 1, М — 1). |
Рис. 5.1. Структурная схема цифрового фильтра в прямой форме
Рис. 5.2. Структурная схема цифрового фильтра в кано нической форме
5.3.3. Оценка по грешности цифрового фильтра. Фильтр может быть описан с помо щью линейной модели, представляющей собой совокуп ность линейного
дискретного фильтра и определенного числа ограниченных по абсолютному значению аддитивных воздействий, учитывающих эффекты квантования.
Пример 5.2. Реализуется в прямой форме рекурсивный фильтр с передаточ ной функцией Н (г) = (z)/A2 (z). Линейная модель погрешностей цифрового фильтра представлена на рис. 5.3 и содержит три источника погрешностей: е0 —
Рис. 5.3. Линейная модель погрешностей цифрового фильтра в прямой форме
шум округления (квантования) входного АЦП; eyj и еу2 — шум цифрового умно жения при реализации нерекурсивной части и рекурсивной соответственно. Как видно из рис. 5.3, шум е0 проходит через фильтр с передаточной функцией Н (г) всего фильтра, а шумы и еу2 — через фильтр, соответствующий только рекур сивной части т. е. с передаточной функцией G (z) — 1М 2 (г).
Пример 5.3. Реализуется в канонической форме рекурсивный фильтр с пере даточной функцией Н (г) = Аг (г)/Л2 (г). Линейная модель погрешностей цифро-
Рис. 5.4. Линейная модель погрешностей цифрового фильтра в канонической форме
вого фильтра представлена иа рис. 5.4. Как видно из рис. 5.4, шумы е0 и вУ2 проходят через весь фильчр, а шум eVi, соответствующий нерекурсивной части, появляется непосредственно на выходе фильтра.
Введем следующие предположения:
1) входной сигнал х (кТ) нормирован в соответствии с неравен ством
max |x(fe7’) j< 1; k>0
2)разрядность входного АЦП после запятой равна пх;
3)разрядности (после запятой) всех умножителей и суммато
ров равны п\ 4) при усечении после умножения используется округление.
Найдем сначала оценку максимальной погрешности цифрового фильтра. Максимальные погрешности АЦП и умножителей:
ео. шах == |
By Шах = <7*/2, |
(5 .54) |
где qx = 2 ~ \ qn = 2~л.
Максимальная погрешность на выходе составляет
Jyi шах '< Г{By max — / i*7n/2>
сумматора фильтра
(5 .55)
где rt — число блоков умножения, подключенных к /-му сумма тору.
Составляющая выходной погрешности, обусловленная кван
тованием в АЦП,
оо
в0. вых = max | е0. вых (kT)| < |
max | е0 (kT)| £ |
IА &Т)\ < |
|
*>о |
|
*>о |
|
< 0 , 5 f c S |
|
\h{kT)\. |
(5 .56) |
k=Q |
|
|
Составляющая выходной погрешности, обусловленная огра ниченностью разрядной сетки блоков умножения, подключенных
к i-му сумматору, |
|
|
00 |
|
|
|
|
ev* вых = max | бу{ ВЫ1 (^^1)| < |
max | В { (kT)\ £ 18i (АТ1)!< |
||
й>о |
|
k>a |
*=о |
< 0 , 5 ? * £ |
|
I f t W I . |
(5-57) |
k=0 |
|
|
|
где gi (kT) — импульсная характеристика |
части фильтра, сле |
дующей после t-ro сумматора.
Полная выходная погрешность цифрового фильтра может
быть оценена |
сверху |
соотношением |
|
|
||
|
|
|
|
00 |
|
|
евых = |
ео. вых |
«> |
e Y* BblX < 0 , 5 ? , 2 |
|А ( * Л | + |
|
|
|
|
|
6-0 |
|
|
|
|
+ 0 |
, 5 ? . S |
r , S l ? i ( M ') |. |
|
(5 .58) |
|
|
|
(») |
k=0 |
|
|
|
Пример 5.4. Цифровой фильтр имеет передаточную функцию Н (г) = |
ft0/(1 — |
—ajZ-1) и реализуется в канонической и в прямой форме. Найти оценку максималь ной погрешности при пх = п = 10, а%= 0,99, Ь0 = 0,01.
В данном случае имеем сочетание нерекурсивного фильтра с передаточной функцией #„ (z) — b0 и рекурсивного с передаточной функцией Я р (z) = 1/(1 —
— ajZ-1), чему соответствуют импульсные характеристики h(kT) = bQa\ ng{kT)=
~ aj. В соответствии с соотношением (5.58) и рис. 5.4 имеем следующую оценку предельной погрешности фильтра для канонической формы:
|
00 |
8пых < |
+ °*54л) £ V * + °>5Чп = 4*0 (* — ° l) _1 + °>5? = |
|
А—О |
— q (I — a»)-1 + 0,5) = 2-м [0,01 (1 — 0.99)-1 + 0,5] « 0,0015.
ев ы х < 0 ’5(?* Е |
bQa i + 2 -°»59л £ °? = |
°»5(7^о ( J |
— e i ) “ l H - ? ( 1 ~ ° l ) " I== |
А=0 |
ft=0 |
|
|
= 0,5-2-10.0,01 (1 — 0.99)-1 + |
2-ю (i _ |
0.99)-1 « 0,10. |
Как видно из примера, прямая форма реализации в данном случае может привести к максимальной погрешности, достигающей 10 % диапазона измерений.
Найдем теперь дисперсию выходной погрешности цифрового фильтра в предположении о некоррелированности соседних шу мовых элементов всех источников шума.
Дисперсии шума квантования входного сигнала и на выходе блоков умножения:
D[e„] = ??/12; D[ev] < ^ /1 2 . |
(5.59) |
Дисперсия эквивалентного шума на выходе сумматора фильтра
D [е?г] = rtD [ev] = г$Ц 12, |
(5.60) |
где Ti — число умножителей, подключенных к f-му сумматору. Дисперсия составляющей выходного шума, обусловленной
квантованием входного сигнала, определяется соотношением
D [8„. „«] = D [е0] S № ФТ)]\ |
(5.61) |
А=0 |
|
а дисперсия составляющей, обусловленной квантованием на вы ходах умножителей, — соотношением
D [ет, „„] = r,D [ev] Б [gi (kT))K |
(5.62) |
Итоговая дисперсия выходного шума цифрового фильтра оце нивается соотношением
D [евых] = D К . >BJ +
+ 4u
Б |
D [ет, вы1] < |
- § S [h фТ)Is + |
(О |
|
lji k—Q |
S Г, s [gi |
(5.63) |
|
(0 |
fe=a |
|
Пример 5.5. Рекурсивный фильтр имеет передаточную функцию Н (г) —
— V(1 — Oiz-1). Найти дисперсию выходной погрешности для прямой и канони ческой формы.
Используя промежуточные результаты, полученные в примере 5.4, получим для канонической формы реализации
Для прямой формы реализации получим следующую оценку дисперсии вы ходного шума:
D [8вых1 ^ |
ЪЪа\к + 2 А12. |
2 * |
12(1 — aj) |
* |
|
к=0 |
А=0 |
|
|
Для конкретных величин, имеющихся в примере 5.4, найдем средние квадра |
||||
тические значения выходного шума: при |
канонической форме |
реализации |
||
~]/D [sDbII] = 0,0003, при прямой форме "|/I) [еВЫ1] = |
0,003. |
|
Часть третья
ПОВЫШЕНИЕ ТОЧНОСТИ И ПОМЕХОУСТОЙЧИВОСТИ ИЗМЕРИТЕЛЬНЫХ ПРЕОБРАЗОВАНИЙ
Глава шестая
ПОВЫШЕНИЕ ТОЧНОСТИ ИЗМЕРИТЕЛЬНЫХ ПРЕОБРАЗОВАНИЙ
Максимизация точности измерений при фиксированных затра тах и обеспечение требуемой точности при минимальных затра
тах |
ресурсов — две основные, тесно связанные между собой за |
|
дачи |
метрологии. |
|
Возможные методы повышения точности измерений можно |
||
разделить на три группы [9, 21, |
28, 36 и др. J: |
|
совершенствование алгоритмов |
измерений; |
|
применение вспомогательных измерений или образцовых сиг |
налов и преобразований для введения поправок в результат изме рения (коррекция);
совершенствование аппаратной части (конструкторско-техно логические решения).
Первая группа методов ориентирована на уменьшение методи ческих погрешностей. Во второй главе охарактеризован общий подход к решению этой проблемы, опирающийся на применение адекватных алгоритмов с выбором рациональных (в идеальном случае — оптимальных) значений управляемых параметров. Иначе говоря, эта группа методов решает задачу приведения вида алгоритма измерений в соответствие с априорной информацией о свойствах объектов измерений, условиях измерений, предъяв ляемых требованиях и наложенных ограничениях.
Вторая группа методов, используя достижения первой, а также дополнительную информацию об условиях измерений и состоянии аппаратуры, получаемую с помощью вспомогательных измере ний или образцовых сигналов и преобразований, позволяет умень шить как методические, так и инструментальные погрешности.
Наконец, третья группа методов, связанная с совершенство ванием конструкторско-технологических решений, обеспечиваю щих повышение стабильности и соответствия реальных характе ристик аппаратуры номинальным, позволяет уменьшить инстру ментальные погрешности.
В предыдущих главах преимущественно рассматривались вопросы, относящиеся к первой группе методов. В данной главе рассматриваются методы второй группы. Третья группа методов в настоящей книге не рассматривается.
6.1. КЛАССИФИКАЦИЯ МЕТОДОВ УМЕНЬШЕНИЯ ПОГРЕШНОСТЕЙ
При синтезе метода (алгоритма) измерений и его реализации приходится решать задачу обеспечения заданных метрологиче ских характеристик, в частности обеспечения требуемой точности. Погрешность измерительного преобразователя (ИП) состоит из многих составляющих с различными характеристиками. В связи с этим необходимо рассмотреть методы уменьшения составляю щих погрешности ИП и взаимосвязь этих методов между собой.
В настоящее время стали очевидными принципиальная огра ниченность возможностей совершенствования ИП и взаимная связь отдельных метрологических характеристик между собой [54]. Поэтому для совершенствования ИП необходимо с помощью доступной технологии добиваться наилучших метрологических и технических характеристик при приемлемых экономических показателях, а затем, используя известные способы обмена одних характеристик ИП на другие, конструировать ИП с заданными характеристиками. В этом смысле можно говорить о принципиаль ной необходимости иметь избыточность по тем или иным харак теристикам ИП для того, чтобы за их счет улучшать другие ха рактеристики.
В практике измерений для целей улучшения точностных ха рактеристик ИП чаще всего используется запас по чувствитель ности, быстродействию и энергообмену с объектом измерения. Однако наличие этого запаса не решает автоматически задачу уменьшения погрешностей. Необходимо практически реализо вать эту возможность по отношению к конкретным составляющим погрешности.
Погрешность преобразования можно записать в виде
by(t> 1> Л» х) = /р (*» fli»***» fl7i)—/ном (*« an)i (6*1)
где /г, (♦) и /ном(-) — реальная и номинальная характеристики измерительного преобразования; а* — параметры характеристики преобразования; t — время; £ — влияющие факторы; rj — неинформативные параметры полезного сигнала х.
В наиболее распространенных линейных измерительных преобразованиях обычно выделяют три составляющие погрешности — погрешность нелинейности Ды (х), аддитивную погрешность Да и мультипликативную Дм:
&У (tf S> x) = AH(x) + Да (/, I» tj) + AM(t>St *1» x). |
(6.2) |
Классификационный |
Способ уменьшения погрешности |
||
|
признак |
|
|
Использование |
апри |
Использование априорной информации; инва |
|
орной |
информации о сиг |
риантный; адаптивный |
|
налах |
и погрешностях |
|
Вид уменьшаемой по грешности
Используемая избыточ ность измер1гтельного преобразования
Спос<)бы |
ре 1лизации |
структур пой |
избыточно |
сти |
|
Уменьшение погрешности: систематической; слу чайной; нелинейности; аддитивной; мультиплика тивной; статической; динамической
Использование избыточности: чувствительности; быстродействия; энергообмена
Стабилизация характеристик; компенсация; адди тивная коррекция; мультипликативная коррекция; введение поправки; самонастройка; способ образцо вых сигналов; способ итераций; фильтрация
Каждая из частных погрешностей в общем случае должна рассматриваться как случайный процесс с определенными характеристиками, которые и опреде ляют эффективность применения различных способов уменьшения погрешностей ИП.
Как уже указывалось, все методы повышения точности измерительных преоб разований делятся на три группы: совершенствование алгоритмов (с целью умень шения методических погрешностей), конструктивно-технологические методы (для уменьшения инструментальных погрешностей), структурные методы (уменьшают инструментальные погрешности и иногда — методические). Детальная классифи кация структурных методов представлена в табл. 6.1.
В основе структурных способов лежит принцип инвариантности (многока нальное™) [64]. Под инвариантностью понимают компенсацию возмущений, т. е. достижение полной или частичной независимости результата измерительных преобразований от дестабилизирующего фактора. В таких ИП помимо основного канала преобразования создается второй канал (рис. 6.1). Выходная величина ИП образуется в результате вычитания соответствующих величин основного ОК и вспомогательного ВК каналов.
Для такого ИП можно записать
Уi 00 = Gi 00 X (s) ф Gf Z (s); К; (s) = G{ (s) X (s) + Gf * (s) Z (s), (6.3)
где Уф (s), У] (s) — операторная запись выходных сигналов в основном и вспомо гательном каналах; Gt (s), G[ (s) — передаточные функции каналов по информа ционному сигналу; Gf (s), Of* (s) — передаточные функции по дестабилизирую щему фактору; X (s), Z (s) — операторная запись
входного и дестабилизирующего сигналов. Выходной сигнал такого ИП определяется
соотношением
У (S) - У, (S) - Y[ (S) = Ог (s) X (s) +
+ G* (?) Z (s) - G{ (s) X (s) - Gt * (s) Z (s).
(6.4)
Рис. 6.1. Структурная схема инвариантного к возмущению двухканального измерительного преобразователя
Если добиться равенства передаточных функций по дестабилизирующему сигналу обоих каналов G* (s) = G** (s) и инвертирования полезного сигнала во втором канале, то
K(s) = [Gx(s) + GHs)]*(s),
т. е. получаем ИП с повышенной чувствительностью и полным отсутствием влия ния дестабилизирующего фактора. Подобный способ используется, например, в дифференциальных преобразователях, в частности во входных каскадах опера ционных усилителей.
Весьма широко для уменьшения погрешностей ИП применяются способы стабилизации, компенсации и коррекции погрешности. Ниже рассмотрены основ ные характеристики указанных способов.
6.2. СТАБИЛИЗАЦИЯ РЕАЛЬНОЙ ХАРАКТЕРИСТИКИ ИЗМЕРИТЕЛЬНОГО ПРЕОБРАЗОВАНИЯ
Основным способом стабилизации реальной характеристики измерительного преобразования, широко применяемым в практике
измерений, |
является способ |
отрицательной |
обратной связи |
(рис. 6.2). |
Здесь ПП и ОП |
измерительные |
преобразователи |
прямой и обратной связи. Этот способ является универсальным в том смысле, что с помощью отрицательной обратной связи умень шается суммарный эффект от действия разнообразных дестаби лизирующих факторов. Однако с помощью только лишь обратной связи невозможно решить задачу повышения точности во всех случаях.
В статическом режиме работы преобразователя с обратной связью предельное значение Ауй аддитивной погрешности, при веденное к выходному сигналу, можно определить из соотношения
+ |
< 6 - 5 > |
где Ах и Аг/0 — предельные значения |
аддитивной погрешности |
на выходе прямой и обратной цепей; К и х — передаточные коэф фициенты преобразователей ПП и ОП. Из выражения (6.5) видно, что если аддитивная погрешность приложена к выходу преобра зователя ПП, то при К% -*• °° она полностью устраняется. Если же аддитивная погрешность действует на входе Л Я , то она не устраняется введением обратной связи. В этом случае относи тельная аддитивная погрешность преобразователя с обратной связью будет
Ayjy = Ах'/х. |
(6.6) |
Предельное значение относительной мультипликативной по грешности преобразователя с отрица тельной обратной связью можно найти из соотношения
Ayly = kcAKIK + ( l - h ) Дх/х, (6.7)
Рис. 6.2. Структурная схема измерительного пре образования с отрицательной обратной связью