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

книги / Цифровая обработка сигналов в измерительной технике

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

AA, < ô t/ ( K l + ft?/Aîr) l l + V r i +

k%

(149)

где 6U = AUJU! — относительная погрешность

действующего

зна­

чения первой гармоники.

Если коэффициент гармоник kT не мал, то оценка (149) вполне эф­ фективна. Однако при малых коэффициентах гармоник эта оценка

получается завышенной из-за малой величины kl в знаменателе. Завы­ шение оценки произошло потому, что при малых коэффициентах гармо­

ник величина VI = Ш 1)\ есть разность двух близких по численно­ му значению величин, каждая из которых зависит от одних и тех же мгновенных значений сигнала,‘так что погрешности в величинах U и Ui фактически в значительной степени компенсируются. Для того чтобы учесть это обстоятельство и получить бодее точную оценку по­ грешности Д£г, справедливую и при малых коэффициентах гармоник, необходимо получить формулу для погрешности, в которой бы явным образом была учтена зависимость от всех мгновенных значений сигна­ ла и их погрешностей. Выраженные через мгновенные значения сигна­ ла все формулы для коэффициента гармоник, приведенные в гл. 4.1, оказываются'эквивалентными, поэтому для определения погрешностей можно использовать любую из них.

Вычислим квадрат действующего значения первой гармоники сигна­

ла

 

 

^1 = “7=Г (^ U +

y)i

 

где

2

m

2

 

и Ух= - ^ - S ‘Wîcosa?l

= — 2 «„sintv

 

Ul

<7=1

m

<7=1

