способ измерения отсчетов характеристической функции

Классы МПК:G01R19/12 измерения скорости изменения 
Автор(ы):
Патентообладатель(и):Омский политехнический институт
Приоритеты:
подача заявки:
1989-04-18
публикация патента:

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

Формула изобретения

СПОСОБ ИЗМЕРЕНИЯ ОТСЧЕТОВ ХАРАКТЕРИСТИЧЕСКОЙ ФУНКЦИИ, основанный на усилении исследуемого сигнала, модуляции по фазе усиленным исследуемым сигналом опорного гармонического колебания, последующем искажении по форме фазомодулированного сигнала, перемножении искаженного фазомодулированного сигнала с опорным, усреднении результатов перемножения и индикации результатов усреднения, отличающийся тем, что, с целью повышения точности, перемножение осуществляют путем последовательного попарного перемножения фазомодулированного сигнала с опорным, после усреднения определяют сумму средних значений результатов перемноженных сигналов.

Описание изобретения к патенту

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

Среднее значение случайной величины ехр(iVx) называется характеристической функцией способ измерения отсчетов характеристической функции, патент № 2015517 (V) случайной величины х (Левин Б.Р. Теоретические основы статистической радиотехники. - М.: Сов. радио, 1969, с. 135):

способ измерения отсчетов характеристической функции, патент № 2015517(V) = Mспособ измерения отсчетов характеристической функции, патент № 2015517exp(iVx)способ измерения отсчетов характеристической функции, патент № 2015517 =способ измерения отсчетов характеристической функции, патент № 2015517 p(x)cosVxdx+iспособ измерения отсчетов характеристической функции, патент № 2015517 p(x)sinVxdx = A(V)+iB(V), где V - параметр характеристической функции, М{ способ измерения отсчетов характеристической функции, патент № 2015517} - знак математического ожидания; p(x) - плотность распределения вероятностей случайной величины х; А(V) и В(V) - соответственно действительная и мнимая части характеристической функции. Известно, Вешкурцев Ю.М. Статистические анализаторы флуктуаций фазы. Омск, 1979, с. 12), что, измерив оценки действительной А(V) и мнимой B(V) частей характеристической функции при различных значениях V, можно определить все одномерные статистические характеристики случайной величины х.

