课设说明书1.doc

课程设计说明书1 引言随着信息时代和数字世界的到来，数字信号处理己成为一门极其重要的学科和技术领域，它在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。数字滤波，是数字信号处理的基本核心内容之一，占有极重要的地位。它通过对采样数据信号进行数学运算处理来达到频域滤波目的，是图像处理、语音处理、模式识别、通信、软件无线电、谱分析等应用中的一个基本处理算法。与模拟滤波器相比，数字滤波器不用考虑器件的噪声、电压漂移、温度漂移等问题，可以容易的实现不同幅度和相位频率等特性指标，能够处理低频信号，频率响应特性可做成非常接近于理想的特性，且精度可以达到很高，容易集成等，这些优势决定了数字滤波器的应用很受欢迎。对于经典滤波器，从频域上可以分为低通、高通、带通、带阻滤波器。从时域特性上来看，数字滤波器还可以分为有限冲激响应数字滤波器(FIR)和无限冲激响应数字滤波器(IIR )。对于有限冲激响应数字滤波器(FIR)其输出y(n)只取决于有限个过去和现在的输入x(n),x(n-1). x(n-m);对于无限冲激响应数字滤波器(IIR )，它的输出不仅取决于过去和现在的输入，还跟过去的输出有关。本课题主要应用MATLAB软件设计FIR和IIR低通数字滤波器，并对所设计的滤波器进行仿真，使用 CCS 5402 Device Simulator 实现正弦波和方波数据输入, 使用FIR和IIR滤波器对输入信号进行滤波，输出相应的波形与频谱。文中采用循环寻址方法在CCS中实现FIR和IIR低通数字滤波器，在MATLAB中采用窗函数法设计FIR数字滤波器和用双线性变换法设计IIR数字滤波器，介绍利用DSP来快速设计数字滤波器的方法，寻找系数的快速传递。2 CCS简介CCS的全称是Code Composer Studio，它是美国德州仪器公司（Texas Instrument, TI）出品的代码开发和调试套件。TI公司的产品线中有一大块业务是数字信号处理器（DSP）和微处理器（MCU），CCS便是供用户开发和调试DSP和MCU程序的集成开发软件。CCS是一种针对TMS320系列DSP的集成开发环境,在Windows操作系统下，采用图形接口界面，提供有环境配置、源文件编辑、程序调试、跟踪和分析等工具。 CCS有两种工作模式，即软件仿真器和硬件在线编程。软件仿真器工作模式可以脱离DSP芯片，在PC上模拟DSP的指令集和工作机制，主要用于前期算法实现和调试。硬件在线编程可以实时运行在DSP芯片上，与硬件开发板相结合进行在线编程和调试应用程序。CCS的功能十分强大，它集成了代码的编辑、编译、链接和调试等诸多功能，而且支持C/C+和汇编的混合编程，其主要功能如下。 （1） 具有集成可视化代码编辑界面，用户可通过其界面直接编写C、汇编、.cmd文件等； （2） 含有集成代码生成工具，包括汇编器、优化C编译器、链接器等，将代码的编辑、编译、链接和调试等诸多功能集成到一个软件环境中； （3） 高性能编辑器支持汇编文件的动态语法加亮显示，使用户很容易阅读代码，发现语法错误； （4） 工程项目管理工具可对用户程序实行项目管理。在生成目标程序和程序库的过程中，建立不同程序的跟踪信息，通过跟踪信息对不同的程序进行分类管理； （5） 基本调试工具具有装入执行代码、查看寄存器、存储器、反汇编、变量窗口等功能，并支持C源代码级调试； （6） 断点工具，能在调试程序的过程中，完成硬件断点、软件断点和条件断点的设置； （7） 探测点工具，可用于算法的仿真，数据的实时监视等； （8） 分析工具，包括模拟器和仿真器分析，可用于模拟和监视硬件的功能、评价代码执行的时钟； （9） 数据的图形显示工具，可以将运算结果用图形显示,包括显示时域/频域波形、眼图、星座图、图像等，并能进行自动刷新； （10） 提供GEL工具。利用GEL扩展语言，用户可以编写自己的控制面板/菜单，设置GEL菜单选项，方便直观地修改变量，配置参数等。3 MATLAB简介MATLAB的全名是Matrix Laboratory，意思是矩阵实验室。早期的MATLAB是建立在DOS操作系统上的，直到20世纪90年代才发展成在Windows操作系统上的，同时它的功能也大为加强，不但可以进行数值计算，还具有从事符号解析运算、逻辑运算、数理统计、控制系统分析、最优化运算、金融的分析、数据的可视化、动力系统的建模和仿真等功能。MATLAB语言是一种高性能、用于科学和技术运算的计算机语言。它使得计算和图像一体化，并且编程容易，它所使用的环境是将问题和解答用大家熟悉的数学标记来表达。与MATLAB组合在一起的软件是Simulink，它是用来对动态系统进行建模、仿真和分析的软件，它支持连续、离散和非线性系统，它不能独立运行，必须在MATLAB环境下运行。MATLAB具有以下特点。（1） MATLAB能执行矩阵运算、符号运算、公式化简、线性和非线性方程式求解、高阶方程求根、线性规划、数理统计、微分和积分运算及自动控制系统的分析等，使得复杂的运算变得简单。（2） MATLAB可以直接在命令窗口进行运算，也可以在M文件窗口中进行编程运算，而无需编译和链接，并且可直接调用MATLAB中包含的数学函数和工程计算函数。（3） MATLAB能够按数据产生高质量的二维和三维数据图形，如散点图、直方图、饼图、树干图、向量图等，给用户提供既直接又精确的图像。（4） MATLAB设置了30多个工具箱，例如生物仪器工具箱、通信工具箱、控制工具箱、最优化工具箱、统计工具箱等。用户可以从工具箱中找到相应的帮助文本、应用举例和演示程序。4总体设计4.1 FIR低通滤波器4.1.1 FIR滤波器的设计原理FIR滤波器无反馈回路，是一种条件稳定系统，可以设计成具有线性相位特性。数字滤波一般具有如下的差分方程形式: （4-1） 式中， 为输入序列，为输出序列，和为滤波器系数，N是滤波器的阶数。若式中所有的均为零，且通常把系数记为,则有： （4-2）上式就是FIR滤波器的差分方程了。FIR滤波器的最主要的特点是没有反馈回路，因此它是无条件稳定系统。它的单位脉冲响应是一个有限长序列。由上面的方程可见，FIR滤波算法实际上是一种乘法累加运算，它不断地输入样本，经延时 ，做乘法累加，再输出滤波结果 。FIR数字滤波器的结构如图4-1所示:Z-1Z-1Z-1Z-1x(n)y(n)h(1)h(0)h(2)h(N-2)h(N-1)图4-1 数字滤波器的结构框图要设计一个FIR滤波器就是要求出它的冲击响应系数设计方法主要有窗函数法和频率抽样法，本设计采用的是Hamming窗函数法。4.1.2 FIR滤波器的设计方法本次设计的基本思路是利用CCS软件完成一个28阶FIR的低通滤波器的设计，通过对方波和不同谐波正弦波的叠加的滤波，完成对高频的滤除，剩下一次基波正弦波。其中在CCS软件中实现的两个不同长度的序列的时域卷积，即实现的是频域中的低通滤波。FIR低通滤波器的时域参数，不同正弦波的叠加所产生的序列，以及方波的序列参数均由MATLAB软件平台实现。通过加载到CCS编程过程中开辟的数据空间，由卷积程序实现滤波过程。实现滤波之后需要对所滤出的波形及数据进行检验，可通过CCS软件所提供的频谱进行分析；滤波之后的数据可通过与MATLAB卷积后的设计数据进行对比检验设计是否正确和符合设计要求。要设计一个FIR滤波器就是要求出它的冲击响应系数设计方法主要有窗函数法和频率抽样法，本设计采用的是Hamming窗函数法。 理想的低通滤波器的频率响应是一个矩形，这意味着它在时域上是无限长的序列，这在实际上是不可能实现的。因此我们要采取某种方法截断 ，可以用一个有限长度的窗函数序列与之相乘。这个窗函数序列的形状和长度都会对最后系统的频率响应特性产生影响，因此对窗函数的分析和选择是设计FIR滤波器的关键问题所在。4.1.3 FIR滤波器的设计流程图置28位空间Hn，BUF；256位空间Xn，Yn 将HDATA值给Hn空间 将Xn空间的首28个数据给BUF；并分别 将AR5，AR4指向Hn和BUF的空间尾地址设置循环缓冲区长度BK为28位，AR0赋值为-1，并将AR3指向Xn空间首地址将Hn和BUF中的数据进行乘累加运算，结果给Yn空间，AR4指向地址自加1将AR3指向的Xn空间的下一个数据给AR4指向的内容结束开始4.1.4 确定FIR滤波器参数用MATLAB确定滤波器系数程序如下。f=0 0.6 0.6 1;m=1 1 0 0;b=fir2(27,f,m)freqz(b,512,1000)由上面程序在MATLAB软件运行的28级系统函数所对应的时域系数为h = Columns 1 through 7 0.0003 -0.0017 -0.0035 -0.0026 0.0040 0.0133 0.0138Columns 8 through 14 -0.0044 -0.0347 -0.0486 -0.0122 0.0834 0.2042 0.2887Columns 15 through 21 0.2887 0.2042 0.0834 -0.0122 -0.0486 -0.0347 -0.0044 Columns 22 through 28 0.0138 0.0133 0.0040 -0.0026 -0.0035 -0.0017 0.0003滤波器的幅频特性和相频特性如图4-2所示图4-2 滤波器的幅频特性和相频特性4.1.5 用MATLAB确定FIR滤波器用MATLAB确定输入信号程序如下。 x=0:2*pi/256:2*pi; y1=square(1*pi*25*x)； x1=sin(2*pi*100*x); x2=sin(2*pi*1.5*x); x3=x2+x1; 在MATLAB中将输入的正弦与方波信号分别与滤波器进行卷积来验证滤波器，卷积的程序如下 x4=conv(x3,b); y2=conv(y1,b); 由上述程序得在MATLAB中用方波验证滤波器，波形如图4-3所示。图4-3 方波分解的正弦基波由图可看出，当输入为方波时，输出为正弦波，确定此滤波器为低通滤波器。在MATLAB中用正弦波验证滤波器，波形如图4-4所示。图4-4混叠正弦波分解的正弦基波 由图可看出，当输入为混叠的正弦波时，输出为低通的正弦信号，滤掉了高频成分，说明此滤波器为低通滤波器。4.1.6 用CCS确定FIR滤波器用CCS确定滤波器的主程序如下。 .title fir.asm .mmregs .global _mainHn: .usect Hn , 28BUF: .usect BUF, 28Xn: .usect Xn, 256 Yn: .usect Yn, 283 .dataHDATA: .word 3*32768/10000, -17*32768/10000 .word -35*32768/10000, -26*32768/10000 .word 40*32768/10000, 133*32768/10000 .word 138*32768/10000, -44*32768/10000 .word -347*32768/10000, -486*32768/10000 .word -122*32768/10000, 834*32768/10000 .word 2042*32768/10000, 2887*32768/10000 .word 2887*32768/10000, 2042*32768/10000 .word 834*32768/10000, -122*32768/10000 .word -486*32768/10000, -347*32768/10000 .word -44*32768/10000, 138*32768/10000 .word 133*32768/10000, 40*32768/10000 .word -26*32768/10000, -35*32768/10000 .word -17*32768/10000, 3*32768/10000 .text_main: SSBX FRCT STM #Hn,AR5 ;AR5指向冲激响应 RPT #27 MVPD HDATA,*AR5+ STM #BUF,AR4 ;AR4指向缓冲区间 RPTZ A,#27 STL A,*AR4+ STM #Xn,AR3 ;AR3指向输入 STM #Yn,AR2 ;AR4指向输出 STM #(Hn+27),AR5 ; STM #(BUF+27),AR4 STM -1,AR0 LD #BUF,DP MVDD *AR3+,*AR4 ;向缓冲去送数 STM #255,BRC RPTB #TABLE ;卷积算法 STM #28,BK RPTZ A,#27 MAC *AR4+0%,*AR5+0%,A STH A,*AR2+ TABLE: MVDD *AR3+,*AR4+0% .end由运行CCS所得对方波的滤波如图4-5所示。图4-5 由方波滤出基波正弦波由于方波输入频率过大，导致方波显示不太明了，但从图可看出，此滤波器和输入的方波基本匹配。由CCS查看方波的输入信号和输出信号频谱图分别如图4-6和图4-7所示。 图4-7 输出信号频谱图图4-6 输入方波频谱图 比较输入和输出信号的频谱可以看出，输入方波的傅里叶级数分解成不同的离散频率的正弦波的叠加，在基频的整数倍处分别存在不同频率的正弦波的频谱。而经过DSP程序滤波后的波形为基频正弦波，即实现了将高频滤除，将低频滤出的目的。再对照MATLAB中用方波验证的结果，可得到结论为：此滤波器为低通滤波器，且能找到和它匹配的方波信号。由运行CCS所得对正弦波的滤波如图4-6所示。图4-8 由正弦波滤出低频正弦波图4-6 由正弦波滤出低频正弦波由图可看出，不同频率正弦波的叠加可分解成不同的离散频率，当输入为混叠的正弦信号时，输出为低频的正弦波，滤掉了高频成分，因此此滤波器和输入的正弦波匹配。由CCS查看正弦波的输入信号和输出信号频谱图分别如图4-9和图4-10所示。图4-9 输入正弦波频谱图图4-10输出信号频谱图如同方波，比较输入和输出信号的频谱可以看出，输入含有低频和高频成分，输出只含有低频成分，说明滤波器滤掉了高频成分，再对照MATLAB中用方波验证的结果，可得到结论为：此滤波器为低通滤波器，能找到和它匹配的正弦波信号。4.2 IIR低通滤波器4.2.1 IIR滤波器的设计原理IIR 数字滤波器的特征是，具有无限持续时间冲激响应，需要用递归模型来实现,其差分方程为： （4-3）系统函数为： （4-4）设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数，使其频率响应满足所希望得到的频域指标，即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。4.2.2 IIR滤波器的设计方法IIR滤波器设计方法有两类,间接法和直接法。这两种方法又包括：双线性变换法和脉冲相应不变法。这两种方法各有优缺点，脉冲响应不变法的优点是频率坐标变换是线性，如果不考虑频率混叠现象，用这种方法设计的数字滤波器会很好的重现原模拟滤波器的频率特性。其主要缺点是会产生混叠现象，使数字滤波器的频响偏移模拟滤波器的频响。双线性变换法的优点是不可能产生频率混叠现象，这是双线性变换法比脉冲响应不变法的最大优点，但是易产生幅度特性和相位特性的失真。这次采用的是先设计Butterworth型IIR低通滤波，再经双线性变换，转化成IIR数字低通滤波器。在进行双线性变换时，将整个频率轴上的频率范围压缩到之间，在用转换到z平面上。为了实现s平面上整个虚轴安全压缩到s1轴之间的转换，利用正切变换实现频率压缩：式中，T是采样间隔，当从经过0变化到时，则由-经过0变化到+，映射了整个虚轴。此时在通过转换到z平面上，得到。从s1平面转换到z平面仍然采用标准转换关系，s1平面的之间水平带的左半部分映射z平面单位圆内部，虚轴映射单位圆，这样因果稳定，转换成的也因果稳定。4.2.3 确定IIR滤波器参数用MATLAB确定滤波器系数程序如下。 fc=4500; fs=30000; fb=3000; wb=fb*2*pi/fs; wc=fc*2*pi/fs; Ap=1; As=15; Wb=2*fs*tan(wb/2); Wc=2*fs*tan(wc/2); N,Wn=buttord(Wb,Wc,Ap,As,s);b,a=butter(N,Wn,s);bz,az=bilinear(b,a,fs);h,w=freqz(bz,az,256);plot(w*fs/(2*pi),20*log10(abs(h);grid on;axis(0 5000 -20 1);由上面程序在MATLAB软件运行后，得到滤波器的阶数为6和滤波器系数，并得到如图4-11所示的滤波器的幅相频特性。 图4-11 滤波器的幅相频特性4.2.4 用MATLAB确定IIR滤波器用MATLAB确定方波和正弦波输入信号程序如下。 x=0:2*pi/256:2*pi; y1=square(1*pi*5.1*x); x1=sin(2*pi*100*x); x2=sin(2*pi*1*x); x3=x2+x1;4.2.5 用CCS确定IIR滤波器用CCS确定滤波器的主程序如下。 .titleiir.asm .mmregs .def _c_int00 x .usect x,7y .usect y,7b .usect b,7a .usect a,7 output .usect output,320 input .usect intput,256 .data table .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 6*16384/10000 .word 35*16384/10000 .word 87*16384/10000 .word 116*16384/10000 .word 87*16384/10000 .word 35 *16384/10000 .word 6*16384/10000 .word 628*16384/10000 .word -5458*16384/10000 .word 20275*16384/10000 .word 41433*16384/10000 .word 49501*16384/10000 .word 33143*16384/10000_c_int00: ssbx frct stm #x,ar1 rpt #5 mvpd #table,*ar1 stm #y,ar1 rpt #5 mvpd #table+6,*ar1+ stm #b,ar1 rpt #6 mvpd table+12,*ar1+ stm #a,ar1 rpt #5 mvpd #table+19,*ar1+ stm #x+6,ar2 stm #a+5,ar3 stm #y+5,ar4 stm #b+6,ar5 stm #7,bk stm #-1,ar0 stm #input,ar6 stm #output,ar7 stm #255,brciir rptb next-1 ld *ar6+,a stl a,-1,*ar2 mpy *ar2+0%,*ar5+0%,a mac *ar2+0%,*ar5+0%,a mac *ar2+0%,*ar5+0%,a mac *ar2+0%,*ar5+0%,a mac *ar2+0%,*ar5+0%,a mac *ar2+0%,*ar5+0%,a mac *ar2,*ar5+0%,a mac *ar4+0%,*ar3+0%,a mac *ar4+0%,*ar3+0%,a mac *ar4+0%,*ar3+0%,a mac *ar4+0%,*ar3+0%,a mac *ar4+0%,*ar3+0%,a mac *ar4+0%,*ar3+0%,a mar *ar3+0% sth a,*ar4 sth a,*ar7+next: b next .end由运行CCS所得对方波的滤波如图4-12所示。图4-12 由方波滤出基波正弦波由图可看出，当输入为方波时，输出基本为正弦波，即此滤波器输入的方波基本匹配。由CCS查看方波的输入信号和输出信号频谱图分别如图4-13和图4-14所示。 图4-14输出信号频谱图图4-13 输入方波频谱图 同FIR的方波，输入方波的傅里叶级数分解成不同的离散频率的正弦波的叠加，在基频的整数倍处分别存在不同频率的正弦波的频谱。而经过DSP程序滤波后的波形为基频正弦波。再对照MATLAB中用方波验证的结果，可得到结论为：此滤波器为低通滤波器，且能找到和它匹配的方波信号。由运行CCS所得对正弦波的滤波如图4-15所示。图4-15 由混叠正弦波滤出低频正弦波由图可看出，当输入为混叠的正弦信号时，输出为低频的正弦波，滤掉了高频成分，因此此滤波器和输入的正弦波匹配。由CCS查看正弦波的输入信号和输出信号频谱图分别如图4-16和图4-17所示。图4-16 输入正弦波频谱图图4-17输出信号频谱图如同方波，比较输入和输出信号的频谱可以看出，输入含有低频和高频成分，输出只含有低频成分，说明滤波器滤掉了高频成分，再对照MATLAB中用方波验证的结果，可得到结论为：此滤波器为低通滤波器，能找到和它匹配的正弦波信号。5 总结两周的课程设计已经接近尾声，首先要感谢我们的学校领导为我们提供了多次动手锻炼的机会，增强了我们的实践能力。感谢此次设计的指导老师的悉心指导和帮助，感谢同学们的互相帮助，在他们的帮助和自己的努力下才使本次此课程设计能够顺利完成。在课设刚刚开始的几天里，我一直弄不懂这次课设到底该从哪里入手，该怎样入手。当从网上搜集了好多资料、从图书馆借了几本书、向老师寻问及与同学时时讨论后，我终于可以着手做课设了。首先是设计FIR低通滤波器，当设计完FIR低通滤波器后，IIR滤波器就简单多了。通过对本次数字低通滤波器的设计，我熟悉了MATLAB的运行环境，初步掌握了MATLAB语言在数字信号处理中一些基本库函数的调用和编写基本程序等应用；对利用循环寻址编程有了更进一步的熟悉；熟悉了滤波器设计的一般原理，对滤波器有了一个理性的认识；学会了数字低通滤波器设计的一般步骤。总之，使理论联系了实际，巩固并深化了对课本基本知识的认识和理解，使理论得以升华。参考文献1 范寿康 DSP技术与芯片M北京：电子工业出版社2 万永格 数字信号处理的MATLAB实现M北京：科学出版社，20073 程佩青 数字信号处理教程M北京：清华大学出版社，20014 高西全 、丁玉美等 数字信号处理M北京：电子工业出版社，20095 范寿康、 伊磊等 DSP原理及应用M北京：电子工业出版社，20096 胡广书 现代数字信号处理教程M北京：清华大学出版社，2005附录 1FIR 程序1 CCS主程序： .title fir.asm .mmregs .global _mainHn: .usect Hn , 28BUF: .usect BUF, 28Xn: .usect Xn, 256 Yn: .usect Yn, 283 .dataHDATA: .word 3*32768/10000,-17*32768/10000,-35*32768/10000, -26*32768/10000 .word 40*32768/10000, 133*32768/10000, 138*32768/10000, -44*32768/10000 .word -347*32768/10000, -486*32768/10000, -122*32768/10000, 834*32768/10000 .word 2042*32768/10000,2887*32768/10000,2887*32768/10000, 2042*32768/10000 .word 834*32768/10000, -122*32768/10000,-486*32768/10000, -347*32768/10000 .word -44*32768/10000,138*32768/10000,133*32768/10000, 40*32768/10000 .word -26*32768/10000, -35*32768/10000, -17*32768/10000, 3*32768/10000 .text_main: SSBX FRCT STM #Hn,AR5 ;AR5指向冲激响应 RPT #27 MVPD HDATA,*AR5+ STM #BUF,AR4 ;AR4指向缓冲区间 RPTZ A,#27 STL A,*AR4+ STM #Xn,AR3 ;AR3指向输入 STM #Yn,AR2 ;AR4指向输出 STM #(Hn+27),AR5 ; STM #(BUF+27),AR4 STM -1,AR0 LD #BUF,DP MVDD *AR3+,*AR4 ;向缓冲去送数 STM #255,BRC RPTB #TABLE ;卷积算法 STM #28,BK RPTZ A,#27 MAC *AR4+0%,*AR5+0%,A STH A,*AR2+TABLE: MVDD *AR3+,*AR4+0% .end中断向量文件：.title vectors.asm .ref _main .sect .vectors B _main .end连接命令文件： vectors.obj example.obj -o example.out -m example.mapMEMORY PAGE 0: EPROM: org=0e000h,len=1000h VECS: org=0ff80h,len=0080h PAGE 1: SPRAM: org=0080h,len=0200h DARAM: org=0281h,len=1400hSECTIONS .text:EPROM PAGE 0 .data:EPROM PAGE 0 Hn : align(32) DARAM PAGE 1 BUF : align(32) DARAM PAGE 1 Xn : DARAM PAGE 1 Yn : DARAM PAGE 1 .vectors:VECS PAGE 02 MATLAB程序b=fir1(72,0.66);freqz(b,1);figure(1);h1,w1=freqz(b,1);figure(2);title(频率特性)plot(w1/pi,20*log10(abs(h1);axis(0,1,-100,20);grid;xlabel(归一化频率/p)ylabel(幅度/dB)%产生测试用正弦信号和噪声信号figure(3)x=0:2*pi/1024:2*pi;x1=4*square(1*pi*40*x);y1=sin(1*x*2*pi+pi/6);y2=sin(100*x*2*pi);y3=y1+y2;y4=conv(y3,b);x2=conv(x1,b);subplot(6,1,1);plot(y1);subplot(6,1,2);plot(y2);subplot(6,1,3);plot(y3);subplot(6,1,4);plot(y4);subplot(6,1,5);plot(x1);subplot(6,1,6);plot(x2);附录 2IIR 程序1 CCS主程序： .titleiir.asm .mmregs .def _c_int00x .usect x,7y .usect y,7b .usect b,7a .usect a,7 output .usect output,320input .usect intput,256 .datatable .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 0 .word 6*16384/10000 .word 35*16384/10000 .word 87*16384/10000 .word 116*16384/10000 .word 87*16384/10000 .word 35 *16384/10000 .word 6*16384/10000 .word 628*16384/10000 .word -5458*16384/10000 .word 20275*16384/10000 .word 41433*16384/10000 .word 49501*16384/10000 .word 33143*16384/10000_c_int00: ssbx frct stm #x,ar1 rpt #5 mvpd #table,*ar1+ stm #y,ar1 rpt #5 mvpd #table+6,*ar1+ stm #b,ar1 rpt #6 mvpd table+12,*ar1+ stm #a,ar1 rpt #5 mvpd #table+19,*ar1+ stm #x+6,ar2 stm #a+5,ar3 stm #y+5,ar4 stm #b+6,ar5 stm #7,bk stm #-1,ar0 stm #input,ar6 stm #output,ar7 stm #255,brciir rptb next-1 ld *ar6+,a stl a,-1,*ar2 mpy *ar2+0%,*ar5+0%,a mac *ar2+0%,*ar5+0%,a mac *ar2+0%,*ar5+0%,a mac *ar2+