75614

Цифровая фильтрация

Лекция

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

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

Русский

2015-01-15

152 KB

13 чел.

ОС. Лекция 8

Цифровая фильтрация

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

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

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

Фильтр скользящего среднего

Пусть мы имеем массив N значений измеренного сигнала, представленный в цифровой форме:

                                      {f1, f2, …fN}, i=1, 2, 3, …N

Для нахождения скользящего среднего в окрестности точки fi берем среднее арифметическое от K предыдущих и K последующих точек, включая и fi  . Таким же образом производим обработку для всех значений i. В результате вычисляем новый массив gi:

или

Точки из интервала 2К+1, используемые при вычислении скользящего среднего,  могут суммироваться с различными весовыми коэффициентами.

Целесообразность весового суммирования может быть вызвана тем, что, чем ближе суммируемая точка к точке с номером i, тем выше ее значимость, а с ее отдалением влияние ее на точку с номером i yменьшается.

Чтобы не исказить масштаб усредняемой функции, весовые коэффициенты нужно выбирать так, чтобы они удовлетворяли условию:

В качестве весовой функции обычно используется функция Гаусса, приведенная на рис. 1.7.

                        Рис. 1.1. Пример весовой функции

Цифровой низкочастотный фильтр

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

          

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

В программах обработки дискретизированных сигналов, представленных в форме числовых массивов, цикл

for k=1:N

x(k) = A*sin(2*pi*KP1*k/N);

end

создает KP1 периодов дискретизированного синусоидального сигнала в числовом массиве, содержащем N значений. Понятие частоты в этом случае отсутствует и появляется только в том случае, если задать шаг дискретности по времени.

Аналогично этому, в результате выполнения быстрого преобразования Фурье

I=1:N;

X1=fft(x,N);

мы получаем числовой массив X1, который будет содержать N элементов, причем информативной будет только первая половина массива, вторая будет зеркально отображать первую половину. В первой половине массива X1, содержащей N/2 элементов, будет представлен «частотный» спектр. Спектр будет отражать не частоту сигнала, установить которую по числовому массиву, представляющему сигнал во временной области, невозможно, а количество периодов. Т.е. в приведенном выше примере «всплеск» в массиве частотной области будет в элементе с номером KP1. Таким образом, массив частотной области укажет на количество периодов сигнала во временной области.

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

for  i=1:N     

H(i)=1/((1+j*i/NC));  

end

Здесь NC - полоса пропускания фильтра по уровню 0,7 амплитуды выражена в количестве отчетов спектра БПФ, пропускаемых фильтром. Остальные отсчеты в массиве частотного спектра будут ослабляться по амплитуде. Таким образом,  понятие постоянной времени фильтра, равно как и полосы пропускания, при дискретизированном представлении линейного фильтра отсутствует.

Ниже приведена программа фильтрации сигналов..

%Низкочастотный фильтр

A=20; %амплитуда сигнала

Q=5; %амплитуда шума

KP=12;% - количество периодов  сигнала

N=1024;%количество точек расчета

NC=12; %NC - полоса пропускания фильтра по уровню 0,7 амплитуды

% выражена в количестве отчетов спектра БПФ, пропускаемых % фильтром. Остальные отсчеты (в частотном спектре!) будут %ослабляться по амплитуде            

for k=1:N % генерация сигнала и шума

s(k) = A*sin(2*pi*KP1*k/N);

q(k)=Q*(randn(size(N))); %шум

x(k)=s(k)+q(k); % суммирование сигнала и шума

end

for i=1:N    

H(i)=1/((1+j*i/NC)); %передаточная функция простого фильтра

%в частотной области

end

i=1:N;

X1=fft(x,N); %частотный спектр сигнала с шумом

Z=ifft(X1.*H);   %свертка зашумленного сигнала с частотной

%характеристикой фильтра

 

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

                   А                                                                    Б

                   

                         В                                                                 Г

Рис.1.8. Исходный сигнал с шумом (А), его частотный спектр (Б), частотная характеристика фильтра в полулогарифмическом масштабе (В) и сигнал после фильтра (Г).


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

Недостатком «простого НЧ-фильтра» является наличие амплитудных, частотных и фазовых искажений сигнала в полосе пропускания.

Более высокое качество фильтрации обеспечивают фильтры Бесселя, Баттерворта, Чебышева.

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

Фильтр Баттерворта

где  - относительная частота среза, - частота среза.

      n – порядок фильтра.

Этим обеспечивается масимально плоская АЧХ в полосе пропускания.

Фильтр Чебышева

- постоянный коэффициент, определяющий степень неравномерности АЧХ в полосе пропускания,

Тn – полином Чебышева 1-го рода n-го порядка.

Этим обеспечивается максимальное подавление шумов вне полосы пропускания (за счет максимальной крутизны спада АЧХ фильтра) при некоторой, задаваемой коэффициентом , неравномерности АЧХ в полосе пропускания.

Фильтрация шумов с использованием прямого и обратного БПФ

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

%Программа БПФ. Сигнал содержит две частотные составляющие.

Fs = 1000;                    % Sampling frequency

T = 1/Fs;                     % Sample time

L = 1000;                     % Length of signal

t = (0:L-1)*T;               % Time vector  

N=1024;

              

% Sum of a 50 Hz sinusoid and a 120 Hz sinusoid

x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);

y = x + 2*randn(size(t));     % Sinusoids plus noise

Y = fft(y,N)/L;

f = Fs/2*linspace(0,1,N/2);

 

% Plot single-sided amplitude spectrum.

figure

plot(f,2*abs(Y(1:N/2)))

title('Single-Sided Amplitude Spectrum of y(t)')

xlabel('Frequency (Hz)')

