CC BY
16
Зависимость текущего такта измерения при адаптивной временной дискретизации в случае восстановления сигнала экстраполяционным
полиномом первой степени
С.В.Кавчук, Г.И.Ткаченко Южный федеральный университет, Таганрог
Аннотация: Эффективность алгоритмов адаптивной временной дискретизации (АВД) определяется сжимаемостью сигнала. В целях теоретической оценки сжимаемости измерительных сигналов рассматривается зависимость текущего такта измерения при АВД от динамических свойств сигнала. Восстановление сигнала в процессе АВД производится экстраполяционным полиномом Тейлора 1-й степени. Качество аппроксимации устанавливается критерием равномерного приближения. Для структурной модели сигнала получена функциональная зависимость текущего такта измерения от структурных свойств сигнала, позволяющая на ее основе определять среднюю длительность такта измерения при АВД для расширенной области допустимой погрешности воспроизведения. В качестве примера приводятся результаты полученные в среде Mathcad, иллюстрирующие восстановление сигнала полиномом первой степени, дискретный фазовый график и соответствующий ему текущий такт измерения на фазовой плоскости.
Ключевые слова: Адаптивная временная дискретизация, такт измерения, структурные свойства сигнала, экстраполяция, погрешность воспроизведения.
Одним из методов сокращения измерительной информации на этапе аналого-цифрового преобразования непрерывных (аналоговых) сигналов является адаптивная временная дискретизация (АВД) [1 - 5]. Эффективность применения алгоритмов АВД по сравнению с равномерной временной дискретизацией (РВД) [3, 4, 6] в цифровых информационно-измерительных системах определяется коэффициентом сокращения числа отсчетов (сжатия), который зависит как от вида алгоритма АВД (инструмента), так и от способности сигнала (материала) к сжатию.
Для априорной оценки сжимаемости аналоговых сигналов x(t) при равномерном критерии приближении нужно прежде всего определять зависимость текущего такта измерения т при идеальной в смысле качества воспроизведения АВД от динамических свойств сигнала [6, 7].
В случае восстановления сигнала при АВД экстраполяционным полиномом Тейлора Tn(t) n-ой степени и величине модуля допустимой погрешности воспроизведения S0 е Д15 принадлежащей области A1 таких значений, при которых на каждом участке аппроксимации длительности At производная сигнала (п+2)-го порядка x(n+2 )(t) = const,t е [tt _1 ,ti _1 +At], зависимость текущего такта измерения т = у э( 50 ,x(Hn+l) ,x(Hn+2}) от структурных (фазовых) свойств сигнала определяется решением двух уравнений [8,9]:
(n +1 )!-x[n+2 Vn+2) + (n + 2)!-x(Hn+1)i(n+1) _ 50 (n +1 )!(n + 2)! = 0 при 5m ^ 0; (1 а) (n +1 )!-x[n+2)x(n+2) + (n + 2)!-x[n+l)x(n+1) + 50(n +1 )!(n + 2)! = 0 при 5m < 0, (1 б) где x(Hn+l) и x(Hn+2) значения (п+1)-ой и (п+2)-ой производных в начале каждого участка аппроксимации; т - текущий такт измерения, при котором текущая погрешность 5(t._ _ t) = x(t._ _ t) _ Tn(ti_1 _ t) достигает своего максимального значения 5m = ±50 в конце каждого участка экстраполяции при t j _ t = т = At.
В общем виде уравнения (1) решить не представляется возможным. Поэтому ограничимся в отличие от [8] вторым распространенным на практике случаем - экстраполяция полиномом 1-й степени. Будем полагать, что любая точка фазовой траектории [x(3}(t),x(2*(t)] структурной модели
сигнала может быть началом участка аппроксимации длительности At = т.
Экстраполяция полшомом 1-й cmeneHU. Пусть на каждом интервале аппроксимации воспроизведение сигнала производится экстраполяционным полиномом Тейлора T1(t'),t' е At первого порядка (n = 1, рис. 5 б), модуль допустимой погрешности воспроизведения соответствует расширенной области ( 50 е A1 ). Требуется найти функциональную зависимость текущего
и
такта измерения от структурных свойств сигнала при линейной экстраполяции (ЛЭ).
При п = 1 из уравнений (1 а, б) имеем
3 3т3 + 3х(2 V - 650 = 0 при 5т > 0; [2 а)
33т3 + 2X23т2 + 650 = 0 при 5т < 0; [2 б)
или, вводя новую переменную 2 = т + х( 2 3 / х( 3 3
23 + 3р1 г + 2q1 = 0 при 5т > 0; [з а)
г3 + 2р2т2 + 2q2 = 0 при 5т -< 0, [з б)
фф X 23 )2 {х( 23 У 35 0 {х( 2 3) 350 где коэффициенты Р1 = р2 = -)-; ql = )-^ - —33; q2 = )-Ьг + -737■
1 2 .3> у2 1 .3> у3 х(3> 2 .3> у х(3>
Дискриминанты уравнений (2 а, б) соответственно равны
О = 4; + Р = 350 (- 2(х( 2' У + 35„ (*'33 ) )■(.'3' )-4;
Э2 = у22 + Р23 = 35 0 (2(.( 23 )3 + 380 X 33 У2\(х( 3' Г Функциональная зависимость такта измерения при АВД от структурных свойств сигнала т = у э1( 5 0 ,х(2 3 ,х( 3 3) определяется действительным,
положительным и наименьшим в случае неоднозначности решением (корнем) уравнений (2 а, б). Известно [10], что независимо от знака дискриминанта областью существования действительных корней вспомогательного кубического уравнения (3 а) или (3 б) является вся фазовая
плоскость Х(33ОХ(23. Однако число действительных корней k зависит от знака дискриминанта О (рис. 1 а, б, где f (х(23) = д/2(х(23 У3 ^350 и
/ (х(2= з(*'23 УУ350 ).
о) б)
Рис. 1. Области и число действительных корней уравнений (2 а, б) Решение каждого вспомогательного уравнения будем искать в наиболее удобной для анализа тригонометрической форме согласно соответствующей нашему случаю р = р2= р < 0) табл. 1 [10], где а - вспомогательная
величина и параметр г = ±^/[р, причем знак г должен выбираться
совпадающим со знаком q, что равносильно условию - ^(а) У 0 при D > О и
cos(а) > 0 при D < 0.
Таблица 1
р < 0
О = д2 + р3 < 0 О = + р3 У 0
cos( а ) = д/г3 ск( а) = д/г3
= -2 г • cos( а/ 3) г2 = -2г • cos(600 - а/3) z3 = +2г • о(600 + а/ 3) z1 = -2 г • сЩ а/3)
Обозначим г = ~лЦр = - х(23/х(33, г2 = +^f\p\ = х(23/х(33 и установим
область существования г\ и г2 на фазовой плоскости соответственно для каждого уравнения. Тем самым, очевидно, будет найдена область существования каждого действительного решения вспомогательных уравнений (3 а, б), поскольку их корни выражаются через параметр г.
В тригонометрической форме уравнение (3 а) представляется в виде:
и
д<0
ДМ)
cos(а3 3 = / г1 = ц -1; сЩа13 = / г1 = ц -1; cos(а4 3 = /г23 = 1 - ц; сЩа2 ^ = /г23 = 1 - ц,
где ц = 350 (х(33 )7(х(23 ^ .
Для случая 01 М 0 при х(23 0 и х(33 > 0 имеем ц 0,сН(а13 < 0 и
ск(а23 М 0, т.е. г = г2. При х(23 > 0 условию 01 М 0, когда х(33 -/2 (х(23)
и
х(33 м /2 (х(23) и х(23 м 0, соответствует условие ц м 2, при котором сИ(а13 м 0 и сИ(а23 < 0, т.е. г = г2.
Условию 01 < 0 (- /2(х(23 )< х(33 < /2(х(23), х(23 м 0) соответствует условие 0 < ц < 2. Отсюда следует:
1) cos(а3^ > 0 и cos(а43 < 0 при 1 <ц< 2, т.е. г = г1 для областей фазовой плоскости /4(х(23 )< х(33 < /2(х(23) и - /2(х(23 )< х(33 <-/4(х(23), где функция
/4(х(233 = р23?/350;
2) cos(а3^ <0 и cos(а4^ >0 при 0 <ц< 1, т.е. г = г2 при
- /4 .23)< X'33 < /4 X23).
Области существования параметров г1 и г2 уравнения (3 а) показаны на рис. 2 а стрелками.
V 4 ' ^
Рис. 2. Области существования действительных корней уравнений (2 а,б)
г
г
и
Уравнение (3 б) подобным образом можно записать в тригонометрической форме, где а' - вспомогательная величина уравнения (3 б),
В2 у0
д<0
ск(а1; = д2 /Г!3 =-( 1 + у); ^(а3) = г3 = -(1 + у); ск( а 2 ) = д2 / г23 = 1 + у; cos( а 4 ) = д2 / г23 = 1 + у.
Аналогичный анализ при О2 < 0 и О2 У 0 дает области существования г1 и г2 уравнения (3 б), представленные на рис. 2 б.
Совмещенная картина действительных и положительных решений х^ и хк уравнений (2 а) и (2 б) соответственно, полученная на основании параметров Г1 и г2 согласно табл. 1, приведена на рис. 3 с учетом их областей существования на фазовой плоскости.
Рис.3. Области действительных и положительных корней
уравнений (2 а, б)
Здесь можно указать четыре области ЫУ, в каждой из которых существует три действительных и положительных корня уравнений (2 а, б), удовлетворяющих очевидным неравенствам: область I - х3 - х1 - х1; область II - х3 - х2 - х1; область III - х3 - х2 - х1; область IV - х'3 — х[ — х1.
Выбирая в областях ЫУ с неоднозначным решением наименьший корень уравнений (2 а, б) получим в конечном итоге функциональную
г
г
и
зависимость уэ1 при АВД такта измерения х = у1э1 + у2э1 из уравнения [2 а) и х' = у1'э1 + у 2'э1 из уравнения [2 б) от структурных свойств сигнала для случая п = 1 и
А) область (х( 2 \х( 3; )еП1
У1э1 =
У1
х
(2)
X
(3;
у 2 =
х
(2)
X
у3 =
X
(3; (2)
X
[2^(а2/3) +1] при х(2) < 0, х(3) у /1 (х(2))-[2ch(а^3) -1] при х(2) > 0,х(3) у /2 (х(2)); [2 а3/3) -1] при х(2) у 0,/2 (х(2))< х(3) < /4 (х(2));
у2
э1
(2)
У 5
У 6
х
(3; (2)
у4 = [260° - а4/3) -1]при х(2) у 0, /4(х(2))< х(3) < 0; х
[2о(600 + а4/3) -1]при х(2) у 0,0 < х(3) < -/4 (х(2))-
х( 3 3 - х(2)
х
(3;
(4 а1)
(4 а2)
[2о(600 + а3/3) +1]при х(2) у 0, - /2 (х(2))< х(3) < -/4 (х(2)).
Б) область (х( 2 3,х( 3; )еП 2 (2)
у1 = ^ [2ch(а2/3) +1] при х(2) < 0, х(3) у /1 (х(2));
У2'э1 =
х1
у 2 =
х
(2)
х
у 3 =
х
(3;
(2)
х
(3;
у2'э1 =
у 4
х
(2)
х
у5 =
х
(3; (2)
[2ch(а1/3) -1] при х(2) > 0, х(3) у /2 (х(2)); [2 а'3/3) -1] при х(2) у 0, /2 (х(2))< х(3) < /4 (х(2))-
[2600 - а4/3) -1]при х(2) у 0, /4 (х(2))< х(3) < 0;
х
(3;
[2600 + а4/3) -1
.(2)
у6 =
х
(2)
(4 б1)
прих-> у 0,0 < х(3)<-/4(х(2)); (4 б2)
х
(3;
[2600 + а'3/3) +1]при х(2) у 0, - /2 (х(2))< х(3) < -/4 (х(2)),
и
где./, (х'2' )=,/- 2 (х'2' У/38 о ; / (х'2 3 М 2 (х'2' )7з80 ; / (х'2 > )=f(^
'38,
/4(х'2')^(х<23У/з8о ; ц = 38о(х'33)7(х'23У5;
а1 = атек'ц -1); а2 = атек' 1 - ц); а3 = arcos'ц -1); а4 = arcos' 1 - ц;
= атек'-ц-1); а2 = атек' 1 + ц); а3 = ат сов(-^-1); а4 = arcos' 1 + ц^. Область задания на фазовой плоскости функциональной зависимости у э1, соответствующей действительному, положительному и наименьшему в случае неоднозначности решению уравнений (2 а, б), показана на рис. 4. Очевидно, в данном случае подобно экстраполяции при п =0 [8] фазовая
плоскость Х(3Ох'23 разделяется кривой Dl0El, имеющей уравнение
^ х3' )=•
/1 (х'2' )^-(х'1' 1^/38
- /2 (х'2 3 )=^(х'2' 1^/38 (
при х
' 2 )
< 0
(5)
' 2 )
/ 2Ух 1= V //38 0 при х> 0,
на две неограниченных замкнутых области - область 0] решений уравнения (2 а) и область 02 решений уравнения (2, 6) (на рис. 4 заштрихована).
Рис.4. Функциональная зависимость у такта измерения т при АВД (п =13 В качестве примера для полигармонического сигнала (рис. 5 а) в среде МаШсаё получены: иллюстрация восстановления сигнала ЛЭ (рис. 5 б), дискретный фазовый портрет (рис. 6 а) и на основании (4) соответствующий ему текущий такт измерения на фазовой плоскости (рис. 6 б).
и
а)
Т1<4)
0 0.2
0.5
б)
Рис. 5. Полигармонический сигнал и его восстановление ЛЭ
а) а) б) б)
Рис. 6. Дискретный фазовый портрет и текущий такт измерения
На рисунке 6 приняты следующие обозначения: Рг 2, Рг 3 - векторы 2-ой и 3-ей производных, И - векторная функция кривой D10E1.
Заключение. Получена зависимость уэ1 (4 а, б) для множества перекрывающихся тактов измерения т(t) при АВД, ориентированной на воспроизведение сигнала экстраполяционным полиномом Тейлора первой степени (п = 1). На основании этих зависимостей можно решать задачу нахождения средней длительности такта измерения при АВД с последующей оценкой сжимаемости аналоговых сигналов.
Полученные соотношения позволяют также оценивать сжимаемость
случайных сигналов, полагая функцию х(0 реализацией дифференцируемого случайного процесса Х(0.
и
Литература
1. Авдеев Б.Я., Антонюк Е.М., Долинов С.Н., Журавин Л.Г., Семенов Е.И., Фремке А.В.Адаптивные телеизмерительные системы / Под ред. Фремке А.В. Л.: Энерго-издат. Ленингр. отд-ние, 1981. 248 с.
2. Дискретизация сигналов по времени (практика, алгоритмы) / Самойлов Л.К., Палазиенко А.А., Сарычев В.В., Ткаченко Г.И., Таганрог: Изд-во ТРТУ, 2000. 81 с.
3. Куревин, В.В., Морозов О.Г., Морозов А.Г. и др. Техническое обеспечение экологической безопасности территориально распределенных систем хранения опасных веществ // Инженерный вестник Дона. -2016.-№3.URL:ivdon.ru/uploads/article/pdf/IVD_79_Sinyutin.pdf_e8c1c28197.p df.
4. Кавчук С.В. Теоретические основы информационно-измерительной техники. Конспект лекций. Часть I. Таганрог: Изд-во ТРТУ, 2001. 125 с.
5. Qaisar, S.M., L.L. Fesquet and M.R. Laurent, 2009. Adaptive Rate Sampling and Filtering Based on Level Crossing Sampling. EURASIP Journal on Advances in Signal Processing, 2009(10.1155/2009/971656), 160 p.
6. Кавчук С.В., Ткаченко Г.И., Ткаченко М.Г. Оценка сжимаемости измерительных сигналов на основании априорных данных об их динамических свойствах // Естественные и технические науки. 2008. №3. С. 15-18.
7. Кавчук, С.В., Ткаченко Г.И., Савченко Я.С. Априорная оценка средней длительности такта измерения и числа отсчетов при адаптивной временной дискретизации // Известия ЮФУ. Технические науки. - 2014. - №4. - С. 147155.
8. Кавчук, С.В. Зависимость текущего такта измерения при адаптивной временной дискретизации экстраполяционного типа от структурных свойств
и
сигнала // Инженерный вестник Дона. - 2016. - №4. URL: ivdon.ru/ru/magazine/archive/n4y2016/3903.
9. Mark, J.W. and T.D. Todd, 1981. A nonuniform sampling approach to data compression. IEEE Transactions on Communications (issue 29), Institute of Electrical and Electronics Engineers, Inc, pp: 24-32.
10. Бронштейн Н.И., Семендяев К.А. Справочник по математике для инженеров и учащихся втузов. М.: Государственное издательство технико-теоретической литературы, 1957. 608 с.
References
1. Avdeev B.Ya., Antonyuk E.M., Dolinov S.N., Zhuravin L.G., Semenov E.I., Fremke A.V. Adaptivnye teleizmeritel'nye sistemy [Adaptive telemeasuring system] Pod red. Fremke A.V. L.: Energo-izdat. Leningr. otd-nie, 1981. 248 p.
2. Samojlov L.K., Palazienko A.A., Sarychev V.V., Tkachenko G.I. Diskretizacija signalov po vremeni (praktika, algoritmy) [Sampling of signals in time (practice, algorithms)]. Taganrog: Izd-vo TRTU, 2000. 81 p.
3. Kurevin V.V., Morozov O.G., Morozov G.A. i dr. Inzenernyj vestnik Dona (Rus). 2016. №3. URL: ivdon.ru/uploads/article/ pdf/IVD _79_ Sinyutin.pdf_e8c 1c28197.pdf.
4. Kavchuk S.V. Teoreticheskie osnovy informatsionno-izmeritel'noy tekhniki. Konspekt lektsiy. Chast' I [Theoretical bases of information-measuring equipment. Lecture notes. Part I]. Taganrog: Izd-vo TRTU, 2001. 125 p.
5. Qaisar, S.M., L.L. Fesquet and M.R. Laurent, 2009. Adaptive Rate Sampling and Filtering Based on Level Crossing Sampling. EURASIP Journal on Advances in Signal Processing, 2009(10.1155/2009/971656), 160 p.
6. Kavchuk S.V., Tkachenko G.I., Tkachenko M.G. Estestvennye i tekhnicheskie nauki. 2008. № 3. pp. 15-18.
7. Kavchuk S.V., Tkachenko G.I., Savchenko Ya.S. Izvestiya YuFU. Tekhnicheskie nauki. 2014. №4. pp. 147-155.
8. Kavchuk, S.V. Inzenernyj vestnik Dona (Rus). 2016. №4. URL: ivdon.ru/ru/magazine/archive/n4y2016/3903.
9. Mark, J.W. and T.D. Todd, 1981. A nonuniform sampling approach to data compression. IEEE Transactions on Communications (issue 29), Institute of Electrical and Electronics Engineers, Inc, pp: 24-32.
10. Bronshtejn N.I., Semendjaev K.A. Spravochnik po matematike dlja inzhenerov i uchashhihsja vtuzov [Handbook of mathematics for engineers and pupils of technical colleges]. M.: Gosudarstvennoe izdatel'stvo tehniko-teoreticheskoj literatury, 1957. 608 p.
