[DSP学习笔记]基于TMS320F28335的IIR滤波实现

已于 2023-08-31 13:58:35 修改
阅读量1k 收藏 25
点赞数 3
文章标签： 学习 笔记 dsp开发 算法
于 2023-08-05 23:31:07 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/luoqicou/article/details/132069652
版权

一、前言

IIR滤波器，递归滤波器，顾名思义，IIR滤波器采用递归型结构，结构上带有反馈回路。IIR数字滤波器相较于FIR数字滤波器具备幅频特性精度高的特点。

二、移植官方库

首先我们先将C2000 wave中我们所需要的函数库给导入到我们的新建工程中（可见于我实现FIR滤波文章中导入函数库的操作）。

工程中可见有以下文件。（仅进行滤波可删除FFT有关文件）

 同时记得添加头文件路径

 至此，对于我们进行滤波所需的函数库已经移植完毕。

 三、借助MATLAB生成测试信号

使用MATLAB生成输入测试信号（相关程序可见于我另两篇文章）

%这里仅给出设定信号相关参数
fs=1e5; %设采样率为100K
f = [1e3,5e3,1e4]; %频率矩阵（1*N）测试信号频率为1K，5K，10K
phase = [0,0,0]; %初相位均为0
amp = [0.3,1.0,0.5];%幅度
N=512; %采样点数

生成信号的时域图与频域图如下：

将刚刚MATLAB生成的record.txt复制到input.c（没有就新建一个）中。

这里笔者给出自己生成好的输入信号

