75599

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

Лекция

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

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

Русский

2015-01-15

189.5 KB

4 чел.

ОС.Лекция 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.

 


 

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

77973. Типизированные файлы с идентичными структурами 30.5 KB
  Типизированный файл – это файл в котором записаны идентичные структуры. Например любой файл может считать файлом байтов – т.е можно читать байт за байтом, можно перейти сразу к любому байту по его номеру, можно сразу узнать сколько байт в файле, можно заменить любой байт на другой не перезаписывая файл.
77974. Условный оператор 28 KB
  Встречаются следующие формы условного оператора: Условный оператор с одной ветвью if условие then команды end При выполнении такого оператора вычисляется условие и если оно истинно то выполняются команды до ключевого слова end в противном случае выполнение программы продолжается со следующей за условным оператором команды. Условный оператор с двумя ветвями if условие then команды else команды end Здесь при истинности условия выполняются команды при ложности команды.
77976. Компоненты ввода-вывода информации 125 KB
  Свойство IsMasked: Boolean доступно только для чтения и содержит True, если строка шаблона задана. Свойство EditText: string содержит текст до наложения на него маски шаблона (т. е. то, что ввел пользователь), а свойство Text: String может (в зависимости от шаблона см. ниже) содержать либо исходный текст, либо результат наложения на него маски шаблона.
77977. Палитра компонентов 271 KB
  Для этого используется специальный редактор окно которого появляется на экране после щелчка правой кнопкой мыши на любой пиктограмме в палитре компонентов и выбора опции properties Свойства. Окно редактора палитры компонентов Данное окно позволяет добавлять или удалять компоненты с палитры компонентов или переименовывать названия вкладок: кнопка dd позволяет добавить новую вкладку палитры компонентов; кнопка Delete предназначена для удаления вкладки палитры компонентов. На нем изображены шесть наиболее важных окон Delphi: главное...
77978. Возможности Delphi для ввода и отображения дат и времен. Таймер 193.5 KB
  Таймер Компонент Delphi Timer очень простой компонент который не виден на экране но тем не менее Delphi Timer выполняет очень важные функции в программе. Delphi Timer позволяет вводить необходимые задержки между выполнением тех или иных действий. Компонент Timer имеет всего четыре свойства и одно событие и работать с компонентом Delphi Timer очень просто. Свойство Назначение Enbled Включение-выключение таймера Intervl Интервал срабатывания в миллисекундах Nme Имя компонента в программе Tg Произвольный числовой параметр Помещаем...
77979. Графические файлы в Delphi 63 KB
  У ряда объектов из библиотеки библиотеки визуальных компонент есть свойство Cnvs канва которое предоставляет простой путь для рисования на них. Cnvs является в свою очередь объектом объединяющим в себе поле для рисования карандаш Pen кисть Brush и шрифт Font. Cnvs обладает также рядом графических методов: Drw TextOut rc Rectngle и др. Используя Cnvs вы можете воспроизводить на форме любые графические объекты – картинки многоугольники текст и т.
77980. Итерационные циклы 47 KB
  Для организации итерационных циклов используются операторы цикла с предусловием цикл ПОКА и цикла с постусловием цикл ДО. Эти операторы не задают закон изменения параметра цикла поэтому необходимо перед циклом задавать начальное значение параметра с помощью оператора присваивания а внутри цикла изменять текущее значение этого параметра. Циклы с предусловием используются тогда когда выполнение цикла связано с некоторым логическим условием. Оператор цикла с предусловием имеет две части: условие выполнения цикла и тело цикла.
77981. Кнопки. Диалоговые окна 67.5 KB
  Виды кнопок Кнопки TButton широко используются для управления программами представляет сабой командную кнопку на странице Stndrd. Определяет цвет стиль размер шрифта прилож Cncel: Boolen; Если имеет значение True событие OnClick кнопки возникает при нажатии клавиши Esc Defult: Boolen; Если имеет значение True событие OnClick кнопки возникает при нажатии клавиши Enter События OnClick Возникает при нажатии на кнопке В отличие от большинства других видимых компонентов кнопка TButton является компонентом самой Windows и...