75598

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

Лекция

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

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

Русский

2015-01-15

140 KB

5 чел.

ОС. Лекция 11

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

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

Результатом БПФ дискретизированного эхо-сигнала спектрометра ЯМР определенной частоты является количество периодов сигнала во временном окне. Если частота отсчетов или интервал дискретности по времени при измерении сигнала известен, то по количеству периодов во временном окне можно установить и частоту измеряемого сигнала. Точность определения частоты в спектре входного сигнала вполне определенна и зависит от количества периодов p сигнала. Если количество периодов целое, то частота с помощью БПФ находится точно (даже при наличии значительной зашумленности сигнала). Если же количество периодов p не является целым, то появляется погрешность определения частоты. Максимальное значение погрешности равняется 1/p. В некоторых, практически важных случаях, например, при  обработке  эхо-сигналов  (рис. 1) импульсных спектрометров ЯМР, количество периодов анализируемого сигнала во временном окне принципиально ограничено величиной около 10 [1]. В этом случае погрешность определения частоты с помощью БПФ достигает 1/10, т. е. 10%.

Непосредственное использование БПФ  не позволяет  получить точное значение количества периодов во временном окне и частоты измеряемого сигнала в том случае, когда количество периодов не является целым, когда анализируемый сигнал занимает только часть временной области, модулирован по амплитуде и зашумлен. В практически важных случаях частотного анализа сигналов спектрометров ЯМР форма и начальная фаза эхо-сигналов известны, неизвестна и требует определения лишь частота. Это дает возможность создать эталонные сигналы, соответствующие ожидаемому эхо-сигналу по форме и начальной фазе, и произвести их корреляционное сравнение. Коэффициент корреляции эхо-сигнала с эталонным сигналом будет равен единице, если частоты эхо-сигнала и эталонного сигнала равны и эхо-сигнал не зашумлен. Поэтому при отсутствии шумов найти частоту эхо-сигнала можно производя корреляционное сравнение с эталонными сигналами, частоту эталонных сигналов подбирать до выполнения условия, когда коэффициент корреляции будет равен единице. Однако коэффициент корреляции уменьшается как при разнице частот эхо и эталонного сигналов, так и при совпадении частот, но из-за наличия шума. Поэтому таким способом определить частоту зашумленного сигнала невозможно.

Повышение точности  определения частоты сигнала может быть достигнуто за счет сочетания положительных качеств корреляционного подхода и БПФ.

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

Рис. 2. Зависимость коэффициента корреляции от частоты при отсутствии шума (А) и при отношении сигнал/шум 1/3 (Б). Точное значение частоты равно 1010. Аппроксимирующие кривые построены с помощью функции сплайн-аппроксимации spaps в MATLAB.

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

Допустим, в исходном массиве положение максимума было на частоте 1010, а после сплайн-аппрроксимации функции кросскорреляции и передискретизации положение максимума стало 1010.3.

Размер окрестности вблизи приближенного значения частоты, найденного, например, с помощью БПФ, должен быть выбран таким, чтобы истинное значение частоты находилось в пределах этой окрестности. Если, например, приближенное значение частоты мы получаем с погрешностью до 10%, то размер окрестности должен составлять не менее 10% от приближенного значения частоты.

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

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

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

Точность определения количества периодов и частоты сигнала описанным выше способом может быть оценена аналитически. Два гармонических сигнала: исходный и эталонный, имеющие частоты f1 и f2 или количество периодов n1 и n2 , отличающиеся по частоте, можно записать в виде:

x=sin(f1*t)    или  sin(2p n1t/T) = sin(2pn1fNt/FsT)=sin(2pn1i/N)

y=sin(f2*t)    или  sin(2p* n2*t/T)

где T – интервал наблюдения, f-частота отсчетов, N – количество отсчетов, Fs- частота дискретизации.

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

Окончательно

Именно эта функция имеет вид параболы, с вершиной, направленной вверх, если выполнено условие Dn<1/n.

     Рис.3. Зависимость коэффициентов корреляции от количества

     периодов сигнала и результат аппроксимации полиномом

     второго порядка. Точное значение количества периодов – 10,2.

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

График производной для частного случая, когда количество периодов сигнала равно 10.2, приведено на рис. 4.

Рис. 4. Зависимость производной  функции коэффициентов корреляции

от количества периодов сигнала. Точное значение количества периодов – 10,2.

Решение этого уравнения Dn=0 означало бы точное определение количества периодов и частоты с помощью корреляционного сравнения исходного сигнала с эталонами и последующей аппроксимации полученных значений коэффициентов кросскорреляции. При большом количестве периодов решением этого уравнения будет действительно Dn=0, т.к. в этом случае последнее уравнения упрощается и приходит к виду:

Если же количество периодов мало (например, n<10) и оно не целое, а сигнал зашумлен, то возникает погрешность определения количества периодов, связанная с погрешностью аппроксимации, которая тем не менее будет в несколько раз меньше, чем при использовании БПФ. В приведенном выше примере рис.3 вычислено функция кросскорреляции построена по 7 значениям коэффициентов корреляции эталонных сигналов с исходным и использована квадратичная аппроксимация. Из приведенного на рис. 3 аналитического выражения аппроксимирующей функции можно вычислить предсказанное значение количества периодов, приравняв производную этой функции нулю:

