Интерполяция — увеличивает частоту дискретизации на целочисленный коэффициент
Создайте синусоидальный сигнал, произведенный на уровне 1 кГц. Интерполируйте его на коэффициент четыре.
t = 0:1/1e3:1; x = sin(2*pi*30*t) + sin(2*pi*60*t); y = interp(x,4);
Постройте исходные и интерполированные сигналы.
subplot(2,1,1) stem(0:30,x(1:31),'filled','MarkerSize',3) grid on xlabel('Sample Number') ylabel('Original') subplot(2,1,2) 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.
