(Биквадратный) БИХ второго порядка цифровая фильтрация
y = sosfilt(sos,x)
y = sosfilt(sos,x,dim)
y = sosfilt(sos,x) применяет раздел второго порядка цифровой фильтр sos к векторному x. Вывод, y, является той же длиной как x.
Если любой вход к sosfilt является одинарной точностью, фильтрация реализована с помощью арифметики с одинарной точностью. Вывод, y, является одинарной точностью.
sos представляет раздел второго порядка цифровой фильтр H (z)
L-by-6 матрица, содержащая коэффициенты каждого раздела второго порядка в его строках.
Если x является матрицей, sosfilt применяет фильтр к каждому столбцу x независимо. Вывод y является матрицей, одного размера, содержа отфильтрованные данные, соответствующие каждому столбцу x.
Если x является многомерным массивом, фильтрами sosfilt по первому неодноэлементному измерению. Вывод y является многомерным массивом, одного размера как x, содержа отфильтрованные данные, соответствующие каждой строке и столбцу x.
Матрица разделов второго порядка, sos, входной сигнал, x, или оба могут быть двойной или одинарной точностью. Если по крайней мере один вход является одинарной точностью, фильтрация сделана с арифметикой одинарной точности.
y = sosfilt(sos,x,dim) управляет по измерению dim.
Загрузите chirp.mat. Файл содержит сигнал, y, который имеет большую часть его силы выше Fs/4 или половину частоты Найквиста. Частота дискретизации составляет 8 192 Гц.
load chirp
t = (0:length(y)-1)/Fs;
Разработайте 7-й порядок Баттерворт highpass фильтр, чтобы ослабить компоненты сигнала ниже Fs/4. Используйте нормированную частоту среза 0.48π рад/выборка. Выразите коэффициенты фильтра с точки зрения разделов второго порядка.
[zhi,phi,khi] = butter(7,0.48,'high');
soshi = zp2sos(zhi,phi,khi);
freqz(soshi)
Отфильтруйте сигнал. Отобразите исходные и highpass-отфильтрованные сигналы. Используйте ту же шкалу оси Y для обоих графиков.
outhi = sosfilt(soshi,y); subplot(2,1,1) plot(t,y) title('Original Signal') ys = ylim; subplot(2,1,2) plot(t,outhi) title('Highpass-Filtered Signal') xlabel('Time (s)') ylim(ys)
Разработайте фильтр lowpass с теми же спецификациями. Отфильтруйте сигнал и сравните результат с оригиналом. Используйте ту же шкалу оси Y для обоих графиков. Результатом является в основном шум.
[zlo,plo,klo] = butter(7,0.48); soslo = zp2sos(zlo,plo,klo); outlo = sosfilt(soslo,y); subplot(2,1,1) plot(t,y) title('Original Signal') ys = ylim; subplot(2,1,2) plot(t,outlo) title('Lowpass-Filtered Signal') xlabel('Time (s)') ylim(ys)
[1] Orfanidis, S.J., введение в обработку сигналов, Prentice Hall, Englewood Cliffs, NJ, 1996.
filter | medfilt1 | sgolayfilt