Спектрально-импульсный метод расширения динамического диапазона фазовой модуляции Текст научной статьи по специальности «Физика»

В.Н.Булатов

СПЕКТРАЛЬНО-ИМПУЛЬСНЫЙ МЕТОД РАСШИРЕНИЯ ДИНАМИЧЕСКОГО ДИАПАЗОНА ФАЗОВОЙ МОДУЛЯЦИИ

В статье приведены основные положения и выводы по линейным зависимостям между установленными участками фазового спектра импульсно-модулированных сигналов и изменениями формы этих сигналов - в области однозначного определения этих зависимостей на всей числовой оси.

Методы и результаты исследований ориентированы на использование их в области синтеза фазовых информационно-измерительных систем.

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

Вместе с тем, цель, обозначенная выше, может достигаться и за счет увеличения диапазона (размаха шкалы) изменения фазы. Анализ немногочисленных работ в этом направлении - главным образом, в области технических решений - наводит на мысль о том, что в данной области существует научная проблема, требующая своего разрешения.

При ближайшем рассмотрении этой проблемы очевидным является недостаток, совершенно несвязанный с природой существования фазовых зависимостей, а

обусловленный формой описания этих зависимостей.

Действительно, если сигнал e(t) имеет спектральную характеристику S(w) , то фазовый спектр определяется как j(w)= arc tg[Im S(w)/Re S(w)], где Im S(w) и Re S(w) соответственно мнимая и вещественная часть спектральной характеристики. При таком определении фазы любой спектральной составляющей имеет смысл говорить только в пределах [-л/2, л/2] . Очевидно, что выход фазы при соответствующем изменении аргумента модуляции за указанный предел вносит неоднозначность определения фазовой характеристики сигнала.

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

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

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

представленном в частотной области (подобный ансамбль далее именуется: импульсный динамический объект — ИДО). Преобразование Фурье, используемое для получения подобных образов, позволяет получать спектральную характеристику, которая не является функцией времени, и отражает поведение амплитудного |£(й)| и фазового д(й) спектров в частотной области только для конкретного конечного отрезка импульсной функции на выбранном интервале времени [^-Т, г,] (или сечении tc). Поэтому, в строгом соответствии с преобразованием Фурье, модуляционную функцию для фазы частотных компонент й, как функцию времени, можно представить только множеством фазо-частотных функций д(й,г1), определенных на интервале \_t -T, г+(к—1)Т] в к сечениях. Дальнейшее восстановление функции фазы как функции времени по полученному множеству (по существу — временных выборок) для отдельных частотных компонент является стандартной процедурой (фильтрация, аппроксимация) и не представляет особых проблем. Основной задачей настоящей работы является установление однозначно определяемых зависимостей д (й) =£ (р, й), где р — параметр ИДО, в потенциально неограниченном диапазоне изменения фазы установленных гармонических составляющих ИДО.

В \1] показано, что методически, для исследований свойств фазовых спектров сигналов типа ИДО, следует из фазового спектра д(й) выделять фазовый спектр дф(й), формируемый формой ИДО, и фазовый спектр дс(й), определяемый смещением г ИДО относительно выравненного ИДО

(выравненный ИДО существует при ^ \<Т/2, где Т — интервал полного определения ИДО). В соответствии с этим и определяются задачи установления модуляционных зависимостей фазы гармонических составляющих для обоих видов фазовых слагаемых.

Фазовый спектр, зависящий от смещения t , определяется как д(й, ) = -jйt и

см 7 * с4 ' см ^ см

однозначно определяется при любых значениях t .

см

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

Пусть функция в^) является сечением некоторой временной функции с

динамическими параметрами, отражающими медленный закон модуляции, и полностью определена на интервале \tc-T,tc] причем t-T=-Х/2, ?с=1/2. Такое расположение функции на оси t соответствует выравненной. И пусть сигнал е(0 удовлетворяет следующим условиям:

1) содержит четную и нечетную составляющие: е^), eнч(t) ;

2) раскладывается в ряд Маклорена.

В таком случае для исследования можно взять решение спектральной характеристики для обобщенного сигнала из [1]. Учитывая, что вещественная часть функции спектра, приведенного в [1], отображает четные параболические составляющие функции сигнала, а мнимая часть — нечетные, произведем в нем следующие замены:

а) для вещественной составляющей спектра: e(t)=eч(t);

б) для мнимой составляющей спектра:

e(t)=eнJt).