//#############################################################################
//! \file   input.c
//! \brief  Input Vector (1024) 
//! \author Vishal Coelho 
//! \date   16-Sep-2016
//! 
//
//  Group:          C2000
//  Target Family:  $DEVICE$
//
//#############################################################################
//
//
// $Copyright: Copyright (C) 2022 Texas Instruments Incorporated -
//             http://www.ti.com/ ALL RIGHTS RESERVED $
//#############################################################################
float test_input[512] = {
                         0.000000,0.621747,1.100913,1.340760,
                         1.319556,1.092705,0.767601,0.461223,
                         0.256783,0.175872,0.176336,0.176103,
                         0.093107,-0.114798,-0.426010,-0.757295,
                         -0.991651,-1.021653,-0.791865,-0.323977,
                         0.285317,0.893485,1.358000,1.582180,
                         1.544357,1.300000,0.956572,0.631123,
                         0.406943,0.305699,0.285317,0.263809,
                         0.159191,-0.070597,-0.403866,-0.757295,
                         -1.013795,-1.065855,-0.857949,-0.411682,
                         0.176336,0.763658,1.207840,1.412279,
                         1.355387,1.092705,0.731771,0.389703,
                         0.149857,0.033962,0.000000,-0.033962,
                         -0.149857,-0.389703,-0.731771,-1.092705,
                         -1.355387,-1.412279,-1.207840,-0.763658,
                         -0.176336,0.411682,0.857949,1.065855,
                         1.013795,0.757295,0.403866,0.070597,
                         -0.159191,-0.263809,-0.285317,-0.305699,
                         -0.406943,-0.631123,-0.956572,-1.300000,
                         -1.544357,-1.582180,-1.358000,-0.893485,
                         -0.285317,0.323977,0.791865,1.021653,
                         0.991651,0.757295,0.426010,0.114798,
                         -0.093107,-0.176103,-0.176336,-0.175872,
                         -0.256783,-0.461223,-0.767601,-1.092705,
                         -1.319556,-1.340760,-1.100913,-0.621747,
                         -0.000000,0.621747,1.100913,1.340760,
                         1.319556,1.092705,0.767601,0.461223,
                         0.256783,0.175872,0.176336,0.176103,
                         0.093107,-0.114798,-0.426010,-0.757295,
                         -0.991651,-1.021653,-0.791865,-0.323977,
                         0.285317,0.893485,1.358000,1.582180,
                         1.544357,1.300000,0.956572,0.631123,
                         0.406943,0.305699,0.285317,0.263809,
                         0.159191,-0.070597,-0.403866,-0.757295,
                         -1.013795,-1.065855,-0.857949,-0.411682,
                         0.176336,0.763658,1.207840,1.412279,
                         1.355387,1.092705,0.731771,0.389703,
                         0.149857,0.033962,0.000000,-0.033962,
                         -0.149857,-0.389703,-0.731771,-1.092705,
                         -1.355387,-1.412279,-1.207840,-0.763658,
                         -0.176336,0.411682,0.857949,1.065855,
                         1.013795,0.757295,0.403866,0.070597,
                         -0.159191,-0.263809,-0.285317,-0.305699,
                         -0.406943,-0.631123,-0.956572,-1.300000,
                         -1.544357,-1.582180,-1.358000,-0.893485,
                         -0.285317,0.323977,0.791865,1.021653,
                         0.991651,0.757295,0.426010,0.114798,
                         -0.093107,-0.176103,-0.176336,-0.175872,
                         -0.256783,-0.461223,-0.767601,-1.092705,
                         -1.319556,-1.340760,-1.100913,-0.621747,
                         -0.000000,0.621747,1.100913,1.340760,
                         1.319556,1.092705,0.767601,0.461223,
                         0.256783,0.175872,0.176336,0.176103,
                         0.093107,-0.114798,-0.426010,-0.757295,
                         -0.991651,-1.021653,-0.791865,-0.323977,
                         0.285317,0.893485,1.358000,1.582180,
                         1.544357,1.300000,0.956572,0.631123,
                         0.406943,0.305699,0.285317,0.263809,
                         0.159191,-0.070597,-0.403866,-0.757295,
                         -1.013795,-1.065855,-0.857949,-0.411682,
                         0.176336,0.763658,1.207840,1.412279,
                         1.355387,1.092705,0.731771,0.389703,
                         0.149857,0.033962,0.000000,-0.033962,
                         -0.149857,-0.389703,-0.731771,-1.092705,
                         -1.355387,-1.412279,-1.207840,-0.763658,
                         -0.176336,0.411682,0.857949,1.065855,
                         1.013795,0.757295,0.403866,0.070597,
                         -0.159191,-0.263809,-0.285317,-0.305699,
                         -0.406943,-0.631123,-0.956572,-1.300000,
                         -1.544357,-1.582180,-1.358000,-0.893485,
                         -0.285317,0.323977,0.791865,1.021653,
                         0.991651,0.757295,0.426010,0.114798,
                         -0.093107,-0.176103,-0.176336,-0.175872,
                         -0.256783,-0.461223,-0.767601,-1.092705,
                         -1.319556,-1.340760,-1.100913,-0.621747,
                         0.000000,0.621747,1.100913,1.340760,
                         1.319556,1.092705,0.767601,0.461223,
                         0.256783,0.175872,0.176336,0.176103,
                         0.093107,-0.114798,-0.426010,-0.757295,
                         -0.991651,-1.021653,-0.791865,-0.323977,
                         0.285317,0.893485,1.358000,1.582180,
                         1.544357,1.300000,0.956572,0.631123,
                         0.406943,0.305699,0.285317,0.263809,
                         0.159191,-0.070597,-0.403866,-0.757295,
                         -1.013795,-1.065855,-0.857949,-0.411682,
                         0.176336,0.763658,1.207840,1.412279,
                         1.355387,1.092705,0.731771,0.389703,
                         0.149857,0.033962,-0.000000,-0.033962,
                         -0.149857,-0.389703,-0.731771,-1.092705,
                         -1.355387,-1.412279,-1.207840,-0.763658,
                         -0.176336,0.411682,0.857949,1.065855,
                         1.013795,0.757295,0.403866,0.070597,
                         -0.159191,-0.263809,-0.285317,-0.305699,
                         -0.406943,-0.631123,-0.956572,-1.300000,
                         -1.544357,-1.582180,-1.358000,-0.893485,
                         -0.285317,0.323977,0.791865,1.021653,
                         0.991651,0.757295,0.426010,0.114798,
                         -0.093107,-0.176103,-0.176336,-0.175872,
                         -0.256783,-0.461223,-0.767601,-1.092705,
                         -1.319556,-1.340760,-1.100913,-0.621747,
                         0.000000,0.621747,1.100913,1.340760,
                         1.319556,1.092705,0.767601,0.461223,
                         0.256783,0.175872,0.176336,0.176103,
                         0.093107,-0.114798,-0.426010,-0.757295,
                         -0.991651,-1.021653,-0.791865,-0.323977,
                         0.285317,0.893485,1.358000,1.582180,
                         1.544357,1.300000,0.956572,0.631123,
                         0.406943,0.305699,0.285317,0.263809,
                         0.159191,-0.070597,-0.403866,-0.757295,
                         -1.013795,-1.065855,-0.857949,-0.411682,
                         0.176336,0.763658,1.207840,1.412279,
                         1.355387,1.092705,0.731771,0.389703,
                         0.149857,0.033962,0.000000,-0.033962,
                         -0.149857,-0.389703,-0.731771,-1.092705,
                         -1.355387,-1.412279,-1.207840,-0.763658,
                         -0.176336,0.411682,0.857949,1.065855,
                         1.013795,0.757295,0.403866,0.070597,
                         -0.159191,-0.263809,-0.285317,-0.305699,
                         -0.406943,-0.631123,-0.956572,-1.300000,
                         -1.544357,-1.582180,-1.358000,-0.893485,
                         -0.285317,0.323977,0.791865,1.021653,
                         0.991651,0.757295,0.426010,0.114798,
                         -0.093107,-0.176103,-0.176336,-0.175872,
                         -0.256783,-0.461223,-0.767601,-1.092705,
                         -1.319556,-1.340760,-1.100913,-0.621747,
                         0.000000,0.621747,1.100913,1.340760,
                         1.319556,1.092705,0.767601,0.461223,
                         0.256783,0.175872,0.176336,0.176103,

     };


