我正在设计一个分数延迟滤波器，我的 5 阶拉格朗日系数 h(n) 在时域中有 6 个抽头。我已经测试过使用 matlab 将 h(n) 与 x(n) 进行卷积，这是 5000 个采样信号，结果似乎还可以。当我尝试使用 FFT 和 IFFT 方法时，输出完全错误。实际上，我的 FFT 是用频域中的 8192 个数据计算的，这是 5000 个信号样本的最接近 2 的幂。对于 IFFT 部分，我将 8192 频域数据转换回时域中的 5000 长度数据。所以，问题是，为什么这个东西在卷积中起作用，但在 FFT 乘法中不起作用。在频域中将我的 6 个抽头 h(n) 转换为 8192 个抽头会导致这个问题吗？
实际上，我已经尝试过使用重叠保存方法，该方法执行 FFT 并与较小的 x(n) 块相乘，并分别执行 5 次。结果似乎比以前稍微好一点，至少我可以看到波形模式，但仍然略有失真。因此，任何想法哪里出了问题，以及解决方案是什么。谢谢你。
我在频域而不是时域实现循环卷积的原因是，我尝试将拉格朗日滤波器与频域中的其他低通滤波器合并，以便实现更高效。当然，我确实相信在频域中实现滤波将比在时域中的卷积快得多。LP 滤波器在时域中有 120 个抽头。由于内存限制，包括填充在内的原始数据长度将限制为 1024，fft bin 也是如此。
因为我的拉格朗日系数只有 6 个抽头，这与 1024 个抽头有很大不同。我怀疑频域中 6 个抽头到 1024 个 bin 的 fft 会导致错误。这是我关于拉格朗日滤波器的 matlab 代码。这只是一个测试代码，不是实现代码。有点乱，抱歉。如果你能给我更多关于这个问题的建议，我真的很感激。谢谢你。
t=1:5000;
fs=2.5*(10^12);
A=70000;
x=A*sin(2*pi*10.*t.*(10^6).*t./fs);
delay=0.4;
N=5;
n = 0:N;
h = ones(1,N+1);
for k = 0:N
index = find(n ~= k);
h(index) = h(index) * (delay-k)./ (n(index)-k);
end
pad=zeros(1,length(h)-1);
out=[];
H=fft(hh,1024);
H=fft([h zeros(1,1024-length(h))]);
for i=0:1:ceil(length(x)/(1024-length(h)+1))-1
if (i ~= ceil(length(x)/(1024-length(h)+1))-1)
a=x(1,i*(1024-length(h)+1)+1:(i+1)*(1024-length(h)+1));
else
temp=x(1,i*(1024-length(h)+1)+1:length(x));
a=[temp zeros(1,1024-length(h)+1-length(temp))];
end
xx=[pad a];
X=fft(xx,1024);
Y=H.*X;
y=abs(ifft(Y,1024));
out=[out y(1,length(h):length(y))];
pad=y(1,length(a)+1:length(y));
end