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

книги / Методы электрических измерений

..pdf
Скачиваний:
8
Добавлен:
19.11.2023
Размер:
53.63 Mб
Скачать

Дискретные и цифровые фильтры принято делить на два класса: нерекурсивные (НФ) и рекурсивные (РФ). Если в выражении (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)| £

&Т)\ <

*>о

 

*>о

 

< 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. Структурная схема измерительного пре­ образования с отрицательной обратной связью

Соседние файлы в папке книги