Известен способ последовательного измерения отсчетов действительной и мнимой частей характеристической функции (авт.св. СССР N 120606, кл. G 01 P 29/00) физических процессов, представленных в виде электрических сигналов. Для этого представляют исследуемый процесс в виде электрического сигнала х(t), усиливают его в число раз V, пропорциональное аргументу характеристической функции, т.е. получают сигнал Vx(t), модулируют сигналом Vx(t) фазу гармонического колебания высокой частоты U(t)=Umsin способ измерения отсчетов характеристической функции, патент № 2015517ot, т.е. получают колебание Uф(t)=Umsin[ способ измерения отсчетов характеристической функции, патент № 2015517ot+Vx(t)], перемножают полученный сигнал Uф(t) с колебаниями Uosin способ измерения отсчетов характеристической функции, патент № 2015517ot и Uocos способ измерения отсчетов характеристической функции, патент № 2015517ot, a продукты перемножения усредняют за время Т, например, с помощью интеграторов, получая при этом постоянные напряжения, пропорциональные действительной A(V)= M{ cosVx(t) и мнимой B(V)=M{ sinVx(t)} частям характеристической функции для аргумента, равного V. Индикаторы постоянного тока (напряжения) показывают измеренные значения мнимой и действительной частей характеристической функции. Повторяя указанную последовательность операций при разных значениях V, последовательно во времени получают ряд отсчетов A(V) и B(V) характеристической функции. Таким образом, для получения N пар отсчетов A(V) и B(V) требуется общее время измерения NT.

Описанный способ измерения отсчетов характеристической функции не позволяет одновременно получить N пар отсчетов A(V) и B(V), что приводит к невысокой скорости измерений и низкой статистической точности измерений. Статистическую точность измерений оценивают величиной среднеквадратического значения способ измерения отсчетов характеристической функции, патент № 2015517 случайной погрешности измерения. Поскольку при измерении отсчеты характеристической функции измеряют как усредненные за время Т значения процессов cosVx(t) и sinVx(t), то среднеквадратическое значение погрешности определения A(V) равно способ измерения отсчетов характеристической функции, патент № 2015517A= A/способ измерения отсчетов характеристической функции, патент № 2015517 , а среднеквадратическое значение погрешности определения B(V) равно способ измерения отсчетов характеристической функции, патент № 2015517B= B/способ измерения отсчетов характеристической функции, патент № 2015517 , где А и В - коэффициенты, зависящие от характеристик усредняемых процессов и усреднителей (Мирский Г.Я. Аппаратурное определение характеристик случайных процессов. М.: Энергия, 1972, с. 456. Следовательно, если задано общее время NT измерения N пар отсчетов A(V) и B(V), то в известном способе измерения получены значения способ измерения отсчетов характеристической функции, патент № 2015517A= A/способ измерения отсчетов характеристической функции, патент № 2015517 и способ измерения отсчетов характеристической функции, патент № 2015517B= B/способ измерения отсчетов характеристической функции, патент № 2015517 . При реализации параллельного принципа обработки информации одновременно определяют все N пар отсчетов A(V) и B(V) и, следовательно, при том же общем времени измерения NT получены значения способ измерения отсчетов характеристической функции, патент № 2015517A= A/способ измерения отсчетов характеристической функции, патент № 2015517 и способ измерения отсчетов характеристической функции, патент № 2015517B= B/способ измерения отсчетов характеристической функции, патент № 2015517 , которые в способ измерения отсчетов характеристической функции, патент № 2015517 раз меньше, чем аналогичные показатели в известном способе.

Наиболее близким по технической сущности к заявляемому является способ (авт. св. СССР N 171448, кл. G 01 P 19/12), позволяющий одновременно получать ряд значений действительной A(V) и мнимой B(V) частей характеристической функции. В указанном способе исследуемый случайный процесс преобразуют в электрический сигнал х(t), усиливают в V раз. Усиленным сигналом Vx(t) модулируют по фазе высокочастотное гармоническое колебание U(t)=Vmsin способ измерения отсчетов характеристической функции, патент № 2015517ot, т. е. получают колебание Uф(t)=Umsin[способ измерения отсчетов характеристической функции, патент № 2015517ot+Vx(t)]. Опорный U(t) и фазомодулированный Uф(t) сигналы в идентичных исказителях искажают по форме, т. е. преобразуют в вид

Uи(t) = Uспособ измерения отсчетов характеристической функции, патент № 2015517L(n)sinспособ измерения отсчетов характеристической функции, патент № 2015517otn и

Uф.и(t) = Uспособ измерения отсчетов характеристической функции, патент № 2015517L(n)sin[nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t).

Из искаженных по форме фазомодулированного и опорного колебаний выделяют с помощью набора узкополосных фильтров ряд гармоник UmL(n)sin[n способ измерения отсчетов характеристической функции, патент № 2015517ot] и UmL(n)sin[n способ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)], которые попарно перемножают. Продукты перемножения типа

способ измерения отсчетов характеристической функции, патент № 2015517 U2mL2(n)[cosnVx(t)-cos[2nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)] усредняют за время Т и получают ряд напряжений, пропорциональных значениям действительной части характеристической функции, т.е. ряд U2mL2(1)A(U); способ измерения отсчетов характеристической функции, патент № 2015517 U2mL2(2)A(2U);...; способ измерения отсчетов характеристической функции, патент № 2015517 U2mL2(n)A(nU); .... Для получения значений мнимой части характеристической функции опорное колебание предварительно сдвигают по фазе на четверть периода, т.е. получают колебание U1(t)= Umcos способ измерения отсчетов характеристической функции, патент № 2015517ot, затем его искажают по форме, т.е. получают колебание

U(t) = Uспособ измерения отсчетов характеристической функции, патент № 2015517L(n)cosnспособ измерения отсчетов характеристической функции, патент № 2015517ot, из которого с помощью набора узкополосных фильтров выделяют ряд гармоник UmL(n)cosn способ измерения отсчетов характеристической функции, патент № 2015517ot и перемножают с соответствующими гармониками искаженного фазомодулированного колебания. Продукты перемножения усредняют во времени и получают ряд напряжений, пропорциональных мнимой части характеристической функции способ измерения отсчетов характеристической функции, патент № 2015517 U2mL2(1)B(U);...; способ измерения отсчетов характеристической функции, патент № 2015517 U2mL2(n)B(nU);. ..Индикаторы постоянного тока (или напряжения) показывают измеренные значения мнимой и действительной частей характеристической функции. Таким образом, за время измерения Т получают N пар отсчетов действительной и мнимой частей характеристической функции, т.е. в N раз быстрее, чем в способе-аналоге.

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

Un(t)= UmL(n)sin[n способ измерения отсчетов характеристической функции, патент № 2015517 ot+nVx(t)] и благодаря этому одновременно получить ряд отсчетов характеристической функции, является фильтрация. Причем по логике известного способа выполнение операции фильтрации должно отвечать двум взаимоисключающим требованиям: в результате фильтрации гармоники искаженного сигнала должны четко разделяться (в противном случае из-за наложения спектров возникает погрешность), т.е. фильтрация должна быть узкополосной; операция фильтрации не должна искажать вероятностные характеристики флуктуаций фазы, а так как колебание, содержащее флуктуации фазы, имеет широкий спектр частот (Баскаков С.И. Радиотехнические цепи и сигналы. М.: Высшая школа, 1988), то фильтрация должна быть широкополосной. В результате противоречивости указанных требований в известном способе всегда присутствует методическая погрешность из-за искажения статистических характеристик выделяемых гармоник Un(t). Искажения вызываются рядом причин: наложением спектров, искажением закона распределения сигналов Un(t), а следовательно, и nUx(t), искажением энергетического спектра колебаний Un(t) после фильтрации.

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

Цель изобретения - повышение точности измерения отсчетов характеристической функции.

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

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

В известном способе (авт. св. СССР N 120606) измерение отсчетов мнимой и действительной частей характеристической функции производят последовательно во времени для разных значений аргумента V характеристической функции. Так как для определения вероятностных характеристик требуется набор отсчетов А(V) и B(V), измеренных при разных V, то этот способ обладает низким быстродействием и низкой статистической точностью измерений. Указанных недостатков лишен способ одновременного измерения ряда значений A(V) и B(V) (авт. св. СССР N 171448). В этом способе одновременное измерение ряда значений A(V), A(2V), A(3V),... обеспечивается тем, что из искаженного по форме фазомодулированного сигнала и из искаженного по форме опорного сигнала выделяют с помощью набора узкополосных фильтров ряд гармоник. Гармоники с одинаковыми номерами перемножают, а продукты перемножения усредняют и получают ряд отсчетов A(V), A(2V), A(3V)... Для получения ряда отсчетов B(V), B(2B), B(3V)... опорное колебание сдвигают по фазе на четверть периода и далее повторяют все вышеназванные операции. Таким образом, в названном способе основной операцией, позволяющей разделить гармоники искаженных по форме колебаний и, благодаря этому одновременно получить ряд отсчетов характеристической функции, является фильтрация, что приводит к появлению методической погрешности из-за наложения спектров, искажения закона распределения и энергетического спектра отдельных гармоник. В предлагаемом способе разделяют не гармоники искаженного по форме сигнала, а напряжения, пропорциональные отсчетам действительной или мнимой частей характеристической функции, так как на выходах усреднителей в предлагаемом способе имеет место суперпозиция напряжений, пропорциональных отсчетам или мнимой, или действительной части характеристической функции при разных значениях аргумента nV, где n=1,2,3,. . .. Операциями, позволяющими провести разделение отсчетов, являются операции последовательного перемножения и вычитания или суммирования. Таким образом, предлагаемый способ позволяет одновременно получить ряд отсчетов характеристической функции без использования операции фильтрации, что повышает точность измерения отсчетов характеристической функции.

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

Заявляемый способ заключается в следующем.

Исследуемый случайный процесс преобразуют в электрический сигнал х(t). Сигнал х(t) усиливают в V раз, т.е. получают сигнал Vx(t). Усиленным сигналом Vx(t) модулируют по фазе высокочастотное колебание U(t)= =Umsin способ измерения отсчетов характеристической функции, патент № 2015517ot, т. е. получают сигнал

Uф(t)=Umsin[ способ измерения отсчетов характеристической функции, патент № 2015517 ot+Vx(t)].

Фазомодулированный сигнал Vф(t) искажают по форме, т.е. преобразуют в сигнал вида

Un(t) = Uспособ измерения отсчетов характеристической функции, патент № 2015517L(n)sin[nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)].

