75599

ЦИФРОВАЯ ОБРАБОТКА КОРОТКИХ СИГНАЛОВ. ОПРЕДЕЛЕНИЕ ВРЕМЕННЫХ ИНТЕРВАЛОВ МЕЖДУ РАДИОИМПУЛЬСАМИ

Лекция

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

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

Русский

2015-01-15

189.5 KB

7 чел.

ОС.Лекция 12

ЦИФРОВАЯ ОБРАБОТКА КОРОТКИХ СИГНАЛОВ. ОПРЕДЕЛЕНИЕ ВРЕМЕННЫХ ИНТЕРВАЛОВ МЕЖДУ РАДИОИМПУЛЬСАМИ

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

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

Для улучшения отношения сигнал/шум в отраженном сигнале могут применяться различные методы, такие как временная фильтрация (накопление), оптимальная частотная фильтрация, прямое и обратное преобразование Фурье, корреляционный анализ.

Влияние шума в регистрируемом отраженном сигнале во временной области можно значительно уменьшать за счет многократного повторения измерения и синхронного накопления отраженных сигналов (см. рис. 9.1Б).

                                 А                                                                 Б

         

Рис. 9.1. Исходный (А) и накопленный (Б)  отраженный сигнал ЭМА толщиномера. Количество накоплений – 16.

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

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

for i=1:N;

if(i<(shift-h)|(i>shift+h))                            

  x(i)= 0;  

end    

if(i>(shift-h)&(i<=shift))

  x(i)= sin(2*3.14*f*(i-shift+h)/sf)*(exp((i-shift+h)/k2)-1);  

end                                          

if ((i>shift)&(i<(shift+h)))    

   x(i)=sin(2*3.14*f*(shift+h-i)/sf)*(exp((shift+h-i)/k2)-1);

end

где

f - частота несущей радиосигнала;

k2 – коэффициент, определяющий масштаб и крутизну фронта и спада сигнала и эталонов;

sf -  частота отсчетов при измерении;

h - полуширина зондирующего и отраженного радиоимпульса. В приведенном примере h=kt/6;

shift - величина сдвига

N - кол-во отсчетов при регистрации отраженного сигнала

          

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

Параметры: f=2660; k2= 300; sf=200000; h=170; shift=180, N=2048;

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

В [1] описан более совершенный способ измерения временных интервалов между двумя радиоимпульсами, в [2] - программа, реализующая этот способ.

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

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

                               А                                                                            Б

                                                                                                  

Рис. 9.3. Зависимость коэффициента корреляции от временного сдвига при отсутствии шума (А) и при отношении сигнал/шум 2/1 (Б). Точное значение времени сдвига равно 500. Вычисленное значение времени сдвига равно 499.998 (А) и 500.327 (Б ). Заданная точность вычисления – 0.001. Для построения графиков использована функция сплайн-аппроксимации spaps в MATLAB.

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

Краткое описание алгоритма.

1) Получение рабочего числового массива отраженного сигнала X[i], i =0,... n.

2) Определение номера элемента jмакс числового массива, соответствующего максимальному значению X . Найденная величина jмакс соответствует приближенному значению центра отраженного сигнала;

3) Создание 2k +1 эталонных числовых массивов по форме соответствующих зондирующему сигналу, центры которых смещены по отношению к центру  jмакс  на величину jмакс *(s/k-1)*b, где s =0, 1, 2,... 2k, b=0..1, значение b задается с учетом возможной погрешности приближенного значения временного сдвига. Например, b=0.1, если погрешность оценки приближенного значения временного сдвига не превышает 10%.  

4) Вычисление коэффициентов корреляции числового массива отраженного сигнала со всеми эталонными и формирование числового массива коэффициентов корреляции KK[m], где m = jmax  +  jmax *(s/k-1)*b.

5) Выполнение аппроксимации для массива KK[m] (Нахождение непрерывной функциональной зависимости F (m), соответствующей массиву KK[m]).

6) Выполнение передискретизации на основе найденной функциональной зависимости F (m) для массива KK[m] с увеличением количества элементов массива в R раз, т. е. формирование массива KK1[m1 ], где m1= jmax  +  jmax *(s1/(k*r)-1)*b, s1 = 0,1,2,…2k*r,  r – коэффициент передискретизации.

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

8) Вычисление разности err = m1макс − jмакс.

9) Если err < errдоп, то jмакс = m1макс и переход к п.3., иначе переход к п. 10.

10) Вывод найденного значения временного сдвига отраженного сигнала относительно зондирующего задержки Т =m1макс*dt  , где dt — шаг дискретности по времени при измерении отраженного сигнала.

