数字信号课件设计.docVIP

电 子 科 技 大 学 实 验 报 告 学生姓名：郑闻尧 学 号： 学生姓名： 樊聪丽 学 号：2012019090026 学生姓名： 王杰 学 号： 指导教师：朱学勇 分工：Matlab仿真程序的撰写 以及报告的修改 郑闻尧 Matlab程序运行，结果分析以及报告撰写 樊聪丽 王杰 一、实验室名称：数字信号处理实验室 二、实验项目名称：题目三 使用Matlab实现抽取和内插 三、实验原理： 1、采样的概念：采样是将连续信号变化为离散信号的过程。 A、理想采样：即将被采样信号与周期脉冲信号相乘 B、实际采样：将被采样信号与周期门信号相乘，当周期门信号的宽度很小，可近似为周期脉冲串。 根据傅里叶变换性质 式中T代表采样间隔， 由上式可知：采样后信号的频谱是原信号频谱以为周期的搬移叠加 结论：时域离散化，频域周期化；频谱周期化可能造成频谱混迭。 C、低通采样和Nyquist采样定理 设且， 即为带限信号。则当采样频率满足时,可以从采样后的信号无失真地恢复。称为奈奎斯特频率，为奈奎斯特间隔。 2.信号内插和抽取的原理 A、降采样率（整数倍抽取）的实现原理，时域和频域的变化情况。 x[n]y[n]降采样率是指每次抽样保留输入序列中的第M个样本，而除去中间的M-1个样本：用框图表示为 x[n] y[n] 可以得到 ，以2倍下抽样器为例，即L=2，可得 ，如下图所示 混叠 混叠 可以知道，在降采样率时，的原形状会丢失，即发生混叠现象。M倍下抽样器的输出和输入之间傅氏变换的关系为： 在下抽样以前，为了避免引起混叠，信号需要通过一个低通滤波器来带限到即： x[n] x[n] y[n] B、升采样率（整数倍内插）的实现原理，时域和频域的变化情况。 升采样率是指通过在对原离散信号的两个连续样本间插入L-1个等距的样本值（不一定为零），亦即抽样因子为L的上抽样。上抽样后的序列长度为原来的L倍： ，框图表示为 x[n] x[n] xu[n] 可以得到：，，对于L=2时，可得下图： 如图，2倍的抽样率扩展导致频谱的2倍重复，表明傅里叶变换以2倍压缩。因此可得输入频谱的一个额外镜像，这个过程也叫做映射。上采样后不必要的镜像必须用一个称为内插滤波器的低通滤波器H(z)来消除，即： x[n] x[n] xu[n] y[n] 五、实验内容： 给定一个连续的低通滤波器，其冲激响应为，频率为1kHz，现利用15kHz的采样率对其采样，得到一个离散的正弦信号x1[n]（取100点长）。要求用MATLAB完成： 1）利用抽取处理将x1[n]采样率变为5kHz，得到x2[n]。绘出x2[n]的时域和频域图形，与之前的x1[n]时域频域进行对比。 2）利用内插处理将x1[n]采样率变为30kHz，得到x3[n]。绘出x3[n]的时域和频域图形，与之前的x1[n]时域频域进行对比。 六、实验器材（设备、元器件）： 安装MATLAB软件的PC机一台 七、实验数据及结果分析： 实验程序 %DSP课程设计 %%%%%%%%%%%%%%%%%%生成连续冲击响应的时域信号%%%%%%%%%%%%%%%%%%%% Nfft=1024; %DTFT仿真长度 N=1000; %仿真长度 fs=10000; %仿真采样率 ts=1/fs; %仿真采样时间间隔 omega0=1000; t=0:ts:(N-1)*ts; %仿真时间范围 x1a=sin(omega0*t)./(omega0*t); %低通滤波器冲激响应的时域表达式 x1a(1)=1; %xla函数在零点处的值 figure(1) plot(t,x1a); %画出低通滤波器的时域图像 title(被采样连续信号);xlabel(t);ylabel(x1a); %%%%%%%%%%%%%%%%%%%生成采样后离散时间信号%%%%%%%%%%%%%%%%%%%%%%% N1=100; FT=1500; %得到离散信号x1[n]的采样频率 T=1/FT; %采样时间间隔 x1n=ones(1,N1); for k=1:N1 x1n(k)=x1a(round((k-1)*T/ts)+1); end n1=0:N1-1; X1=fftshif