C语言 低通滤波器\带通滤波器\高通滤波器

[知识交流]当信号和干扰信号同时存在时,不会发生叠加吗?要如何滤掉干扰信号?

文章发表于:2007-10-21 18:16

这是对滤波的理解有误,滤波不是指滤掉干扰信号的幅值,而是滤掉它的频率。有用信号和干扰信号的频率不同,如果干扰信号的频率高,则用一个低通滤波器滤掉干扰信号,如果信号频率高,则用一个高通滤波器滤掉干扰信号。

MATLAB:

   y(m)=-(a(2)*y(m-1)+a(3)*y(m-2)+a(4)*y(m-3)+a(5)*y(m-4)+a(6)*y(m-5)+a(7)*y(m-6)+a(8)*y(m-7)+a(9)*y(m-8))+(b(1)*x(m)+b(2)*x(m-1)+b(3)*x(m-2)+b(4)*x(m-3)+b(5)*x(m-4)+b(6)*x(m-5)+b(7)*x(m-6)+b(8)*x(m-7)+b(9)*x(m-8));

C语言
以下一个结果:通带 5Hz,1db, 阻带8Hz,40db

b=[0.0117   -0.0329    0.0451   -0.0329    0.0117]
a=[1.0000   -3.6048    4.9787   -3.1156    0.7447]

int        filterBegin=5;
double xBuf[5];
double yBuf[5];
doub

  • 0
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
float DigFil(invar, setic) float invar; int setic; /******************************************************************************/ /* Filter Solutions Version 2009 Nuhertz Technologies, L.L.C. */ /* www.nuhertz.com */ /* +1 602-279-2448 */ /* 3rd Order Band Pass Butterworth */ /* Bilinear Transformation with Prewarping */ /* Sample Frequency = 5.000 KHz */ /* Standard Form */ /* Arithmetic Precision = 4 Digits */ /* */ /* Center Frequency = 300.0 Rad/Sec */ /* Pass Band Width = 20.00 Rad/Sec */ /* */ /******************************************************************************/ /* */ /* Input Variable Definitions: */ /* Inputs: */ /* invar float The input to the filter */ /* setic int 1 to initialize the filter to zero */ /* */ /* Option Selections: */ /* Standard C; Initializable; Internal States; Not Optimized; */ /* */ /* There is no requirement to ever initialize the filter. */ /* The default initialization is zero when the filter is first called */ /* */ /******************************************************************************/ /* */ /* This software is automatically generated by Filter Solutions */ /* no restrictions from Nuhertz Technologies, L.L.C. regarding the use and */ /* distributions of this software. */ /* */ /******************************************************************************/ { float sumnum=0.0, sumden=0.0; int i=0; static float states[6] = {0.0,0.0,0.0,0.0,0.0,0.0}; static float znum[7] = { -7.968e-09, 0.0, 2.39e-08, 0.0, -2.39e-08, 0.0, 7.968e-09 }; static float zden[6] = { .992, -5.949, 14.88, -19.86, 14.92, -5.981 }; if (setic==1){ for (i=0;i<6;i++) states[i] = [i] = [i]*invar; return 0.0; } else{ sumnum = sumden = 0.0; for (i=0;i<6;i++){ sumden += states[i]*zden[i]; sumnum += states[i]*znum[i]; if (i<5) states[i] = states[i+1]; } states[5] = invar-sumden; sumnum += states[5]*znum[6]; return sumnum; } }
数字滤波器的C语言实现通常分为三种类型:高通滤波器低通滤波器带通滤波器。其中,高通滤波器可以通过滤除低频信号来保留高频信号,低通滤波器可以滤除高频信号来保留低频信号,带通滤波器则是同时滤除高、低频信号,只保留一个中间的频率段。 实现数字滤波器可以采用一些经典的滤波算法,如FIR(有限脉冲响应滤波器)和IIR(无限脉冲响应滤波器)。以FIR滤波器为例,其C语言实现需要定义一组系数,可以通过差分方程或直接计算获得。一般而言,FIR滤波器的算法步骤包括采样、计算、存储等过程,其代码实现如下: ```c #define FILTER_LEN 51 // 滤波器长度 float filter[FILTER_LEN] = { 0.0 }; // 滤波器系数 float input_data = 0.0; // 输入数据 float output_data = 0.0; // 输出数据 void fir_filter(float data) { int i; float tmp = 0.0; // 滑动数组 for (i = (FILTER_LEN - 1); i >= 1; i--) { filter[i] = filter[i - 1]; } // 存储新数据 filter[0] = data; // 计算输出数据 for (i = 0; i < FILTER_LEN; i++) { tmp += filter[i] * filter[FILTER_LEN - 1 - i]; } output_data = tmp; } ``` 对于高通、低通和带通滤波器实现,可以根据其特点来定义不同的滤波算法和系数。例如,对于低通滤波器而言,滤波器系数应该反映出其去除高频信号的特性,可以通过对一组低通滤波器的系数进行加权和以获得不同的截止频率。而对于带通滤波器,则需要同时考虑其阻带和通带的特性,进而确定其滤波器系数和算法。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值