两份2016心电信号的采集与处理

1、中北大学信息商务学院 课程设计任务书 2015-2016 学年第 二 学期 学 院： 中北大学信息商务学院 专 业： 电子信息工程 学生姓名 ： 苏慧敏 学 号： 1305034211 学生姓名 ： 王晓腾 学 号： 1305034217 学生姓名 ： 李 康 学 号： 1305034243 课程设计题目： 心电信号的采集与处理 起 迄 日 期 ： 2016年6 月13日2016年7月1 日 课程设计地点： 系专业实验室 指 导 教 师 ： 王浩全 系 主 任： 王浩全 下达任务书日期: 2016年6 月 9日课 程 设 计 任 务 书1设计目的：掌握信号的采集、存储和处理方法2设计内容和要求

2、（包括原始数据、技术参数、条件、设计要求等）：1、掌握PCI总线的基本结构，了解基于PCI总线A/D卡的通用结构。写出关于基于PCI总线发展趋势的报告。2、通过A/D卡，编写检测信号的采集、存储程序。3、对心电信号进行滤波、分析处理。设计要求：1双线性变换法设计一个butterworth滤波器。设计指标参数为：通带内频率低于 0.2时，最大衰减小于dB;阻带内 0.3,频率区间上，最小衰减大于15dB； 2 为采样间隔，打印出数字滤波器在频率区间 上的频率响应曲线。3用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，滤除其中的高频干扰。观察总结滤波作用与效果4 已知心电图数字信号序列如

3、下：x=-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,.-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,.-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0;S3 设计工作任务及工作量的要求包括课程设计计算说明书(论文)、图纸、实物样品等：程序演示、毕业设计说明书一份课 程 设 计 任 务 书4主要参考文献：l 要求按国标GB 771487文后参考文献著录规则书写，例：1 傅承义，陈运泰，祁贵中.地球物理学基础.北京：科学出版社，1

4、985 （5篇以上）5设计成果形式及要求：程序6工作计划及进度：2016年6 月13 日 2016年6 月20 日：查资料，写PCI 总线的发展现状2016年6 月21 日 2016年6 月29 日：对检测信号进行分析处理2016年6月 30 日 2016年7 月 1 日：完成课程设计说明书系主任审查意见： 签字： 年 月 日设计说明书应包括以下主要内容： （1）封面：课程设计题目、班级、姓名、指导教师、时间 （2）设计任务书 （3）目录 （4）设计方案简介 （5）设计条件及主要参数表 （6）设计主要参数计算 （7）设计结果 （8）设计评述,设计者对本设计的评述及通过设计的收获体会 （9）参考

5、文献中北大学信息商务学院2016信息处理综合实践目 录一、 基于PCI总线A/D卡的报告11 基于PCI总线的基本结构11.1 PCI总线11.2 PCI总线的基本含义12 基于PCI的A/D卡的通用结构23 基于PCI总线发展趋势24 PCI总线的特点35 PCI-8310数据采集卡硬件说明3二、 心电信号的采集与处理101 设计方案简介102 巴特沃斯（butterworth）模拟滤波器102.1 特性102.2 传递函数103 用双线性变换法设计低通滤波器步骤114 设计主要参数表125 程序内容及结果13三、 设计评述19四、 参考文献201、 基于PCI总线A/D卡的报告1 基于PC