// End of File

 四、借助MATLAB计算IIR滤波器系数

首先进入MATLAB的filterDesigner。

滤波器截止频率选择2000Hz，滤除高频噪声，指定阶数为50阶。

 将生成的滤波器系数导入到工作区

由于系数在SOS矩阵的行中排序为（b0,b1,b2,1,-a1,-a2），于是我们需要对其第五列与第六列进行取反操作。 

这点在MATLAB提供的文档以及TI的例程注释中有提出

MATLAB提供的文档：

TI给出的注释：

 

G矩阵的元素为每一个对应的二阶节的增益。

编写MATLAB程序，处理导出到工作区的滤波器系数，并保存到本地便于我们后续使用

%将通过filterDesigner生成的SOS、G矩阵处理，保存到本地的coeffs_B.txt、coeffs_A.txt和scaleFactors.txt

SOS_num=size(SOS,1);%获取SOS矩阵的行数
G_num=size(G,1);%获取G矩阵的行数
%对生成的系数A取反
SOS(:,5)=-SOS(:,5);
SOS(:,6)=-SOS(:,6);

%将滤波器设计器生成的滤波器系数打印到coeffs_B,coeffs_A.txt文件中
fileID = fopen('coeffs_B.txt','w');%写文件
for i=1:1:SOS_num
    fprintf(fileID,'%6f,%6f,%6f,\n',SOS(i,1),SOS(i,2),SOS(i,3));
end
fclose(fileID);%将生成的初始信号保存到.m文件路径下的coeffs_B.txt

fileID = fopen('coeffs_A.txt','w');%写文件
for i=1:1:SOS_num
    fprintf(fileID,'%6f,%6f,%6f,\n',SOS(i,4),SOS(i,5),SOS(i,6));
end
fclose(fileID);%将生成的初始信号保存到.m文件路径下的coeffs_A.txt

%将滤波器设计器生成的比例系数打印到scaleFactors.txt文件中
fileID = fopen('scaleFactors.txt','w');%写文件

TEMP=G(1:G_num-1,1);%取G前G_num-1个数
count=0;
for i=1:1:(G_num-1)
    fprintf(fileID,'%6f,',TEMP(i,1));
    count=count+1;
    if count==3
        fprintf(fileID,'\n');%每三个元素后就换行
        count=0;
    end
end

生成的滤波器系数及对应增益如下：

/*
 * Filter Coefficients (C Source) generated by the Filter Design and Analysis Tool
 * Generated by MATLAB(R) 9.11 and Signal Processing Toolbox 8.7.
 * Generated on: 24-Jul-2023 16:04:36
 */

/*
 * 离散时间 IIR 滤波器(实数)
 * ----------------
 * 滤波器结构  : 直接 II 型，二阶节
 * 节数     : 25
 * 稳定     : 是
 * 线性相位   : 否
 */

/* General type conversion for MATLAB generated C-code  */

#define SOS_NUM         (25U)    // Number of Second Order Stages (biquad)