Формируют ряд произведений искаженного по форме фазомодулированного сигнала и опорного сигнала U1(t)=Uи(t)U(t); U2(t)=Uи(t)U(t)U(t); U3(t)= Uи(t)U(t)U(t)U(t);...Un(t) = Uи(t) способ измерения отсчетов характеристической функции, патент № 2015517 ;...

Продукты каждого из произведений подвергают низкочастотной фильтрации, в результате чего определяют напряжения, пропорциональные средним значениям каждого n-го произведения Vn(t). При n четном указанные средние равны

Ucn= Unm+1 способ измерения отсчетов характеристической функции, патент № 2015517(-1)способ измерения отсчетов характеристической функции, патент № 2015517 L(n-2k)Bспособ измерения отсчетов характеристической функции, патент № 2015517V(n-2k)способ измерения отсчетов характеристической функции, патент № 2015517, а при n нечетном средние равны

Ucn= Unm+1 способ измерения отсчетов характеристической функции, патент № 2015517-1способ измерения отсчетов характеристической функции, патент № 2015517L(n-2k)Aспособ измерения отсчетов характеристической функции, патент № 2015517V(n-2k)способ измерения отсчетов характеристической функции, патент № 2015517.

Определяют разности напряжений: для n четного

Upn= способ измерения отсчетов характеристической функции, патент № 2015517Ucn-Unm+1 способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517(-1)способ измерения отсчетов характеристической функции, патент № 2015517L(n-2k)Bспособ измерения отсчетов характеристической функции, патент № 2015517V(n-2k)способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517= B(Vn), для n нечетного

