设为首页收藏本站
开启辅助访问 切换到宽版

 找回密码
 立即注册

QQ登录

只需一步，快速开始

热门： 51单片机 | 24小时必答区 | 单片机教程 | 单片机DIY制作 | STM32 | Cortex M3 | 模数电子 | 电子DIY制作 | 音响/功放 | 拆机乐园 | Arduino | 嵌入式OS | 程序设计

搜索
»论坛 › 嵌入式/单片机论坛 › DSP论坛 › DSP C语言的IIR滤波器程序源码
返回列表 发新帖
查看: 5884|回复: 8
收起左侧

DSP C语言的IIR滤波器程序源码

  [复制链接]
ID:209315 发表于 2018-8-8 16:32 | 显示全部楼层 |阅读模式
分享个DSP C语言的IIR滤波器

单片机源程序如下:
  1. #include <math.h>
  2. #include <stdio.h>
  3. #include "fdacoefs.h"

  5. //*********************************************
  6. //fs=1000000,fpass=10000,fstop=50000(-80db)
  7. //butterworth iir-df1
  8. //*********************************************
  9. #define fs 1000000
  10. #define f1 4000//测试正弦频率1
  11. #define f2 50000//测试正弦频率2
  12. #define f_rec 5000//测试方波频率

  14. #define Impulse 1
  15. #define Sine    2
  16. #define Rect    3

  18. #define pi 3.1415926

  20. #define STG_NUM (MWSPT_NSEC-1)/2
  21. #define DAT_NUM 800


  24. float IIR_out[DAT_NUM];


  27. float input_step[STG_NUM+1][3];// n; n-1; n-2

  29. void iir_filter_init(float *ptr,int num)
  30. {
  31.     int i;       
  32.         for(i=0;i<num;i++)
  33.         {
  34.                 *ptr++=0;
  35.         }
  36.        
  37. }

  39. float iir_filter(float in)
  40. {
  41.       
  42.       int i;
  43.           float temp;
  44.      
  45.       input_step[0][0] = in;
  46.       for(i=0;i<STG_NUM;i++)
  47.       {
  48.                 
  49.           temp =NUM[i*2][0]*(input_step[i][0]*NUM[i*2+1][0]+input_step[i][1]*NUM[i*2+1][1]+input_step[i][2]*NUM[i*2+1][2])-
  50.                                             input_step[i+1][1]*DEN[i*2+1][1]-input_step[i+1][2]*DEN[i*2+1][2];
  51.                                             
  52.           temp *= DEN[i*2+1][0];
  53.          
  54.           input_step[i][2]=input_step[i][1];             
  55.                 input_step[i][1]=input_step[i][0];                           
  56.           input_step[i+1][0]=temp;
  57.          
  58.       
  59.       }
  60.                                          
  61.       return  temp*NUM[MWSPT_NSEC-1][0];   
  62.    
  63. }

  65. float x;
  66. int main()
  67. {
  68.        
  69.         int i,Test_type;

  71.         printf("请选择测试类型：1、冲击测试 2、正弦测试 3、方波测试\n");
  72.         scanf("%d",&Test_type);
  73.                
  74.         while(1){
  75.                 iir_filter_init(input_step[0],sizeof(input_step)/sizeof(input_step[0][0]));
  76.                        
  77.                 switch(Test_type)
  78.                 {
  79.                         case Rect://方波输入
  80.                                 for(i=0;i<DAT_NUM;i++)
  81.                                 {                               
  82.                                         if(((2*f_rec*i/fs)%2)==0)
  83.                                                 x=5000;
  84.                                         else
  85.                                                 x=-5000;                                                               
  86.                                         IIR_out[i]=iir_filter(x);                                                                                                                                       
  87.                                 }               
  88.                         break;
  89.                        
  90.                         case Impulse://冲激输入
  91.                                 IIR_out[0]=iir_filter(5000);
  92.                                 for(i=1;i<DAT_NUM;i++)                                                               
  93.                                         IIR_out[i]=iir_filter(0);                                                                                                                                                       
  94.                         break;
  95.                        
  96.                         case Sine://正弦输入                       
  97.                                 for(i=0;i<DAT_NUM;i++)
  98.                                 {
  99.                                         //***************************************************************
  100.                                         x=500*(cos(2*pi*f1*i/fs)+cos(2*pi*f2*i/fs)); //正弦测试                                                                                       
  101.                                         IIR_out[i]=iir_filter(x);                                                                                                                                       
  102.                                 }                       
  103.                         break;       
  104.                 }
  105.         //请在printf前设置断点，以便停止运行 查看输出结果        否则，看不到输出
  106.         printf("请选择测试类型：1、冲击测试 2、正弦测试 3、方波测试\n");
  107.         scanf("%d",&Test_type);       
  108.                                
  109.         }
  110.        
  111.         while(1);

  113. }
复制代码

所有资料51hei提供下载:
IIR.7z (23.26 KB, 下载次数: 78)

评分

参与人数 1黑币 +50 收起 理由
admin + 50 共享资料的黑币奖励！

查看全部评分

回复

使用道具 举报

ID:95821 发表于 2018-8-10 11:55 | 显示全部楼层
學習了,正是在找的
回复

使用道具 举报

ID:209315 发表于 2018-8-10 16:14 | 显示全部楼层
alai318 发表于 2018-8-10 11:55
學習了,正是在找的

努力加油
回复

使用道具 举报

ID:95821 发表于 2018-8-12 12:36 | 显示全部楼层
剛學了FFT基2算法，FFT還沒學透，學好再來學濾波方面的
回复

使用道具 举报

ID:209315 发表于 2018-8-14 14:05 | 显示全部楼层
alai318 发表于 2018-8-12 12:36
剛學了FFT基2算法，FFT還沒學透，學好再來學濾波方面的

老哥厉害了
回复

使用道具 举报

ID:216962 发表于 2018-8-14 14:16 | 显示全部楼层
谢谢分享，论坛有您更精彩！
回复

使用道具 举报

ID:424128 发表于 2018-11-10 18:29 | 显示全部楼层
老哥，这个程序适合用于DSP C54x做IIR滤波器的实物吗，就是IIR数字滤波器的MATLAB仿真及DSP实现这种毕业设计题，能做出实物不？
回复

使用道具 举报

ID:392405 发表于 2019-7-27 10:49 | 显示全部楼层
感谢楼主分享
回复

使用道具 举报

ID:276779 发表于 2019-10-22 17:59 | 显示全部楼层
看一下正需要呢
回复

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

手机版|小黑屋|51黑电子论坛 |51黑电子论坛6群 QQ 管理员QQ:125739409;技术交流QQ群281945664

Powered by 单片机教程网

快速回复 返回顶部 返回列表