各位大神好
我想要在频率域中设计一个有Delay的filter
但是实际delay的秒数并非取样点能够取到的地方
例如以下code，取样率为16k，需要delay的时间为0.0017秒
这样的话我的sample点为16k * 0.0017 = 27.2点 (非整数点0.2无法被取到)
导致画出时域的filter很毛 (正常应该是只有一个点为1其余为零)
想请问大神们若要处理像这样的fractional delay问题该使用什么方法较好呢?
ps.可能不能直接近似sample的点数到整数点，因为之后会设计多个等间的delay channel
若是直接近似的话可能会使channel跟channel间的delay 间隔不一致。
谢谢!
fs = 16000;
H_pts = 1024;
NN=2^nextpow2(H_pts);
df=fs/NN;
freqs=fs/2*linspace(0,1,NN/2+1);%0:df:(NFFT/2)*df;
delay_second = 0.0017;
delay = zeros(length(freqs), 1);
for dd = 1:length(freqs)
delay(dd,1) = exp(-1j*2*pi*freqs(dd)*delay_second); %delay time term
end
delay_symm = [delay;conj(flipud(delay(2:end-1)))];
delay_timed = real(ifft(delay_symm)); % Impulse respons
delay_tt = 0:1/fs:(length(delay_timed)-1)/fs;% time axis
figure;plot(delay_tt, delay_timed); % Plot the impulse response in time domain