const float scaleFactors[SOS_NUM]= {0.891251,1.000000,1.000000,
                                    1.000000,1.000000,1.000000,
                                    0.999999,0.999998,0.999996,
                                    0.999990,0.999976,0.999942,
                                    0.999861,0.999668,0.999202,
                                    0.998086,0.995413,0.989026,
                                    0.973869,0.938459,0.858783,
                                    0.694654,0.418945,0.129278,
                                    0.004035,
};
const float coeffs_B[SOS_NUM*3] = {1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984229,1.000000,
                                   1.000000,-1.984228,1.000000,
                                   1.000000,-1.984227,1.000000,
                                   1.000000,-1.984223,1.000000,
                                   1.000000,-1.984214,1.000000,
                                   1.000000,-1.984193,1.000000,
                                   1.000000,-1.984141,1.000000,
                                   1.000000,-1.984016,1.000000,
                                   1.000000,-1.983712,1.000000,
                                   1.000000,-1.982959,1.000000,
                                   1.000000,-1.981005,1.000000,
                                   1.000000,-1.975321,1.000000,
                                   1.000000,-1.952864,1.000000,
                                   1.000000,-1.685740,1.000000,
                                   1.000000,-1.984229,1.000000,
};
const float coeffs_A[SOS_NUM*3] = {1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-1.000000,
                                   1.000000,1.984229,-0.999999,
                                   1.000000,1.984228,-0.999999,
                                   1.000000,1.984227,-0.999997,
                                   1.000000,1.984224,-0.999992,
                                   1.000000,1.984216,-0.999980,
                                   1.000000,1.984197,-0.999953,
                                   1.000000,1.984150,-0.999886,
                                   1.000000,1.984040,-0.999727,
                                   1.000000,1.983777,-0.999348,
                                   1.000000,1.983154,-0.998450,
                                   1.000000,1.981703,-0.996359,
                                   1.000000,1.978466,-0.991696,
                                   1.000000,1.971931,-0.982279,
                                   1.000000,1.961429,-0.967146,
                                   1.000000,1.951510,-0.952853,
                                   1.000000,1.984229,-1.000000,
};

至此已完成了大部分准备工作，接下来开始进行IIR滤波操作。

五、IIR滤波的代码实现

仿照TI给出的例程，完成了以下文件的编写

IIR.h文件

/*
 * IIR.h
 *

 */


#ifndef APP_IIR_H_
#define APP_IIR_H_

#include "DSP28x_Project.h"
#include <stdio.h>
#include <string.h>


void IIR_filter_run(void);


#endif /* APP_IIR_H_ */

IIR.c文件

/*
 * FIR.c
 *
 *  Created on: 2023年7月28日
 *      Author: yang
 */


#include <IIR.h>
#include "dsp.h"
#include "fpu_filter.h"
#include <math.h>
#include "complex.h"
#include "fdacoefs.h"

#define pi 3.14159
#define IIR_SIZE       (512U)
#define NUM_SOS         (25U)    // Number of Second Order Stages (biquad)
#define FILTER_ORDER    (NUM_SOS<<1)


#pragma DATA_SECTION(coeffs_B, "IIR_buffer0");
#pragma DATA_SECTION(coeffs_A, "IIR_buffer0");
#pragma DATA_SECTION(IIR_output, "IIR_buffer1");

// IIR_f32 object
IIR_f32 iir;
// Handle to the IIR_f32 object
IIR_f32_Handle hnd_iir = &iir;


// The filter coefficients and scale factors are tacked on to the end of the
// MATLAB generated input
extern const float coeffs_B[];//系数数组大小在fdacoefs.h中定义，大小为滤波器节数*3
extern const float coeffs_A[];
extern const float scaleFactors[];

// The delay line buffer
float delayLine[NUM_SOS*4U];