6、I总线的基本结构1.1 PCI总线PCI是由Intel公司1991年推出的一种局部总线。从结构上看，PCI是在CPU和原来的系统总线之间插入的一级总线，具体由一个桥接电路实现对这一层的管理，并实现上下之间的接口以协调数据的传送。管理器提供了信号缓冲，使之能支持10种外设，并能在高时钟频率下保持高性能，它为显卡、声卡、网卡、MODEM等设备提供了连接接口，它的工作频率为33MHz/66MHz。PCI是Peripheral Component Interconnect(外设部件互连标准）的缩写，它是目前个人电脑中使用最为广泛的接口，几乎所有的主板产品上都带有这种插槽。PCI插槽也是主板带有最多数量

7、的插槽类型，在目前流行的台式机主板上，ATX结构的主板一般带有56个PCI插槽，而小一点的MATX主板也都带有23个PCI插槽，可见其应用的广泛性。PCI总线是一种不依附于某个具体处理器的局部总线。管理器提供了信号缓冲，使之能支持10种外设，并能在高时钟频率下保持高性能。PCI总线也支持总线主控技术，允许智能设备在需要时取得总线控制权，以加速数据传送。图1-1 典型的PCI系统总线构成 1.2 PCI总线的基本含义不同于ISA总线，PCI总线的地址总线与数据总线是分时复用的。这样做的好处是，一方面可以节省接插件的管脚数，另一方面便于实现突发数据传输。在做数据传输时，由一个PCI设备做发起者（主

8、控，Initiator或Master），而另一个PCI设备做目标（从设备，Target或Slave）。总线上的所有时序的产生与控制，都由Master来发起。PCI总线在同一时刻只能供一对设备完成传输，这就要求有一个仲裁机构（Arbiter），来决定在谁有权力拿到总线的主控权。当PCI总线进行操作时，发起者（Master）先置REQ#，当得到仲裁器（Arbiter）的许可时（GNT#），会将FRAME#置低，并在AD总线上放置Slave地址，同时C/BE#放置命令信号，说明接下来的传输类型。所有PCI总线上设备都需对此地址译码，被选中的设备要置DEVSEL#以声明自己被选中。然后当IRDY#与T

9、RDY#都置低时，可以传输数据。当Master数据传输结束前，将FRAME#置高以标明只剩最后一组数据要传输，并在传完数据后放开IRDY#以释放总线控制权。2 基于PCI的A/D卡的通用结构信号连接器有通道信号，PCI总线接口，外触发，多卡扩展同步接口等4种，以下分别介绍。通道信号接口：位于采集卡的尾部，信号用同轴电缆接入。一张采集卡最多只有4个通道，也可能少于4个通道，视用户的需求而定。PCI总线接口：采集卡与上位PC机的数据和控制信号通信接口。其采用32位PCI总线。外触发接口：用于接外触发信号用。多卡扩展同步接口：用于多卡扩展。其传输的信号包括主卡的同步时钟源、触发以及状态信号。一般用1

10、6针带缆连接。3 基于PCI总线发展趋势 从1992年创立规范到如今，PCI总线已成为了计算机的一种标准总线。PCI总线取代了早先的ISA总线。当然与在PCI总线后面出现专门用于显卡的AGP总线，与现在PCI Express总线，但是PCI能从1992用到现在，说明他有许多优点，比如即插即用(Plug and Play)、中断共享等。在这里我们对PCI总线做一个深入的介绍。从数据宽度上看，PCI总线有32bit、64bit之分；从总线速度上分，有33MHz、66MHz两种。目前流行的是32bit 33MHz，而64bit系统正在普及中。改良的PCI系统，PCI-X，最高可以达到64bit 13

11、3MHz，这样就可以得到超过1GB/s的数据传输速率。如果没有特殊说明，以下的讨论以32bit 33MHz为例。1991 年下半年，Intel 公司首先提出了PCI 的概念，并联合IBM、Compaq、AST、HP、DEC 等100 多家公司成立了PCI 集团，其英文全称为：Peripheral Component Interconnect Special Interest Group(外围部件互连专业组)，简称PCISIG。PCI 有32 位和64 位两种，32位PCI 有124 引脚，64 位有188 引脚，目前常用的是32 位PCI。32 位PCI 的数据传输率为133MBs，大大高于I

12、SA。PCI 总线的主要性能（1） 支持10 台外设（2） 总线时钟频率33.3MHz/66MHz（3） 最大数据传输速率133MB/s（4） 时钟同步方式（5） 与CPU 及时钟频率无关（6） 总线宽度 32 位（5V）/64 位（3.3V）（7） 能自动识别外设PCI （Peripheral Component Interconnect）总线是一种高性能局部总线，是为了满足外设间以及外设与主机间高速数据传输而提出来的。在数字图形、图像和语音处理，以及高速实时数据采集与处理等对数据传输率要求较高的应用中，采用PCI总线来进行数据传输，可以解决原有的标准总线数据传输率低带来的瓶颈问题。4 PC

13、I总线的特点 1.数据总线32位，可扩充到64位。2.可进行突发（burst)式传输。3.总线操作与处理器-存储器子系统操作并行。4.总线时钟频率33MHZ或66MHZ，最高传输率可达528MB/S。5.中央集中式总线仲裁6.全自动配置、资源分配、PCI卡内有设备信息寄存器组为系统提供卡的信息，可实现即插即用（PNP）。7.PCI总线规范独立于微处理器，通用性好。8.PCI设备可以完全作为主控设备控制总线。9.PCI总线引线：高密度接插件，分基本插座（32位）及扩充插座（64位）。5 PCI-8310数据采集卡硬件说明5.1 概述 PCI-8310 模入接口卡适用于提供了PCI 总线插槽的PC

