采样定理及MATLAB实现报告,采样定理及matlab0.ppt

采样定理及matlab0

2014-8-12 §1.4 线性时不变系统的描述 n阶前向差分方程 §1.4 线性时不变系统的描述 差分方程的求解 递推法 经典解法 时域解法 Z域分析法 例 2 经典解法 §1.5 连续时间信号的数字处理 对连续(模拟)信号实施数字处理的典型框图 *1.5信号的抽样(难点及重点) 连续信号频谱与抽样信号频谱的关系 时域抽样定理 抗混叠滤波 信号的重建 连续信号的离散处理 1.5.1抽样过程的数学模型 上式中δ(t)是单位冲激信号，在上式中只有当t=nT时，才可能有非零值，因此抽样信号写成下式： 在傅里叶变换中，两信号在时域相乘的傅里叶变换等于两个信号分别的傅里叶变换的卷积，按照(1.5.2)式，推导如下： 设 式中，Ωs=2π/T，称为采样角频率，单位是弧度/秒， 上式表明采样信号的频谱是原模拟信号的频谱沿频率轴，每间隔采样角频率Ωs重复出现一次，或者说采样信号的频谱是原模拟信号的频谱以Ωs为周期，进行周期性延拓而成的。 在图1.5.3中，设xa(t)是带限信号，最高截止频率为Ωc，其频谱Xa(jΩ)如图1.5.3(a)所示。 连续信号频谱X(jw)与理想抽样信号频谱Xs(jw)的关系 点抽样信号频谱X(ejW)与理想抽样信号频谱Xs(jw)的关系 连续信号频谱X(jw)与点抽样信号频谱X (ejW )的关系 带限(band limit)信号 周期序列及其产生的算例 例 设x=[2,3,4,5]，求将它延拓6个周期所得的序列。 解：x=[2,3,4,5]; N=length(x); K=6; ny=0:K*N－1； y=x(mod(ny,N)＋1) 序列的对称性分解算例 例 设x(n)=u(n+2)-u(n-7), 将它分解为偶和奇分量. 解：n = [-2:10]; x = stepseq(-2,-2,10)-stepseq(7,-2,10); [x1,n1] = seqfold(x,n); [xe,nxe] = seqadd(0.5*x,n,0.5*x1,n1); [xo,nxo] = seqadd(0.5*x,n,-0.5*x1,n1); 序列的对称性分解算例 设n1=2:8,x1=sin(n1);n2=-1:3;x2= x2=ones(1,5) (a) 求x1与x2之和。 (b) 求x1与x2之积。 (c) n3=0:4;x3=[6,5,4,3,2]; 求它对n=0点折叠后的序列y3。 (d) 利用(c)的结果，构成一个对称序列y4， (e) 利用(c)的结果，构成一个反对称序列y5 解: (a) 的核心语句如下，(b)相仿 n1=2:8,x1=sin(n1);n2=-1:3; x2= ones(1,5); [y1,ny1]=seqadd(x1,n1,x2,n2); 序列的对称性分解算例 (d)和(e)的核心语句： n3=0:4;x3=[6,5,4,3,2]; % 将n扩展到负区域 ny4=-4:4; %对称序列,把n=0处去一点 y4=[y3(1:4),x3]; %反对称序列,把n=0处置0。 y5=[-y3(1:4),0,x3(2:5)] 序列卷积用MATLAB运算算例 对每一个m值，求y(m)可用下列语句 [h1,nh1]= seqfold(h,nh); ％ h1(nh1)=h(-nh) [h2,nh2]= seqshift(h1,nh1,m); ％h2(nh2)=h1(-nh1-m) [y1,ny1]= seqmult(x,nx,h2,nh2); ％ y1=h2.*x y=sum(y1) 用MATLAB函数conv可计算所有m值的卷积 y=conv(x,h) 不过它没有给出y的位置向量n. 序列卷积用MATLAB运算算例 自编带有位置向量的卷积子程序convwthn.m function [y,ny] = convwthn (x,nx,h,nh) nys = nx(1)+nh(1); nyf = nx(end) + nh(end); % MATLAB允许用end表示最后一个下标 y = conv(x,h); ny = [nys:nyf]; 将例2.3.4中的 x,nx,h,nh代入，执行 [y,ny] = convwthn (x,nx,h,nh) 得到 y= 6 3 5 6 19 -31 30 18 -27 -1 9 2 ny=-5 -4 -3 -2 -1 0 1 2 3 4 5 6 序列的相关 序列互相关的定义： 它提供了两