Отсюда x=77.676/7.6146=10.2009 (точное значение равно 10.2000)

С помощью БПФ было бы получено значение x=10.


Текст базовой программы

%Комбинированное использование ключевых операций ЦОС

%Для повышения точности определения количества периодов и частоты

%"короткого" сигнала используется комбинация

%БПФ, кросскорреляции, сплайн-аппроксимации,передискретизации

kt=1024; % количество отсчетов

f=12; %частота сигнала

dt=2;%шаг дискретности по времени при измерении

k2=500; %декремент затухания огибающей сигнала

shum=0;%шум

kp=11.2;%количество периодов сигнала

%1. ГЕНЕРАЦИЯ МОДЕЛЬНОГО СИГНАЛА

for i=1:kt %обнуление массива сигнала

   y(i)=0;

end    

for i=1:kt %генерация модельного сигнала

   if(i>0)&(i<=kt/2)             

      y(i)=sin(2*3.14*kp*i/kt)*i*exp(i/k2)/1200;

   end

   if(i>kt/2)&(i<(kt))             

      y(i)=sin(2*3.14*kp*i/kt)*(kt-i)*exp((kt-i)/k2)/1200;

   end

    y(i)=y(i)+shum*(2*rand(1)-1);                                               

end

i=1:kt; %отображение модельного сигнала во временной области

figure

plot(i,y);

axis tight;

title('Time domain')

xlabel('Sample number')

%2. ФУНКЦИОНАЛЬНОЕ ПРЕОБРАЗОВАНИЕ (БПФ)

bpfy=fft(y,kt);%БПФ

bpf=bpfy.*conj(bpfy)/kt;%БПФ

%f=1000*(0:256)/512;

figure

plot(i(1:257),bpf(1:257));

axis tight;

title('Frequency domain')

xlabel('frequency')

%нахождение макс. знач. функции БПФ для массива Y

C=max(bpf);

for i=1:kt %поиск количества периодов, соответствующих максимуму БПФ

   if (bpf(i)==C)          

       kpbpf=(i-1);         

       break

   end

end

kp_bpf=kpbpf

%3. СОЗДАНИЕ ЭТАЛОНОВ И КРОССКОРРЕЛЯЦИЯ

fr=kpbpf;

   procch=0.1;%область поиска в процентах относит. kp_bpf    

shagkor=fr*procch/3;%шаг поиска

       k=0;

 for iii=fr-fr*procch:shagkor:fr+fr*procch %цикл для создания 6 эталонов в окрестности приближенного

           %значения количества периодов, определенных с помощью БПФ.           

           k=k+1;

           xkor(k)=iii;

           kor(k)=0;

           for i=1:kt           

               x(i)=0;    

           end

%Вычисление массивов эталонных сигналов              

       for i=1:kt                                                

           if(i>0)&(i<=kt/2)                       

               x(i)=sin(2*3.14*iii*i/kt)*i*exp(i/k2)/1200;

           end

           if(i>kt/2)&(i<(kt))                        

               x(i)=sin(2*3.14*iii*i/kt)*(kt-i)*exp((kt-i)/k2)/1200;

           end

       end    

%вычисление средних значений модельного и эталонных сигналов

           x_sr=mean(x);

           y_sr=mean(y);

           x_sko=0;

           y_sko=0;

%вычисление СКО и коэф. корреляции модельного и эталонных сигналов

           for i=1:kt

               x_sko=x_sko+(x(i)-x_sr)*(x(i)-x_sr);

               y_sko=y_sko+(y(i)-y_sr)*(y(i)-y_sr);

               kor(k)=kor(k)+(x(i)-x_sr)*(y(i)-y_sr);

           end

           kor(k)=kor(k)/(sqrt(x_sko*y_sko));

 end %конец цикла создания эталонов и вычисления массива коэф. корр.

 %СПЛАЙН-АППРОКСИМАЦИЯ И ПЕРЕДИСКРЕТИЗАЦИЯ

       xx=1:k;

       xi=1:0.1:k;

       r1=sin(xx); %только для тестирования сплайн-аппроксимации       

       yint=interp1(xx,kor,xi,'spline');% сплайн-аппроксимация коэф корреляции               

       r1=kor;

       %%apr=csaps(xx,r1);

       apr=spaps(xkor,kor,0.000001);

        fnplt(apr)

        hold on

        plot(xkor,r1,'ro')

        hold off  

%НАХОЖДЕНИЕ УТОЧНЕННОГО ЗНАЧЕНИЯ КОЛИЧЕСТВА ПЕРИОДОВ СИГНАЛА

       cmax=max(yint); %нахождение максимума коэф. корр.       

       for i=1:round((k-1)/0.1+1)

           if (yint(i)==cmax)     

               kp_int=fr-fr*procch+(i-1)*shagkor/10 %уточненное значение частоты по МАХ функции коэф. корр.

           end           

       end