Upn= способ измерения отсчетов характеристической функции, патент № 2015517Ucn-Unm+1 способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517(-1)способ измерения отсчетов характеристической функции, патент № 2015517L(n-2k)Aспособ измерения отсчетов характеристической функции, патент № 2015517V(n-2k)способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 = A(Vn), в результате чего имеют ряд напряжений, пропорциональных искомым отсчетам действительной А(Vn) и мнимой B(Vn) частей характеристической функции.

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

Структурная схема устройства, реализующего способ, содержит усилитель 1, вход которого является входом устройства, а выход подключен к управляющему входу фазового модулятора 2. Выход модулятора 2 присоединен к входу исказителя 3, выход которого подключен к первому входу первого перемножителя 4= 1, начинающего цепочку 4 последовательно включенных перемножителей. Выход каждого из перемножителей подключен к первому входу последующего перемножителя, а к вторым входам перемножителей цепочки 4 и второму входу фазового модулятора 2 подключен генератор 5. К выходам перемножителей присоединены усреднители 6, к выходам которых подключены суммирующие усилители 7, к выходу которых подключены индикаторы 8 постоянного тока (или напряжения). Причем выход i-го перемножителя через последовательно включенный i-й усреднитель соединен с первым входом i-го суммирующего усилителя, второй вход которого соединен с выходом первого усреднителя 6-1, если i - нечетное число, или с выходом второго усреднителя 6-2, если i - четное число, третий вход соединен с выходом третьего усреднителя 6-3, если i - нечетное число, или с выходом четвертого усреднителя 6-4, если i - четное число и т.д., последний 0(i/2)-й вход i-го суммирующего усилителя соединен с выходом (i-2)-го усреднителя (0(.) - функция округления числа, стоящего в скобках).

Фазовый модулятор может быть выполнен на базе дискретных элементов или микросхем (Колтик Е.Д. Фазосдвигающие устройства. М., 1981, с.164). Суммирующие усилители могут быть выполнены на базе дискретных элементов или микросхем (Применение прецизионных аналоговых микросхем./А.Г.Алексенко, Е.А. Коломбет, Г. И.Стародуб. - М.: Радио и связь, 1985). Исказитель может быть выполнен, например, как усилитель с отсечкой (Баскаков С.И. Радиотехнические цепи и сигналы. - М.: Высшая школа, 1988, с. 44). Усреднитель может быть выполнен на базе дискретных элементов (например, RC-цепочка) или микросхем (Мирский Г.Я. Аппаратурное определение характеристик случайных процессов. - М.: Энергия, 1972, с. 456). Перемножители могут быть реализованы на базе дискретных элементов или микросхем, например, серий 140, 525, 526 и т. д. (Тимонтеев В.Н., Величко Л.М., Ткаченко В.А. Аналоговые перемножители сигналов в радиоэлектронной аппаратуре. - М.: Радио и связь, 1982).

