dsp.MedianFilter
Средний фильтр
Описание
dsp.MedianFilter Система object™ вычисляет движущуюся медиану входного сигнала вдоль каждого канала независимо в зависимости от времени. Объект использует метод раздвижного окна, чтобы вычислить движущуюся медиану. В этом методе окно заданной длины отодвинуто каждый канал, выборка выборкой, и объект вычисляет медиану данных в окне. Для получения дополнительной информации см. Алгоритмы.
Вычислить движущуюся медиану входа:
Создайте
dsp.MedianFilterобъект и набор его свойства.Вызовите объект с аргументами, как будто это была функция.
Чтобы узнать больше, как Системные объекты работают, смотрите то, Что Системные объекты? MATLAB.
Создание
Синтаксис
Описание
medFilt = dsp.MedianFiltermedFilt, использование свойств по умолчанию.
medFilt = dsp.MedianFilter(Len)WindowLength свойство к Len.
medFilt = dsp.MedianFilter(Name,Value) задает WindowLength свойство с помощью Name,Value пара.
Пример:
movMin = dsp.MedianFilter('WindowLength',5);Свойства
Использование
Синтаксис
Описание
Входные параметры
Выходные аргументы
Функции объекта
Чтобы использовать объектную функцию, задайте Системный объект как первый входной параметр. Например, чтобы выпустить системные ресурсы Системного объекта под названием obj, используйте этот синтаксис:
release(obj)
Примеры
Удалите высокочастотный шум Используя средний фильтр
Отфильтруйте высокочастотный шум от шумного синусоидального сигнала с помощью среднего фильтра. Сравните производительность среднего фильтра с фильтром усреднения.
Инициализация
Настройте dsp.MedianFilter объект, medFilt, и dsp.MovingAverage объект, movavgWin. Эти объекты используют метод раздвижного окна с длиной окна 7. Создайте осциллограф времени для просмотра выхода.
Fs = 1000; medFilt = dsp.MedianFilter(7); movavgWin = dsp.MovingAverage(7); scope = dsp.TimeScope('SampleRate',Fs,... 'TimeSpanOverrunAction','Scroll',... 'TimeSpan',1,'ShowGrid',true,... 'YLimits',[-3 3],... 'LayoutDimensions',[3 1],... 'NumInputPorts',3); scope.ActiveDisplay = 1; scope.Title = 'Signal + Noise'; scope.ActiveDisplay = 2; scope.Title = 'Moving Average Output (Window Length = 7)'; scope.ActiveDisplay = 3; scope.Title = 'Median Filter Output (Window Length = 7)'; FrameLength = 256; count = 1; sine = dsp.SineWave('SampleRate',Fs,'Frequency',10,... 'SamplesPerFrame',FrameLength);
Отфильтруйте шумный Sine wave
Сгенерируйте шумный синусоидальный сигнал с частотой 10 Гц. Примените средний фильтр и объект скользящего среднего значения к сигналу. Просмотрите выход на осциллографе времени.
for i = 1:500 hfn = 3 * (rand(FrameLength,1) < 0.02); x = sine() + 1e-2 * randn(FrameLength,1) + hfn; y1 = movavgWin(x); y2 = medFilt(x); scope(x,y1,y2); end
Средний фильтр удаляет высокочастотный шум эффективнее, чем объект скользящего среднего значения.
Удалите высокочастотный шум из данных о гироскопе
В этом примере показано, как удалить высокочастотные выбросы из сигнала потоковой передачи использование dsp.MedianFilter Системный объект?.
Используйте dsp.MatFileReader Системный объект, чтобы считать файл MAT гироскопа. Файл MAT гироскопа содержит 3 столбца данных с каждым столбцом, содержащим 7 140 выборок. Эти три столбца представляют Ось X, Ось Y и данные оси Z из датчика движения гироскопа. Выберите формат кадра 714 выборок так, чтобы каждый столбец данных содержал 10 систем координат. dsp.MedianFilter Системный объект использует длину окна 10. Создайте dsp.TimeScope возразите, чтобы просмотреть отфильтрованный выход.
reader = dsp.MatFileReader('SamplesPerFrame',714,'Filename','LSM9DS1gyroData73.mat',... 'VariableName','data'); medFilt = dsp.MedianFilter(10); scope = dsp.TimeScope('NumInputPorts',1,'SampleRate',119,'YLimits',[-300 300],... 'ChannelNames',{'Input','Filtered Output'},'TimeSpan',60,'ShowLegend',true);
Отфильтруйте данные о гироскопе с помощью dsp.MedianFilter Системный объект. Просмотрите отфильтрованные данные оси Z в осциллографе времени.
for i = 1:10 gyroData = reader(); filteredData = medFilt(gyroData); scope([gyroData(:,3),filteredData(:,3)]); end
Исходные данные содержат несколько выбросов. Увеличьте масштаб данных, чтобы подтвердить, что средний фильтр удаляет все выбросы.
Алгоритмы
Метод раздвижного окна
В методе раздвижного окна выход для каждой входной выборки является медианой текущей выборки и Len - 1 предыдущая выборка. Len является длиной окна. Чтобы вычислить первый Len - 1, выходные параметры, когда окно еще не имеет достаточного количества данных, алгоритм, заполняют окно нулями. Как пример, чтобы вычислить среднее значение, когда вторая входная выборка входит, алгоритм заполняет окно Len - 2 нуля. Вектор данных, x, является затем двумя выборками данных, сопровождаемыми Len - 2 нуля. Этот объект выполняет медианную фильтрацию на входных данных в зависимости от времени.
Рассмотрите пример вычисления движущейся медианы потоковой передачи входные данные с помощью метода раздвижного окна. Алгоритм использует длину окна 4. С каждой входной выборкой, которая входит, окно длины 4 проходит данные.
Ссылки
[1] Боденхэм, декан. “Адаптивное обнаружение фильтрации и изменения для потоковой передачи данных”. PH.D. Тезис. Имперский колледж, Лондон, 2012.
Расширенные возможности
Смотрите также
Объекты
dsp.MovingAverage|dsp.MovingMaximum|dsp.MovingMinimum|dsp.MovingRMS|dsp.MovingStandardDeviation|dsp.MovingVariance