float IIR_output[IIR_SIZE];
extern float test_input[];
//*****************************************************************************
// the function definitions
//*****************************************************************************
void IIR_filter_run(void)
{
    // Locals
    uint16_t i;
    float32u_t in, out;

    // Configure the object
    IIR_f32_setCoefficientsBPtr(hnd_iir, coeffs_B);
    IIR_f32_setCoefficientsAPtr(hnd_iir, coeffs_A);
    IIR_f32_setDelayLinePtr(hnd_iir, delayLine);
    IIR_f32_setInputPtr(hnd_iir, (float *)&in);
    IIR_f32_setOutputPtr(hnd_iir, (float *)&out);
    IIR_f32_setScalePtr(hnd_iir, scaleFactors);
    IIR_f32_setOrder(hnd_iir, FILTER_ORDER);
    IIR_f32_setInitFunction(hnd_iir, (v_pfn_v)IIR_f32_init);
    IIR_f32_setCalcFunction(hnd_iir, (v_pfn_v)IIR_f32_calc);

    // Run the initialization function
    hnd_iir->init(hnd_iir);

    for(i = 0U; i < IIR_SIZE; i++)
    {
        out.f32  = FLT_MAX;
        in.f32   = test_input[i];

        // Call the calculation routine
        hnd_iir->calc(hnd_iir);

        IIR_output[i] = out.f32;
     }

}


// End of File

由于笔者使用了自定义段，故要在cmd文件中分配其对于的存储空间。Cmd文件如下：

/*
// TI File $Revision: /main/3 $
// Checkin $Date: July 9, 2008   14:12:45 $
//###########################################################################
//
// FILE:    28335_RAM_lnk.cmd
//
// TITLE:   Linker Command File For IQmath examples that run out of RAM
//
// NOTE; The example project uses memory protected by the
//       Code Security Module (CSM).  Make sure the CSM is
//       unlocked before you load the project.  One quick way
//       to do this on an erased device is to open a memory
//       window to the CSM password locations.  If these locations
//       read back 0xFFFF (or non-zero), then the CSM is unlocked:
//
//      Device    Password locations
//      28335:    0x33FFF8 - 0x33FFFF
//
//###########################################################################
//
//
// $Copyright: Copyright (C) 2014-2022 Texas Instruments Incorporated -
//             http://www.ti.com/ ALL RIGHTS RESERVED $
//###########################################################################
*/




MEMORY
{
PAGE 0 :
   BEGIN      : origin = 0x000000, length = 0x000002     /* Boot to M0 will go here                      */
   BOOT_RSVD  : origin = 0x000002, length = 0x00004E     /* Part of M0, BOOT rom will use this for stack */

   RAML0      : origin = 0x008000, length = 0x000800
   RAML1L2    : origin = 0x008800, length = 0x004800
   ZONE7A     : origin = 0x200000, length = 0x00FC00    /* XINTF zone 7 - program space */

   FLASHH     : origin = 0x300000, length = 0x008000     /* on-chip FLASH */
   FLASHG     : origin = 0x308000, length = 0x008000     /* on-chip FLASH */
   FLASHF     : origin = 0x310000, length = 0x008000     /* on-chip FLASH */
   FLASHE     : origin = 0x318000, length = 0x008000     /* on-chip FLASH */
   FLASHD     : origin = 0x320000, length = 0x008000     /* on-chip FLASH */
   FLASHC     : origin = 0x328000, length = 0x008000     /* on-chip FLASH */
   FLASHA     : origin = 0x338000, length = 0x007F80     /* on-chip FLASH */


   CSM_RSVD   : origin = 0x33FF80, length = 0x000076     /* Part of FLASHA.  Program with all 0x0000 when CSM is in use. */
   CSM_PWL    : origin = 0x33FFF8, length = 0x000008     /* Part of FLASHA.  CSM password locations in FLASHA            */
   OTP         : origin = 0x380400, length = 0x000400     /* on-chip OTP */
   ADC_CAL    : origin = 0x380080, length = 0x000009

   IQTABLES    : origin = 0x3FE000, length = 0x000b50     /* IQ Math Tables in Boot ROM */
   IQTABLES2   : origin = 0x3FEB50, length = 0x00008c     /* IQ Math Tables in Boot ROM */
   FPUTABLES   : origin = 0x3FEBDC, length = 0x0006A0     /* FPU Tables in Boot ROM */
   ROM         : origin = 0x3FF27C, length = 0x000D44     /* Boot ROM */

   RESET      : origin = 0x3FFFC0, length = 0x000002
   VECTORS     : origin = 0x3FFFC2, length = 0x00003E     /* part of boot ROM  */







PAGE 1 :
   BOOT_RSVD  : origin = 0x000002, length = 0x00004E     /* Part of M0, BOOT rom will use this for stack */
   RAMM0      : origin = 0x000050, length = 0x0003B0     /* on-chip RAM block M0 */
   RAMM1      : origin = 0x000400, length = 0x000400     /* on-chip RAM block M1 */
   RAML3      : origin = 0x00D000, length = 0x001000
   RAML4      : origin = 0x00E000, length = 0x001000
   RAML5      : origin = 0x00F000, length = 0x001000
   ZONE7B     : origin = 0x20FC00, length = 0x000400     /* XINTF zone 7 - data space */
   FLASHB     : origin = 0x330000, length = 0x008000     /* on-chip FLASH */

}
 