Устройство работает следующим образом.

Исследуемый процесс, представленный в виде электрического сигнала х(t), усиливают в V раз в усилителе 1. Выходным сигналом усилителя Vx(t) с помощью фазового модулятора 2 модулируют фазу гармонического колебания высокой частоты U(t)=Umsin способ измерения отсчетов характеристической функции, патент № 2015517ot, поступающего с генератора 5. Выходной сигнал фазового модулятора 2 Uф(t)=Umsin[способ измерения отсчетов характеристической функции, патент № 2015517ot+ Vx(t)] в исказителе 3 преобразуют в вид

Uф.и(t) = Uспособ измерения отсчетов характеристической функции, патент № 2015517L(n)sin[nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)].

Для получения отсчета мнимой части характеристической функции сигнала х(t) при аргументе, равном Um, где m - четное число, т.е. m=2,4,6..., перемножают m раз искаженный фазомодулированный сигнал Uфи(t) с сигналом генератора 5 U(t). Продукты перемножения

Unm(t) = Uфн(t)способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 = Uспособ измерения отсчетов характеристической функции, патент № 2015517L(n)sin[nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)]способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517Vmm способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517(-1)m/2-k Cmk cosспособ измерения отсчетов характеристической функции, патент № 2015517(m-2k)способ измерения отсчетов характеристической функции, патент № 2015517otспособ измерения отсчетов характеристической функции, патент № 2015517+Cmm/2/ 2способ измерения отсчетов характеристической функции, патент № 2015517, где Cmk = способ измерения отсчетов характеристической функции, патент № 2015517, подвергаются низкочастотной фильтрации, в результате чего из сигнала Unm(t) формируется его среднее значение, образуемое произведениями при n= m-2k, т.е. получаемое напряжение - среднее значение есть результат операции типа

Uусm= способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517 = способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 L(m-2k)CmkBспособ измерения отсчетов характеристической функции, патент № 2015517V(m-2k)способ измерения отсчетов характеристической функции, патент № 2015517(-1)m/2-k, где черта сверху обозначает усреднение; Кm - коэффициент передачи m-го усреднителя. Выбирая Um= 1B; Km=2Km-1(-1)m-1; K1=2, получают на выходе второго усреднителя 6=2 сигнал

Uyc2= L(2)B(2V), на выходе усилителя 2 с коэффициентом усиления 1/L (2) имеют сигнал

Uy2=B(2V).

На выходе четвертого усреднителя 6-4 получают сигнал

Uyc4= L(4)B(4V)-4L(2)B(2V), а на выходе суммирующего усилителя 7-4, выполняющего преобразование

(Uyc4+4Uyc2)/L(4), формируется сигнал B(4V),

На выходе шестого усреднителя 6-6 получают сигнал

Uyc6= L(6)B(6V)-6L(4)B(4V)+15L(2)B(2V), а на выходе суммирующего усилителя 7-6, выполняющего преобразование

(Uyc6+6Uyc4+9Uyc2)/L(6), формируется сигнал B(6V).

На выходе восьмого усреднителя 6-8 получают сигнал

Uyc8=L(8)B(8V)-8L(6)B(6V)+

+ 28L(4)B(4V)-56L(2)B(2B), а на выходе суммирующего усилителя 7-8, выполняющего преобразование

(Uyc8+8Uyc6+20Uyc4+16Uyc2)/L(8), имеют сигнал B(8V).

На выходе десятого усреднителя 6-10 получают сигнал

Uyc10=L(10)B(10V)-10L(8)B(8V)+

+ 45L(6)B(6V)-120L(4)B(4V)+

+ 210L(2)B(2V), а на выходе усилителя 7-10, выполняющего преобразование

(Uyc10+10Uyc8+35Uyc6+50Uyc4+

+ 25Uyc2)/L(10), формируется сигнал B(10V).

