《数字信号处理实验3抽样.docxVIP

信号处理实验实验三：抽样一、实验题目：抽样二、实验原理：本部分说明抽样过程有基本原理：混叠和重建，涉及正弦波和线性调频信号的混叠。然后使用DTFT在频域中展开混叠过程。最后实践几种可从其抽样样本把信号恢复的不同手段。对连续时间信号抽样时，因为频域以抽样频率延拓，故其频谱显现出混叠效应。在MATLAB中，只能仿真这一效应。仿真包括抽样运算，D/A转换（包括一个重建滤波器）。为了仿真模拟信号，必须用非常高的采样率—至少是任何模拟信号所允许的具有的最高频率的5倍。三、实验内容与结果分析：1.抽样引起的混叠由于在MATLAB中不能产生模拟信号，所以需要做实轴t轴的仿真。因此，把仿真时的△t与所研究的抽样周期明确的放开是很重要的对连续时间正弦信号考虑下面的表达式：可以按抽样频率对x(t)抽样来获得离散时间信号以不同组合的和 绘出x[n], =8kHz, =300Hz,在10ms长间隔上采样，相位任意指定。（这里为方便，相位统一取零）使用stem绘出产生的离散信号程序：nn=linspace(1,80,80);a=sin(2*pi*300/8000*nn);stem(nn,a)使用plot绘图程序：nn=linspace(1,80,80);a=sin(2*pi*300/8000*nn);plot(nn,a)把正弦频率从100Hz变至475Hz，每次增加125Hz,用subplot指令把四个图放在一个屏上。程序：nn=linspace(1,80,80);a1=sin(2*pi*100/8000*nn);a2=sin(2*pi*225/8000*nn);a3=sin(2*pi*350/8000*nn);a4=sin(2*pi*475/8000*nn);subplot(411),plot(nn,a1)subplot(412),plot(nn,a2)subplot(413),plot(nn,a3)subplot(414),plot(nn,a4)把正弦频率从7525Hz变至7900Hz，每次增加125Hz,用subplot指令把四个图放在一个屏上。程序：nn=linspace(1,80,80);a1=sin(2*pi*7525/8000*nn);a2=sin(2*pi*7650/8000*nn);a3=sin(2*pi*7775/8000*nn);a4=sin(2*pi*7900/8000*nn);subplot(411),plot(nn,a1)subplot(412),plot(nn,a2)subplot(413),plot(nn,a3)subplot(414),plot(nn,a4)把正弦频率从32,100Hz变至32，475Hz，每次增加125Hz,用subplot指令把四个图放在一个屏上。提前预测显现的频率将会增加还是减少。预测：由于频谱混叠，显现的频率会增加。程序：nn=linspace(1,80,80);a1=sin(2*pi*32100/8000*nn);a2=sin(2*pi*32225/8000*nn);a3=sin(2*pi*32350/8000*nn);a4=sin(2*pi*32475/8000*nn);subplot(411),plot(nn,a1)subplot(412),plot(nn,a2)subplot(413),plot(nn,a3)subplot(414),plot(nn,a4)结果分析：从c、d、e中可以看出，当采样频率过小时，就会发生频谱混叠现象。只有当采样频率大于正弦波频率的2倍时，才能不发生频谱混叠，进一步验证了采样定理。2. 抽样的频域视图：（1）实验原理：为了仿真模拟信号，必须用非常高的抽样率——至少是任何模拟信号所允许具有的最高频率的5倍。因此，在这个问题中将会有两个“抽样率“，一个用于所研究的实际抽样中，另一个用于仿真连续时间信号。第二个问题是如何显示连续时间信号的傅里叶变换。同样，这也只能仿真实现。（2）实验内容：编写M文件用来绘出模拟频谱：程序：function fmagplot(xa,dt)L=length(xa);Nfft=round(2.^round(log2(5*L)));Xa=fft(xa,Nfft);range=0:(Nfft/4);ff=range/Nfft/dt;plot(ff/1000,abs(Xa(1:length(range))))title(CONT-TIME FOURIER TRANSFORM(MAG))xlabel(FREQUENCY(kHz)),gridpause3.a.产生一个模拟频率为的余弦信号作为仿真的模拟信号在长度为T的时间间隔上，以速率产生抽样样本。选择仿真模拟信号长度T以取得大约900到1000个样本。b.用plot指令将绘出的时间信号样本连接起来，