Сверхдискретизируйте, примените КИХ-фильтр и субдискретизируйте
yout = upfirdn(xin,h)yout = upfirdn(xin,h,p)yout = upfirdn(xin,h,p,q)Измените частоту дискретизации сигнала рациональным коэффициентом преобразования от уровня DAT 48 кГц к частоте дискретизации CD 44,1 кГц. Используйте функцию rat, чтобы найти числитель L и знаменателем M рационального фактора.
Fdat = 48e3; Fcd = 44.1e3; [L,M] = rat(Fcd/Fdat)
L = 147
M = 160
Сгенерируйте синусоиду на 1,5 кГц, выбранную в в течение 0,25 секунд. Постройте первую миллисекунду сигнала.
t = 0:1/Fdat:0.25-1/Fdat;
x = sin(2*pi*1.5e3*t);
stem(t,x)
xlim([0 0.001])
hold on
Разработайте сглаживание lowpass фильтр с помощью окна Kaiser. Установите ребра полосы фильтра как 90% и 110% частоты среза, . Задайте пульсацию полосы пропускания 5 дБ и затухание полосы задерживания 40 дБ. Установите усиление полосы пропускания на L.
f = (Fdat/2)*min(1/L,1/M); d = designfilt('lowpassfir', ... 'PassbandFrequency',0.9*f,'StopbandFrequency',1.1*f, ... 'PassbandRipple',5,'StopbandAttenuation',40, ... 'DesignMethod','kaiserwin','SampleRate',48e3); h = L*tf(d);
Используйте upfirdn с фильтром h, чтобы передискретизировать синусоиду. Вычислите и компенсируйте задержку, введенную фильтром. Сгенерируйте соответствующий передискретизируемый временной вектор.
y = upfirdn(x,h,L,M); delay = floor(((filtord(d)-1)/2-(L-1))/L); y = y(delay+1:end); t_res = (0:(length(y)-1))/Fcd;
Наложите передискретизируемый сигнал на графике.
stem(t_res,y,'*') legend('Original','Resampled','Location','southeast') hold off
Допустимые комбинации размеров xin и h:
xin является вектором, и h является вектором.
Входные параметры являются одним фильтром и одним сигналом, таким образом, функция применяет операцию свертки к xin с h. yout выходного сигнала является вектором - строкой, если xin является вектором - строкой; в противном случае yout является вектор-столбцом.
xin является матрицей, и h является вектором.
Входные параметры являются одним фильтром и многими сигналами, таким образом, функция применяет операцию свертки к h с каждым столбцом xin. Получившийся yout является матрицей с одинаковым числом столбцов как xin.
xin является вектором, и h является матрицей.
Входные параметры являются несколькими фильтрами и одним сигналом, таким образом, функция применяет операцию свертки к каждому столбцу h с xin. Получившийся yout является матрицей с одинаковым числом столбцов как h.
xin является матрицей, и h является матрицей, обоими с одинаковым числом столбцов.
Входные параметры являются несколькими фильтрами и несколькими сигналами, таким образом, функция применяет операцию свертки к соответствующим столбцам xin и h. Получившийся yout является матрицей с одинаковым числом столбцов как xin и h.
upfirdn использует многофазную структуру интерполяции. Количество умножается - добавляют, что операции в многофазной структуре приблизительно (Lh Lx – pLx)/q, где Lh и Lx являются длинами h (n) и x (n), соответственно.
Более точное количество флопс вычисляется в программе, но точное количество является все еще аппроксимированным. Для длинных сигналов x (n) формула часто точна.
upfirdn выполняет каскад трех операций:
Сверхдискретизируйте входные данные в матричном xin фактором целочисленного p (вставка нулей)
КИХ-фильтр сверхдискретизированные данные сигнала с импульсной последовательностью ответа, данной в векторном или матричном h
Субдискретизируйте результат фактором целочисленного q (выбрасывающий выборки)
КИХ-фильтр обычно является фильтром lowpass, который необходимо разработать использование другой функции, такой как firpm или fir1.
[1] Crochiere, R. E. и Лоуренс Р. Рэбинер. Многоскоростная Цифровая обработка сигналов. Englewood Cliffs, NJ: Prentice Hall, 1983, стр 88–91.
[2] Crochiere, R. E. “Общая Программа, чтобы Выполнить Преобразование Уровня Выборки Данных Рациональными Отношениями”. Программы для Цифровой обработки сигналов (Комитет по цифровой обработке сигналов Акустики IEEE, Речи, и Общества Обработки сигналов, редакторов). Нью-Йорк: Нажатие IEEE, 1979, Программы 8.2-1-8.2-7.