Здесь для сокращения записей обозначено и (tq) = uq. и (ûtq = <xq. Тогда

2 т т

У |= - г г Е S “Л cos(aç— as).

т<7 = 1 s = 1

Таким образом, для квадрата действующего значения высших гармоник получаем

1

т

т т

и \ = и2 — и\ = —

2 <4 —

2 « Л cos (aq— a s).

m

<7=1

m <7=1 S=1

Найдем^ погрешности Д (i/o) и Д (Uf) и выразим их через погрешности

£а в мгновенных значениях сигнала: lq =

uQuqi где uq — прибли­

женное, а tig — точное мгновенные значения сигнала:

А (и Ь = -^г 2

<&«cos (a, — a s) + -^ г 2 M icos К “ °0'- (150)

<7 ,s

m

<7,s

A(£/o) = 4 - ? u^ + 4 T Ç ^ — ^ 5 ^ C0S K - « s ) , (151)

где u0BII = tig Utxsin a , — U% cos a„ — мгновенное значение суммы высших гармоник сигнала.

Формулы (150), (151) дают явную зависимость погрешностей Д (t/î)

и Д (ÛI) от мгновенных значений сигнала и° и их погрешностей %, и являются исходными для вычисления погрешности коэффициента

гармоник по формуле (148).

 

величины

ДАГ. Обозначим

Оценим максимальную погрешность

максимальную погрешность величин

через £ =

макс | £в |. Оце­

ним входящие в формулы (150),

(151) члены, используя соотношение

т2 S U,£sCOS (®,

as)

и? |< 4 |( / ,

Q.S

 

 

 

записанное исходя из того, что среднее

значение величины | и° | не

превышает ее среднеквадратичного значения:

 

Принимая во внимание, что

 

 

 

тл• S 1,1, cos (a, — as) = -К-

S l / “* Г <

2ES,

находим J Д ((/?) I < 4| U + 2£2.

 

 

 

Получаем также

 

 

 

^ S | u ? , | < 2 £t/B—

 

Тогда |А((/“) | < 2 а д + 5а.

 

 

 

Подставляя полученные оценки Д (Ul) и Д ((/?)

в формулу (148)

и обозначая ô£0= | Ц}г — максимальную относительную погрешность

мгновенного значения, находим оценку абсолютной погрешности ко­ эффициента гармоник:

IДА,К Ôio +

+ [2^1 + А? 6Е0+ (Sio)*] А,-

(152)

Из анализа формул (152) и (149) видно, что при kr ~ 1

оценки ДА,

дают величины одного порядка, а при Ar

1 оценка по формуле (152)

дает значительно меньшую величину погрешности и является поэтому значительно точнее оценки по формуле (149). При малых коэффициентах гармоник (kf 1) третьим и четвертым членом в формуле (152) можно пренебречь по сравнению с первыми двумя. Это соответствует тому, что при малых коэффициентах гармоник составляющей результирую­

щей погрешности, вносимой величиной Д ((/?), можно пренебречь и учи­

тывать только составляющую погрешности Д ((/в).

Для измерения коэффициента гармоник с высокой точностью, т. е. с малой относительной погрешностью | ДАrlkr | <£ 1, необходимо вы­

полнение условия 8Çjkt <£ 1. При таком условии можно пренебречь и вторым членом в формуле (152) и получить оценку относительной по­ грешности

|A £ ,K ô |0.

(153)

В тех случаях, когда погрешности £? являются случайными (на­ пример, за счет погрешностей квантования либо наличия помех), целесообразна оценка погрешности Akv по статистическим парамет­ рам — математическому ожиданию и дисперсии.

При вычислении математического ожидания и дисперсии относи­ тельной погрешности Ak„ будем считать, что величины образуют стационарную, некоррелированную последовательность с нулевым средним значением (если среднее значение величии 19 не равно нулю, то при необходимости можно внести поправку -на систематическую погрешность в величины и9).

Вычислим математическое ожидание относительной погрешности коэффициента гармоник Akv. Используя формулы (150) и (151), нахо­ дим

 

М[Д ( i / ? ) ] = ~ о2; М[А<У\)\

от—2

 

от

 

 

где о* — дисперсия величин £,.

 

Из соотношения (148) получаем

 

При т

1 эта формула упрощается!

 

где 61г = о!иг — относительное среднеквадратичное значение по­ грешности мгновенного значения.

При вычислении дисперсии погрешности коэффициента гармоник D (Akr) можно ограничиться первыми, линейными по членами в фор­ мулах (150) и (151), так как учет следующих членов даст составляющие более высокого порядка по величине дисперсии (при вычислении мате­ матического ожидания ограничиться линейными членами нельзя* так как математическое ожидание этих членов равно нулю). Для вычис­ ления дисперсии абсолютной погрешности величину àkp приближенно запишем так:.

26/? ( k2 т 2

JJ UQ^Scos (а7 a À .

В этой формуле величина |0заменена на

потому, что среднее значе­

ние величин равно нулю. Тогда дисперсия абсолютной погрешности

D (Айг) = <А/гг) = - ^ 4- 1 ^2 S и°щим +

+ ~ r S M!W»I (A|sM st) cos (a„ — a s) cos (ce0i — a si) —

4.s

çl.sl

S « <A?,Ais> COS K t — a s) J .

Для вычисления сумм, входящих в D 6г),вослользуемся стацио­ нарностью и некоррелированностью величин iQ:

 

 

(о2,

q = s;

 

 

 

( W

- | 0i

Яфь.

 

 

"Тогда для первого члена в фигурных скобках получаем

 

 

7 Г Г s « X

<Д£^Ь> -

kbn

 

 

 

«r/7t4 q,s

 

 

 

 

для второго члена

 

 

 

 

 

16

S “X i <AIjA?si) COS (а, —css) cos (o„i — a si) =

 

•^г

m

 

Q,S

 

 

 

 

 

£?1. Si

 

 

 

 

 

•а третий член обращается в нуль.

 

 

 

 

Тогда дисперсия погрешности

 

 

 

 

 

D {№<)■■

кт

_ т *

,и*

1).

 

mU\

 

 

( « +

 

mVj

 

 

 

 

При

кг ^ \

 

 

 

 

 

 

D(Aftr)=»(6y 2/m.

 

 

 

В качестве результирующей статистической

погрешности, учиты­

вающей как смещение (математическое ожидание), так и дисперсию

•относительной погрешности, можно

принять

М (Д£г) +

\(D (ДАГ).

Учитывая полученные

соотношения,

находим

при m ^

Г, kv 1

М (Aft,) +

У Щ Щ с*

 

 

(154)

В зависимости от относительной дисперсии б |х, числа точек отсчета m и значения коэффициента гармоник kr может оказаться/определяющим либо первый, либо второй член формулы (154). Так, при ô |x

kjV m

второй

член значительно превосходит первый и

 

 

M (ДА,) + У Щ Щ ы à y ^ r n .

При

krtY т

6| х

М(ДА,) + V П(ААг) гьг (« д а ? .

Ив том, и в другом случаях для получения высокой точности измерения коэффициента гармоник необходимо обеспечить выполнение условия

«&L « К

Оценим погрешность коэффициента гармоник, вносимую арифме­ тическими узлами.

Если считать, что все вычисления в арифметических узлах выполня­ ются без сбоев, то погрешность этих узлов определяется только по-

грешностью округления результатов, т. е. в конечном йтог’е определя­ ется разрядностью счетчиков. При этом следует руководствоваться сле­ дующим соображением: вычисления должны проводиться с такой точностью, чтобы не загрублять точности измерения, т. е. погрешность вычисления должна быть не больше (желательно меньше) погрешности измерения мгновенных значений исследуемого сигнала. Однако пого­ ня за очень высокой точностью вычислений приводит к чрезмерному увеличению разрядности арифметических узлов и объема перерабаты­ ваемой информации без существенного повышения точности окончатель­ ного результата. Поэтому при вычислении следует иметь в запасе один, максимум два разряда, которые в окончательном результате следует округлить.

При известной относительной погрешности измерения мгновенных значений сигнала (в качестве 0£ следует взять либо 660при оценке по максимуму, либо б£г при статистической оценке погрешности) разрядность входного десятичного счетчика должна выбираться иэ

условия = 1g (1/66) + 1. Это соответствует объему счетчика 10*' =

=10/ 6£.

Точность, с которой следует возводить в квадрат, и объем счетчика,

Тв котором накапливаются квадраты, будут разными для алгоритмов (142) и (145). Алгоритм (143) не рассматривается, поскольку в нем нужно точное вычисление и задание фазы i|)lf что затруднительно. При реализа­ ции алгоритма (145) объем счетчика квадрата равен (10£r/ô|)2, а его разрядность /2 = 2— 21g (1//гг), так как величины (и„ U\xsin aq

UiyCùsaq) имеют порядок U1Мак<Л- При реализации алгоритма (142) разрядность счетчиков должна быть = 2/х. При не очень малых

kpразличие в объеме счетчиков /2и 12не очень существенно, но для ма­ лых kp может стать значительным.

Например, при kr = 0,1 % это дает шесть лишних десятичных разря­ дов.

Таким образом, применяя алгоритм (145) при малых ftr, можно существенно уменьшить объем счетчиков по сравнению с алгоритмом (142) либо при фиксированном объеме счетчиков значительно повысить

точность опредёления kr.

 

Заключительные операции: вычисление величины U\ =

(t/î* +

+ U\y)I2, деление на U\ и извлечение корня следует делать с

погреш­

ностью, меньшей (в' несколько раз или на порядок) относительной по­ грешности 6£г, т. е. при малых kt здесь можно допустить меньшую точнось вычислений.

Таким образом, измерение коэффициента гармоник kr с малой от­ носительной погрешностью ставит (при малых kr) исключительно вы­ сокие требования к АЦП. Необходимую точность АЦП можно опреде­ лить по формулам (152) или (153). Вычислительные узлы и их объемы, разрядность задания sin а и cos а следует выбирать так, чтобы на всех этапах вычислений вносимая ими погрешность была меньше погреш­ ности квантования. Число точек квантования т следует выбирать из условия, чтобы при измерении коэффициента kr влиянием гармоник,, начиная с номера 2), можно было пренебречь.

3.СПОСОБ ИЗМЕРЕНИЯ

МАЛЫХ КОЭФФИЦИЕНТОВ ГАРМОНИК

Простой способ измерения малых коэффициентов гармоник, поз­ воляющий уменьшить разрядность арифметических узлов, получим, представив формулу (НО) в виде

Аг = V V -U ù V + UùfVv

(155)

Поскольку при малых kPвеличина AU = U Ux мала по сравне­ нию с Uu то при перемножении величины AU на (U + Ut) необхо­ димая разрядность арифметических узлов, в которых происходит пе-

_2, Зависимость методических погрешностей от измеряемого значения коэффициента гармоник

fer, %

1,0

5,0

10,0

20,0

30,0

40,0

50,0

ДЙг, %

0,00001

0,0015

0,01

0,1

0,33

0,8

1,5

V, %

10

10

10

0,02

0,1

0,3

0,8

ремножение и фиксация результата, оказывается меньше,

чем при не­

посредственном использовании формул (142) и

(140).

 

 

При малых коэффициентах гармоник возможно упрощение принци­

па построения ИНИ за счет исключения операций сложения

(U +

иг)

и перемножения AU (U + £Д) [50]. Для этого выражение

(155)

пре­

образуем к виду

 

 

k„= V AU {2Ut + AU)IU1= V 2 VàUJÜ[ V i + AU/2Ult

 

где ДU ~ i/ — Ui — разность действующих значений исследуемого

‘Сигнала и его первой гармоники. _________

Используя разложение функции V 1 + Д{//2[/х в степенной ряд, получаем

К = V2 К Д Щ " (1 + AU/4UJ.

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

k, = V2AU(UX= V2(U UJ/UV

(156)

Получим оценку методической погрешности, допускаемой при этом.

Абсолютное, значение погрешности.

 

и

Ди

Ut

4Ut '

ИЛИ \Akr\ = (AU/AU1)kr

(157)

Из формулы (156) находим AU!U1 —~kl/2. Подставив

это равенство

в формулу (157), получаем

 

Akt = *?/8.

(158)

S6

В табл. 2 приведены методические погрешности для нескольких значений Аг.

Из этой таблицы видно, что методическая погрешность не оказыва­ ет заметного влияния на точность измерения в области малых зна­ чений коэффициента гармоник. Особенно незначительна эта погреш­ ность в пределах Аг до 10 %. Она может оказаться значительно мень­ ше, чем погрешность, вносимая квадраторами при определении Аг по формуле (140). Покажем это. Обозначим относительные погрешнос­

ти вычисления величин U2 и U\ через у и уг соответственно, тогда из

формулы (140) получаем приближенное значение коэффициента гар­ моник

/ ^ 2 Q + ï ) - ^ ( l + V i )

1

+С/?(1 + Ъ)

 

YU ' - У г Щ

V

V

+

t/f (1 + Yx)

Пренебрегая в знаменателе

величиной у2 (уг

1),

получаем

kr = к У \ + {уи * - ъ и ])1 Ш ).

Разложив выражение под квадратным корнем в степенной ряд до линейного члена, найдем

kr = К [1 + (yU2yxUi)/çzuikrfl

Абсолютная погрешность измерений, вносимая неточностью вы­

полнения операции возведения в квадрат,

 

д*; = £ _ й г = (уЦ2 -

ytuW u \kr).

Ее максимальное значение

 

 

ААр.макс = {\yU2\ +

\4iU\\)l2U\kr.

 

Принимая в этом выражении U ~

Ux и | у | =

| уг |, окончатель­

но получим

 

 

А&г = у/Аг.

(159)

Из сравнения выражений (158) и (159) следует, что методическая погрешность ДАР и погрешность ДАг, вносимая операцией возведения

в квадрат действующих значений, равны при условии у = А*/8. По этой формуле определим погрешность вычисления квадратов,

при которых обеспечивается погрешность измерения Аг, равная мето­ дической погрешности для значений Аг, приведенных в табл. 2. Полу­ ченные значения погрешности у запишем в эту же таблицу. Из таблицы следует, что для обеспечения одной и той же погрешности измерения малых коэффициентов гармоник по точной (140) и приближенной (156) формулам в первом случае необходимо обеспечить исключительно вы­ сокую точность выполнения операции возведения в квадрат (напри­ мер, при Ар = 20 % необходимо у < 0,02 %). Такие точности при реа­ лизации способа на аналоговых элементах вообще не достижимы. При реализации на дискретных элементах указанные точности реальны,

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

ГЛАВА 5

ЦИФРОВЫЕ ИЗМЕРИТЕЛИ ПАРАМЕТРОВ АМПЛИТУДНО-МОДУЛИРОВАННЫХ СИГНАЛОВ

1.ЦИФРОВЫЕ АМПЛИТУДНЫЕ МОДУЛОМЕТРЫ

Цифровые амплитудные модулометры основаны на измерении максимального 1/шкс и минимального U ma значений модулированного сигнала и вычисления по этим значениям глубины (коэффициента) модуляции по формуле

Р — (^макс

^мин)/(^макс "Ь Цмин)*

(160)

Максимальное и минимальное значения амцдитудно-модулированного (AM) ‘сигнала могут быть измерены по «выборочным» мгновенным зна­ чениям и способом сканирования.

Измерители коэффициента модуляции по «выборочным» мгновенным значениям

Запишем выражение для AM сигнала

 

U(t) = £/„ес.макс(1 + P SÎIl Qt) SÎn (ùt,

(161)

где UHeс.макс — амплитуда

несущей

в

отсутствие модуляции;

Q —

=

2лF —- модулирующая

частота;

о = 2л/ — несущая

частота.

 

Условие равенства нулю начальных фаз модулирующего и несу­

щего сигналов не меняет существа рассматриваемого

способа, но

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

L =

=

Мгновенные значения AM сигнала и (t) измеряют в точках

1(2д — 1)/2] (л/©), (q =

1, 2, 3,

 

2я),

где п = со/й.

 

 

В эти моменты

 

 

 

 

 

 

 

 

« iU) = Ûнес.макс 1 1 + Ц SÎn

2?~

'■

(— 1)"

 

 

 

преобразуют в цифровой код Ng — ktlu (tg),

где ku — коэффициент

пропорциональности.

 

 

 

 

 

 

 

 

Цифровые коды соседних мгновенных значений N„

и A/^+i)

по­

очередно сравнивают друг с другом до момента выделения их макси­ мального и минимального значений:

Л^иакс ~

МаКС [ | — kuU нес.макс [1 "f* P- Sifl (îl/2 -j~ ô j)]:

Nиод =

МИН | Nq| = АцС/нес.макс 11 — JJ.SÎn (îl/2 -(- 02)],

где ôj и 02—фазовые погрешности выделения моментов', соответствую­ щих максимальному и минимальному значениям напряжения огибаю­ щей.

Суммируя и вычитая М*акс и NMmt находим

=

Л/мако Ч" ЭДмн = 2^(7„ес.макСО

(162)

AL-. ~

Л^Мака Л^мин — 2 ku\xUнес.макс (1 + &У_).

(163)

В выражениях (162) и (163) ÔJV+ и 6ЛГ_ — относительные методи­ ческие погрешности определения кодов N+ и N -, причем

| W + 1= -J- р (cos Sj — cos S21~ p 16? — б*I /4;

|Ô A 4 = |l - 4 - ( c o s ô 1 + cos ô2) | ~ (Ô? -f- ô|)/4.

Поскольку

при

переходе от

 

одной точки квантования, <ùtq,

 

к другой,

©VH ),

аргумент

 

синуса ^модулирующей

изме­

 

няется на

величину ля,

т, е.

 

(<Ц?+ц /л — a>tJri) =

п/п, то

Рис. 25. Структурная схема цифрового ам-

] бх . <

п /2 я ,

I 6 ,

I <

я /2 л .

Тогда, относительные

методи-

плитудного модулометра

ческие

погрешности

 

 

 

 

 

| ÔJV+ к

 

(Vie)и(л2//г2); ! 6Л^_I < (Ve)(я*/л*).

Обычно л величина достаточно большая (гг ^ 1), так что | ÛN+ | <£ 1, | ÔAL. | <£ 1 . Поэтому методическими погрешностями в определе­

нии кодов N+ и

практически можно пренебречь:

 

N+ = 2kuUнес.макс» Af—= 2ku[lUнес.макс*

Амплитуда несущей (в отсутствие модуляции) и коэффициент мо­ дуляции определяются по найденным кодам N+ и ЛГ*,:

Uиес.макс в /v+/2A0;

(164)

р =

У _ /У + .

(165)

Их относительные методические

погрешности

ÔU = 6N+; 0 р < |0 У _ | +

| 6У+ |-

В основу построения амплитудного цифрового модулометра поло­ жены соотношения (162) — (165). Структурна^ схема такого модуло­

метра показана на рис. 25 [29J.

Амплитудно-модулированный сигнал и (/) через ВУ поступает на БУ и на потенциальный вход АЦП. Блок управления выделяет момен­ ты времени, соответствующие максимальным значениям напряжения несущей и (tg) и запускает в эти моменты АЦП. Последний преобразу­ ет мгновенные значения напряжения в цифровые коды Мд, которые поступают в ЦК, где они поочередно сравниваются друг с другом до момента выделения их максимального Nmкс и минимального Nma значений. Эти значения вводятся в ЛУ, где они подвергаются цифровой обработке в соответствии с выражениями (162), (163) й (165).

Нижний диапазон измерения коэффициента модуляции р определя­ ется чувствительностью и дискретностью АЦП, причем, чем меньше

амплитуда несущей, тем выше должна быть чувствительность преобра­ зователя Л£/п, которая выбирается из условия

K U n ^ Р (^макс)мни,

где (имакс)мив— минимально возможная амплитуда несущей.

Так, для (х = 0,001 и (U„аКС)ыин = 10 мВ получим Лt/„ ^ 0,01 мВ. Существующие АЦП такой чувствительности не обеспечивают [58, 63], поэтому исследуемый сигнал и ((] необходимо усилить.

Коэффициент усиления усилителя ky следует выбирать с учетом верхнего предела работы АЦП, который у АЦП на полупроводников

вых приборах может быть обеспечен и п.т«с = 10 В,

 

 

k, - Uп.макс/(^макс)иин Ш 10/(10 •

К )'3) -

10».

Тогда

чувствительность преобразователя

ДUn <

р£у (f/MaKC)MHH~

~ 0,001

• 10 = 0,01 В.

 

 

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

вателей. Таким образом, при использовании высокоточных

АЦП

и ЦК амплитудный цифровой модулометр позволяет обеспечить

высо­

кие метрологические характеристики.

Рассмотренному измерителю присущ и ряд недостатков: ограни­ ченный верхний предел частоты несущей; ограниченный нижний пре­ дел коэффициента модуляции и амплитуды несущей; наличие методи­ ческих погрешностей за счет нелинейных искажений модулирующего сигнала; сложность схемной реализации. Кроме того, он не обеспечи­

вает измерение фазового сдвига AM сигнала. Устранить в

какой-то

степени указанные недостатки позволяет цифровой способ

измерения

параметров AM сигналов их сканированием, причем возможно равно­ мерное и неравномерное сканирование.

Цифровой модулометр с равномерным шагом сканирования AM сигналов

Способ сканирования основан на измерении максимального и мини­ мального значений AM сигнала с последующим вычислением коэффи­ циента модуляции. В отличие от известных способов это измерение осуществляется сканированием ÀM сигнала на определенных калибро­ ванных уровнях. Суть способа поясняется временными диаграммами (рис. 26) [47].

Пусть .несущая модулируется несинусоидальным напряжением (рис. 26, а). Как известно, глубина модуляции вверх рв и вниз. рн определяется из выражений

Рв “

U u J U c , ~ Щмакс U CV) I U ср;

Pu =

Ц и.н/U ср — (Ц ср 1Ли!н)/1^СР»

где и к.ъ и Uи.в — максимальные положительное и отрицательное значения модулирующего напряжения ии (/); Ucp— среднее значение напряжения.

Напряжения £/Иако и UmB определяют следующим образом. Формируют калиброванное напряжение, диапазон изменения которр*