实验一离散信号的时域分析实用技术.docxVIP

实验一 离散系统时域分析 一、实验目的 1．学习MATLAB语言的编程和调试技巧； 2．差分方程的求解； 3．掌握笔算离散卷积方法和MATLAB语言实现。 二、 实验内容 时域中，离散时间系统对输入信号或延迟信号进行运算处理，生成具有所需特性的输出信号。本实验通过MATLAB仿真一些简单的离散时间信号和系统，并研究其时域特性。涉及到离散时间信号、离散时间系统、系统性质及线性卷积等知识点。 三、实验原理与方法和手段 一个离散时间系统，输入信号为x(n)，输出信号为y(n)，运算关系用T[﹒]表示，则输入与输出的关系可表示为y(n)＝T[x(n)]。 在《信号与系统》和《数字信号处理》课程中，我们知道描述线性移不变离散时间系统的数学模型是常系数差分方程，它与系统的结构流图之间可以互相推导。迭代解法（也称递推解法）是求解差分方程的最简单也最适用的方法，也是实现数字滤波器的一种基本方法。 差分方程通式为： x(n)与y(n)分别为系统的激励和响应。 MATLAB以函数filter(num , den , x)，来计算在给定输入和差分方程系数时求差分方程的数值解。num,den分别为系统方程的系数向量。X是输入序列。 线性时不变系统的输入输出关系可通过单位脉冲响应h(n)表示： 式中*表示卷积运算。 可物理实现的线性时不变系统是稳定的、因果的。这种系统的单位脉冲响应是因果的（单边）且绝对可和的，即：h(n)=0,n0； 。在MATLAB语言中采用conv实现卷积运算即：y=conv(x,h)，它默认从n=0开始。 四、实验组织运行要求 1．学生在进行实验前必须进行充分的预习，熟悉实验内容； 2．学生根据实验要求，读懂并理解相应的程序； 3．教师在学生实验过程中予以必要的辅导，独立完成实验； 4．采用集中授课形式。 五、实验条件 1. 具有WINDOWS 操作系统的笔记本电脑； 2. MATLAB编程软件。 六、实验步骤 在“开始--程序”菜单中，找到MATLAB程序，运行启动；进入MATLAB后 ，首先熟悉界面；在Command Window中输入参考程序，并执行；记录运行结果图形，并与笔算结果对照。 具体步骤如下： 已知某一系统方程为： y[n]+3y[n-1]+2y[n-2]=x[n] 计算并画出脉冲响应h(n),n=(-10,100) 具体程序： n=[-10:100]; num=[1]; den=[1 3 2]; x=impseq(-10,100,0); h=filter(num,den,x); stem(n,h); xlabel(时间序号N); ylabel(脉冲响应h); title(脉冲响应); 若x（n）=u（n），计算并画出系统的零状态响应y（n）。 具体程序： n=[-10:100]; num=[1]; den=[1 3 2]; x=stepseq(-10,100,0); h=filter(num,den,x) stem(n,h); xlabel(时间序号); ylabel(零状态响应); title(阶跃响应); 2．设某LTI的单位脉冲响应和h(n)=, (1)判断此系统是否可实现； (2)当输入为时，用MATLAB计算此LTI的输出y(n)，并画出图形(h(n)的n取值范围自定)。 具体程序： n1=0:4; x=[ones(1,5)];%x=u(n)-u(n-10) N1=length(x); n2=0:20; h=0.5.^n2; N2=length(h); y=conv(x,h); N=N1+N2-1; n=0:N-1; subplot(3,1,1);stem(n1,x);xlabel(时间序号N1);ylabel(系统输入x);title(系统输入); subplot(3,1,2);stem(n2,h);xlabel(时间序号N2);ylabel(系统响应x);title(系统输入); subplot(3,1,3);stem(n,y);xlabel(时间序号N);ylabel(零状态响应h);title(阶跃响应); 3．x(n)=[2,3,1,1,2,-1,0,3]，-2≤n≤5；h(n)=[2,4,1,-2,0,-1]，-3≤n≤2，手工用图解法计算卷积y(n)=x(n)*h(n),然后用MATLAB计算卷积,验证两次计算的结果。 具体程序： function[y,ny]=conv_m(x,nx,h,nh) nyb=nx(1)+nh(1); nye=nx(length(x))+nh(length(h)); ny=[nyb,nye]; y=conv(x,h); x=[2,3,1,1,2,-1,0,3];nx=[-2:5]; h=[2,