Файл: Ермолов Р.С. Цифровые частотомеры.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.06.2024

Просмотров: 116

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

и н т е г р а л ь н ых микросхем) позволяют р е а л и з о в а т ь

высоко

точные

и н а д е ж н ы е цифровые

л и н е а р и з а т о р ы .

 

 

 

Существуют три метода цифровой

л и н е а р и з а ц и и :

 

 

1. Использование

специального

функционального

генератора,

который

инвертирует

выходную характеристику

П И П

цифровым

способом.

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

ограни­

ченного класса специальных функций и является относительно до­ рогим.

2. Л и н е а р и з а ц и я выходной характеристики П И П .

П о с л е д н я я

заменяется по участкам отрезками п р я м ы х

(кусочно-линейная ап­

п р о к с и м а ц и я ) , различный наклон которых

реализуется

с помощью

двоичных умножителей . Недостатком метода является сравни­ тельно невысокая точность аппроксимации .

3.

Л и н е а р и з а ц и я кривой ошибок,

п р е д с т а в л я ю щ е й собой

раз ­

ность

м е ж д у линеаризируемой кривой

и ж е л а е м о й линейной

зави ­

с и м о с т ь ю . Этот метод позволяет получить высокую точность ли­

неаризации д л я широкого класса

функций.

 

 

 

 

 

Н и ж е р а с с м а т р и в а ю т с я

д в а способа реализации

третьего

ме­

тода цифровой

линеаризации .

 

 

 

 

 

 

 

 

В

частотных

П И П

и з м е р я е м а я

величина

воспринимается

эле­

ментами

частотозависимой

цепи

автоколебательной системы. Вы ­

х о д н а я частота

П И П о к а з ы в а е т с я

однозначной функцией

измеря ­

емой величины. Эта функция может быть

к а к линейной, т а к и не­

линейной.

Анализируя

в ы р а ж е н и я

д л я

функций,

с в я з ы в а ю щ и х

частоту с

измеряемой

величиной,

д л я

П И П , различных

по прин­

ципу

действия,

м о ж н о сделать вывод,

что большинство из них сво­

д и т с я

к следующему ряду [41]:

 

 

 

 

 

 

 

 

 

 

 

1)

/ = f t i

 

Ух,

5)

І = к х

У

X—m

 

 

 

 

 

 

2)

/ =КХ,

 

6)

f = k x

y

m — X

 

 

 

 

 

 

3)

/ =

kJX,

7)

f = k r

y

X + m

 

 

 

 

 

 

4)

/ =КіУх,

 

 

 

 

 

 

 

 

 

где f—выходная

 

частота

П И П ;

X—измеряемая

величина;

ти

к{ — постоянные.

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

ботку результата, чтобы была

воспроизведена функция,

о б р а т н а я

выходной характеристике П И П .

 

 

 

 

Д л я

выходных

характеристик

П И П

(5-42) м о ж н о

з а п и с а т ь

•обратные

функции:

X

=kf\

5)

X = m + kf\

 

 

1)

 

 

2)

X

=kf,

6)

X =

m—kp,

(5-43)

 

 

X

 

7)

X =

kf*—m.

 

3)

=klf,

 

 

4)

X

=kff*,

 

 

 

 


 

П р о с т е й ш ей формой аппроксимации

является

кусочно-линей­

ная. Л ю б у ю

функцию

