Лабораторная работа №3 «Линейная фильтрация» Выполнена студентом гр. 3085/2 Рахманкулов А.. Цель

Лабораторная работа №3
«Линейная фильтрация»
Выполнена студентом гр. 3085/2 Рахманкулов А..
Цель
Изучить воздействие ФНЧ на тестовый сигнал.
Постановка задачи
Сгенерировать гармонический сигнал с шумом и синтезировать ФНЧ. Получить сигнал во временной и
частотной областях до и после фильтрации. Сделать выводы о воздействии ФНЧ на спектр сигнала.
Теоретическое обоснование
Линейный фильтр — динамическая система, применяющая некий линейный оператор ко
входному сигналу для выделения или подавления определённых частот сигнала и других функций по
обработке входного сигнала. Линейные фильтры широко применяются в электронике, цифровой обработке
сигналов и изображений, в оптике, теории управления и других областях.
Наиболее часто они используются для того, чтобы подавить нежелательные частоты входного сигнала или
для того чтобы выделить нужную полосу частот в сигнале. Существует большое количество различных
типов и модификаций линейных фильтров.
Фильтр ни́жних часто́т (ФНЧ) — один из видов аналоговых или электронных фильтров,
эффективно пропускающий частотный спектр сигнала ниже некоторой частоты (частоты среза), и
уменьшающий (подавляющий) частоты сигнала выше этой частоты. Степень подавления каждой частоты
зависит от вида фильтра.
Фильтр Баттерво́рта — один из типов электронных фильтров. Фильтры этого класса отличаются
от других методом проектирования. Фильтр Баттерворта проектируется так, чтобы его амплитудночастотная характеристика была максимально гладкой на частотах полосы пропускания.
Ход работы
1. Выполнение задачи в командном окне MATLAB
Построим низкочастотный гармонический сигнал y=sin(2*pi*f), (частота сигнала 1 Гц). Частоту
дискретизации установим равной 100 Гц. Будем использовать встроенную функцию BUTTER,
синтезирующую ФНЧ Баттерворта. В нашей работе будем использовать фильтр 16 порядка и
нормализированную относительно частоты Найквиста частоту среза равную 0.98 (частота Найквиста равна
половине частоты дискретизации, т.е.100Гц/2=50 Гц, тогда нормализованная частота среза равна 1(1Гц/50Гц)=0.98.
Построим сигнал и его спектр до и после фильтрации.
x = 0:0.01:4*pi;
f=100*(0:255)/512;
%исходный сигнал с шумом
figure
noise=rand(size(x));
yy = sin(2*pi*x); %гармонический сигнал
y = yy+0.3*noise; %шум
plot(x(1:200),y(1:200))
grid
%моделирование ФНЧ
[B,A] = BUTTER(16,0.98);
B=B./sum(B);
A=A./sum(A);
%обработка сигнала ФНЧ
figure
y2 = conv(y,[B,A]);
plot(x(1:200),y2(1:200))
grid
%спектр исходного сигнала с шумом
figure
s = fft(y,512);
ss = s.*conj(s)/512;
plot(f,ss(1:256))
axis([0 max(f) 0 2])
grid
%спектр отфильтрованного сигнала
figure
s2 = fft(y2,512);
ss2 = s2.*conj(s2)/512;
plot(f,ss2(1:256))
axis([0 max(f) 0 2])
grid
Рис.1 Гармонический сигнал с шумом y=sin(2*pi*f)
Рис. 3. Спектр сигнала y=sin(2*pi*f) с шумом.
Рис. 2. Сигнал с шумом после фильтрации.
Рис. 4. Спектр сигнала после фильтрации.
2. Выполнение задачи в среде Simulink
Аналогично п.1, смоделируем гармонический сигнал с шумом. Для синтеза ФНЧ будем
использовать стандартный блок Digital Filter Design. Построим сигнал и его спектр до и после фильтрации
по следующей схеме.
Рис.3 Схема simulink для визуализации исходного сигнала и его спектра до и после фильтрации
Результат работы
Рис. 5. Исходный сигнал с шумом в среде Simulink.
Рис. 7. Спектр сигнала с шумом в Simulink.
Рис. 6. Сигнал с шумом после фильтрации в Simulink.
Рис. 8. Спектр сигнала после фильтрации в Simulink.
Вывод
Несомненно, фильтрация позволяет убрать часть шума из сигнала, улучшив таким образом его
«читаемость». Однако, полностью устранить шум она не в состоянии, это обусловлено тем, что фильтр не
может разделить шум и сигнал в одной частотной области.