2016年北航数字信号处理上机实验一实验报告

离散时间信号处理实验报告实验一信号的采样与重构班级学号姓名同组者日期实验一信号的采样与重构1实验介绍连续时间信号采样是获得离散时间信号的一种重要方式，但是时域上的离散化会带来信号在频域上发生相应的变化。在本实验中，我们将分别看到低通信号和带通信号在不同的采样率下得到的离散信号波形与连续信号波形在时域和频域上的对应关系。同时，离散信号的二次采样在实际的应用中可能是必须的，有时甚至是非常重要的。在实验的最后，我们也会看到离散信号的抽取和内插所带来的频谱变化。由于matlab语言无法表达连续信号，实验中我们采用足够密的采样点来模拟连续信号（远大于奈奎斯特采样的要求），即：t=0:Ts:T（Ts=1/fs奈奎斯特采样频率）实验中，为了分析离散信号与连续信号之间的频谱关系，加深对采样定理的理解，了解模拟频谱、数字频谱、以及离散信号被加窗后各自的频谱，从而直观的理解采样频率对频谱的影响和加窗后对频谱的影响。由此可以掌握数字处理方法对模拟信号进行频谱分析的基本原则，即：如何选择合适的信号长度、采样周期以使得对模拟信号的频谱分析的误差达到分析的要求。在该实验中，用到的Matlab函数有:plot(x,y)，其作用是在坐标中以x为横坐标、y为纵坐标的曲线，注意x和y都是长度相同的离散向量；xlabel(‘xxx’)，其作用是对x轴加上坐标轴说明“xxx”；ylabel(‘yyy’)，其作用是对y轴加上坐标轴说明“yyy”；title(‘ttt’)，其作用是对坐标系加上坐标轴说明“ttt”；subplot(m,n,w)，其作用是当需要在同一显示面板中显示多个不同的坐标系时，m、n分别指明每行和每列的坐标系个数，w为当前显示坐标系的流水号(1到m*n之间)。在实验中我们需要画出信号的频谱，对于连续信号频谱的逼近需要你自己编写，原理如下：连续时间非周期信号()xt的傅里叶变换对为：()()jtXjxtedt用DFT方法对该变换逼近的方法如下：1、将)(tx在t轴上等间隔（宽度为T）分段，每一段用一个矩形脉冲代替，脉冲的幅度为其起始点的抽样值)(()(nxnTxtxnTt），然后把所有矩形脉冲的面积相加。该方法实际为平顶处理，利用采样和零阶保持器就可以完成,则有：nnTjTenTxjX)()(2、将序列)()(nTxnx截断成从t=0开始长度为0T的有限长序列，包含有N个采样，即时域加矩形窗，则上式又可以进一步近似为：实验一信号的采样与重构210()()NjnTnXjTxnTe由于时域采样，采样频率为1sfT，则频域产生以sf为周期的周期延拓，如果是限带信号，则有可能不产生混叠，成为连续周期频谱序列，频域周期为sf。3、由于数值计算的限制，在频域上也只能计算离散点（频域抽样）上的数值。我们将频域的一个周期)(sf中也分成N段，即0NFfs。每个频域采样点间的间隔为0F。则上式可以进一步化简为：001002211()()00()()()()[()]|sNjknTnFNNjnkjnkfNxnxnTnnXjkTxnTeTxnTeTxnTeTDFTxn由此我们看到了DFT与CTFT之间的近似关系。如果将T和0F取得尽量小，则我们可以得到模拟信号的较精确的时频特性。实验解答一、低通采样考虑以下模拟信号xa(t)，xa(t)=cos(2*pi*f0*t)+2*cos(2*pi*10*f0*t),令f0=100Hz。显然该信号的最高频率1kHz。1．采样频率fs大于信号最高频率10f0的2倍是信号的恢复。1）对信号以fs1=2.2kHz的频率进行采样，得到抽样信号x(n)=xa(t)|t=nT；2）利用内插公式xr(t)=x(n)*[sin(pi*(t-n*Ts)/Ts)/(pi*(t-n*Ts)/Ts)]恢复出原始信号x^a(t)；3）绘出Δx(t)=xa(t)-x^a(t)的时域波形；4）利用前面介绍的画频谱的方法画出xa(t)和x(n)的频谱，对应3）中的误差曲线和频谱图说明是否有新的频率分量产生，是何原因？2．采样频率fs低于信号最高频率10f0的2倍是信号的恢复。分别取fs=6*f0和fs=14*f0，重复1中的操作，注意观察时域波形和频谱的对应变化。实验一信号的采样与重构3（一）、MATLAB源代码clear;f0=100;Fs=3000;Ts=1/Fs;ts=0:Ts:0.1;X=xa(ts);t=0:Ts/10:0.1;L=length(t);Xa=xa(t);Xr=zeros(1,L);forn=1:LXr(n)=xr(X,Ts,t(n));endsubplot(6,1,1);plot(t,Xa);title('原始信号');subplot(6,1,2);stem(ts,X);title('采样信号');subplot(6,1,3);plot(t,Xr);title('恢复信号');subplot(6,1,4);plot(t,Xa-Xr);title('原信号与恢复信号之差');fs=20000;y1=fft(Xa,fs);mag=abs(y1);f=(0:length(y1)-1)'*(10*Fs)/length(y1);subplot(6,1,5)plot(f(1:fs/2),mag(1:fs/2));axis([0120003000]);fs=3000;y2=fft(X,fs);title('原信号频谱');mag=abs(y2);f=(0:length(y2)-1)'*(Fs)/length(y2);subplot(6,1,6)plot(f(1:fs/2),mag(1:fs/2));title('抽样信号频谱');functionf=xa(t)实验一信号的采样与重构4f0=100;f=cos(2*pi*f0*t)+2*cos(2*pi*10*f0*t);endfunctionx=xr(X,Ts,t)x=0;forn=1:length(X)x=x+X(n)*sinc((t-(n-1)*Ts)/Ts);end（二）、运行结果1、fs=2200Hz2、fs=600Hz实验一信号的采样与重构53、fs=1400Hz（三）、分析及结论1（4）：当fs=2200Hz时，基本可认为没有新的频率分量产生，当fs=600Hz或fs=1400Hz时，有新的频率分量产生。可能原因之一是当不满足fs2f0时，采样信号不能完全恢复原始信号，因此在频域上会产生新的频率分量。2、当不满足fs2f0时，可以看到时域恢复信号出现失真，且fs越小，恢复信号与原信号差距越大。在频域上表现为产生了新的频率分量（200/400Hz）。综上所述，欲使采样信号能无失真地恢复原信号，要求是fs=2f0，若不满足该条件，会出现恢复波形失真的情况。二、带通采样通过该题目进一步加深对带通信号采样定理的认识，即不需要象低通采样那样需要信号最高频率的2倍才能恢复出原始信号。带通信号的采样定理可描述如下：如果模拟信号f(t)为带通信号，其角频率限制在fL和fm之间，则必须的最低采样频率fs=2（fm－fL）。且当采样率满足fs＝4f0/(2n+1)（f0＝（fm＋fL）/2为信号的中心频率，n=0,1,2,……），则可以无失真的从采样信号中恢复原始信号。注意当fs的取值大于2fm时，与低通采样定理意义相同。考虑模拟信号xa(t)，xa(t)=cos(2*pi*f0*t)+2*cos(2*pi*f1*t),令f0=1.6kHz，f1=2kHz。实验一信号的采样与重构6该信号带宽0.4kf0=2kHz，显然为带通信号。1．请你根据带通采样定理，选取适当的采样速率以得到抽样信号x(n)，然后同样用内插公式恢复成模拟信号，并思考如何由该模拟信号得到原来的带通信号。2．请绘出带通信号的频谱及x(n)的频谱，观察并比较它们。（一）、MATLAB源代码1、fc=400;%信号带宽400HZfmin=2*fc;fc0=(1600+2000)/2;fc1=1600;%信号一的频率fc2=2000;%信号二的频率fs1=400;%欠采样的采样频率fs2=800;%临界采样采样频率n=0;fs3=4*fc0/(2*n+1);%n取1时的采样速率%fs3=5000;%n取1时的采样速率f0=80000;%用以模拟连续信号的离散信号采样速率%（大nyquist率）t=0:1/f0:1;N=1*f0;xt0=cos(2*pi*fc1*t)+2*cos(2*pi*fc2*t);figure(1);subplot(3,1,1);plot(t,xt0);axis([00.03-44]);title('待采样信号波形');xlabel('x');ylabel('xt');yjw0=fft(xt0,N);absy0=abs(yjw0);f=(0:N-1)*f0/N;实验一信号的采样与重构7%subplot(2,1,2);%plot(f,absy0);%axis([08000-410000]);%title('正常采样信号频谱')%xlabel('数字频率');ylabel('幅度谱');t=0:1/fs3:1;N=1*fs3;xts=cos(2*pi*fc1*t)+2*cos(2*pi*fc2*t);figure(1);subplot(3,1,2);plot(t,xts);axis([00.03-44]);title('采样信号波形');xlabel('x');ylabel('xt');yjw0=fft(xts,N);absy0=abs(yjw0);f=(0:N-1)*f0/N;fs=fs3;%采样频率Ts=1/fs3;%采样周期T0=0.01;n=0:(3*T0)/Ts;t1=0:Ts:3*T0;%disp('t1='+t1);%x1=sin(2*pi*f0*n/fs)+1/3*sin(6*pi*f0*n/fs);%采样信号x1=cos(2*pi*fc1*n/fs)+2*cos(2*pi*fc2*n/fs);T_N=ones(length(n),1)*t1-n'*Ts*ones(1,length(t1));%表示t-nTxa=x1*sinc(fs*T_N);subplot(3,1,3);plot(t1,xa);axis([00.03-44]);title('恢复信号波形');实验一信号的采样与重构82、fc=400;%信号带宽400HZfmin=2*fc;fc0=(1600+2000)/2;fc1=1600;%信号一的频率fc2=2000;%信号二的频率fs1=400;%欠采样的采样频率fs2=800;%临界采样采样频率n=0;fs3=4*fc0/(2*n+1);%n取1时的采样速率f0=80000;%用以模拟连续信号的离散信号采样速率t=0:1/f0:1;N=1*f0;xt0=cos(2*pi*fc1*t)+2*cos(2*pi*fc2*t);figure(1);yjw0=fft(xt0,N);absy0=abs(yjw0);f=(0:N-1)*f0/N;subplot(2,1,1);plot(f,absy0);%axis([08000-410000]);title('待采样信号频谱')xlabel('数字频率');ylabel('幅度谱');%（大nyquist率）t=0:1/fs3:1;N=1*fs3;xt0=cos(2*pi*fc1*t)+2*cos(2*pi*fc2*t);figure(1);yjw0=fft(xt0,N);absy0=abs(yjw0);f=(0:N-1)*f0/N;subplot(2,1,2);plot(f,absy0);axis([080000-410000]);title('正常采样信号频谱')实验一信号的采样与重构9xlabel('数字频率');ylabel('幅度谱');（二）、运行结果1、2、实验一信号的采样与重构10（三）、分析及结论1、对采样信号及频谱