Интерполяция — увеличивает частоту дискретизации целочисленным фактором
y = interp(x,r)y = interp(x,r,n,cutoff)[y,b] = interp(x,r,n,cutoff)Создайте синусоидальный сигнал, выбранный на уровне 1 кГц. Интерполируйте его фактором четыре.
t = 0:0.001:1; x = sin(2*pi*30*t) + sin(2*pi*60*t); y = interp(x,4);
Постройте исходные и интерполированные сигналы.
subplot 211 stem(0:30,x(1:31),'filled','markersize',3) grid on xlabel 'Sample number',ylabel Original subplot 212 stem(0:120,y(1:121),'filled','markersize',3) grid on xlabel 'Sample number',ylabel Interpolated
Интерполяция увеличивает исходную частоту дискретизации последовательности к более высокому уровню. Это - противоположность десятикратного уменьшения. interp вставляет нули в исходный сигнал и затем применяет lowpass, интерполирующий фильтр к расширенной последовательности. Функция использует lowpass алгоритм интерполяции 8,1 описанных в [1]:
Расширьте входной вектор до правильной длины путем вставки 0s между исходными значениями данных.
Разработайте специальный симметричный КИХ-фильтр, который позволяет исходным данным проходить неизменный и интерполирует, чтобы минимизировать среднеквадратичную погрешность между интерполированными точками и их идеальными значениями. Фильтр, используемый interp, совпадает с фильтром, возвращенным intfilt.
Примените фильтр к расширенному входному вектору, чтобы произвести вывод.
[1] Комитет по Цифровой обработке сигналов Акустики IEEE, Речи, и Общества Обработки сигналов, программ редакторов для Цифровой обработки сигналов. Нью-Йорк: Нажатие IEEE, 1979.
[2] Oetken, G., Томас В. Парки и Х. В. Шюсслер. “Новые результаты в проекте цифровых интерполяторов”. IEEE® Transactions на Акустике, Речи и Обработке сигналов. Издание ASSP-23, № 3, июнь 1975, стр 301–309.