Для получения отсчета действительной части характеристической функции процесса х(t) при аргументе, равном pV, где p - нечетное число, т.е. p= 1,3,5. .., перемножают p раз искаженный сигнал Uфи(t) с сигналом генератора 5 V(t). Продукты перемножения

Unp(t) = Uфи(t)способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 = Um способ измерения отсчетов характеристической функции, патент № 2015517L(n)sin[nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)]способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517Upmспособ измерения отсчетов характеристической функции, патент № 2015517sinpспособ измерения отсчетов характеристической функции, патент № 2015517ot = Uспособ измерения отсчетов характеристической функции, патент № 2015517L(n)sin[nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)]способ измерения отсчетов характеристической функции, патент № 2015517Upm способ измерения отсчетов характеристической функции, патент № 2015517-1способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517Cpk способ измерения отсчетов характеристической функции, патент № 2015517 sin[(p-2k)способ измерения отсчетов характеристической функции, патент № 2015517ot] усредняются, в результате чего все высокочастотные слагаемые в сигнале Uпр(t) обращаются в ноль, а остаются произведения при n= p-2k. Поэтому на выходе p-го усреднителя получают сигнал (черта сверху обозначает усреднение)

Uуср= KpUpm+1 способ измерения отсчетов характеристической функции, патент № 2015517(p-2k)sin[(p-2k)способ измерения отсчетов характеристической функции, патент № 2015517ot+(p-2k)Vx(t)]способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517(-1)(p-1)/2-kспособ измерения отсчетов характеристической функции, патент № 2015517Cpk sin[(p-2k)способ измерения отсчетов характеристической функции, патент № 2015517ot] = способ измерения отсчетов характеристической функции, патент № 2015517-1способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517способ измерения отсчетов характеристической функции, патент № 2015517 L(p-2k)Cpk способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517Aспособ измерения отсчетов характеристической функции, патент № 2015517V(p-2k)способ измерения отсчетов характеристической функции, патент № 2015517, где Kp - коэффициент передачи p-го усреднителя. Выбирая Um=1B; Kp=2Kp-1(-1)p-1; K1=2, на выходе первого усреднителя 6-1 получают сигнал

Uyc1= L(1)A(V), а на выходе суммирующего усилителя 7-1 с коэффициентом усиления 1/L (1) получают сигнал

Uy1=A(V).

На выходе третьего усреднителя 6-3 получают сигнал

Uyc3= L(3)A(3V)-3L(1)A(V), а на выходе суммирующего усилителя 7-3, выполняющего преобразование

(Uyc3+Uyc1 3)/L(3), получают сигнал А(3V), являющийся отсчетом действительной части характеристической функции процесса х(t) при аргументе 3V.

На выходе пятого усреднителя 6-5 получают сигнал

Uyc5=L(5)A(5V)-5L(3)A(3V)+10L(1)A(V), а на выходе усилителя 7-5, выполняющего преобразование

(Uyc5+5Uyc3+5Uyc1/L(5), получают сигнал A(5V).

На выходе седьмого усреднителя 6-7 получают сигнал

Uyc7=L(7)A(7V)-7L(5)A(5V)+

+ 21L(3)A(3V)-35L(1)A(V), а на выходе усилителя 7-7, выполняющего преобразование

(Uyc7+7Uyc5+14Uyc3+7Uyc1)/L(7), получают сигнал A(7V).

На выходе девятого усреднителя 6=9 получают сигнал

Uyc9=L(9)A(9V)-9L(7)A(7V)+

+ 36L(5)A(5V)-84L(3)A(3V)+

+ 126L(1)A(V), а на выходе усилителя 7-9, выполняющего преобразование

(Uyc9+9Uyc7+27Uyc5+30Uyc3+

+ 9Uyc1)/L(9), получают сигнал A(9V).

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

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

Погрешность от наложения спектров. Фазомодулированный сигнал Uф(t)= Umsin[ способ измерения отсчетов характеристической функции, патент № 2015517ot+Vx(t)] в исказителе сигнала преобразуется в вид Uфн(t) = способ измерения отсчетов характеристической функции, патент № 2015517UmL(n)sin[nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)] , где L(n)Um - амплитуда n-й гармоники искаженного сигнала. Допустим определяются отсчеты характеристической функции для аргумента, равного V. Для этого из сигнала Vи(t) c помощью фильтра с центральной частотой способ измерения отсчетов характеристической функции, патент № 2015517о выделяют гармонику с номером n=1. В результате такой фильтрации получают колебание

