连续信号采样和重构.docVIP

数字信号处理实验（综合） 实验题目：连续信号采样和重构 一、实验目的 通过利用MATLAB实现对信号采样、求频谱、滤波以及时域，域重构熟悉通信系统的整个过程。 二、实验原理 奈奎斯特采样定理，连续信号傅立叶变换（CTFT）、连续信号傅立叶逆变换、sample函数时域重构原理、巴特沃兹低通滤波器的设计、时域卷积定理等。 三、实验内容 （1）绘制原信号及其频谱，采样信号及其频谱 图A 连续信号及其采样信号对应频谱图 图1 为y= 3*cos(3*pi*t)+2*sin(2*pi*t)+cos(5*pi*t)的信号，时间间隔为0.01秒。 因为CTFT公式只适用于求连续信号，但本实验中采用的是MATLAB数值计算方法，所以将上面的积分式变成以下的求和式为：，在程序中采用For循环和sub函数实现求解，最后用abs求出其模值输出。 从原信号时域表达式可以看出，信号角频率为5pi,若要应用奈奎斯特采样定理，则采样角频率必须大于2*5pi,于是我们采用15pi的采样角频率。而，所以对应到时域，采样周期为2/15秒。于是在绘制图3时，我们的时间间隔为2/15秒，于是得到许多离散点。同样，利用公式可求的采样信号的频谱图。从图4可以看出，频谱得到了搬移，又由于满足奈奎斯特采样定理，没有出现混频的现象。 （2）离散信号时域重构 图B 离散信号时域重构过程图 重构原理为生成大量自变量点，在每个采样点处，生成一个以该采样点的幅值为中央最大值、为采样时间间隔的sample函数，最后把所有sample函数自变量点的函数值相加，及得到了原信号在这些点处的值，从而重构出原信号。图5中，每个采样点用蓝色圆圈标记出，蓝色虚线是在这些点处生成的sample曲线，红色虚线为自变量点的函数值累加后的值，及为原信号的包络。图6即重构出的原信号。 （3）设计滤波器，频域重构 图C 滤波及频域重构信号 图7即图4，表示采样后的信号频谱图。图8是针对图7频谱设计的低通滤波器，采用的是巴特沃兹型，由于原信号最高频率为2.5赫兹，所以在此我们设计的滤波器通带边界频率为3赫兹，阻带边界频率为4赫兹，通带波纹为0.5dB，阻带波纹为30 dB,以得到测试，满足给定的指标。根据时域卷积性质，将图7与图8 两图中的数值对应相乘，滤去高频分量，得到所需要的低频分量，如图9.然后根据时域傅立叶反变换公式：，同样将其变成求和式，从而得到频域的重构信号，如图10. 四、实验小结 本次实验为一次综合性实验，涉及的内容较多，难度也较大。清晰的设计思路很重要，但编程的过程中遇到了很多麻烦，最关键的是数组相乘中的长度一致问题，所以一开始就应该站在全局、系统的角度来进行设计。通过本次实验，增进了对信号处理的认识，提高了编程能力。 附录 图A程序代码 t=0:0.01:10; %原信号及频谱 采样信号频谱 y1=3*cos(3*pi*t)+2*sin(2*pi*t)+cos(5*pi*t); subplot(2,2,1) plot(t,y1) ylabel(幅度) title((1) 原信号) grid on t2=0:2/15:10; y2=3*cos(3*pi*t2)+2*sin(2*pi*t2)+cos(5*pi*t2); subplot(2,2,3) plot(t2,y2,r:);hold on plot(t2,y2,. ) xlabel(时间 (秒)),ylabel(幅度) title((3)采样后信号 ) grid on t=0:0.01:30; y1=3*cos(3*pi*t)+2*sin(2*pi*t)+cos(5*pi*t); w=-8*2*pi:15*pi/511: 8*2*pi; for i=1:length(w) a(i)=abs(sum(y1.*exp(-j*w(i)*t)))*0.01; end subplot(2,2,2) plot(w/2/pi,a); ylabel(幅度 ) title((2) 原信号频谱); grid on t2=0:2/15:30; w=-8*2*pi:15*pi/511:8*2*pi; y2=3*cos(3*pi*t2)+2*sin(2*pi*t2)+cos(5*pi*t2); for i=1:length(w) m(i)=abs(sum(y2.*exp(-j*w(i)*t2)))*2/15; end; subplot(2,2,4) plot(w/2/pi,m); axis([-8 8 -0.5 60]); ylabel(幅度),xlabel(频率（赫兹）) title((4) 采样后频谱搬移) grid on