Тогда спектральная плотность сечения tc

на основании материалов в [1]:

5 (ю,*с) = — ю

- ■ ют , ют

А (ю.г --------------+ В(юл )сов

с 2 с 2

+ ]

- ■ ют , ют

С (ю,1с )вт —— + О (ю,1с) со^ —2—

(1)

где

А (ю,гс) = £ е, (2‘

г =0

В (ю,^) = Ё е, <2‘ +2)(0)£

2(г -к)

к =о [2(г - к)]!ю"

(т/2)2(г -к )+1

С (ю,гс) = -£е нч <2‘ +>)£

к=0 [2(г - к) + 1]!ю ’

(т/2)2(г -к)

к =0 [2(г - к )]!ю2*+1 ’

О (ю,гс) = £ е„ч (2‘ +>)£

(т/2)

2(г -к )+1

к=0 [2(г - к) + 1]!ю ■

Здесь eч(p) и e нч(ч) производные соответственно четной и нечетной

составляющих временной функции e(t) .

Если принять

А&) =-

ют ют

А (ю,гс ^т-^-+В (ю,гс )со^~2~

А 2(/с) =-ю

ют ют

С (ю,1с + О (ю,£с )со^-^

то выражение (1) для некоторой

+

фиксированного значения w=w I можно записать так:

S(wI,tc) = Aj(t) + jA2(t) .

Отсюда получаем выражение фазовой характеристики для частотной компоненты wI в традиционном ее представлении для рассматриваемого сечения tc:

j(WlJj = arc tg [A2(tJ/Ai(tJl (2)

которая может быть определена в смысле главного значения аргумента на интервале [-л/ 2,л/2] .

В общем случае функция фазы j=arctg(y) является ограниченной функцией. Первая ее производная — непрерывна и также ограничена:

dj 1

(3)

dt 1 + y

2 .

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

С этой целью определим частные производные (для упрощения записей формул опустим обозначение переменной t ):

1

- A 2

9A1 1 + (A 2/A1 )2 A12 A2 + A

dj 11 A

-= P (A„ A 2);

9A2 1 + (A2/Aj)2 A1 A2 + A22

= Q (A„ A 2);

. (4)

(5)

Соответственно полный дифференциал

«— =-- «А 1 +---- ¿А 2 = Р (А 1, А 2) ■ ¿А 1 + е (А >, А 2) ■ «А ^

ЭА1 ЭА 2

Теперь, чтобы найти - , необходимо вычислить:

j = Jdj = J P(A1; A 2) ■ dAj + Q (Aj, A 2)-dA2

(6)

Произведем проверку на независимость пути интегрирования по замкнутому контуру Ь по известному условию [2]. Для этого вычислим частные производные:

ЭР/ЭА2 = (А/-А/;/( А/-А// ,

Эе/ЭА, = (А22-А;2;/( а22-а;2;2 ,

откуда следует, что ЭР/ЭА2 = 3Q/9A; , то есть проверка дает положительный результат. Но этот же результат является признаком того, что интеграл по замкнутому контуру может быть равен нулю [2]. Это возможно, если в односвязной области D зависимости P(A1,A2) u Q (A1,A2) всюду будут непрерывны.

Исследование позволяет определить единственную точку в области D , где P(AfAJ u Q(Ai,A2) разрывны (особая точка):

— A

lim P (Ai,A 2) = lim-712 =—¥;

A2 ®0 A2®0 a 2

Ai®0

—A

lim Q (Ai,A 2) = lim ~at = ¥.

A2 ®0 Ai®0 Ai

Следовательно, если точка (A =0, A2=0) будет лежать внутри контура интегрирования, то интеграл (6) будет определять зависимость фазовой характеристики от A1 и А2согласно поставленной задаче. Для этого достаточно задаться нужной зависимостью А]=/(А2) или

Arf(Ai).

Примечательно, что если в качестве критерия брать линейность функциональной зависимости j от вновь вводимой переменной, то она реализуется следующей зависимостью:

A j = yjR2 - A 22 , R = const ,

где A] и A2 можно выразить, например, через x :

A} = R cos x,

(7)

A2 = R sin x .

(8)

Тогда, подставляя (7), (8) в (4) и (5), получим:

dj = cos x _ dj = sin x

,

R

3Aj R

Далее, на основании (7), (8) определяем: dA} = —Rxsin x.dx ; dA2 = Rxcos x.dx . Теперь находим зависимость j(x) :

j(x) = Jp(A1,A2) ■ dA1 + Q(A1,A2) ■ dA2 =

2

=j

- sin x 1/ 4 (cos x I/ -

(- R sin x dx ) +1------------ ( cos x dx

I R

Поскольку все шкалы в измерительных преобразованиях конечны, то достаточно взять х1 = 0, х2=х. Тогда, учитывая, что

Фг

рассматривается фазовая спектральная составляющая для частного сечения tc, окончательно можно записать:

Ф(Ч^ = X.

Таким образом, получена зависимость фазовой характеристики —(t¡) от параметра х^) , связывающего зависимость между А(с) и А2^), позволяющая однозначно определять фазу спектральной составляющей на всей числовой оси.

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

вышеизложенного анализа, которые обязательно или рекомендательно надо учитывать при разработке необходимых и достаточных ИДО.

1 Обязательные условия

1.1 Необходимо, чтобы выражение спектральной плотности Б (и, t¡) содержало вещественную А1 (и,tJ и мнимую уА2(ц,tc) части; иначе все комплексные амплитуды будут лежать или только в вещественной, или только в мнимой плоскостях. Это, в свою очередь, требует наличия в ИДО четной и нечетной динамической составляющих функции времени для любого сечения t .

1.2 В связи с тем, что функции Р(А1,А2) и Q(A1,A2) имеют разрывы при А1=А2=0 , и эта точка должна лежать внутри контура интегрирования, чтобы

(А1?А 2) ■ «А 1 + е(А1?А 2) ■ «А2 Ф 0,

то диапазоны изменения Яе Б (и,, tо), 1т Б (и.Л) от переменной t должны формироваться таким образом, чтобы эта точка была расположена внутри каждого из этих диапазонов.

1.3 Из частотного диапазона должны быть исключены участки, где спектральная плотность ИДО обращается в ноль, иначе пункт 1.2 не выполним.

2 Рекомендательные условия

Эти условия диктуются практической целесообразностью.

2.1 Ориентируясь на возможность высокоточной реализации управления функциями А^/^ ^) и А =/^,0 техническими средствами, составляющие ИДО аппликаты должны быть на интервале о их определения

элементарными, не имеющими производных.

2.2 Ансамбль импульсов, составляющих ИДО, должен иметь минимально возможное число аппликат; очевидно, что меньше двух аппликат (например, четная и нечетная) быть не может.

Исходя из составленных выше условий, было установлено [3], что аппликатой для построения ИДО должен быть импульс прямоугольной формы длительностью ?, а общий вид спектральной функции S(u,tc) на основании должен выглядеть следующим образом:

S (w,tc) =Wsin Wr[i(ffl,ic) + jA 2(a,tc)], (9)

где сомножители перед скобкой отражают характер спектра одиночного прямоугольного импульса с единичной амплитудой, а выражение в скобках представляет динамику соотношения четной и нечетной составляющих ИДО, представляющего собой ансамбль импульсов прямоугольной формы с вариацией их амплитуд, посредством которой осуществляется фазовая модуляция отдельных частотных компонент.

Очевидно, что выражение (9) поставленной цели не удовлетворяет при u=uo=±2kd/d (условие пункта 1.3).

Для обеспечения линейной зависимости фаз гармонических составляющих ИДО с частотами й. от некоторой модулирующего сигнала x(t--T) (выборка x по моменту начала сечения) в соответствии с (7), (8) должно выполняться

Ai(ojt,tc) = R(at) cos[kM.x(tc-T)], (10)

A2(rn,tc) = R(wt) sin[kM.x(tc-T)], (11)

RЦ) = VAi2(W ,tc) + A22(W ,tc) = const для w, (12)

где км —коэффициент масштаба. Тогда выражение (9) преобразуется в следующий вид:

2R(m i) mix

S(mi,tc)=--------sin - - exp [/(p§(mi,tc)] , 13)

mi 2

где

0i>(<$utc)=ku.x(tc-T)±{n/2-(n/2)sgn\sin(mT;/2)M<m}. (14) Примечание — Если в пределах сечения [tc-T, tj с допустимой погрешностью можно считать, что функция x не зависит от t, а

значения А1 и А2 при управлении меняются только от сечения к сечению: А1.=Г1(1с+1Т) и А2.=Г2(1с+1Т) — такой процесс модуляции классифицируется как процесс с медленным законом модуляции (МЗМ).

Уравнение (14) является главным уравнением фазовой модуляции. Но для того, чтобы реализовать это уравнение, нужно найти решения системы уравнений для функций амплитуд Е .(7 о) импульсов прямоугольной формы, входящих в выражения А1 (и), А2 (и) конкретных ИДО, с использованием (10)-(12). В результате анализа свойств ИДО, подходящих под условия

пунктов 1.1 и 1.2, были выявлены четыре базисных варианта построения ИДО, на основе изменения формы которых можно сформировать фазовые модуляционные уравнения в области определения фазового спектра [-¥ ,+ ¥ ]. Эти варианты схематично представлены таблицей 1, а их временные

диаграммы для частного сечения г, - на рисунках 1-4.

Для установленного некоторого множества { щк } были выявлены оптимальные уравнения модуляции фазы гармонических составляющих ИДО, которые в конечном виде сведены в таблицу 2.

Таблниа 2

Вари- ант Частоты вдспрон> б*денил *4 Гшн<4 и [гриыгжу зонные уранэеяня фазовой модуляции с нргун^кттн ^дутцчи *

А £ . = К со.чрсм.х}; Е; = Н(-¡) Д„.г]/2; К{'

Б {я/2±Ьг}/т £¡-9, {- 1/с <м [кн. х]-,ч с ?г{км. х ?}; Ез-О, У +4-гв[кыл}}{ К (а^) = ,а + Е\ = согШ

Б (1г±21ис)/1 <рс=-ки-*<-!? ± {тп2-ШН-1}'}: Е] = В. сох{к*.х]\ £з = Л(св(>=^1£|1+£1! = сею1

Г Е^0^Я{&£пШсо.1[кягх]-фф}ят[кь.х]}, Е^0У1Я{зёГ1(о^со1[кяг.х]+.-1^ф)зт[км.х]}, где &=Ф1{л/4±Ы2)> 0=см(п?4±кк/2); Л[а^) =т}Е / +Е] =<хяШ

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

а) реализацию нелинейного периодического преобразования типа Е=Ясо8(км.х) необходимо осуществлять посредством программируемых матриц и аналого-цифровых преобразователей; следовательно, величина х должна задаваться дискретным способом (например, числом в двоичном коде);

б) коэффициент к выбора частотной компоненты для модуляции по заданной модуляционной зависимости при некотором своем множестве вносит отрицательную зависимость в главное уравнение модуляции и смещение (в виде методической систематической ошибки); оба явления просто могут быть устранены перекоммутация-ми в зависимости от вариантов: Е1®Е2, Е2®Е1; Е1®-Е1, Е2®-Е2; в итоге всегда можно реализовать главное уравнение фазовой модуляции в виде:

^ф(юк,?с)=км.х(?о-Т) ; (15)

в) реализацию моментов 1 начала и конца

прямоугольных импульсов, входящих в ИДО, осуществлять в соответствии с отношением эквивалентности этих событий стробированным событиям периодического повторения меры времени, равной ? .

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

линейно в зависимости от х определяться (или воспроизводиться) на всей числовой оси.

Обозначенные свойства особенно важны при использовании указанного выше узла модуляции в том или ином виде, в основе которого реализуются установленные в таблице 2 зависимости, в фазовых системах управления и регулирования [4], а также в широкодиапазонных информационно-измерительных фазовых системах - с медленным законом изменения фазы (МЗМ).

Список использованной литературы

1 Булатов В.Н. Спектральная характеристика для обобщенного сигнала с динамическими параметрами// Анализ структур электронной и вычислительной техники: Меж-вуз. сб. научн. тр.- Оренбург, ОГТУ.- 1995.- С. 25-30.

2 Бермант А.Ф., Араманович И.Г. Краткий курс математического анализа. - М.:

Наука, Гл. ред. ф.-мат. лит-ры, 1971 - С. 483 - 487.

3 Булатов В.Н. Формирование фазовых модуляционных уравнений на основе импульсных динамических объектов // Методы и средства управления технологическими процессами: Сб. тр. III Междунар. научн. конф. - Саранск: Изд-во Мордов. ун-та, 1999. - С. 56 - 60.

4 Широкодиапазонный фазовый модулятор //Методы и средства управления технологическими процессами: Вторая Междунар. научн. конф.- Саранск: Изд-во Мордов. университета, 1997.- С. 61 - 62.

Статья поступила в редакцию 20.12.99г.