decimate
Децимация — уменьшает частоту дискретизации целочисленным коэффициентом
Описание
Примеры
Десятикратно уменьшите сигнал
Создайте синусоидальный сигнал, произведенный на уровне 4 кГц. Десятикратно уменьшите его фактором четыре.
t = 0:.00025:1; x = sin(2*pi*30*t) + sin(2*pi*60*t); y = decimate(x,4);
Постройте исходные и подкошенные сигналы.
subplot 211 stem(0:120,x(1:121),'filled','markersize',3) grid on xlabel 'Sample number',ylabel 'Original' subplot 212 stem(0:30,y(1:31),'filled','markersize',3) grid on xlabel 'Sample number',ylabel 'Decimated'
Десятикратно уменьшите сигнал Используя чебышевский фильтр
Создайте сигнал с двумя синусоидами. Десятикратно уменьшите его фактором 13 использований Чебышевского БИХ-фильтра порядка 5. Постройте исходные и подкошенные сигналы.
r = 13; n = 16:365; lx = length(n); x = sin(2*pi*n/153) + cos(2*pi*n/127); plot(0:lx-1,x,'o') hold on y = decimate(x,r,5); stem(lx-1:-r:0,fliplr(y),'ro','filled','markersize',4) legend('Original','Decimated','Location','south') xlabel('Sample number') ylabel('Signal')
Исходные и подкошенные сигналы имеют соответствие с последними элементами.
Десятикратно уменьшите сигнал Используя КИХ-фильтр
Создайте сигнал с двумя синусоидами. Десятикратно уменьшите его фактором 13 использований КИХ-фильтра порядка 82. Постройте исходные и подкошенные сигналы.
r = 13; n = 16:365; lx = length(n); x = sin(2*pi*n/153) + cos(2*pi*n/127); plot(0:lx-1,x,'o') hold on y = decimate(x,r,82,'fir'); stem(0:r:lx-1,y,'ro','filled','markersize',4) legend('Original','Decimated','Location','south') xlabel('Sample number') ylabel('Signal')
Исходные и подкошенные сигналы имеют соответствие с первыми элементами.
Входные параметры
Выходные аргументы
Алгоритмы
Децимация уменьшает исходную частоту дискретизации последовательности к более низкому уровню. Это - противоположность интерполяции. decimate lowpass фильтрует вход, чтобы принять меры против искажения и прореживает результат. Функция использует алгоритмы децимации 8.2 и 8.3 от [1].
decimateсоздает фильтр lowpass. Значением по умолчанию является Чебышевский Тип, я фильтрую спроектированное использованиеcheby1. Этот фильтр имеет нормированную частоту среза0.8/rи неравномерность в полосе пропускания 0,05 дБ. Иногда, заданный порядок фильтра производит искажение полосы пропускания из-за ошибок округления, накопленных от сверток, должен был создать передаточную функцию.decimateавтоматически уменьшает порядка фильтра, когда искажение заставляет ответ величины в частоте среза отличаться от пульсации больше чем 10-6.Когда
'fir'опция выбрана,decimateиспользованиеfir1спроектировать КИХ-фильтр lowpass с частотой среза1/r.При использовании КИХ-фильтра,
decimateфильтрует входную последовательность только в одном направлении. Это сохраняет память и полезно для работы с длинными последовательностями. В БИХ-случае,decimateприменяет фильтр в прямых и обратных направлениях с помощьюfiltfiltудалить искажение фазы. В действительности этот процесс удваивает порядка фильтра. В обоих случаях функция минимизирует переходные эффекты в обоих концах сигнала путем соответствия с условиями конечной точки.Наконец,
decimateпередискретизирует данные путем выбора каждогоrth указывают от внутренней части отфильтрованного сигнала. В передискретизируемой последовательности (y),y(end)соответствияx(end)когда БИХ-фильтр используется, иy(1)соответствияx(1)когда КИХ-фильтр используется.
Ссылки
[1] Комитет по Цифровой обработке сигналов IEEE® Acoustics, Речи, и Общества Обработки сигналов, программ редакторов для Цифровой обработки сигналов. Нью-Йорк: Нажатие IEEE, 1979.