Точность определения основной частоты при использовании предложенного алгоритма зависит от значений k и r и тем выше, чем они больше, однако, если анализируемый сигнал зашумлен, существенное ее увеличение происходит при росте k и r лишь до некоторого предела. В частности, при соотношении сигнал/шум >2/1, k =3 и r = 10, как можно установить путем моделирования, оказывается наилучшим выбором по критерию точность/время анализа.

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

На рис. 9.4 приведена экспериментально полученная в режиме компьютерного моделирования зависимость погрешности определения временного сдвига зашумленных эхо-сигналов  при частоте радиолокационного сигнала 2660 кГц и частоте отсчетов 200000 кГц. Генерация нормально распределенного шума выполнена с помощью функции rand в среде MATLAB с амплитудой, равной 0.5, при амплитуде эхо-сигнала, равной 1. Усреднение произведено по 20 реализациям, значение временного сдвига изменялось с шагом 1.

                         А                                                                  Б                                       

Рис. 9.4. Среднеквадратическое отклонение от среднего значения вычисленного временного сдвига при использовании традиционного (А) и описанного выше (Б) алгоритмов.

Максимальное среднеквадратическое значение погрешности определения временного сдвига при использовании описанного выше алгоритма, как следует из приведенного рис. 1.28, в этих условиях не превышает 0,65 ед. шкалы, а при использовании традиционного алгоритма - 4,5 ед. шкалы. С учетом диапазона шкалы измерения от 200 до 800 относительная среднеквадратическая погрешность измерения составит не более 0.3%  при использовании описанного выше алгоритма и 2.25% - при использовании традиционного алгоритма. Эффективность описанного выше алгоритма возрастает с увеличением зашумленности эхо-сигнала. Временной сдвиг эхо-сигнала может быть вычислен с высокой точностью даже в условиях, когда амплитуда шума превышает амплитуду эхо-сигнала в 1,5-2 раза. Обработка таких значительно зашумленных эхо-сигналов должна производиться при  b>>0.1 и k>>3.  При                         погрешность определения приближенного значения временного сдвига не имеет значения. На рис. 9.5 приведен пример зависимости коэффициента корреляции от временного сдвига при b=0.9, k=50, на этапе первой итерации (А) и последней (Б). Отношение сигнал/шум в эхо-сигнале – 1/2.

                  А                                                                Б

 

Рис.9.5. Зависимость коэффициента корреляции от временного сдвига при b=0.9, k=50, отношении сигнал/шум 1/2 на этапе первой итерации (А) и последней (Б).

Описанный алгоритм исследован в режиме компьютерного моделирования с нормальным (функция randn MATLAB) и   белым  шумом (функция awgn MATLAB). Программа обработки данных была реализована в среде MATLAB.

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

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

Коэффициент корреляции r, как известно, вычисляется как отношение коэффициента ковариации kxy к нормирующему множителю:

Коэффициент ковариации сигналов x и y вычислится по формуле:

Введем обозначение:, тогда окончательно получим:

Для получения нормирующего множителя вычислим

Введем обозначение в последнем выражении. В результате получим

Окончательно выражение для коэффициента корреляции будет иметь вид:

Именно эта зависимость имеет вид параболы, обращенной вершиной вверх, принимающей максимальное значение при  при любых (целых и нецелых) значениях n (рис. 9.6А).   

При знаменатель будет приблизительно равен ½ и выражение приближенно может быть записано как

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

Решение последнего тригонометрического уравнения также приводит к  при любых (целых и нецелых) значениях n (рис 9.6Б), что доказывает получение точного значения временного сдвига радиолокационных эхо-сигналов предлагаемым способом.   

                         А                                                                           Б

Рис. 9.6. Пример аналитической зависимости коэффициента корреляции (А) и первой производной (Б) от относительного временного сдвига при n=4.2.

Литература

1. В.С.Тутыгин Способ измерения временных интервалов между радиоимпульсами. Патент РФ на изобретение №2456632.

2. Тутыгин В.С.  Цифровая обработка коротких сигналов./В.С.Тутыгин - СПб.: Изд-во Политехн. ун-та, 2012. – 164с. ISBN 978-5-7422-3723-5.

 


 

А также другие работы, которые могут Вас заинтересовать

