2021.09.10

最新推荐文章于 2022-11-20 21:16:41 发布

嘀嘀叮叮当当

最新推荐文章于 2022-11-20 21:16:41 发布

阅读量153

点赞数

文章标签：算法

本文链接：https://blog.csdn.net/jiangzecui/article/details/120226774

版权

Matlab设计巴特沃斯滤波器

命令行窗口输入fdatool
设计2阶、低通、IIR、巴特沃斯、采样频率100Hz、截止频率10Hz
file -> Export -> export
命令行输入[b,a]=sos2tf(Hd.sosMatrix,Hd.ScaleValues)
得到设计结果
- 与代码计算得到的基本一致


/******************************************************************************
 * @brief	计算二阶低通巴特沃斯滤波器参数
 * @param	sample_frequent 待处理数据的采样频率
 * @param cutoff_frequent 截止频率
 * @param *LPF 滤波器参数结构体
 * @return 
 *****************************************************************************/
void Set_Cutoff_Frequency(float sample_frequent, float cutoff_frequent,Butter_Parameter *LPF)
{
  float fr = sample_frequent / cutoff_frequent;
  float ohm = tanf(M_PI_F / fr);
  float c = 1.0f + 2.0f * cosf(M_PI_F / 4.0f) * ohm + ohm * ohm;
  if (cutoff_frequent <= 0.0f) {
    // no filtering
    return;
  }
  LPF->b[0] = ohm * ohm / c;
  LPF->b[1] = 2.0f * LPF->b[0];
  LPF->b[2] = LPF->b[0];
  LPF->a[0] = 1.0f;
  LPF->a[1] = 2.0f * (ohm * ohm - 1.0f) / c;
  LPF->a[2] = (1.0f - 2.0f * cosf(M_PI_F / 4.0f) * ohm + ohm * ohm) / c;
}