Uфn(t) = [K(способ измерения отсчетов характеристической функции, патент № 2015517o)+способ измерения отсчетов характеристической функции, патент № 2015517K1(t)]UmL(1)sin[способ измерения отсчетов характеристической функции, патент № 2015517ot+Vx(t)]+способ измерения отсчетов характеристической функции, патент № 2015517UmL(n)[K(nспособ измерения отсчетов характеристической функции, патент № 2015517o)+способ измерения отсчетов характеристической функции, патент № 2015517

способ измерения отсчетов характеристической функции, патент № 2015517+ Kn(t)]sin(nспособ измерения отсчетов характеристической функции, патент № 2015517ot+nVx(t)], где K(способ измерения отсчетов характеристической функции, патент № 2015517 ) - коэффициент передачи фильтра на частоте способ измерения отсчетов характеристической функции, патент № 2015517; способ измерения отсчетов характеристической функции, патент № 2015517 Kn(t) - функция, описывающая паразитную амплитудную модуляцию, возникающую на n-ой гармонике фильтруемого колебания из-за широкого спектра фазомодулированного колебания и зависимости ослабления сигнала при фильтрации от частоты. Для получения искомых отсчетов характеристической функции продукты фильтрации Uфn(t) перемножают с колебаниями, полученными в результате фильтрации искаженных по форме опорного и сдвинутого на четверть периода опорного колебаний, т.е. с сигналами вида Us(t) = K(способ измерения отсчетов характеристической функции, патент № 2015517o)UmL(1)sinспособ измерения отсчетов характеристической функции, патент № 2015517ot+способ измерения отсчетов характеристической функции, патент № 2015517UmL(n)K(nспособ измерения отсчетов характеристической функции, патент № 2015517o)sinnспособ измерения отсчетов характеристической функции, патент № 2015517ot и

Uc(t) = K(способ измерения отсчетов характеристической функции, патент № 2015517o)UmL(1)cosспособ измерения отсчетов характеристической функции, патент № 2015517ot+способ измерения отсчетов характеристической функции, патент № 2015517UmL(n)K(nспособ измерения отсчетов характеристической функции, патент № 2015517o)cosnспособ измерения отсчетов характеристической функции, патент № 2015517ot.

Продукты перемножения усредняют во времени, в результате чего получают оценки действительной способ измерения отсчетов характеристической функции, патент № 2015517(V) и мнимой способ измерения отсчетов характеристической функции, патент № 2015517(V) частей характеристической функции:

способ измерения отсчетов характеристической функции, патент № 2015517(V) способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517K2(способ измерения отсчетов характеристической функции, патент № 2015517o)L2(1)U2mA(V)способ измерения отсчетов характеристической функции, патент № 2015517K2(nспособ измерения отсчетов характеристической функции, патент № 2015517o)L2(n)U2mA(nV)способ измерения отсчетов характеристической функции, патент № 2015517;

способ измерения отсчетов характеристической функции, патент № 2015517(V) способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517K2(способ измерения отсчетов характеристической функции, патент № 2015517o)L2(1)U2mB(V)способ измерения отсчетов характеристической функции, патент № 2015517K2(nспособ измерения отсчетов характеристической функции, патент № 2015517o)L2(n)U2mB(nV)способ измерения отсчетов характеристической функции, патент № 2015517.

Таким образом, в способе по авт.св. N 171448 за счет наложения спектров имеют относительную погрешность способ измерения отсчетов характеристической функции, патент № 2015517А определения А(V) (при условии K2( способ измерения отсчетов характеристической функции, патент № 2015517o)L2(1)Um2=2):

способ измерения отсчетов характеристической функции, патент № 2015517A= способ измерения отсчетов характеристической функции, патент № 2015517 и соответственно относительную погрешность определения B(V)

способ измерения отсчетов характеристической функции, патент № 2015517B= способ измерения отсчетов характеристической функции, патент № 2015517.

