#简单的数字信号插值、抽取及成型,MMSE误差分析
主题:
对于一个信号,奈奎斯特采样后的序列进行9倍插值,之后2倍抽取,选择不同的升余弦滚降因子插值滤波器,分析4.5倍插值后序列和实际采样后的对应序列误差情况。
此处的插值与抽取倍数均可调整,以实现不同分数倍的信号插值。使用升余弦滚降滤波器成型。
算法流程图如下:
需要注意虽然理论上奈奎斯特采样频率大于2倍频即可,但实际仿真为了追求效果,最好大于10倍。
A = 5;
f0 = 9;
x0 = 0:0.001:6/9;
y0 = A * sin(2 * pi * f0 .* x0);
plot(x0,y0);
hold on
fN = 360;
dt = 1/fN;
T = 0:dt:6/9;
yN = A * sin(2 * pi * f0 .* T);
stem(T,yN);
title('奈奎斯特采样后信号图');
xlabel('时间t');
ylabel('幅度A');
len = length(yN);
y1 = zeros(1,9 * len);
a = 1;
for n1 = 1:9 * len
if mod(n1 - 1,9) == 0
y1(n1) = yN(a);
a = a + 1;
end
end
figure
x1 = 1:9 * len;
stem(x1,y1,'fill','r--');