SECTIONS
{
   FPUmathTables    : > FPUTABLES, PAGE = 0, TYPE = NOLOAD
   /* Setup for "boot to SARAM" mode:
      The codestart section (found in DSP28_CodeStartBranch.asm)
      re-directs execution to the start of user code.  */

   codestart        : > BEGIN,     PAGE = 0
 
#ifdef __TI_COMPILER_VERSION__
    #if __TI_COMPILER_VERSION__ >= 15009000
        .TI.ramfunc : {} LOAD = RAML0,
                         RUN = RAML0,
                         LOAD_START(_RamfuncsLoadStart),
                         LOAD_END(_RamfuncsLoadEnd),
                         RUN_START(_RamfuncsRunStart),
                         LOAD_SIZE(_RamfuncsLoadSize),
                         PAGE = 0
    #else
        ramfuncs       : LOAD = RAML0,
                         RUN = RAML0,
                         LOAD_START(_RamfuncsLoadStart),
                         LOAD_END(_RamfuncsLoadEnd),
                         RUN_START(_RamfuncsRunStart),
                         LOAD_SIZE(_RamfuncsLoadSize),
                         PAGE = 0
    #endif
#endif

	.text           : {*(.text)}>> RAML1L2|RAML0 PAGE = 0
	.cinit          : > RAML0,     PAGE = 0
    .pinit           : > RAML0,     PAGE = 0
    .switch          : > RAML0,     PAGE = 0

    .stack           : > RAMM1,     PAGE = 1
    .ebss            : > RAML3,     PAGE = 1
    .econst          : > RAML3,     PAGE = 1
    .sysmem          : > RAML3,     PAGE = 1
    .esysmem         : > RAML3,     PAGE = 1

    .sysmem          : > RAML3,     PAGE = 1

	IIR_buffer0		: > RAML4,     PAGE = 1
	IIR_buffer1		: > RAML5,     PAGE = 1


   .cio             : > RAML3,     PAGE = 1
   ZONE7DATA        : > ZONE7B,    PAGE = 1

   DMARAML4         : > RAML4,     PAGE = 1
   DMARAML5         : > RAML5,     PAGE = 1

   .reset           : > RESET,     PAGE = 0, TYPE = DSECT /* not used                    */
   csm_rsvd         : > CSM_RSVD   PAGE = 0, TYPE = DSECT /* not used for SARAM examples */
   csmpasswds       : > CSM_PWL    PAGE = 0, TYPE = DSECT /* not used for SARAM examples */

   /* Allocate ADC_cal function (pre-programmed by factory into TI reserved memory) */
   .adc_cal     : load = ADC_CAL,   PAGE = 0, TYPE = NOLOAD

}

/*
//===========================================================================
// End of file.
//===========================================================================
*/

在主函数中，我们调用IIR_filter_run();函数即可完成滤波操作。

借助CCS的Graph画图工具测试滤波效果如下：

原始信号：

滤波后的1KHz信号：

可以看到IIR数字滤波器滤除了5KHz与10KHz的信号，仅保留了1KHz的信号（图中信号1个完整周期的采样点数为100点，滤波后信号幅度约为0.3，与理论值相符），从而验证滤波结果的准确性。 