X = (f(f)

м о ж н о

представить

в виде

л о м а н о й

кривой (рис. 5-10), которая получается делением

в о з м о ж н ы х зна ­

чений

аргумента

/

на

поддиапазоны

/ — 5,

внутри которых

функ­

ция

ф(/ )

з а м е н я е т с я

прямой .

Че м

у ж е

поддиапазон,

тем

точнее

аппроксимация .

Д л я

к а ж д о г о

из поддиапазонов

м о ж н о

записат ь

уравнения

прямой:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х в

= Х м + а Л ,

h<f<f3,

 

\

 

 

 

(5-44)

 

 

 

 

 

 

X3

=

XS0

+ a3f,

 

 

f3<f<h-

 

 

 

 

 

 

Определи м

 

необходимое

 

 

 

 

 

 

 

 

 

 

 

число

поддиапазонов,

 

обеспе­

 

 

 

 

 

 

 

 

 

 

 

ч и в а ю щ е е

з а д а н н у ю

точность

 

 

 

 

 

 

 

 

 

 

 

линеаризации .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

м а т е м а т и к е при

 

аппрок­

 

 

 

 

 

 

 

 

 

 

 

симации

функций

одним

из

 

 

 

 

 

 

 

 

 

 

 

наиболее

 

употребительных

 

 

 

 

 

 

 

 

 

 

 

способов

является

отыскание

 

 

 

 

 

 

 

 

 

 

 

аппроксимирующей

 

функции

 

 

 

 

 

 

 

 

 

 

 

по методу наименьших квад ­

 

 

 

 

 

 

 

 

 

 

 

ратов

[47]. П о

этому

 

методу

 

 

 

 

 

 

 

 

 

 

 

аппроксимирующу ю

функцию

 

 

 

 

 

 

 

 

 

 

 

находят из соотношений:

 

 

 

 

 

 

 

 

 

 

 

 

М = / [ /

(х)—ср

(x)f dx,

(5-45)

 

 

 

 

 

 

 

 

 

fm f

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о = Ум/(Ь—а),

 

 

(5-46)

 

Рис. 5-10. Кусочно-линейная

аппрок­

 

 

 

 

 

 

 

 

 

 

 

симация

нелинейной

выходной

ха­

где

интервалг

аппроксимации

 

 

 

 

рактеристики ПИ П

 

 

з а д а н значениями

аргумента

 

 

 

 

 

 

 

 

 

 

 

a^x^b;

 

f(x)

— аппроксимируемая

функция;

с р ( х ) — а п п р о к с и м и ­

р у ю щ а я функция;

о — среднеквадратическа я

погрешность

аппрок­

симации .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В наше м случае а п п р о к с и м и р у ю щ а я функция представляет со­

бой

пряму ю

вида

(5-44). Если

по

методу

наименьших

к в а д р а т о в

а п п р о к с и м и р у ю щ а я

функция отыскивается

в

виде

ср(%) = аофо(*)

+

+ аіц>і(х),

где Фо=1 , фі = х, то

д л я

определения коэффициентов

а0

и а,\ получается система

линейных уравнений:

 

 

 

 

 

 

1

дМ

п

= 2

l

oak

і = 0

ь

ь

щ (х) dx = 0,

 

at J

q>,(х) ф А (х) dx—П(х)

(5-47)

о

о

 

 

где

k = 0;

1.

 

 

Р а с с м о т р и м первую функцию

из р я д а (5-43), имеющую вид:

f(X)

=kx2,

и пусть Ь = ас. Тогда,

воспользовавшись в ы р а ж е н и е м


(5-47), дл я коэффициентов

аппроксимирующей

функции

получаем

 

 

 

3 с4 1

+ 1

 

 

 

 

 

 

 

 

 

 

 

 

4 б 3 — 1

 

 

•ка\

 

 

 

 

 

 

 

 

 

3

( с » - 1 ) ( с + 1 )

 

 

 

 

 

 

 

 

 

4

с 3 — 1

 

 

 

 

 

 

 

(5-48)

 

 

 

 

 

 

 

3

с4 1

с + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 (с — 1)

2

 

 

З ( С 2 - 1 ) (

£ а 2

 

 

 

 

 

 

t

 

С + 1 )

 

 

 

 

С учетом (5-46)

м о ж н о получить формулу:

 

 

 

 

 

б =

 

1

a 1 0 ( c 4 - l )

 

( 2 a o o - a i o ) ( с

3

- ^

 

 

 

 

Г

5

Гс —

2 ( с — 1)

 

 

3 (с — 1)

 

 

 

 

 

 

 

( с - 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

а

я 1 0 ( с + 1 ) +

а2

(5-49)

где

6 = а/(/га 2 ), a00 = aQ/(ka2),

aw = aj

(ka).

 

 

 

 

 

 

 

Втора я

функция

р я д а (5-43)

является

линейной.

 

 

 

 

Третью и четвертую м о ж н о

 

представить

так : X = kT;

X =

kTz,

где

T=l/f.

 

Третья

функция

является линейной

от периода,

а

чет­

вертая аппроксимируется прямой, с коэффициентам и и среднеква -

дратической

погрешностью

аппроксимации,

определяемым и

 

выра ­

ж е н и я м и

(5-48) и

(5-49).

 

 

 

 

 

 

 

 

 

 

 

 

 

П я т а я функция

отличается

от первой только

наличием

началь ­

ного

значения

т,

которое будет

д о б а в л я т ь с я к коэффициенту

а0

в аппроксимирующей

функции. В ы р а ж е н и е

д л я

среднеквадратичен

ской погрешности аппроксимации совпадает с в ы р а ж е н и е м

(5-49).

Совершенно очевидно, что первая и пятая функции при одина­

ковых интервала х

аппроксимации

(при одинаковых

с) будут

иметь

одинаковую

погрешность аппроксимации .

 

 

 

 

 

 

 

 

В

шестой и седьмой

функциях

нелинейная с о с т а в л я ю щ а я

совпа­

д а е т с первой функцией. Поэтому

и в ы р а ж е н и я

д л я коэффициен ­

тов аппроксимирующей функции

и погрешности

б совпадаю т с вы­

р а ж е н и я м и

(5-48)

и

(5-49). Отличие

з а к л ю ч а е т с я

только

в том,

что

при

расчете

аппроксимирующей

функции

вместо

 

коэффици ­

ента

с 0 , определяемого

в ы р а ж е н и е м

(5-48),

необходимо

принимать

коэффициент

ao' = ao — гп. Д л я шестой

функции

такое

утверждени е

не к а ж е т с я

очевидным. О д н а к о м о ж н о показать, что эт а функция

сводится

к

виду седьмой функции,

д л я которой

приведенные

рас ­

с у ж д е н и я очевидны. Шестая функция в числе

импульсов

може т

быть представлена

в виде Nx

= Nm

— Nf.

 

 

 

 

 

 

 

 

Если перед измерением в счетчике установить число N, а счет­

чик включить в р е ж и м вычитания, то результат

измерения

примет

вид:

NX

= N + Nf — Nm.

З а д а в

N = 0,

получаем

Nx

= Nf

— Nm,

что

эквивалентно

функции

вида

X = kf2

т. Если

включить

счетчик

в р е ж и м

суммирования, то

результатом

измерения

будет

число


Nx = Nm — Nf.

Считывая

этот

результат

в

дополнительном

коДб,

получаем

NX

= NC4 + Nj — Nm,

где iV C 4 — объем счетчика.

Изме ­

нив а0

так, чтобы N'm = Nm

+ NC4, м о ж н о

получить

NX = NC4 + Nf —

— NC4

— Nm = Nf — Nm,

что

эквивалентно

 

седьмой

функции.

Т а к и м

образом, д л я

всех

рассмотренных типичных

выходных

характеристик

частотных преобразователей

с помощью

в ы р а ж е н и й

(5-48) и (5-49) м о ж н о рассчитать зависимость погрешности ап­

проксимации

от относительного

интервала аппроксимации,

т. е.

зависимость

б = гр(Аа ), где

Д а = ( Ь —

а)/а

= с—1. Боле е того,

при

одинаковых

относительных

интер­

 

 

 

в а л а х погрешности

аппроксимации

 

 

 

д л я всех преобразователей будут

 

 

 

одинаковы.

Н а рис.

5-11

приведена

г

8, г)

 

зависимость

8 = t p ( A a ) ,

рассчитан­

10'

 

 

ная по в ы р а ж е н и я м

(5-48)

и

 

 

 

(5-49).

 

 

 

 

 

 

 

 

 

К а к

видно

 

из рис. 5-11,

при

от­

 

 

 

?— /

 

 

 

 

носительном

 

интервале

аппрокси­

 

 

 

TU

/

/

 

 

 

 

мации,

не

п р е в ы ш а ю щ е м 4%, по­

10'

 

 

'

 

,4

 

 

 

грешность аппроксимации не пре­

 

 

і '

J

 

 

 

 

 

вышает

 

0,012%.

 

Следовательно,

 

 

 

 

 

 

 

 

при 20% - ной девиации частоты на

 

 

/

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

выходе преобразователя, что харак ­

 

и

 

 

 

 

 

 

 

 

терно

д л я

большинства

частотных

10 —t

 

 

 

 

 

 

 

 

преобразователей,

достаточно

 

пяти

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

интервалов

аппроксимации

 

выход­

 

 

 

 

 

 

 

 

 

 

ной

 

характеристики

преобразова ­

 

 

 

 

 

 

 

 

 

 

теля.

Пр и

этом

погрешность

ап-

 

 

 

 

 

 

 

 

Аа

роксимации

не

превысит

0,012%.

10 0

 

 

 

 

 

 

 

 

Пусть

функция,

обратная

 

выход­

1

2

3

 

4-5

6

1 8

9

%

ной

 

характеристике

П И П ,

имеет

 

 

 

 

 

 

 

 

 

 

вид, приведенный на рис. 5-10

(кри­

Рис.

5-11. Зависимость погрешно­

вая

 

/ ) ,

 

а

нелинейность

выходной

сти

от

интервала

аппроксимации

характеристики

описывается

функ­

 

 

 

 

 

 

 

 

 

 

цией

f = kX2.

 

Пр и этом частота на

 

 

 

 

 

 

 

 

 

 

выходе

П И П

может

изменяться

в пределах от /о до fm.

З а д а в ш и с ь

допустимой

погрешностью

аппроксимации,

с

помощью

в ы р а ж е н и й

(5-48)

и

(5-49)

определяем

интервалы

аппроксимации,

внутри

к а ж д о г о

 

из

которых

з а м е н я е м

реальную

функцию

прямой

вида

(5-44).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выходная

 

характеристика

П И П

будет

линейной

(кривая

/ /

на

рис. 5-10),

если

о б р а т н а я

функция

ее имеет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

'

= a'o + a'ifx>

 

 

 

 

 

 

 

( 5

- 5 0 )

где

а'о и

а\

определяются в ы р а ж е н и я м и

(5-48)

при

c =

/ i / / 0 .

 

 

 

В

нашем

случае

о б р а т н а я

характеристика

нелинейна. Н а

пер­

вом участке

аппроксимации

эта характеристика описывается урав ­

нением:

Xi — aoi + On/jc, где

Сої и an

определяются в ы р а ж е н и я м и

(5-48)

при

c = fi/f0. Таким

образом,

д л я линеаризаци и выходной


х а р а к т е р и с т и ки датчика в результат измерения необходимо ввести

поправку,

равную

д л я

первого участка

аппроксимации

разности

 

 

{Х'~ХХ)

 

 

= АХ,

= a'Q-aol

+

{а[-ап)

 

fx,

 

 

или в числе импульсов:

 

 

 

 

 

 

 

 

 

 

 

 

 

Mi

 

=

{f'o-rm)T0+(T[-Tn)fx,

 

 

 

 

(5-51)

где f'0T0=a'0;

foiT0=a0i;

 

Т\ =

а\\

Тн—-аи.

 

 

 

 

Р а с с у ж д а я

аналогично, м о ж н о

показать, что

в общем случае

поправка

в числе

импульсов

определяется в ы р а ж е н и е м (48]:

 

 

 

 

 

Мс

=

(Го-ЦТ0+(Т-Ти)[х>

 

 

 

 

(5-52)

где і — порядковый

номер участка аппроксимации

(і = 0, 1,2,. .. , 5);

Т'і — образцовый интервал, тождественный коэффициенту а\;

Тц—

образцовый интервал, тождественный коэффициенту ац.

В ы р а ж е ­

ние (5-52) представляет собой аппроксимированную кривую

оши­

бок. Коэффициенты айг и ац рассчитываются с помощью

в ы р а ж е ­

ний (5-48), где

 

c=fi+i/fi.

 

 

 

 

 

 

 

 

Н а

первом

участке

при

fx = fo

имеем

A/Vi = 0.

Тогда

получаем

уравнение: (f0

— / ш ) Т 0 + (Т\

— Tn)f0

 

= 0,

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

' 01

' 0

 

 

 

 

П о д с т а в л я я

в ы р а ж е н и е (5-53) в

(5-52), получаем

 

 

 

 

AN і =

^If-(Т[-Ти)

/0+

[Т[-Ти)

 

f,

 

(5-54)

 

 

 

 

 

/оі

' 0

 

 

 

 

 

 

 

 

Частоту fx

можно

представить

в

виде суммы: fx = fo+(їх

— fo).

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f0^(T-Tn)fo+(T-Tu){i-f0).

 

 

 

 

(5-55)

К а к

видно

из

в ы р а ж е н и я

(5-54),

поправка

является

линейной

функцией

частоты.

Следовательно,

 

в счетчик

результата

необхо­

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

времени,

пропорционального измеренной

частоте.

Пусть

время

введения

поправки

 

 

 

 

 

 

TK = (Nx~N0)/fOK,

 

 

(5-56)

где

Тк — время

введения поправки или

время коррекции

резуль ­

тата

измерения;

/ок — некоторая о б р а з ц о в а я частота;

No, Nx

— по­

к а з а н и я

счетчика,

соответствующие частотам

f0,

fx

соответственно,

причем

NQ = T[fQ;

Nх = T[fx.

Тогда частота

импульсов

коррек­

ции результата измерения на t'-ом участке (fx

= fi)

д о л ж н а

быть

 

 

1 _|_

( ^ J - J j O i o ...

 

ТІТІ

/ок- (5-57)

 

 

 

 

 

 

 

 

 

fm-f'oj

T'i{h~fo)

 

 

T[

 

 

136