Интерполяция - увеличение частоты дискретизации на целочисленный коэффициент
Создайте синусоидальный сигнал, дискретизированный с частотой дискретизации 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, eds. Программы для цифровой обработки сигналов. Нью-Йорк: IEEE Press, 1979.
[2] Oetken, G., Thomas W. Parks, and H. W. Schüssler. «Новые результаты в проекте цифровых интерполяторов». IEEE® Транзакции по акустике, речи и обработке сигналов. Том ASSP-23, № 3, июнь 1975, стр. 301-309.