14、系列微机，具有即插即用（PnP）的功能。其操作系统可选用目前流行的 Windows 系列、高稳定性的Unix等多种操作系统以及专业数据采集分析系统LabVIEW 等软件环境。在硬件的安装上也非常简单，使用时只需将接口卡插入机内任何一个PCI总线插槽中并用螺丝固定，信号电缆从机箱外部直接接入。 PCI-8310 模入接口卡允许采用32路单端输入方式或16路双端输入方式。用户可根据需要选择测量单极性信号或双极性信号。其输入的模拟信号由卡前端的37芯D型插头直接接入。本卡还提供了TTL电平的16路输入和16路输出信号通道，这些信号通道由卡后端的40芯扁平电缆转换为37芯D型插头提供给用户。5.2 主

15、要技术参数5.2.1 模入部分(标*为出厂标准状态，下同)（1）输入通道数：单端32路*；双端16路（2）输入信号范围： 0V10V*；5V5V； -10V10V（3）输入阻抗：10M（4）AD转换分辨率：12位（5）AD转换速率：10S（6）AD启动方式：程序启动（7）AD转换非线性误差：1LSB（8）AD转换输出码制：单极性原码*双极性偏移码（9）系统综合误差：0.1 F.S5.2.2 开关量部分（1）输入路数：16路TTL电平（2） 输出路数：16路TTL电平（3） 电源功耗：5V(10) 500mA（4）环境要求：工作温度： 1040相对湿度： 4080存贮温度：-5585 （5） 外

16、型尺寸(不含档板)：长高164.8mm106.7mm ( 6.5英寸4.2英寸)5.3 工作原理PCI-8310 模入接口卡主要由模拟多路开关选通电路、差分放大器电路、模数转换电路、开关量输入输出电路和接口控制逻辑电路组成。5.3.1 模拟多路开关选通电路 模拟通道开关由4片MPC508（或同类产品）及跨接器KJ1、KJ2组成，可以从32路单端信号或16路双端信号中任选一路，送入差分放大器。5.3.2 差分放大器电路 差分放大器由3个运算放大器以及相关的电阻、电容组成一个标准的仪用差分放大器，用以对通道开关选中的模拟信号进行变换处理。5.3.3 模数转换电路12位逐次逼近式AD转换器ADS77

17、4片内自带精密基准源，并经激光修调，具有较高的转换速率和转换精度，其转换时间仅为 10S。AD转换器由程序启动，其转换状态的结束可由程序查询读出或产生结束中断申请。AD转换器的模拟输入信号幅度由跨接器KJ3 选择， AD 转换后的输出代码形式由跨接器KJ4选择，可分别输出二进制原码或双极性偏移二进制码。电位器W1 用于零点调节，W2 用于双极性偏移调节，W3 用于满量程增益调节。5.3.4 开关量输入输出电路本卡还提供了各16路的开关量输入输出信号通道。使用中需注意对这些信号的要求应严格符合TTL 电平规范。5.3.5 接口控制逻辑电路接口控制逻辑电路用来将PCI总线控制逻辑转换成与各种操作相