ylabel('|Y(f)|')

На рис. 1.11 приведен результат работы программы.

 

Рис. 1.11. Зашумленный сигнал, представляющий сумму двух синусоидальных сигналов разных частот (50Гц и 120 Гц) (А) и частотный спектр, полученный с помощью БПФ.

Если теперь обнулить участки спектра от 0 до 40, от 55 до 110 и от 125 до 500, а затем выполнить обратное преобразование БПФ (ifft), то получим спектр и сигнал, представленные на рис. 1.12.                                  А                                                                     Б

 

Рис. 1.12. Частотный спектр после удаления спектральных составляющих шума (А) и восстановленный с помощью обратного преобразования Фурье сигнал (Б)

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

%Программа фильтрации помех с использованием БПФ-ОБПФ

Fs = 1000;                    % Sampling frequency

T = 1/Fs;                     % Sample time

L = 1000;                     % Length of signal

t = (0:L-1)*T;                % Time vector

% Sum of a 50 Hz sinusoid and a 120 Hz sinusoid

x = 1*sin(2*pi*50*t) + sin(2*pi*120*t);

y = x + 0.2*randn(size(t));     % Синусоиды+шум

%N = 2^nextpow2(L); % Next power of 2 from length of y

N=1024;

Y = fft(y,N)/L;%вычисление спектра сигнала с шумом

for i=1:N %"вырезание" спектральных составляющих шума

if (i>125)|(i<40)|((i>55)&(i<110))

   Y(i)=0;

end 

end

f = Fs/2*linspace(0,1,N/2);

%восстановление исходного сигнала после "вырезания"

%спектральных составляющих шума

z=ifft(Y);

f = Fs/2*linspace(0,1,N);

PAGE  7


 

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

79260. Значение совершенствования организации основного производства. Направления совершенствования организации основного производства 16.21 KB
  В условиях рыночных отношений обеспечить промышленному предприятию эффективное функционирование и конкурентные преимущества может только эффективная система управления его производственной деятельностью.Без четкой научно обоснованной организации процесса вряд ли можно добиться высоких показателей системы управления.Поэтому организация труда управленческого персонала во многом определяет эффективность всей системы управления производством как в текущем так и в оперативном режиме её функционирования.Формирование рациональной...
79261. Структура и взаимосвязь элементов системы организации производства. Предприятие как организационная система. Особенности, функции предприятия 17.23 KB
  Предприятие как организационная система. Особенности функции предприятия Предприятие как организационная система Системаобъективное единство закономерно связанных друг с другом предметов или целое состоящее из частей упорядоченных по определенному закону или принципу или Взаимосвязь взаимообусловленных сопокупность входящих в ее состав элемнтов Системный подход инструмент исследования систематизированный способ мышления согласно которому процесс обоснования решений базируется на определении общей цели и подчинении ей подсистем...
79263. Отраслевой состав национальной экономики. Особенности черной металлургии 14.01 KB
  Отраслевой состав национальной экономики. Особенности черной металлургии Отраслевая структура национальной экономики заключается в группировке хозяйствующих субъектов в однородные по своему составу группы связанные однородными функциональными характеристиками отрасли национальной экономики. Отраслевая структура национальной экономики проходит следующие этапы своего развития: 1 первый связан с активным развитием и преобладанием первичных отраслей экономики таких как сельское хозяйство добыча полезных ископаемых; 2 второй связан с...
79264. Характерные признаки предприятий, образующих металлургическую отрасль 14.74 KB
  На территориальную организацию черной металлургии влияет ряд факторов. Например в черной металлургии ФРГ этот показатель равен всего 22. Так на 1 т готового проката у нас приходится 85 человекочасов что в 15 2 раза больше чем в странах с развитой черной металлургией Южной Корее Бразилии Китае Тайвань. Низкая производительность труда на предприятиях черной металлургии ведет к тому что конкурентоспособность отечественной металлопродукции может быть обеспечена только при сохранении самого низкого из существующих уровня заработной...
79265. Основы организации и управления производством коммерческого предприятия 15.88 KB
  Основы организации и управления производством коммерческого предприятия Финансы коммерческих предприятий представляют собой экономические отношения возникающие в процессе формирования производственных фондов производства и реализации продукции образования собственных ресурсов привлечения внешних источников финансирования их распределения и использование. Финансами коммерческих организаций и предприятий присущи те же функции что и общегосударственным финансам распределительная и контрольная. Положительный финансовый результат...
79266. Методы организации производства и их технико-экономическая характеристика 15.09 KB
  Различают три основных метода организации производства: индивидуальный единичный партионный и поточный. При индивидуальном единичном методе организации производства выпуск продукции осуществляется единичными экземплярами или очень мелкими партиями изделий широкой номенклатуры. Равномерная работа на всех участках производства при партионном методе достигается за счет определения оптимальных размеров партий заделов длительности цикла графиков запускавыпуска.
79267. Понятие и особенности производственного процесса. Основные принципы организации производственных процессов 14.73 KB
  По операциям ведется увязка объемов работ на участках: нормирование оперативное планирование и учет производства Основные операции операции технологического цикла состоящие из операционных циклов в результате которых изменяются форма размеры свойства и взаимное расположение деталей Технологический процесс совокупность основных операций Вспомогательные операции операции обеспечивающие бесперебойность производственного процесса связанные с перемещением предметов труда с одного рабочего места на другое снятием и установкой...
79268. Классификация производственных процессов 53.5 KB
  Ручные процессы – осуществляются рабочими без помощи механизмов. Машинно-ручные процессы – выполняются машинами или механизмами при непосредственном участии в них рабочих – вскрытие чугунной летки с помощью ВЧП