標(biāo)題: MATLAB聯(lián)合ccs做數(shù)字濾波器程序設(shè)計(jì)，適合純小白 [打印本頁]
作者: DJN    時(shí)間: 2020-6-24 09:34
標(biāo)題: MATLAB聯(lián)合ccs做數(shù)字濾波器程序設(shè)計(jì)，適合純小白
之前學(xué)習(xí)的內(nèi)容，因?yàn)闆]有開發(fā)板，所以只能ccs聯(lián)合MATLAB一起做數(shù)字濾波器的設(shè)計(jì)。這個(gè)對(duì)于剛剛接觸的新手還是非常有難度的，像我就熬了一個(gè)通宵才做出來，當(dāng)然還是有點(diǎn)菜。詳細(xì)的程序都在安裝包里�？梢詤⒖�。而且ccs這個(gè)軟件有點(diǎn)迷，有時(shí)候沒有圖或者一些錯(cuò)誤，可能重啟一下軟件就好了。


源程序如下:
  1. /*
  2. * main.c
  3. */

  5. #include "math.h"
  6. #include "stdio.h"
  7. #include <stdint.h>
  8. #include "stdlib.h"
  9. #include "string.h"
  10. #include "limits.h"
  11. //#include"fdacoefsiir900.h"
  12. //#include "fdacoefsiir600.h"
  13. //#include "fdacoefsiir50+900.h"
  14. //#include "fdacoefsiir50.h"
  15. #include "fdacoefs50001.h"

  17. #define length 2048
  18. #define pi 3.1415926

  20. long fs=10000;
  21. int f1=50;
  22. int f2=200;
  23. int f3=600;
  24. int f4=900;

  26. #define w1 2*pi*f1/fs
  27. #define w2 2*pi*f2/fs
  28. #define w3 2*pi*f3/fs
  29. #define w4 2*pi*f4/fs

  31. double input[length];
  32. double output[length];

  34. static double xlast[2];
  35. static double mlast[2];

  37. static double IIR_DR2(double x,double *plast,const double (*A)[3],const double (*B)[3])
  38. {
  39.         double tmp,last;

  41.         tmp = x*B[0][0];

  43.         last = tmp - (A[1][1]*plast[0] + A[1][2]*plast[1]);
  44.         tmp = last + (B[1][1] * plast[0] + B[1][2]*plast[1]);

  46.         plast[1] = plast[0];
  47.         plast[0] = last;

  49.         return tmp;
  50. }
  51. double IIR_Filter(double x)
  52. {
  53.         double mid,y;

  55.         mid = IIR_DR2(x,xlast,DEN,NUM);

  57.         y   = IIR_DR2(mid,mlast,&DEN[2],&NUM[2]);//二階濾波器組合成更高階數(shù)的濾波器

  59.         //更多階數(shù)...

  61.         return y;
  62. }
  63. void Init_Filter(void)//初始化中間數(shù)值
  64. {
  65.         xlast[0] = 0;
  66.         xlast[1] = 0;
  67.         mlast[0] = 0;
  68.         mlast[1] = 0;
  69. }

  71. int main(void)
  72. {
  73. unsigned int i,n;

  75. Init_Filter();

  77. for(i=0;i<length;i++)
  78. {

  80.         input[i]=2048*sin(w1*i)+2048*sin(w2*i)+2048*sin(w3*i)+2048*sin(w4*i);
  81. }


  84. for(n=0;n<length;n++)
  85. {
  86.         output[n]=IIR_Filter(input[n]);


  89. }


  92. for(;;);
  93. }

復(fù)制代碼

所有資料51hei提供下載:
7.zip (101.16 KB, 下載次數(shù): 50)

作者: wzg-235    時(shí)間: 2023-9-19 10:35
你好加個(gè)好友交流一下可以嗎？



歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1