18、关的控制信号。 5.4 安装本卡的安装十分简便，只要将主机机壳打开，在关电情况下，将本卡插入主机的任何一个空余PCI扩展槽中，再将档板固定螺丝压紧即可。本卡采用的模拟开关是COMS 电路，容易因静电击穿或过流造成损坏，所以在安装或用手触摸本卡时，应事先将人体所带静电荷对地放掉，同时应避免直接用手接触器件管脚，以免损坏器件。禁止带电插拔本接口卡。本卡跨接选择器较多，使用中应严格按照说明书进行设置操作。设置接口卡开关、跨接套和安装接口带缆时均应在关电状态下进行。当模入通道不全部使用时，应将不使用的通道就近对地短接，不要使其悬空，以避免造成通道间串扰和损坏通道。为保证安全及采集精度，应确保系统地线（

19、计算机及外接仪器机壳）接地良好。特别是使用双端输入方式时，为防止外界较大的共模干扰，应注意对信号线进行屏蔽处理。当本卡使用的信号环境较为恶劣时，为保护本卡和主机，用户可以在本卡前端的预留位置加装双向TVS舜态电压保护管。但加装TVS管后，TVS管的特性将使本卡的输入阻抗下降，同时对信号源的驱动能力有一定的要求，否则将降低本卡的采样精度。5.5 使用注意事项5.5.1 模入部分本卡前端37芯D型插座（CZ1）的信号定义见表1 ，用户可根据需要选择连接信号线(单端)或信号线组(双端)。为减少信号杂波串扰和保护通道开关，凡不使用的信号端应就近与模拟地短接。 表1 CZ1模拟输入信号端口定义(括号内表

20、示双端方式)插座引脚号信 号 定 义插座引脚号信 号 定 义1模拟地20模拟地2CH1(CH1+)21CH17(CH1-)3CH2(CH2+)22CH18(CH2-)4CH3(CH3+)23CH19(CH3-)5CH4(CH4+)24CH20(CH4-)6CH5(CH5+)25CH21(CH5-)7CH6(CH6+)26CH22(CH6-)8CH7(CH7+)27CH23(CH7-)9CH8(CH8+)28CH24(CH8-)10CH9(CH9+)29CH25(CH9-)11CH10(CH10+)30CH26(CH10-)12CH11(CH11+)31CH27(CH11-)13CH12(CH1

21、2+)32CH28(CH12-)14CH13(CH13+)33CH29(CH13-)15CH14(CH14+)34CH30(CH14-)16CH15(CH15+)35CH31(CH15-)17CH16(CH16+)36CH32(CH16-)18NC（空脚）37模拟地19模拟地5.5.2 开关量部分本卡后端40芯扁平线插座(CZ2)的信号定义见表2 。 表2 CZ2开关量输入输出信号端口定义插座引脚号信 号 定 义插座引脚号信 号 定 义1+5V电源输出2+5V电源输出3DI14DI25DI36DI47DI58DI69DI710DI811DI912DI1013DI1114DI1215DI1316

22、DI1417DI1518DI1619DO120DO221DO322DO423DO524DO625DO726DO827DO928DO1029DO1130DO1231DO1332DO1433DO1534DO1635数字地36数字地37数字地38数字地39数字地40数字地 40芯扁平电缆转换为37芯D型插头后的信号定义见表3 。 表3 转换为37芯D型插头时开关量输入输出信号端口定义插座引脚号信 号 定 义插座引脚号信 号 定 义1+5V电源输出20+5V电源输出2DI121DI23DI322DI44DI523DI65DI724DI86DI925DI107DI1126DI128DI1327DI149

23、DI1528DI1610DO129DO211DO330DO412DO531DO613DO732DO814DO933DO1015DO1134DO1216DO1335DO1417DO1536DO1618数字地37数字地19数字地5.5.3 控制端口与数据格式(1)各控制端口的地址与功能见表4： 表4 端口地址与功能表端口地址操作命令(字操作)功 能基地址+0写写通道代码基地址+2写启动AD转换，清除中断申请基地址+2读读A/D转换标志和结果基地址+4读读DI开关量数据基地址+6写写DO开关量数据注：所有操作均为16位。(2)模入通道代码数据格式见表5(端口地址为基地址+0)： 表5 模入通道代码数

