[数学]数字信号处理实验指导书打印.doc

实验要求 在进入实验室前，按要求仔细阅读实验内容和相关资料，并熟悉或编写上机程序。 在实验过程中，不可大声喧哗、走动，出现问题及时向指导老师提问或者向附近同学请教。 保持实验室的整洁卫生，个人物品注意保管； 凡调试成功的程序必须由实验辅导教师认可后方可离开实验室，注意先关电脑，再关电源； 实验报告要求 实验原理和目的； 程序的功能介绍及主要流程图； 完成思考题，记录实验过程中遇到的问题及处理的方法； 对实验提出改进意见和想法； 实验报告以书面形式提交，程序及实验绘图可以打印； 每次实验报告在下一次实验课前交给老师。  目 录 实验要求 1 目录 2 实验一、离散时间信号与系统时域分析 3 实验二、 离散时间信号与系统频域分析 6 实验三、离散傅里叶变换及其快速算法 9 实验四、 脉冲响应不变法和双线性变换法设计IIR滤波器 13 实验五、 IIR滤波器的设计方法 16 实验六、窗函数法设计FIR数字滤波器 23 实验七、频率采样法及优化设计法设计FIR 26 实验八、离散系统的频率响应分析和零、极点分布 30 附录1、 MATLAB应用简介 35 附录2、 MATLAB信号处理工具箱函数简介 41 实验一 离散时间信号与系统时域分析 实验学时：2学时 实验类型：验证 一 实验目的 1 学习MATLAB语言编程和调试技巧 2 学会简单的矩阵输入和图形表示法 3 掌握简单的绘图命令 二、实验原理 本实验主要为了熟悉MATLAB环境，重点掌握简单的矩阵（信号）输入和绘图命令,特别是绘图命令stem()和plot()。 实验内容中涉及到信号的无失真采样、离散卷积运算和差分方程求解三个主要的问题。其基本原理分别如下： 对一个模拟信号进行采样离散化，为了不失真地从采样信号中恢复原始信号，采样时必须满足采样定理，即采样频率必须大于等于模拟信号中最高频率分量的2倍。 一个离散时间系统，输入信号为x(n)，输出信号为y(n)，运算关系用T[﹒]表示，则输入与输出的关系可表示为y(n)＝T[x(n)]。 式中*表示卷积运算。 （2）LTI系统的实现 可物理实现的线性时不变系统是稳定的、因果的。这种系统的单位脉冲响应是 因果的（单边）且绝对可和的，即： 在MATLAB语言中采用conv实现卷积运算，即：Y=conv(x,h),它默认从n=0开始。 常系数差分方程可以描述一个LTI系统，通过它可以获得系统的结构，也可以求 信号的瞬态解。利用MATLAB自带的filter（），可以代替手工迭代运算求解系统的 差分方程，求解的过程类似于对输入信号进行滤波处理。 三、实验内容 1、试画出如下序列的波形 （1） （2） 解：用MATLAB描述波形。 （1）程序 x=[3 1 2 0 -4 2 -3]; % 矩阵输入x n=-3:1:3; %输入自变量n,以间隔为1从-3到3变化 stem(n,x); % stem()函数绘制火柴杆图，注意n,x元素个数必须相等 xlabel(n); %横坐标显示n ylabel(x(n)); %纵坐标显示x(n) grid; %绘制网格 (2) 程序n=0:9; x=0.5.^n; stem(n,x); xlabel(n); ylabel(x(n)); grid; 2、用MATLAB计算序列{-2? 0? 1? –1? 3}和序列{1? 2? 0? -1}的离散卷积 解：用MATLAB描述波形。 程序 x=[-2 0 1 -1 3]; h=[1 2 0 -1]; ??????c=conv(x, h); % conv()函数计算序列卷积 ??????M=length(c)-1; %计算序列c的有效长度 ??????n=0:1:M; stem(n,c); ?xlabel(n); ylabel(幅度); ,进行理想采样。 （1）由MATLAB仿真程序画出x(t)的波形。 （2）为了不失真地从采样信号中恢复原始信号，求奈奎斯特采样频率。 （3）分别以采样周期T=0.2s和0.1s对模拟信号进行采样，求采样输出序列 解：（1）x(t)波形仿真程序如下： t=0:1/256:2; % t的取值在0~2s之间，且每隔1/256s进行采样 x=3*cos(2*pi*t)-2*cos(6*pi*t)+ cos(10*pi*t); %输入x plot(t,x); grid; （2）信号x(t)的最高频率分量fmax=5Hz，由此奈奎斯特采样频率=10Hz =2fmax。 （3）由MATLAB仿真程序画出波形及其采样点， 下面可取T=0.2或0.1,可分别获得所要求的信号。 T=0.2; t=0