Пусть в качестве фильтра используется колебательный контур с центральной частотой способ измерения отсчетов характеристической функции, патент № 2015517о/2 способ измерения отсчетов характеристической функции, патент № 2015517 =128 кГц, К( способ измерения отсчетов характеристической функции, патент № 2015517о)=1, ширина энергетического спектра процесса х(t) способ измерения отсчетов характеристической функции, патент № 2015517 =3,5 кГц, половина полосы пропускания фильтра способ измерения отсчетов характеристической функции, патент № 2015517 =35 кГц. Тогда K2(n способ измерения отсчетов характеристической функции, патент № 2015517o)способ измерения отсчетов характеристической функции, патент № 2015517 1/(1+(((n-1) способ измерения отсчетов характеристической функции, патент № 2015517o-способ измерения отсчетов характеристической функции, патент № 2015517 )/ /способ измерения отсчетов характеристической функции, патент № 2015517)2).

Пусть в качестве исказителя используется усилитель с отсечкой, у которого

L(1) = способ измерения отсчетов характеристической функции, патент № 2015517 (способ измерения отсчетов характеристической функции, патент № 20155171-sinспособ измерения отсчетов характеристической функции, патент № 20155171cosспособ измерения отсчетов характеристической функции, патент № 20155171),

L(n) = способ измерения отсчетов характеристической функции, патент № 2015517 способ измерения отсчетов характеристической функции, патент № 2015517 где n=2,3..., способ измерения отсчетов характеристической функции, патент № 2015517 1 - угол отсечки (Баскаков С. И. Радиотехнические цепи и сигналы. - М.: Высшая школа, 1988), процесс х(t) гауссовский с А(V)=exp(- способ измерения отсчетов характеристической функции, патент № 2015517x2V2/2). Для таких характеристик фильтра при способ измерения отсчетов характеристической функции, патент № 2015517 V=1 и способ измерения отсчетов характеристической функции, патент № 20155171=90ополучают способ измерения отсчетов характеристической функции, патент № 2015517А=0,5%.

Искажение энергетического спектра колебаний в результате фильтрации.

Известно, что при прохождении случайных процессов через линейные системы (фильтры) происходит изменение вида энергетического спектра и уменьшение энергии процесса. Допустим, что в качестве разделяющих фильтров используются колебательные контуры, а корреляционная функция процесса Vx(t) экспоненциальная: Rx(способ измерения отсчетов характеристической функции, патент № 2015517 )=exp(- способ измерения отсчетов характеристической функции, патент № 2015517/ способ измерения отсчетов характеристической функции, патент № 2015517) способ измерения отсчетов характеристической функции, патент № 2015517x2 V2. В результате фильтрации дисперсия способ измерения отсчетов характеристической функции, патент № 2015517x2 V2 фазовых флуктуаций Vx(t) уменьшается и становится равной (Вешкурцев Ю. М. ). О влиянии коэффициента взаимной корреляции сигналов на точность оценок статистических характеристик фазовых флуктуаций.//Вопросы радиоэлектроники, сер. Техника радиосвязи, 1969, вып. 9, с. 12-23)

способ измерения отсчетов характеристической функции, патент № 2015517Vспособ измерения отсчетов характеристической функции, патент № 2015517= V2способ измерения отсчетов характеристической функции, патент № 20155172x/(1+способ измерения отсчетов характеристической функции, патент № 2015517/способ измерения отсчетов характеристической функции, патент № 2015517), где 2 способ измерения отсчетов характеристической функции, патент № 2015517 - полоса пропускания колебательного контура по уровню 3 дБ. Поэтому при измерении отсчетов характеристической функции гауссовского случайного процесса х(t) вместо истинной характеристической функции способ измерения отсчетов характеристической функции, патент № 2015517 (V) = exp (-способ измерения отсчетов характеристической функции, патент № 2015517x2V2/2) получают ее оценку способ измерения отсчетов характеристической функции, патент № 2015517 (V) = exp (-способ измерения отсчетов характеристической функции, патент № 2015517V2/2). При способ измерения отсчетов характеристической функции, патент № 2015517хV= 1, способ измерения отсчетов характеристической функции, патент № 2015517/способ измерения отсчетов характеристической функции, патент № 2015517 = 0,1 имеют способ измерения отсчетов характеристической функции, патент № 2015517 (V)=0,60653; способ измерения отсчетов характеристической функции, патент № 2015517 (V) = 0,66151 и, следовательно, относительная погрешность равна 9%.

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

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

Наверх