24、据格式通道号十进制代码十六进制代码输入方式通道号十进制代码十六进制代码输入方式1000H单双171610H单2101H单双181711H单3202H单双191812H单4303H单双201913H单5404H单双212014H单6505H单双222115H单7606H单双232216H单8707H单双242317H单9808H单双252418H单10909H单双262519H单11100AH单双27261AH单12110BH单双28271BH单1312OCH单双29281CH单1413ODH单双30291DH单15140EH单双31301EH单16150FH单双32311FH单(3)读AD转换

25、标志和结果，数据格式见表6(端口地址为基地址+2)： 表6 读AD转换标志和结果数据格式(x表示任意)命令D15D14D13D12D11D10D9D8D7D6D5D4D3D2D1D0AD转换状态读1xxxxxxxxxxxxxxx正在转换读0000DB11DB10DB9DB8DB7DB6DB5DB4DB3DB2DB1DB0转换结束注：双极性时DB11为符号位(4)开关量输入输出信号的数据格式： 开关量输入输出信号的数据格式采用的是位方式，即一个字中的任意一位对应一路输入输出信号（端口地址为基地址+4；+6）。 表7 开关量输入输出信号数据格式端口地址操作命令D15D14D13D12D11D10D

26、9D8D7D6D5D4D3D2D1D0基地址+4读DI16DI15DI14DI13DI12DI11DI10DI9DI8DI7DI6DI5DI4DI3DI2DI1基地址+6写DO16DO15DO14DO13DO12DO11DO10DO9DO8DO7DO6DO5DO4DO3DO2DO12、 心电信号的采集与处理1 设计方案简介双性变换法可以实现从s平面到z平面的单值映射关系，可以使数字滤波器的频率响应模仿模拟滤波器的频率响应。1.先设计模拟滤波器，再转化数字滤波器 2.将模拟指标转变成数字指标3.选择滤波器的最小阶数 4.创建butterworth模拟滤波器5.用双线性变换法实现模拟滤波器到数字滤

27、波器的转换 6.绘制频率响应曲线2 巴特沃斯（butterworth）模拟滤波器 巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。这种滤波器最先由英国工程师斯替芬巴特沃斯（Stephen Butterworth）在1930年发表在英国无线电工程期刊的一篇论文中提出的。2.1 特性 巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。 一阶巴特沃斯滤波器的衰减率为每倍频6分贝，每十倍频20分贝。二阶巴特沃斯滤波器的衰减

28、率为每倍频12分贝、三阶巴特沃斯滤波器的衰减率为每倍频18分贝、如此类推。巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高，在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。2.2 传递函数巴特沃斯低通滤波器可用如下振幅的平方对频率的公式表示：其中, n= 滤波器的阶数 c= 截止频率 = 振幅下降为 -3分贝时的 频率 p= 通频带边缘频率 在通频带边缘的数值图 2-1 巴斯沃特滤波器3 用双线性变换法设计低通滤波器步骤设计数字滤波器的具体步骤如下：（1）确定数字低通滤波器的

29、技术指标：通带截止频率、通带最大衰减系数、阻带截止频率、阻带最小衰减系数。（2）将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。采用双线性变换法，频率的转换关系为 （3）根据转换后的性能指标，确定滤波器最小阶数n和固有频率wn。（4）由最小阶数n得到巴特沃斯模拟低通滤波器原型。（5）将模拟滤波器系统函数Ha(S)从s平面转换到z平面，得到数字低通滤波器系统函数H(z)。H(z)=Ha（S）|4 设计主要参数表1、Buttord 巴特沃思模拟滤波器阶数的获得N，Wnbuttord(Wp，Ws，Rp，Rs，s)s表示获取模拟滤波器的阶数Wp通带截止频率，Ws阻带截止频率Rp通带最大衰减，