pause;

close all;

                          

Алгоритм определения частоты

Литература 

    

  1.  В.И.Тарханов, В.С.Тутыгин Приборный комплекс для поиска и исследования сигналов ЯМР в магнитоупорядоченных веществах. // «Научное приборостроение», 2003, том 13, №1, с.58-63.
  2.  В.С.Тутыгин, А.В.Дебелова. Новый метод частотного анализа. // Вестник СПбГУ, Серия 10 «Прикладная математика, информатика и процессы управления», №1, 2009.
  3.  В.С.Тутыгин. Способ измерения частоты сигнала. Патент РФ на изобретение №2478213.
  4.  Тутыгин В.С.  Цифровая обработка коротких сигналов./В.С.Тутыгин - СПб.: Изд-во Политехн. ун-та, 2012. – 164с. ISBN 978-5-7422-3723-5.


 

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

45931. Типы зажимных устройств приспособлений. Краткая характеристика по составу, типу производства 12.18 KB
  По составу зажимные устройства делят на группы. 1Зажимные устройства состоящие из силового механизма и привода который обеспечивает перемещение контактного элемента и создаёт исходное усилие преобразуемое силовым механизмом в зажимное усилие. 2Зажимные устройства в котором силовой механизм приводится в действие рабочим прилагающим исходное усилие на орпеделёное плечё.Такие зажимные устройства с ручным приводом.
45932. Правила определения силы зажима заготовок в приспособлении 2.1 MB
  Для этого составляют расчетную схему где изображают все действующие силы и моменты резания зажимного усилия реакции опор и силы трения в местах контакта заготовки с опорными элементами и зажимными устройствами. По этому уравнению выводят формулу для расчета силы зажима Пример: расчетная схема на фрезерные операции. условий применительно к которым рассчитывались силы и моменты резания то их надо увеличить введением коэффициента запаса надежности закрепления согласно требованиям безопасности.
45933. Приводы зажимных устройств 1.73 MB
  Недостатки: незначительная плавность перемещения рабочих органов особенно при переменой нагрузке; низкое давление воздуха 04 мПа обуславливающие большие размеры приводов для приложения значительных усилий. на всех производственных участках подаётся воздушная среда давлением до 1МПа. Пневмоприводы рассчитываются на прочность при Р=06мПа а исходное усилие определяется при р=04МПа. Испытания их осуществляют при р не менее 09МПа.
45934. Цели, принципы, функции и основные задачи стандартизации 16.4 KB
  В соответствии с Федеральным Законом О техническом регулировании стандартизация осуществляется в целях: повышения уровня безопасности жизни или здоровья граждан имущества физических или юридических лиц государственного или муниципального имущества экологической безопасности безопасности жизни или здоровья животных и растений и содействия соблюдению требований технических регламентов; повышения уровня безопасности объектов с учетом риска возникновения чрезвычайных ситуаций природного и техногенного характера; обеспечения...
45935. Основные понятия в области метрологии. Метрология. Измерение. Погрешности измерения. Средство измерения. Единство измерений. Проверка средств измерений 18.03 KB
  Единство измерений. Проверка средств измерений.Рассматривает общие теоретические проблемы разработка теории и проблем измерений физических величин их единиц методов измерений.Устанавливает обязательные технические и юридические требования по применению единиц физической величины методов и средств измерений.
45936. Погрешности средств измерений. Систематическая погрешность средств измерений. Случайная погрешность средств измерений. Абсолютная, относительная погрешность. Точность средств измерений. Класс точности средств измерений 12.85 KB
  Погрешности средств измерений. Систематическая погрешность средств измерений. Случайная погрешность средств измерений. Точность средств измерений.
45937. Эталоны единиц физической величины. Эталон еденицы физической величины. Поверочная схема для средств измерений. Рабочий эталон. Вторичный эталон. Международный эталон 12.86 KB
  Эталоны единиц физической величины. Эталон еденицы физической величины. Рабочий эталон. Вторичный эталон.
45938. Средства измерительной техники. Средство измерений. Автоматичесое средство измерений. Автоматизированное средство измерений 12.24 KB
  Средство измерений. Автоматичесое средство измерений. Автоматизированное средство измерений. Средства измерительной техники измерительная техника – обобщающее понятие охватывающее технические средства специально предназначенные для измерений.
45939. Классификация размерных цепей. Основные термины и определения. Метод расчета размерных цепей, обеспечивающие полную взаимозаменяемость 35.97 KB
  Размерные цепи отражают объективные размерные связи в конструкции машины технологических процессах изготовления ее детали и сборки при измерении возникающие в соответствии с условиями решаемых задач. Обозначаются размерные цепи прописными буквами русского алфавита и строчными буквами греческого алфавита кроме . Размеры образующие размерную цепь называют звеньями размерной цепи. Одно звено в размерной цепи замыкающее исходное а остальные – составляющие.