Создание КИХ-фильтра интерполяции
b = intfilt(l,p,alpha)2* p ненулевые выборки, при использовании в последовательности, перемежающейся с l-1 последовательные нули каждый l выборки, принимая исходную полосу alpha умножение частоты Найквиста. Возвращенный фильтр b идентичен используемому interp.
b = intfilt(l,n,'Lagrange')nполиномиальная интерполяция Лагранжа I-го порядка на последовательности, перемеженной с l-1 последовательные нули каждый l выборки.
Создайте цифровой интерполяционный фильтр, чтобы увеличить частоту сигнала на семь, используя метод с ограничением полосы пропускания. Задайте коэффициент «ограниченности полосы пропускания» 0,5 и используйте выборки в интерполяции.
upfac = 7; alpha = 0.5; h1 = intfilt(upfac,2,alpha);
Фильтр работает лучше всего, когда исходный сигнал ограничен alpha умножение частоты Найквиста. Создайте полосно-ограниченный сигнал шума путем генерации 200 Гауссовых случайных чисел и фильтрации последовательности с помощью конечной импульсной характеристики lowpass фильтра 40-го порядка. Сбросьте генератор случайных чисел для воспроизводимых результатов.
lowp = fir1(40,alpha);
rng('default')
x = filter(lowp,1,randn(200,1));Увеличьте частоту дискретизации сигнала, вставив нули между каждой парой выборок x.
xr = upsample(x,upfac);
Используйте filter функция для формирования интерполированного сигнала.
y = filter(h1,1,xr);
Компенсируйте задержку, введенную фильтром. Постройте график исходных и интерполированных сигналов.
delay = mean(grpdelay(h1));
y(1:delay) = [];
stem(1:upfac:upfac*length(x),x)
hold on
plot(y)
xlim([400 700])
intfilt также выполняет полиномиальную интерполяцию Лагранжа.
Полиномиальная интерполяция первого порядка является просто линейной интерполяцией, которая выполняется с помощью треугольного фильтра.
Интерполяция порядка нулей выполняется с помощью фильтра скользящего среднего и напоминает выход отображения sample-and-hold.
Интерполируйте исходный сигнал и наложите результат.
h2 = intfilt(upfac,1,'Lagrange'); y2 = filter(h2,1,xr); y2(1:floor(mean(grpdelay(h2)))) = []; plot(y2) hold off
Бандлимитированный метод использует firls для разработки интерполяционной конечной импульсной характеристики. Полиномиальный метод использует формулу полиномиальной интерполяции Лагранжа на одинаково разнесенных выборках, чтобы создать соответствующий фильтр. Оба типа фильтров являются в основном lowpass и имеют коэффициент усиления l в полосе пропускания.
decimate | downsample | interp | resample | upsample