30、Rs阻带最小衰减；N符合要求的滤波器最小阶数Wn为Butterworth滤波器固有频率(3dB)。2、buttap巴特沃思模拟滤波器的设计Z,P,K = BUTTAP(N)N为阶数，Z零点，P极点，K为增益3、zp2tf零极点增益模型到传递函数模型的转换B，Azp2tf(Z，P，K)；输人参数：Z，P，K分别表示零极点增益模型的零点、极点和增益；输出参数：B，A分别为传递函数分子和分母的多项式系数。4、Lp2lp 低通到低通b,a=lp2lp(B,A,Wn); B，A分别为截止频率为1的模拟滤波器传递函数分子和分母的多项式系数b，a分别为截止频率为Wn的模拟滤波器传递函数分子和分母的多项式系数

31、，5、Bilinear双线性变换法设计数字滤波器bz，a2二bilinear(b，a，Fs)；b，a分别为模拟滤波器传递函数分子和分母的多项式系数，Fs是采样频率bz，az分别为数字滤波器传递函数分子和分母的多项式系数6、Freqz 数字滤波器的频响特性H,W=freqz(bz,az); H为幅度，W为相位7、filter 滤波Y = filter(bz,az,X)bz，az分别为数字滤波器传递函数分子和分母的多项式系数X 为输入信号，Y为输出信号5 程序内容及结果方法一：%先设计模拟滤波器，再转化数字滤波器 wp=0.2*pi;%通带截止频率ws=0.3*pi;%阻带截止频率Rp=1;%通带

32、最大衰减Rs=15;%阻带最小衰减Ts=0.02*pi;Fs=1/Ts;%采样频率wp1=2/Ts*tan(wp/2);%将模拟指标转变成数字指标ws1=2/Ts*tan(ws/2);N,Wn=buttord(wp1,ws1,Rp,Rs,s); %选择滤波器的最小阶数 Z,P ,+K=buttap(N);%创建butterworth模拟滤波器Bap,Aap=zp2tf(Z,P,K);%零极点增益模型到传递函数模型的转换b,a=lp2lp(Bap,Aap,Wn);%低通到高通 bz,az=bilinear(b,a,Fs);%用双线性变换法实现模拟滤波器到数字滤波器的转换（具体原理） H,W=fr

33、eqz(bz,az,50);%绘制频率响应曲线L=length(W)/2+1;figure(1),plot(W(1:L)/pi,abs(H(1:L),grid,xlabel(角频率(pi),ylabel(频率响应幅度); x=-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,.-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,.-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0;y=filter(bz,az,x); %滤波figure(2

34、),subplot(2,1,1),plot(x),title(原始信号);subplot(2,1,2),plot(y),title(滤波后信号);figure(3)subplot(211)%原始信号频谱plot(abs(fft(x);title(原始信号频谱);xlabel(频率/Hz);ylabel(幅度);subplot(212)%滤波后信号频谱plot(abs(fft(y);title(滤波后信号频谱);xlabel(频率/Hz);ylabel(幅度);结果：图2-1 滤波图2-2 滤波前后信号图2-3 滤波前后信号频谱图方法二： %直接设计数字滤波器 wp=0.2*pi;%通带截止频率

35、 ws=0.3*pi;%阻带截止频率 Rp=1;%通带最大衰减 Rs=15;%阻带最小衰减 Ts=0.02*pi; Fs=1/Ts;%采样频率 wp1=0.2; %归一化指标 ws1=0.3; N,Wn=buttord(wp1,ws1,Rp,Rs); %选择滤波器的最小阶数 bz,az=butter(N,Wn);%Bilinear双线性变换法设计数字滤波器 H,W=freqz(bz,az);%绘制频率响应曲线 figure(1),plot(W*Fs/(2*pi),abs(H),grid,xlabel(频率Hz),ylabel(频率响应幅度); x=-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,.-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,.-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0;y=filter(bz,az,x); %滤波figure(2),subplot(2,1,1),plot(x),title(原始信号);subplot(2,1