43161. ИСПОЛНИТЕЛЬНЫЙ МЕХАНИЗМ 4.44 MB
  Провести расчет элементов и параметров конструкции исполнительного механизма прибора комплекса ЛА с учетом указанных в задании системных особенностей. Введение 4 Расчет кинематических параметров 5 Выбор двигателя 5 Расчет мощности двигателя 5 Кинематический расчет редуктора 6 Определение передаточного числа 6 Выбор кинематической схемы и типа используемых зп 7 Расчет числа зубьев 7 Ошибка по скорости 8 Расчет КПД...
43162. Проектирование технических расчетов зон ТО, диагностики и ТР на примере подвижного состава автотранспортных предприятий 273.5 KB
  Автомобильный транспорт является наиболее массовым и удобным видом транспорта обладающим большой манёвренностью хорошей проходимостью и приспособленностью для работы в различных климатических и географических условиях. Техническое обслуживание ТО является профилактическим мероприятием проводимым в плановом порядке через определенные длительность пробега или срок работы подвижного состава. ТО1 и ТО2 включают контрольнодиагностические крепёжные...
43163. Водный транспорт леса 2.26 MB
  В данном курсовом проекте рассмотрен пример организации первоначального лесосплава, представляющий собой комплекс производственных и подготовительных работ, связанных с перемещением лесных грузов по водным путям. В проекте рассматриваются наиболее распространенные виды водной транспортировки леса - молевой лесосплав, сплав леса в пучках, плотах и в баржах. Также необходимо оптимальным образом подобрать технику и оборудование на технологических участках, что, в свою очередь, обеспечивало бы беспрерывность работы и снижало простой данного оборудования.
43164. Восстановление детали оси пульта управления автокрана К-64 и разработка технологической планировки кабино-жестяницкого участка завода по ремонту тракторов Т-130 374.5 KB
  Курсовой проект является завершающим этапом изучения дисциплины ремонт машин и оборудования позволяющим в ходе работы над ним углубить и закрепить умение и навыки более детально изучить вопросы восстановления детали в частности оси пульта управления автокрана К64 углубить и закрепить умение и навыки в разработке технологической планировки медницкорадиаторного участка завода по ремонту тракторов Т130. В настоящее время ремонт детали достаточно широко применяется в практике эксплуатации строительных машин что и делает тему...
43165. Тепловой расчет конвективной туннельной сушильной установки для зимнего и летнего режимов 1.72 MB
  Определяем по заданным температурам tол=20.4 Определяем влагосодержание do г кгс.5 Определяем энтальпию ho кДж кгс.6 Определяем плотность природного газа при нормальных условиях.
43166. Тепловой расчет конвективной туннельной сушильной установки для зимнего (январь) и летнего (июль) периода 1.57 MB
  Выполнить тепловой расчет конвективной туннельной сушильной установки, определить длительность сушки, размеры установки, выбрать вентилятор для подачи наружного воздуха, дымосос, циклон и сожигательное устройство, на основании следующих данных.
43167. ОСКОРБЛЕНИЕ КАК ИЛЛОКУТИВНЫЙ ЛИНГВОКУЛЬТУРНЫЙ КОНЦЕПТ 194 KB
  Научная новизна данной работы заключается в применении концептологического подхода к рассмотрению лингвистических проблем права и в историко-этимологическом описании социальных явлений, которые стали основой современного толкования концепта «оскорбление». В работе была исследована дискурсная реализация этого концепта и выделена типовая базовая структура иллокутивных концептов, объясняющая прагматическую природу лингвосоциальных явлений
43168. ОБРАЗ ШЕРЛОКА ХОЛМСА В ПРОИЗВЕДЕНИЯХ СЭРА АРТУРА КОНАН ДОЙЛА 248.5 KB
  При этом мироощущение и fin de siècle и неоромантизма было подчеркнуто инаким особенным что не устраивало консервативное викторианское общество и образ Шерлока Холмса начал меняться под воздействием такого экстралитературного фактора как цензура: образ редактировался и упрощался чтобы умещаться в строгие рамки жанра семейного чтения. Таким образом образ Шерлока Холмса прошел в процессе своего формирование через влияние fin de siècle и неоромантизма чтобы прийти к викторианским традиционным ценностям. В данной работе...
43169. Дистанционное зондирование Земли из космоса 412.04 KB
  Система правового регулирования ДЗЗ в России и в мире 9 Глава 1. Общие понятия ДЗЗ 9 Глава 2. Международноправовые акты регулирующие ДЗЗ 10 Глава 3. Ранее являясь исключительно прерогативой военных структур сегодня ДЗЗ решает множество гражданских задач и является крайне важной для обеспечения защиты окружающей природной среды разведки полезных ископаемых кадастрового учета и иных направлений деятельности.