luoqicou
关注
  • 3
    点赞
  • 踩
  • 25
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
FFT在DSP28335中的使用
sysu_liuyang的博客
01-27 8665
引言： 在做谐波提取的过程中，用到了FFT，由于开关频率的限制，最开始考虑的是DFT以及滑窗DFT，之后又考虑过FFT。过程中走了许多弯路，记录之！ 1.傅里叶变换理论知识 1.1傅里叶级数 引用矩阵里面基的概念，一个非周期信号，可以用泰勒级数去展开，所选取的幂次信号就是一组基底，系数则是权重；于此对应，一个周期信号，就可以选取不同频率的正余弦信号作为基底，进行表示，从而得到傅里叶级数
基于DSP的iir数字滤波器设计
qq_37182918的博客
06-14 2261
IIR滤波器的设计可以利用模拟滤波器原型,借鉴成熟的模拟滤波器的设计结果进行双线性变换,将模拟滤波器变换成满足预定指标的数字滤波器,即根据模拟设计理论设计出满足要求的传递函数H(s),然后将H(s)变换成数字滤波器的传递函数H(z).IIR滤波器的设计可以利用模拟滤波器原型,借鉴成熟的模拟滤波器的设计结果进行双线性变换,将模拟滤波器变换成满足预定指标的数字滤波器,即根据模拟设计理论设计出满足要求的传递函数H(s),然后将H(s)变换成数字滤波器的传递函数H(z).
IIR滤波器的DSP实现
12-10
IIR滤波器的DSP实现
基于DSP的音频信号滤波技术研究
最新发布
QQ_2193276455的博客
04-27 770
*单片机设计介绍，基于DSP的音频信号滤波技术研究。
IIR滤波器的DSP实现 (2013年)
05-09
IIR滤波器用较少阶数就获得较高的频率选择特性。根据滤波器的技术指标，采用FDATool工具来设计阶数最小的滤波器，导出滤波器的抽头系数；然后用DSP汇编语言编程，实现IIR算法；用MATLAB产生合成输入信号，导入滤波系统；最后，用CCS软件观察滤波前后的波形变化。波形仿真结果和理论值相吻合。
几种软件滤波算法的比较
裸奔人生
01-04 1238
第1种方法限幅滤波法（又称程序判断滤波法）A方法根据经验判断，确定两次采样允许的最大偏差值（设为A）每次检测到新值时判断：如果本次值与上次值之差则本次值有效如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值B优点能有效克服因偶然因素引起的脉冲干扰C缺点无法抑制那种周期性的干扰平滑度差第2种方法中位值滤波法A方法连续采样N次（N取奇数）把N次采样值按大小排列取中间值为本次有效值
基于TMS320F28335的FIR滤波实现
05-29
【基于TMS320F28335的FIR滤波实现】 TMS320F28335是一款由Texas Instruments（TI）公司生产的高性能浮点数字信号处理器（DSP），常用于实时信号处理应用，如滤波、控制、通信等领域。在本项目中，我们将探讨如何...
手把手教你学DSP：基于TMS320F28335
10-16
手把手教你学DSP：基于TMS320F28335 手把手教你学DSP：基于TMS320F28335 手把手教你学DSP：基于TMS320F28335 手把手教你学DSP：基于TMS320F28335 手把手教你学DSP：基于TMS320F28335 手把手教你学DSP：基于TMS320F...
手把手教你学DSP—基于TMS320F28335的应用开发及实战1-5章课件.zip
01-07
《手把手教你学DSP—基于TMS320F28335的应用开发及实战》是一本深入浅出的教程，主要针对Texas Instruments（TI）的TMS320F28335数字信号处理器（DSP）进行讲解。这门课程涵盖的内容丰富，从基础到实践，为读者提供...
嵌入式DSP的原理与应用 基于TMS320F28335
04-24
为了深入理解和熟练使用TMS320F28335，可以参考TI提供的技术文档、应用笔记、在线课程以及“嵌入式DSP的原理与应用 基于TMS320F28335.pdf”这本书，这些资料将帮助开发者掌握从基本概念到复杂应用的完整知识体系。...
dsp 28335的实验报告 带滤波的AD采样
01-14
内面附带完整源码和注释，这是我们学校的28335的研究生实验课的报道，关于计算采样周期，采样信号的周期，幅值等
基于DSP的FIR滤波器设计与实现
11-11
用可编程DSP芯片实现数字滤波可通过修改滤波器的参数十分方便地改变滤波器的特性。因此，我们有必要对滤波器的设计方法进行研究，理解其工作原理，优化设计方法，设计开发稳定性好的滤波器系统。我们将通过DSP设计平台，实现较为重要的FIR和自适应滤波器系统。从而通过本课题的研究，掌握滤波器的设计技术，为通信、信号处理等领域实用化数字滤波器设计提供技术准备。
基于TMS320F28335的温控系统设计DSP.zip
04-21
“基于TMS320F28335的温控系统设计DSP课设报告.docx”和“基于TMS320F28335的温控系统设计DSP课设PPT.pptx”提供了项目的设计思路、实施过程、结果分析和结论，是学习者理解系统设计的详细指南。 综上所述，这个...
DSP28335实现FFT的工程实现方法
weixin_52940637的博客
05-13 1615
对于含有控制器的DSP来实现FFT，FFT的结果反映了控制的效果，FFT的结果不是控制的目标。因此，为了保证控制效果的可观测，在不使用谐波测量仪器的前提下，直接使用控制器的一些资源，计算出控制目标的谐波含量还是一个不错的方案。在FFT的工程应用的痛点就是它的蝶形运算的循环。对于采用FPGA+DSP的控制器架构实现FFT来看，利用FPGA的特性，使用FPGA实现高精度采样，DSP仅完成FFT的解算也是一个不错的方案，这将解决DSP控制中断与高速数据采样中断的冲突，使控制器资源最大化。
三电平有源电力滤波器 软硬件介绍 基于DSP28335
2301_78835236的博客
06-22 301
5. 数字信号处理（DSP）：这是一种处理数字信号的技术，涉及到数字信号的获取、处理、分析和控制等方面。6. 电力电子：这是研究电力系统中电力的转换和控制的学科领域。电力电子技术广泛应用于电力系统中的各种设备和系统，包括电力滤波器、变频器、逆变器等。三电平有源电力滤波器是一种用于减小电力系统中谐波和电压波动的电力滤波器。它采用有源电力电子器件来实现滤波功能，能够有效地改善电力质量。三电平有源电力滤波器主要由一个有源电力电子器件（如IGBT或MOSFET）、控制电路和滤波电路组成。
DSP-IIR滤波器设计
vizio.blog.csdn.net
11-27 2463
DSP-IIR滤波器设计学习笔记
DSP28335声音降噪(未完成)
weixin_30652897的博客
10-23 503
1. 确定使用的模块是Webrtc-NS，采集声音的芯片TLV32AIC23，实际测试发现Webrtc-NS无法使用，所以改成FIR滤波器。 从时域特性上来看，数字滤波器还可以分为有限冲激响应数字滤波器(FIR)和无限冲激响应数字滤波器(IIR )。对于有限冲激响应数字滤波器(FIR)其输出y(n)只取决于有限个过去和现在的输入x(n),x(n-1)... x(n-m);对于无限冲激响应数字滤波...
有源电力滤波器APF仿真模型，实验平台，基于dsp28335
m0_67349476的博客
02-25 490
有源电力滤波器APF仿真模型，实验平台，基于dsp28335 编号:1650664625138188花开花落永不凋零！！
基于DSP的IIR滤波器的设计
SwiftlyFly
10-28 3687
一、基本要求 1 设计FIR滤波器，其中 N=10、Fs = 8000Hz、Fc = 2500Hz; 2 定点实现，不允许使用float、double; 二、FIR滤波器基本知识介绍
手把手教你学dsp:基于tms320f28335 pdf
05-02
学习DSP，首先要对DSP的基本原理有基本的了解，如采样、变换、滤波等基础知识。其次，需要掌握DSP芯片的硬件架构和指令集，以及编写DSP程序所需的软件开发工具。 基于TMS320F28335 PDF的学习过程可以分为以下几个步骤： 一、学习DSP原理基础知识： 1、理解采样定理及信号的采样方法； 2、学习变换——FFT、DFT、傅里叶级数等； 3、学习滤波——数字滤波器的设计与实现等。 二、学习TMS320F28335 DSP芯片的硬件架构和指令集： 1、了解TMS320F28335的硬件架构及内部组成； 2、学习TMS320F28335的指令集和寄存器的使用方法。 三、学习DSP程序的编写： 1、学习并掌握Keil C51、Code Composer Studio等开发工具； 2、了解DSP程序设计的基本开发流程和方法； 3、学习并掌握DSP程序的调试和测试技术。 四、应用实例： 1、根据实际需求编写DSP程序； 2、将编写好的程序在硬件上进行测试和调试。 需要注意的是，学习DSP需要有较扎实的电子电路基础和一定的程序设计经验。在学习过程中，应尽量多实践自己的代码，参考书籍及互联网上的开源代码，并结合具体应用场景进行学习。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值