dsp滤波要点
数字信号处理DSP数字滤波器的基本概念及一些特殊滤波器课件
件
目录
• 数字滤波器的基本概念 • 数字滤波器的实现方法 • 特殊数字滤波器介绍 • 数字滤波器的优化与设计 • 数字滤波器的应用实例
01
数字滤波器的基本概念
数字滤波器的定义与分类
数字滤波器定义为能够实现离散时间信号处理的算法或 方程，通常在数字信号处理系统中用于改善信号的质量 和特征。
03
特殊数字滤波器介绍
梳状滤波器
总结词
减小频率范围
详细描述
梳状滤波器是一种特殊类型的数字滤波器，其频率响应类似于“梳子”，在一 定频率范围内减小了信号的传递，而在这个范围之外则允许信号通过。这种滤 波器通常用于减小信号中的高频噪声。
陷波滤波器
总结词
抑制特定频率
详细描述
陷波滤波器是一种特殊的数字滤波器，其频率响应类似于“陷坑”，在某一特定 频率处完全抑制信号的传递，而在这个频率之外则允许信号通过。这种滤波器通 常用于消除信号中的特定频率成分。
数字滤波器的应用场景与优势
数字滤波器广泛应用于图像处理、语音信号处理、雷达信号处理等领域。
数字滤波器的优势在于能够克服模拟滤波器的一些缺点，如易受干扰、精度低、不易复制等， 同时具有处理速度快、精度高、稳定性好等优点。
02
数字滤波器的实现方法
IIR数字滤波器的实现方法
直接形式
通过串联、并联或反馈连接的方式将基本运算单元（如 加法器、乘法器和延迟器）组合起来，构成IIR数字滤 波器的系统函数。
在图像处理中的应用实例
图像去噪
数字滤波器可以用于图像信号的去噪。例如，可以使用适应性滤波器来消除图像中的噪声 和干扰，或者使用形态学滤波器来填补图像中的空洞和去除小的噪声点。
DSP课程设计滤波器
DSP课程设计--滤波器目录摘要...................................... 错误!未定义书签。
前言.. (1)1 方案设计与论证 (2)1.1 设计方案概论 (2)1.2 设计方案详论 (2)1.3 设计工具CCS及SEED-DTK2812 实验系统简介 (3)2 系统设计 (4)2.1 IIR数字滤波器的设计方法及原理 (4)2.2 程序设计流程图 (6)2.3 系统设计步骤 (7)4 总结 (10)参考文献 (12)致谢 (13)附录 (14)前言本文介绍了滤波器的滤波原理以及模拟滤波器、数字滤波器的设计方法。
重点介绍了IIR数字滤波器的设计方法。
即脉冲响应不变法和双线性变换法。
在此基础上，用DSP虚拟实现任意阶IIR滤波器。
此设计扩展性好，便于调节滤波器的性能，可以根据不同的要求在DSP上加以实现。
数字信号处理是利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理，以得到符合人们需要的信号形式。
数字信号处理是围绕着数字信号处理的理论、实现和应用等几个方面发展起来的。
数字信号处理在理论上的发展推动了数字信号处理应用的发展。
反过来，数字信号处理的应用又促进了数字信号处理理论的提高。
而数字信号处理的实现则是理论和应用之间的桥梁。
数字信号处理是以众多学科为理论基础的，它所涉及的范围极其广泛。
例如，在数学领域，微积分、概率统计、随机过程、数值分析等都是数字信号处理的基本工具，与网络理论、信号与系统、控制论、通信理论、故障诊断等也密切相关。
近来新兴的一些学科，如人工智能、模式识别、神经网络等，都与数字信号处理密不可分。
可以说，数字信号处理是把许多经典的理论体系作为自己的理论基础，同时又使自己成为一系列新兴学科的理论基础。
数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。
DSP几种滤波算法比较
第1种方法限幅滤波法（又称程序判断滤波法）A 方法根据经验判断，确定两次采样允许的最大偏差值（设为A）每次检测到新值时判断：如果本次值与上次值之差<=A,则本次值有效如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值B优点能有效克服因偶然因素引起的脉冲干扰C 缺点无法抑制那种周期性的干扰平滑度差第2种方法中位值滤波法A方法连续采样N次（N取奇数）把N次采样值按大小排列取中间值为本次有效值B优点能有效克服因偶然因素引起的波动干扰对温度、液位的变化缓慢的被测参数有良好的滤波效果C缺点对流量、速度等快速变化的参数不宜第3种方法算术平均滤波法A方法连续取N个采样值进行算术平均运算N值较大时：信号平滑度较高，但灵敏度较低N值较小时：信号平滑度较低，但灵敏度较高N值的选取：一般流量，N=12；压力：N=4B优点适用于对一般具有随机干扰的信号进行滤波这样信号的特点是有一个平均值，信号在某一数值范围附近上下波动C缺点对于测量速度较慢或要求数据计算速度较快的实时控制不适用比较浪费RAM第4种方法递推平均滤波法（又称滑动平均滤波法）A方法把连续取N个采样值看成一个队列队列的长度固定为N每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据.(先进先出原则)把队列中的N个数据进行算术平均运算,就可获得新的滤波结果N值的选取：流量，N=12；压力：N=4；液面，N=4~12；温度，N=1~4B优点对周期性干扰有良好的抑制作用，平滑度高适用于高频振荡的系统C缺点灵敏度低对偶然出现的脉冲性干扰的抑制作用较差不易消除由于脉冲干扰所引起的采样值偏差不适用于脉冲干扰比较严重的场合比较浪费RAM第5种方法中位值平均滤波法（又称防脉冲干扰平均滤波法）A方法相当于“中位值滤波法”+“算术平均滤波法”连续采样N个数据，去掉一个最大值和一个最小值然后计算N-2个数据的算术平均值N值的选取：3~14B优点融合了两种滤波法的优点对于偶然出现的脉冲性干扰，可消除由于脉冲干扰所引起的采样值偏差C缺点测量速度较慢，和算术平均滤波法一样比较浪费RAM第6种方法限幅平均滤波法相当于“限幅滤波法”+“递推平均滤波法” 每次采样到的新数据先进行限幅处理再送入队列进行递推平均滤波处理B 优点融合了两种滤波法的优点对于偶然出现的脉冲性干扰，可消除由于脉冲干扰所引起的采样值偏差C 缺点比较浪费RAM第7种方法一阶滞后滤波法A方法取a=0~1本次滤波结果=（1-a）*本次采样值+a*上次滤波结果B优点对周期性干扰具有良好的抑制作用适用于波动频率较高的场合C缺点相位滞后，灵敏度低滞后程度取决于a值大小不能消除滤波频率高于采样频率的1/2的干扰信号第8种方法加权递推平均滤波法A方法是对递推平均滤波法的改进，即不同时刻的数据加以不同的权通常是，越接近现时刻的资料，权取得越大给予新采样值的权系数越大，则灵敏度越高，但信号平滑度越低B优点适用于有较大纯滞后时间常数的对象和采样周期较短的系统C缺点对于纯滞后时间常数较小，采样周期较长，变化缓慢的信号不能迅速反应系统当前所受干扰的严重程度，滤波效果差第9种方法消抖滤波法A方法设置一个滤波计数器将每次采样值与当前有效值比较：如果采样值＝当前有效值，则计数器清零如果采样值<>当前有效值，则计数器+1，并判断计数器是否>=上限N(溢出)如果计数器溢出,则将本次值替换当前有效值,并清计数器B优点对于变化缓慢的被测参数有较好的滤波效果,可避免在临界值附近控制器的反复开/关跳动或显示器上数值抖动C缺点对于快速变化的参数不宜如果在计数器溢出的那一次采样到的值恰好是干扰值,则会将干扰值当作有效值导入系统第10种方法限幅消抖滤波法A方法相当于“限幅滤波法”+“消抖滤波法”先限幅后消抖B优点继承了“限幅”和“消抖”的优点改进了“消抖滤波法”中的某些缺陷,避免将干扰值导入系统C缺点对于快速变化的参数不宜第11种方法IIR数字滤波器A方法确定信号带宽，滤之。
DSP应用课程设计课件第4讲利用DSP实现信号滤波
实现FIR滤波器的方法
FIR滤波器是非递归滤波器，可以通过对输入信号的加权平均来滤波信号。常 用的实现方法有直接形式、线性相位形式和FFT快速算法。
滤波器的性能评估指标
滤波器的性能可以通过频率响应、相位响应、群延迟、截止频率、通带波动 等指标进行评估。这些指标可以帮助我们选择适合的滤波器。
滤波器的设计流程概述
滤波器设计工具的使用
现代DSP平台提供了一系列工具和算法，用于滤波器设计、优化和实现。MATLAB等软件还提供滤波器 设计函数和可视化工具，简化了设计流程。
MATLAB中滤波器设计方法的演示
MATLAB提供了多种滤波器设计函数和工具箱，可以方便地实现各种滤波器的设计和验证。通过演示这 些方法，可以更好地理解滤波器设计的过程。
常见的滤波器类型介绍
低通滤波器
允许低于截止频率的信号通过。
带通滤波器
只允许特定频率范围内的信号通过。
高通滤波器
允许高于截止频率的信号通过。
带阻滤波器
阻止特定频率范围内的信号通过。
实现IIR滤波器的方法
IIR滤波器基于递归算法，通过使用过去的输出作为当前输入的一部分来滤波 信号。常用的实现方法有直接形式、级联形式和并联形式。
滤波器的实现方式及优化方法
滤波器可以通过离散时间系统、差分方程、传输函数等多种方式进行实现。 优化方法包括结构优化、编程优化等，以提高滤波器的性能和效率。
离散时间滤波器与连续时间滤 波器的关系
离散时间滤波器是连续时间滤波器的一种离散采样版本。通过离散化连续时 间滤波器的传输函数，可以获得离散时间滤波器的差分方程和频率响应。
带通滤波器的应用实例
带通滤波器可以用于音频信号音乐分析、图像纹理提取、数据频段分析等应 用。通过选择合适的通带范围和滤波器参数，可以实现对特定频率范围内的 信号处理。
数字信号处理中的频域滤波方法
数字信号处理中的频域滤波方法数字信号处理（Digital Signal Processing，DSP）是一门研究如何对数字信号进行变换、操作和分析的学科。
其中，频域滤波方法是一种常用的信号处理技术，用于去除信号中的噪声或改善信号质量。
本文将介绍数字信号处理中的频域滤波方法，包括傅里叶变换、傅里叶变换的性质以及滤波器设计。
一、傅里叶变换傅里叶变换是一种将信号从时域（时序）转换到频域（频率）的方法，它将信号表示为正弦和余弦函数的线性组合。
傅里叶变换可以将信号分解为不同频率成分的和，通过分析这些频率成分可以实现频域滤波。
在数字信号处理中，傅里叶变换通常使用离散傅里叶变换（Discrete Fourier Transform，DFT）来实现。
DFT将连续时域信号离散化为一系列离散频率，从而可以在计算机上进行处理。
二、傅里叶变换的性质1. 线性性质：傅里叶变换具有线性性质，即信号的线性组合的傅里叶变换等于信号各自的傅里叶变换的线性组合。
2. 积移性质：信号在时域上的平移会导致其在频域上的相位变化，即频谱随时间的平移而变化。
3. 对称性质：实信号的傅里叶变换具有共轭对称性，即其频谱是一个关于零频率对称的函数。
三、频域滤波器设计频域滤波器是根据信号在频域的特性来选择和调整信号成分的方法。
常见的频域滤波器有低通滤波器、高通滤波器、带通滤波器和带阻滤波器。
1. 低通滤波器：低通滤波器用于去除高频成分，只保留低频成分。
在频域上，低通滤波器会在截止频率以下的频率范围内透传，而在截止频率以上的频率范围内抑制信号。
2. 高通滤波器：高通滤波器用于去除低频成分，只保留高频成分。
高通滤波器在截止频率以下的频率范围内抑制信号，而在截止频率以上的频率范围内透传。
3. 带通滤波器：带通滤波器用于滤除不在指定频率范围内的信号。
它可以让指定范围的频率通过，而将其他频率抑制。
4. 带阻滤波器：带阻滤波器用于滤除指定频率范围内的信号。
它可以让指定范围外的频率通过，而将指定范围内的频率抑制。
在示波器上使用DSP滤波技术的探讨
在示波器上使用DSP滤波技术的探讨在今天的科技遍地开花的社会中，数字信号处理技术的应用越来越广泛。
其中，DSP滤波技术在很多领域得到了广泛的应用，而在示波器上也不例外。
在本文中，我们将探讨在示波器上使用DSP滤波技术的相关问题。
一、DSP滤波技术概述DSP滤波技术是数字信号处理中的一项重要技术，它是指对数字信号进行滤波的处理方法。
它可以通过数字滤波器来实现，数字滤波器将输入信号分成不同的频率段，然后对不同的频率段进行滤波，以去除或衰减滤波器不希望的频率成分。
数字信号处理技术能够对数字信号进行高速处理和现实时间处理，不易受到噪声和失真的干扰。
与传统的模拟滤波器相比，它的设计和调试过程更加简单，能够快速地进行实现。
二、在示波器上使用DSP滤波技术示波器是一种常用的测试仪器，用于观察和分析电信号的波形和参数。
它能够显示信号的幅度、频率、相位和脉冲等相关信息，对于电子行业的研发、测试和维修都有着重要的作用。
在示波器上使用DSP滤波技术，可以对输入的电信号进行数字滤波处理，以去除或衰减噪声、干扰等对于信号质量的负面影响。
可以有效提高示波器信号的准确度和可靠性，保证测试结果的正确性。
三、DSP滤波器的分类在示波器上使用的DSP滤波器可以分为一阶数字滤波器、二阶数字滤波器、三阶数字滤波器等几种类型。
1、一阶数字滤波器一阶数字滤波器是指只有一个延迟器和一个加权器的数字滤波器，其部件数最少。
它可实现低通滤波、高通滤波、带通滤波和带阻滤波等基本型式的滤波。
2、二阶数字滤波器二阶数字滤波器通常由两个一阶数字滤波器级联而成，可以实现更多的滤波器类型，如butterworth滤波器、elliptical滤波器等。
3、三阶数字滤波器三阶数字滤波器通常由三个一阶数字滤波器级联而成，还可以实现更多滤波器类型，如Chebyshev滤波器、Bessel滤波器等。
四、DSP滤波技术的应用在示波器上使用DSP滤波技术主要应用于以下几个方面：1、信号处理在电子设计中，信号必须在其范围内进行滤波，以达到电子设计的成功。
dsp滤波要点
dsp滤波要点数字信号处理综合设计性实验报告学院：电子信息工程学院班级：自动化0706班电子0701班指导教师：高海林学生：张越07213056陈冠宇07214004北京交通大学电工电子教学基地2010年1月20日利用DSP实现信号滤波一：实验目的实现离散信号滤波是DSP的基本功能，本实验中我们尝试实现分别了设计FIR 和IIR滤波器实现低通，高通，带通，带阻四种滤波器对正弦离散信号进行滤波。
（1）学会编写滤波程序和输入信号程序。
（2）熟悉CCS集成开发环境，熟悉DSK板的使用。
（3）通过实验比较FIR和IIR在设计上和滤波效果上的区别。
二：FIR滤波器的设计例：设计一个采样频率Fs为8000Hz，输入信号频率为1000HZ、2500HZ与4000HZ 的合成信号，通过设计FIR滤波器分别实现低通，高通，带通，带阻的滤波功能。
一：实验原理一个线性位移不变系统的输出序列y(n)和输入序列x(n)之间的关系，应满足常系数线性差分方程：x(n): 输入序列，y(n): 输出序列，ai、bi : 滤波器系数，N: 滤波器的阶数。
在式上式中，若所有的ai 均为0，则得FIR 滤波器的差分方程：对上式进行z 变换，可得FIR 滤波器的传递函数：FIR 滤波器的结构#include <stdio.h>#include <math.h>void main(){int i;double f[256];FILE *fp;if((fp=fopen("in.inc","wt"))==NULL){printf("can't open file! \n");return;}for(i=0;i<=255;i++){f[i]=sin(2*3.14159*i*1000/8000)+sin(2*3.14159*i*2500/8000)+sin(2* 3.14159*4000/8000);fprintf(fp," .word %1d\n",(long)(f[i]*32768/3));}fclose(fp);}该文件会生成一个.inc的文件，由于比较冗长，在此不再赘述。
利用DSP实现自适应滤波
序列值，M 是串联的总级数。格式结构的优点是按阶递归，故增加或减少级数不会影响存在 的阶数设计。
DSP 课程设计
设计报告
设 计 题 目 ： 利 用 DSP 实 现 自 适 应 滤 波
2013 年 7 月 12 日
DSP 系统课程设计
目
录
1 设计任务要求………………………………………………………1
2 理论分析…………………………………………………………
2.1 TMS320C55x 概况……………………………………………………… 2.2 自适应滤波器原理…………………………………………………… 2.3 自适应滤波器结构…………………………………………………… 2.4 LMS 算法与 RLS 算法………………………………………………… 2.5 自适应滤波器的应用…………………………………………………
尽管 LMS 算法的应用广泛，但是它也有缺点。当环境噪声不是平稳随机信号时，LMS 算 法很难自适应的跟踪统计特性变化的外部噪声干扰，因而其收敛效果一般。而基于 RLS 的自 适应滤波算法克服了上述缺点，能在非平稳环境下取得较满意的滤波效果。RLS 算法是在 LMS 算法的基础上而来的，所不同的是在求均方误差时观测数据的长度是变化的。且随着观测数 据的时间先后顺序分别乘了加权因子。即 RLS 算法的均方误差变为：
k
2
ξ（k）= β( k,n) (n)
n1
式中：β（k，n）是加权因子，满足 0＜β（k，n）≤1，n=1，2，…，k，也称作遗忘 因子。这样会使多次迭代之前的信号被遗弃掉，当滤波器工作于非平稳环境时，观测数据仍
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
.mmregs
.global start
.def start,_c_int00
KS .set 256 ;输入样本数据个数
N .set 46 ;FIR滤波器阶数
COEF_FIR .sect "COEF_FIR" ;FIR滤波器系数
fprintf(fp," .word %1d\n",(long)(f[i]*32768/3));
}
fclose(fp);
}
该文件会生成一个.inc的文件，由于比较冗长，在此不再赘述。然后，在DSP汇编语言中通过.copy汇编命令将生成的数据文件iirin.inc复制到汇编程序中，作为IIR滤波器的输入数据。此过程的作用相当于把信号通过了AD转换。
BOS .usect "STACK",0Fh
TOS .usect "STACK",1
.text
.asg AR4,DATA_P ;数据x(n)缓冲区指针
.asg AR6,INBUF_P ;模拟输入数据指针
.asg AR7,OUTBUF_P ;FIR滤波器输出数据指针
_c_int00
b start
nop
.word 193, 108, -74, -635, -620, 112, 784, 773, 94
.word -533, -431, -18, -243, -946, -741, 1064, 2933, 2307
.word -1241, -4633, -4091, 551, 5088, 5088, 551, -4091, -4633
RPTZ A,N-1
MACD *DATA_P-,COEF_FIR,A ;FIR滤波运算
MAR *DATA_P+ ;调整DATA_P指针指向DATABUF第一个单元
STH A,*OUTBUF_P+
LOOP:
EEND B EEND
.end
（4）带阻滤波器
*********一个FIR带阻滤波器程序bandstop.asm************
DATABUF .usect "FIR_BFR",N
BOS .usect "STACK",0Fh
TOS .usect "STACK",1
.text
.asg AR4,DATA_P ;数据x(n)缓冲区指针
.asg AR6,INBUF_P ;模拟输入数据指针
.asg AR7,OUTBUF_P ;FIR滤波器输出数据指针
.word -1241, 2307, 2933, 1064, -741, -946, -243, -18, -431
.word -533, 94, 773, 784, 112, -620, -635, -74, 108, 193
.data
INPUT .copy "bandpass.inc" ;输入数据在数据区0x2400
程序算法实现采用循环缓冲区法。
二：算法原理：
①在数据存储器中开辟一个N个单元的缓冲区(滑窗)，用来存放最新的N个输入样本；
②从最新样本开始取数；
③读完最老样本后，输入最新样本来代替最老样本，而其他数据位置不变；
④用BK寄存器对缓冲区进行间接寻址，使缓冲区地址首尾相邻。
三：实验程序设计步骤
1.FIR滤波器设计
.word -263, 396, 572, -142, -806, -319, 840, 949, -519
.word -1650, -355, 2294, 2303, -2746, -10014, 19293, -10014, -2746
.word 2303, 2294, -355, -1650, -519, 949, 840, -319, -806
.mmregs
.global start
.def start,_c_int00
KS .set 256 ;输入样本数据个数
N .set 46 ;FIR滤波器阶数
COEF_FIR .sect "COEF_FIR" ;FIR滤波器系数
.word 297, 226, 180, 24, -213, -451, -585, -532, -271
**************使用线性缓冲区法************************
由此即可得到16位二进制数据表示的离散数据文件，即离散的三频谱叠加的正弦信号。
3.编写滤波器的汇编源程序
（1）低通滤波器
**********一个FIR低通滤波器程序lowpass.asm************
**************使用线性缓冲区法************************
FIR滤波器的单位冲激响应h(n)为有限长序列。
若h(n)为实数。
在数字滤波器中，FIR滤波器具有如下几个主要特点：
①FIR滤波器无反馈回路，是一种无条件稳定系统；
②FIR滤波器可以设计成具有线性相位特性。
本实验程序设计的就是一种偶对称的线性相位滤波器。
.asg AR7,OUTBUF_P ;FIR滤波器输出数据指针
_c_int00
b start
nop
nop
start: SSBx FRCT ;小数乘法编程时，设置FRCT（小数方式）位
STM #DATABUF,DATA_P ;数据缓冲区清零
RPTZ A,N-1
STL A,*DATA_P+
FIR_TASK:
（1）学会编写滤波程序和输入信号程序。
（2）熟悉CCS集成开发环境，熟悉DSK板的使用。
（3）通过实验比较FIR和IIR在设计上和滤波效果上的区别。
二：FIR滤波器的设计
例：设计一个采样频率Fs为8000Hz，输入信号频率为1000HZ、2500HZ与4000HZ的合成信号，通过设计FIR滤波器分别实现低通，高通，带通，带阻的滤波功能。
一：实验原理
一个线性位移不变系统的输出序列y(n)和输入序列x(n)之间的关系，应满
足常系数线性差分方程：
x(n):输入序列，y(n):输出序列，ai、bi:滤波器系数，N:滤波器的阶数。
在式上式中，若所有的ai均为0，则得FIR滤波器的差分方程：
对上式进行z变换，可得FIR滤波器的传递函数：
FIR滤波器的结构
_c_int00
b start
nop
nop
start: SSBx FRCT ;小数乘法编程时，设置FRCT（小数方式）位
STM #DATABUF,DATA_P ;数据缓冲区清零
RPTZ A,N-1
STL A,*DATA_P+
FIR_TASK:
STM #DATABUF,DATA_P
STM #INPUT,INBUF_P
利用MATLAB中的FDA工具导出正确的参数。
2.产生滤波器输入信号文件
以下是一个产生输入信号的C语言程序，信号是频率为1000Hz、2500Hz和4000Hz的正弦波合成的波形：
#include <stdio.h>
#include <math.h>
void main()
{
int i;
double f[256];
数字信号处理
综合设计性实验报告
学 院：电子信息工程学院
班 级：自动化0706班
电子0701班
指导教师：高海林
学 生：张越 07213056
陈冠宇07214004
北京交通大学电工电子教学基地
2010年1月20日
利用DSP实现信号滤波
一：实验目的
实现离散信号滤波是DSP的基本功能，本实验中我们尝试实现分别了设计FIR和IIR滤波器实现低通，高通，带通，带阻四种滤波器对正弦离散信号进行滤波。
.word -142, 572, 396, -263, -448, -6, 371, 211, -184
.word -284, 30, 632, -358
.data
INPUT .copy "highpass.inc" ;输入数据在数据区0x2400
OUTPUT .space 1024 ;输出数据在数据区0x2500
nop
start: SSBx FRCT ;小数乘法编程时，设置FRCT（小数方式）位
STM #DATABUF,DATA_P ;数据缓冲区清零
RPTZ A,N-1
STL A,*DATA_P+
FIR_TASK:
STM #DATABUF,DATA_P
STM #INPUT,INBUF_P
STM #OUTPUT,OUTBUF_P
.word 125, 505, 686, 535, 40, -651, -1260, -1456, -979
.word 248, 2064, 4092, 5839, 6848, 6848, 5839, 4092, 2064
.word 248, -979, -1456, -1260, -651, 40, 535,686, 505
STM #KS-1,BRC
RPTB LOOP-1
LD *INBUF_P+,A ;装载输入数据
FIR_FILTER: ;FIR滤波运算
STL A,*DATA_P ;输入样本值x(n)
STM #(DATABUF+N-1),DATA_P ;数据缓冲区指针指向x[n-(N-1)]
RPTZ A,N-1
MACD *DATA_P-,COEF_FIR,A ;FIR滤波运算
OUTPUT .space 1024 ;输出数据在数据区0x2500