Интерполяция - увеличение частоты выборки на целочисленный коэффициент
Создайте синусоидальный сигнал, дискретизированный на частоте 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 вставляет нули в исходный сигнал и затем применяет фильтр интерполяции нижних частот к расширенной последовательности. Функция использует алгоритм 8.1 интерполяции нижних частот, описанный в [1]:
Разверните входной вектор до правильной длины, вставив 0s между исходными значениями данных.
Спроектируйте специальный симметричный КИХ-фильтр, который позволяет исходным данным проходить через неизмененные и интерполяции, чтобы минимизировать среднеквадратическую ошибку между интерполированными точками и их идеальными значениями. Фильтр, используемый interp совпадает с фильтром, возвращенным intfilt.
Примените фильтр к развернутому входному вектору для получения выходных данных.
[1] Комитет по цифровой обработке сигналов Общества по акустической, речевой и сигнальной обработке IEEE, eds. Программы для цифровой обработки сигналов. Нью-Йорк: IEEE Press, 1979.
[2] Эткен, Г., Томас В. Паркс и Х. В. Шюсслер. «Новые результаты в разработке цифровых интерполяторов». Транзакции IEEE ® для обработки акустики, речи и сигналов. т. ASSP-23, № 3, июнь 1975, с. 301-309.