DSP学习心得
在做DSP实验的内插部分时,突然发现了如下现象:
用sinc函数做内插,具体代码如下:
T=0.1;
n=0:10; x=cos(5*pi*n*T); stem(n,x);
dt=0.001; t=ones(11,1)*[0:dt:1]; n=n'*ones(1,1/dt+1);
y=x*sinc((t-n*T)/T); hold on; plot(t/T,y,'r')
输出图像如下:
图 1
倘若改变采样频率,令 T=0.2,则所得图像如下:
图 2
理论上,上述算法并无错误,重构后的图像应很大程度上与离散的图像重合(如图 1
所示),然而在图 2 中,重构后的图像在时域上的长度仅为离散图像的一半,也即上述重
构代码在时间轴上其实并不能完全的重构出原信号,虽然理论上可行。那么,问题究竟出
在哪里呢?
经仔细研究代码,发现为使重构信号与原信号重合,必须让 y 函数的横轴量度归一化到
与 x 函数一致,即将 t 变为 t/T,而此时 t 的范围由(0-1)变为了(0-1/T),倘若 1/T 不
评论1
最新资源