EDA365电子论坛网

标题: DSP28 fir低通滤波器的设计源码 [打印本页]

作者: kiygb    时间: 2021-9-28 13:45
标题: DSP28 fir低通滤波器的设计源码
#include"stdio.h"
#include"math.h"
#include"DSP28_Device.h"
#defineorder 26
#defineN  128         
#definepi 3.1415926
floatFIR_LPF[order] ={-1.52864422149721e-05,-0.000133359311426970,-0.000544636220343477,-0.00137835208978255,-0.00218948983490142,-0.00124340086554490,0.00478134865608317,0.0198838342801837,0.0464361812966399,0.0825102845162328,0.120763491355644,0.150360916315540,0.161536936687780,0.150360916315540,0.120763491355644,0.0825102845162328,0.0464361812966399,0.0198838342801837,0.00478134865608317,-0.00124340086554490,-0.00218948983490142,-0.00137835208978255,-0.000544636220343477,-0.000133359311426970,-1.52864422149721e-05};
floats;
inti,n;
floatinput[N],output[N];  
intfs=10000;   
intf1=600;   
intf2=1500;
#definefa 2*pi*f1/fs
#definefb 2*pi*f2/fs
void inputwave()
{
    for(n=0;n<N;n++ )
input[n]=sin(fa*n)+sin(fb*n);
}
voidmain()
{   int m=0,n=0;
    inputwave();           
for(n=0;n<N+order;n++)
{
  for(s=0,m=0;(m<order)&&(m<n);m++)
{
   s=FIR_LPF[m]*input[n-m]+s;
}
        output[n]=s;}
while(1);  
5 o$ g5 [- ~8 U' }, w

作者: jspij1    时间: 2021-9-28 14:20
学习了  感谢分享
作者: grand    时间: 2021-9-28 16:43
建议下次加个注释+ r$ u/ ]4 p; k0 }





欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2