下载文档原格式
第三章DSP芯片系统实验实验3.1 ：数据存取实验一．实验目的1.了解TMS320F2812A的内部存储器空间的分配及指令寻址方式。
2.了解ICETEK-F2812-A评估板扩展存储器空间寻址方法，及其应用。
3.了解ICETEK-F2812-EDU实验箱扩展存储器空间寻址方法，及其应用。
4.学习用Code Composer Studio修改、填充DSP内存单元的方法。
5.学习操作TMS32028xx内存空间的指令。
二．实验设备计算机，ICETEK-F2812-A-EDU实验箱（或ICETEK仿真器+ICETEK-F2812-A评估板+相关连线及电源）。
三．实验内容在外部SARAM的0x80000~0x8000f单元置数0～0xf，将该单元块存储的数据复制到0x80100~0x8010f处，最后通过“Memory”查看窗口观察各存储区中的数据。
四．实验原理TMS32028xx DSP内部存储器资源介绍：TMS32028xx系列DSP基于增强的哈佛结构，可以通过三组并行总线访问多个存储空间。
它们分别是：程序地址总线（PAB）、数据读地址总线（DRAB）和数据写地址总线（DW AB）。
由于总线工作是独立的，所以可以同时访问程序和数据空间。
TMS32028xx系列DSP的地址映象请参考第一章1.2.4节ICETEK-F2812-A评估板的存储空间定义及寄存器映射说明中的介绍。
五．实验步骤1.实验准备连接实验设备。
参见第一章1．3．1节中的“硬件连接方法”。
连接仿真器USB口接线，打开实验箱电源开关，接通评估板电源（关闭实验箱上的扩展模块和信号源电源开关）。
2.设置Code Composer Studio 2.21在硬件仿真(Emulator)方式下运行。
参见第一章1．4．2节中的“设置CCS工作在硬件仿真环境”。
3.启动Code Composer Studio 2.21选择菜单Debug→Reset CPU。
DSP实验报告（第二次实验）实验四、正弦信号发生器学院：信息工程学院班级：11级通信工程3班姓名：李慧学号：2011551309指导老师：姚志强完成日期：2014.4.12验四、正弦信号发生器一、 实验目的1. 掌握利用DSP 产生正弦信号的原理2. 熟悉子程序调用的程序结构以及堆栈的使用3. 掌握CCS 的图形输出操作二、 实验设备1. 集成开发环境Code Composer Studio （简称CCS ）2. 实验代码Sin.s54、Lab.cmd 和Lab.gel三、实验内容1. 阅读理解多项式逼近正弦的文档2. 阅读和理解Sin.s543. 调试正弦波发生器4. 加入断点，并选取图形观测，利用动画及时更新5. 试利用迭代的方法来实现正弦信号发生器四、实验结果和提示1.2345sin()= 3.140625 + 0.02026367 - 5.325196 + 0.5446778 + 1.800293x x x x x x ，x 为第一象限内的弧度值。
因为sin()sin(),sin()sin()x x x x π-=-=-，所以只需将第二，三，四象限内的弧度值转换到第一象限即可计算出相应的正弦函数值。
由于有限精度，规定弧度值从~ππ-，其中π=0x7FFF ，π/2=0x4000，π-=0x8000。
利用级数展开产生正弦波，必须在调用计算子程序之前备份好累加器A 中的当前弧度值，以便计算结束后实现x 增量。
正弦波的频率可以通过增幅的大小来进行控制，如果假定程序循环一次为一个时间单位，则正弦波的周期为65536/步长，频率为周期倒数。
x 自动增长时要注意当x 超过π后必须调整到~ππ-的范围内才能调用计算子程序，即若,2x x x ππ>=-则。
2. 需要使用临时数据时，必须用frame 语句留出所需空间，使用结束后要将堆栈指针还原以防堆栈内存泄漏。
要注意的是frame 的下一条指令不能使用直接寻址。
实验三 IIR 数字滤波器的设计04013222 张嘉俊一、实验目的（1）掌握双线形变换法及脉冲响应不变法设计IIR 数字滤波器的具体设计方法及其原理，熟悉用双线形变换法及脉冲响应不变法设计低通，高通和带通IIR 数字滤波器的计算机编程。
（2）观察双线形变换法及脉冲响应不变法设计的滤波器的频域特性，了解双线形变换法及脉冲响应不变法的特点。
（3）熟悉巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器的频率特性。
（3）实验中有关变量的定义cr sf f At f T 通带边界频率阻带边界频率通带波动最小阻带衰减采样频率采样周期（4）设计一个数字滤波器一般包括以下两步a. 按照任务要求，确定滤波器性能指标b. 用一个因果稳定的离散时不变系统的系统函数去逼近这一性能要求（5）数字滤波器的实现对于IIR 滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。
先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。
用双线形变换法设计IIR 数字滤波器的过程：a. 将设计性能指标中的关键频率点进行“预畸”b. 利用“预畸”得到的频率点设计一个模拟滤波器。
c. 双线形变换，确定系统函数二、实验内容（1）fc=0.3kHz，δ=0.8dB，fr=0.2kHz，At=20dB，T=1ms；设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。
源程序：clc,clear;wc=2000*tan(2*pi*300*0.001/2);wr=2000*tan(2*pi*200*0.001/2);[N,wn]=cheb1ord(wc,wr,0.8,20,'s');[B,A]=cheby1(N,0.8,wn,'high','s');[num,den]=bilinear(B,A,1000);[h,w]=freqz(num,den);f=w/(2*pi)*1000;plot(f,20*log10(abs(h)));axis([0,500,-80,10]);grid;xlabel('Frequency/Hz');ylabel('Amplitude/dB');title('Chebyshev High-pass Filter');实验结果：分析与结论：观察图形可知，δ趋近于0；f=200Hz时，幅度约为-30dB，满足At=20dB的要求，故其通带损耗和阻带衰减满足要求。
DSP技术及课程设计实验报告二(精)东南大学自动化学院实验报告课程名称： D SP 原理及C 程序开发第二次实验实验名称：基于DSP 系统的实验——指示灯、拨码开关和定时器院（系）：自动化专业：自动化姓名：学号：实验室：实验组别：同组人员：实验时间：2012 年 4 月 18日评定成绩：审阅教师：第一部分实验：基于DSP 系统的实验——指示灯和拨码开关一．实验目的1. 了解ICETEK –F28335-A 评估板在TMS320F28335DSP外部扩展存储空间上的扩展。
2. 了解ICETEK –F28335-A 评估板上指示灯和拨码开关扩展原理。
3. 学习在C 语言中使用扩展的控制寄存器的方法。
二．实验设备计算机，ICETEK –F28335-A 实验箱（或ICETEK仿真器+ICETEK–F28335-A 评估板+相关连线及电源）。
三．实验原理1．TMS320F28335DSP 的存储器扩展接口存储器扩展接口是DSP扩展片外资源的主要接口，它提供了一组控制信号和地址、数据线，可以扩展各类存储器和存储器、寄存器映射的外设。
-ICETEK –F28335-A 评估板在扩展接口上除了扩展了片外SRAM外，还扩展了指示灯、DIP 开关和D/A 设备。
具体扩展地址如下：0x180004- 0x180005：D/A 转换控制寄存器0x180001：板上DIP 开关控制寄存器0x180000：板上指示灯控制寄存器-与ICETEK –F28335-A 评估板连接的ICETEK-CTR显示控制模块也使用扩展空间控制主要设备：208000-208004h ：读-键盘扫描值，写-液晶控制寄存器208002-208002h ：液晶辅助控制寄存器208003-208004h ：液晶显示数据寄存器2．指示灯与拨码开关扩展原理图1指示灯扩展原理图2拨码开关扩展原理四．实验步骤LED 程序如下：#define LED (*(unsigned short int *0x180000 for(;;{LED=0x01; Delay(1000; LED=0x02;Delay(1000;LED=0x04;Delay(1000;LED=0x08;Delay(1000;}开关程序如下;#define SW (*(unsigned short int *0x180001 for(;;{LED=SW;}五．实验结果可知：映射在扩展存储器空间地址上的指示灯寄存器在设置时是低4位有效的，数据的最低位对应指示灯D1，次低位对应D2，... 依次类推。
DSP实验报告(二)实验二应用FFT对信号进行频谱分析一、实验目的1、在理论学习的基础上，通过本次实验，加深对快速傅里叶变换的理解，熟悉FFT算法及其程序的编写。
2、熟悉应用FFT对典型信号进行频谱分析的方法。
3、了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。
二、实验原理与方法①一个连续信号的频谱可以用它的傅立叶变换表示为+ Xa(jW)=-jWtx(t)edtòa-如果对该信号进行理想采样，可以得到采样序列x(n)=xa(nT)同样可以对该序列进行z变换，其中T为采样周期X(z)=+ x(n)z-n+ -令z为ejw，则序列的傅立叶变换X(ejw)=x(n)ejwn-其中ω为数字频率，它和模拟域频率的关系为w=WT=W/fs式中的是采样频率。
上式说明数字频率是模拟频率对采样率的归一化。
同模拟域的情况相似。
数字频率代表了序列值变化的速率，而序列的傅立叶变换称为序列的频谱。
序列的傅立叶变换和对应的采样信号频谱具有下式的对应关系。
1X(e)=Tjw+ - w-2pXa(j)T即序列的频谱是采样信号频谱的周期延拓。
从式可以看出，只要分析采样序列的谱，就可以得到相应的连续信号的频谱。
注意：这里的信号必须是带限信号，采样也必须满足Nyquist定理。
在各种信号序列中，有限长序列在数字信号处理中占有很重要的地位。
无限长的序列也往往可以用有限长序列来逼近。
有限长的序列可以使用离散傅立叶变换。
当序列的长度是N时，定义离散傅立叶变换为：X(k)=DFT[x(n)]=其中W=e2pj-NN-1n=0WNkn它的反变换定义为：1x(n)=IDFT[X(k)]=N根据式和，则有N-1n=0X(k)WNknX(z)|z=Wnk=NN-1n=0x(n)WNnk=DFT[x(n)]j2pN可以得到X(k)2pk的点，就NN是将单位圆进行N等分以后第k个点。
所以，X(k)是z变换在单位圆上的等距采样，或者说是序列傅立叶变换的等距采样。
DSP课程设计实验报告语音噪声滤波院（系）：电子信息工程学院设计人员：刘聪学号：08213013叶鸣08213023成绩：工程设计50 报告20 答辩30 总分评语：指导教师签字：日期：目录一、设计任务书 (1)二、设计内容 (2)三、设计方案、算法原理说明 (3)四、程序设计、调试与结果分析 (11)五、设计（安装）与调试的体会 (37)六、参考文献 (39)语音通信的目的是传递声音信息。
位于甲地的通信者发出的声音经语音传感器变换成为电信号，经发送端设备变换为适合传输的形式，通过传输信道传输到乙地。
在乙地经接收端设备恢复出原来的语音信号，经耳机或者喇叭转换为接收者可以听到的声音信号。
这就是最基本的语音通信系统，但是由于周围环境的原因，我们采集到语音信号经常含有不同程度的噪声。
噪声是由于发生体作无规则振动产生的。
典型语音通信系统中的噪声来自三个方面：①信号处理设备产生的电噪声及传输信道中的电噪声；②信号发送端空间环境中的音频噪声信号经麦克风变换为电信号之后，与有用信号共同传递到接收端；③信号接收端空间环境中的音频噪声对信号接收者的影响。
在很多情况下，环境中的背景噪声是通信系统中噪声干扰的主要来源。
当语音信号受到背景噪声干扰时语音通信质量变得不可接受，因此要对语音信号中的噪声滤除。
DSP利用直接存储器访问方式DMA（Direct Memory Access）采集数据时不打扰CPU，因此利用DMA方式工作时，CPU可以对语音信号进行实时地滤波。
本设计要求利用DSP的DMA方式进行信号采集和信号输出，对语音信号进行数字编码，滤波后进行解码。
自适应滤波不仅能够选择信号，而且能够控制信号的特性。
自适应滤波器具有跟踪信号和噪声变化的能力，它的系数能够被一种自适应算法所修改。
利用DSP可以实时地对信号进行自适应滤波。
DSP利用直接存储器访问方式DMA采集数据时不打扰CPU，因此CPU可以对信号进行实时地滤波。
评定成绩:审阅教师:东南大学自动化学院实验报告第三次实验实验名称：基于DSP 系统的实验——液晶屏、键盘外设控制实验姓 名:院（系八自动化 专 业: 自动化同组人员: 实验时间: 2017年3月30日课程名称:DSP 原理及C 程序开发 学 号:《DSP技术及课程设计》实验报書学1；・08014102一-实验目的二-实验S备三-实验原理四.基本实验内容五.提高要求六.实验小结12《DSP技术及课程设计》实验报倂学号08014102一.实验目的通过实验学习使用F28335ADSP的扩展端I I控制外W设备的方法• 了解（1）发光一•极tf的控制编程方法（2）液晶显示器编程方法（3）键盘外设控制编程方法。
二.实验设备计算机，ICETEK・F28335・EDU实验箱(或ICETEK仿真器+ICETEK - F28335-A系统板+相关连线及电源）。
三.实验原理1.发光二极管显示阵列（交通灯）控制TMS32OF28335DSP的扩展存储器接I l（5MnO用来与大多数外W设备进行连接.典型应用如连接片外扩展存储器等。
这一接II提供地址连线、数据连线和一组控制线。
ICETEK-F28335-A将这些扩展线引到了板上的扩展插座上供扩展使用。
实验箱中ICETEK.CTR板上的发光二极管显示阵列（交通灯）的显示是由扩展端I 1控制,扩展在EMIF接口的两个寄存器提供具体控制。
CTRGR为全W控制寄存器，地址为0x208000, 作用为对液晶控制板上的设备进行初始化。
CTRLR为辅助控制寄存器，地址为0x208007,作用为対交通灯进行控制。
交通灯一共12个，使其顺序亮灭的例程（labHC2）如下所示:#defuie LEDS * (iiit *)0xc0000#defiiie CTRGR ♦ (int *)0x208000 //page 29 expeTiinent instructor #defiiie CTRLCDCMDR * (iiit *)0x208001 #defiiie CTRKEY » (iiit *)0x208001 ftdefiiie CTRLCDCR * (int *)0x208002#defiiie CTRCLKEY * (int *)0x208002 ftdefiiie CTRLCDLCR ♦(int *)0x208003 #d€fine CTRLCDRCR * (iiit *)0x208004 #defiiie CTRLA • (iiit *)0x208005 f tdefiiie CTRLR ♦ (int *)0x208007 luisigned nituLed[12][2] = {{0x01,0x40},{0x02,0x40},{0x04,0x40},{0x00,0x48}, {0x00,0x50},{0x00,0x60}, {0x08,0x40},{0x10,0x40 )40x20,0x40}.{0x00,0x41}, {0x00,0x42}, {0x00,0x44},void inaui(void)intnCoiintInitSysCtrlQ;InitXintfl6GpioO：CTRGR = 0x80;// 初始化ICETEK 一CTR《DSP技术及课程设计＞实验报書学1；・08014102CTRGR = 0x0;CTRGR = 0x80;CTRLR = 0;//关闭东酋方向的交通灯CTRLR = 0x40; //关闭南北方向的交通灯11 Count = 0;* (iiit *)0x208007 = uLed[nCoimt][0];〃设置指示灯状态» (iiit *)0x208007 = uLed[nCoimt][l];nCount + +； nCoiinl% =12;Dehy(512);全局控W寄存S CTRGROS：全局控制标志位:BUZZE:蜂鸣器使能:PWME： PWM控制使能:lOPE：通用I/O端「lCPA・LED)n接控制交通灯北方向红灯使能:DCME：直流电机使能;例如需要使能II 流电机•町以用以卞C 语言语句：CTRGR=1：CTRLR 寄存器对应端【I 地址为:0x208007＞需要连续两次写入8位数据。
dsp信号处理实验报告DSP信号处理实验报告一、引言数字信号处理（DSP）是一种将连续信号转换为离散信号，并对其进行处理和分析的技术。
在现代通信、音频处理、图像处理等领域中，DSP技术被广泛应用。
本实验旨在通过对DSP信号处理的实践，加深对该技术的理解与应用。
二、实验目的本实验旨在通过对DSP信号处理的实践，掌握以下内容：1. 学习使用DSP芯片进行信号采集和处理；2. 理解离散信号的采样和重构过程；3. 掌握常见的DSP信号处理算法和方法。
三、实验原理1. 信号采集与重构在DSP信号处理中，首先需要对模拟信号进行采样，将连续信号转换为离散信号。
采样过程中需要注意采样频率的选择，以避免混叠现象的发生。
采样完成后，需要对离散信号进行重构，恢复为连续信号。
2. DSP信号处理算法DSP信号处理涉及到多种算法和方法，如滤波、频谱分析、时域分析等。
其中，滤波是一种常见的信号处理方法，可以通过滤波器对信号进行去噪、增强等处理。
频谱分析可以将信号在频域上进行分析，了解信号的频率成分和能量分布。
时域分析则关注信号的时序特征，如幅值、相位等。
四、实验步骤1. 信号采集与重构在实验中，我们使用DSP芯片进行信号采集与重构。
将模拟信号输入DSP芯片的模拟输入端口，通过ADC（模数转换器）将模拟信号转换为数字信号。
然后，通过DAC（数模转换器）将数字信号转换为模拟信号输出。
2. 滤波处理为了演示滤波处理的效果，我们选择了一个含有噪声的信号进行处理。
首先，使用FIR滤波器对信号进行低通滤波，去除高频噪声。
然后，使用IIR滤波器对信号进行高通滤波，增强低频成分。
3. 频谱分析为了对信号的频率成分和能量分布进行分析，我们使用FFT（快速傅里叶变换）算法对信号进行频谱分析。
通过观察频谱图，可以了解信号的频率特性。
4. 时域分析为了对信号的时序特征进行分析，我们使用时域分析方法对信号进行处理。
通过计算信号的均值、方差、峰值等指标，可以了解信号的幅值、相位等特性。
数字信号处理第二次实验报告--------IIR 数字滤波器的设计姓名：印友进 学号：04012540一、实验目的（1）掌握双线形变换法及脉冲响应不变法设计IIR 数字滤波器的具体设计方法及其原理，熟悉用双线形变换法及脉冲响应不变法设计低通，高通和带通IIR 数字滤波器的计算机编程。
（2）观察双线形变换法及脉冲响应不变法设计的滤波器的频域特性，了解双线形变换法及脉冲响应不变法的特点。
（3）熟悉巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器的频率特性。
二、实验原理（1）脉冲响应不变法用数字滤波器的单位脉冲响应序列()h n 模仿模拟滤波器的冲激响应()a h n ,让()h n 正好等于()a h n 的采样值，即()()a h n h n =，其中T 为采样间隔，如果以()a H n 及()H z 分别表示()a h n 的拉式变换及()h n 的Z 变换，则有)2(1)(m T j s H T z H m a e z sT ∑∞-∞==+=π（2）双线性变换法S 平面与z 平面之间满足以下映射关系：111212,112T sz s z T T z s --+-=⋅=+- s 平面的虚轴单值地映射于z 平面的单位圆上，s 平面的左半平面完全映射到z 平面的单位圆内。
双线性变换不存在混叠问题。
双线性变换是一种非线性变换（tan22TωΩ=），这种非线性引起的幅频特性畸变可通过预畸而得到校正。
IIR 低通、高通、带通数字滤波器设计采用双线性原型变换公式：（3）实验中有关变量的定义c rs f f At f T δ通带边界频率阻带边界频率通带波动最小阻带衰减采样频率采样周期（4）设计一个数字滤波器一般包括以下两步 a. 按照任务要求，确定滤波器性能指标b. 用一个因果稳定的离散时不变系统的系统函数去逼近这一性能要求（5）数字滤波器的实现 对于IIR 滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。
先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。
用双线形变换法设计IIR 数字滤波器的过程： a. 将设计性能指标中的关键频率点进行“预畸”b. 利用“预畸”得到的频率点设计一个模拟滤波器。
c. 双线形变换，确定系统函数三、实验内容（1）fc=0.3kHz ，δ=0.8dB ，fr=0.2kHz ，At=20dB ，T=1ms ；设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。
源程序：clc,clear;% 预畸变，W=2/T*tan(w/2) % w 为数字角频率，w=2*pi*fc/fswc=2000*tan(2*pi*300*0.001/2);wr=2000*tan(2*pi*200*0.001/2);[N,wn]=cheb1ord(wc,wr,0.8,20,'s');[B,A]=cheby1(N,0.8,wn,'high','s');[num,den]=bilinear(B,A,1000);[h,w]=freqz(num,den);f=w/(2*pi)*1000;plot(f,20*log10(abs(h)));axis([0,500,-80,10]);grid;xlabel('Frequency/Hz');ylabel('Amplitude/dB');title('Chebyshev High-pass Filter');实验结果：分析与结论：观察图形可知，δ趋近于0；f=200Hz时，幅度约为-30dB，满足At=20dB的要求，故其通带损耗和阻带衰减满足要求。
（2）fc=0.2kHz，δ=1dB，fr=0.3kHz，At=25dB，T=1ms；分别用脉冲响应不变法及双线形变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。
比较这两种方法的优缺点。
源程序：clc,clear;%脉冲响应不变法对应的频率wc=2*pi*200;wr=2*pi*300;%双线性变换法对应的频率（预畸变）wc2=2*1000*tan((2*pi)*200*0.001/2);wr2=2*1000*tan((2*pi)*300*0.001/2);%脉冲响应不变法[N1,wn]=buttord(wc,wr,1,25,'s');[B1,A1]=butter(N1,wn,'s');[num1,den1]=impinvar(B1,A1,1000);[h1,w1]=freqz(num1,den1);%双线性变换法[N2,wn2]=buttord(wc2,wr2,1,25,'s');[B2,A2]=butter(N2,wn2,'s');[num2,den2]=bilinear(B2,A2,1000);[h2,w2]=freqz(num2,den2);%作图f1=w1/pi*500;f2=w2/pi*500;plot(f1,20*log10(abs(h1)),'-',f2,20*log10(abs(h2)),'-.'); axis([0,500,-80,10]);grid;xlabel('Frequency/Hz');ylabel('Amplitude/dB');title('Butterworth Lowpass Filter');程序结果：分析与结论：从图中可见，通带边界和阻带边界分别为200hz,300hz,衰减量也满足为25Db。
脉冲响应不变法的一个重要特点是频率坐标的变换是线性的，ω＝ΩΤ，ω与Ω是线性关系：在某些场合，要求数字滤波器在时域上能模仿模拟滤波器的功能时，如要实现时域冲激响应的模仿，一般使用脉冲响应不变法。
脉冲响应不变法的最大缺点：有频谱周期延拓效应，因此只能用于带限的频响特性，如衰减特性很好的低通或带通,而高频衰减越大，频响的混淆效应越小，至于高通和带阻滤波器,由于它们在高频部分不衰减，因此将完全混淆在低频响应中，此时可增加一保护滤波器，滤掉高于Ωs/2的频带，再用脉冲响应不变法转换为数字滤波器，这会增加设计的复杂性和滤波器阶数，只有在一定要满足频率线性关系或保持网络瞬态响应时才采用。
双线性变换法的主要优点是S平面与Z平面一一单值对应，s平面的虚轴(整个jΩ)对应于Z平面单位圆的一周，S平面的Ω=0处对应于Z平面的ω=0处，Ω= ∞处对应于Z平面的ω= π处,即数字滤波器的频率响应终止于折叠频率处，所以双线性变换不存在混迭效应。
双线性变换缺点: Ω与ω成非线性关系，导致：a. 数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。
b. 线性相位模拟滤波器经双线性变换后，得到的数字滤波器为非线性相位。
c.要求模拟滤波器的幅频响应必须是分段恒定的，故双线性变换只能用于设计低通、高通、带通、带阻等选频滤波器（3）利用双线性变换法分别设计满足下列指标的巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器，并作图验证设计结果：fc=1.2kHz，δ<=0.5dB，fr=2kHz，At>=40dB，fs=8kHz。
比较这三种滤波器的阶数。
源程序：clc,clear;wr=2*8000*tan(2*pi*2000/(2*8000));wc=2*8000*tan(2*pi*1200/(2*8000));[N1 wn1]=buttord(wc,wr,0.5,40,'s');%巴特沃思[B1 A1]=butter(N1,wn1,'low','s');[num1 den1]=bilinear(B1,A1,8000);[h1 w]=freqz(num1,den1);[N2 wn2]=cheb1ord(wc,wr,0.5,40,'s');%切比雪夫[B2 A2]=cheby1(N2,0.5,wn2,'low','s');[num2 den2]=bilinear(B2,A2,8000);[h2 w]=freqz(num2,den2);[N3 wn3]=ellipord(wc,wr,0.5,40,'s');%椭圆[B3 A3]=ellip(N3,0.5,40,wn3,'low','s');[num3 den3]=bilinear(B3,A3,8000);[h3 w]=freqz(num3,den3);f=w/(2*pi)*8000;plot(f,20*log10(abs(h1)),'-',f,20*log10(abs(h2)),'-.',f,20*log10(abs(h3)),'.');axis([0 4000 -100 10]);grid;xlabel('Frequency/Hz');ylabel('Amplitude/dB');title('Three Lowpass Filter');实验结果：分析与结论：在相同指标下：Butterworth型滤波器需要9阶，Chebyshev型需要5阶，椭圆型需要4阶。
（4）分别用脉冲响应不变法及双线形变换法设计一巴特沃思数字带通滤波器，已知fs=3kHz，其等效的模拟滤波器指标为δ<3kHz，2kHz<f<=3kHz，At>=5dB，f>=6kHz，At>=20dB，f<=1.5kHz。
源程序：clc,clear;%脉冲响应不变法w1=2*pi*2000;w2=2*pi*3000;wr1=2*pi*1500;wr2=2*pi*6000;[N1 wn1]=buttord([w1 w2],[wr1 wr2],3,20,'s');[B A]=butter(N1,wn1,'s');[num1 den1]=impinvar(B,A,30000);[h1 w]=freqz(num1,den1);%双线性变换法w1=2*30000*tan(2*pi*2000/(2*30000));w2=2*30000*tan(2*pi*3000/(2*30000));wr1=2*30000*tan(2*pi*1500/(2*30000));[N2 wn2]=buttord([w1 w2],[wr1 wr2],3,20,'s');[B A]=butter(N2,wn2,'s');[num2 den2]=bilinear(B,A,30000);[h2 w]=freqz(num2,den2);f=w/pi*15000;plot(f,20*log10(abs(h1)),'-.',f,20*log10(abs(h2)),'-');axis([0 15000 -80 10]);grid;xlabel('Frequency/Hz');ylabel('Amplitude/dB');title('Butterworth Bandpass Filter');实验结果：分析与结论：由上图可以看出，用脉冲响应不变法由于滤波器的混叠作用在过度带和阻带都衰减的较双线性变换法慢。
