数字信号处理实验报告.pdfVIP

....

数字信号处理实验报告

实验一：混叠现象的时域与频域表现

实验原理：当采样频率Fs不满足采样定理，会在0.5Fs附近引起频谱混叠，造成频谱分析误差。

实验过程：考虑频率分别为3Hz，7Hz，13Hz的三个余弦信号，即：g1(t)=cos(6πt),g2(t)=cos(14πt),

g3(t)=cos(26πt),当采样频率为10Hz时，即采样间隔为0.1秒，则产生的序列分别为：g1[n]=cos(0.6

πn),g2[n]=cos(1.4πn),g3[n]=cos(2.6πn)

对g2[n]，g3[n]稍加变换可得：

g2[n]=cos(1.4πn)＝cos((2π-0.6π)n)=cos(0.6πn)

g3[n]=cos(2.6πn)=cos((2π+0.6π)n)=cos(0.6πn)

利用Matlab进行编程：

n=1:300;

t=(n-1)*1/300;

g1=cos(6*pi*t);

g2=cos(14*pi*t);

g3=cos(26*pi*t);

plot(t,g1,t,g2,t,g3);

k=1:100;

s=k*0.1;

q1=cos(6*pi*s);

q2=cos(14*pi*s);

q3=cos(26*pi*s);

holdon;plot(s(1:10),q1(1:10),bd);

figure

subplot(2,2,1);plot(k/10,abs(fft(q1)))

subplot(2,2,2);plot(k/10,abs(fft(q2)))

subplot(2,2,3);plot(k/10,abs(fft(q3)))

通过Matlab软件的图像如图所示：

....

....

如果将采样频率改为30Hz，则三信号采样后不会发生频率混叠，可运行以下的程序，观察序列的频谱。

程序编程改动如下：

k=1:300;

q=cos(6*pi*k/30);

q1=cos(14*pi*k/30);

q2=cos(26*pi*k/30);

subplot(2,2,1);plot(k/10,abs(fft(q)))

subplot(2,2,2);plot(k/10,abs(fft(q1)))

subplot(2,2,3);plot(k/10,abs(fft(q2)