Butterworth数字滤波器设计

MATLAB语句: [b,a]=butter(n,wn,’ftype’)

设计一个阶数为n,截止频率为wn的数字低通滤波器。若参数ftype取‘high’或‘stop’,则用来设计相应的高通或带阻滤波器。返回值a和b为系统函数的分子和分母的系数。

例 : 设计一个阶数为9阶、截止频率为300Hz的高通Butterworth滤波器,采样频率为1000。并画出它的频率特性。

程序:       n=9;

             wn=300/500;

             [b,a]=butter(n,wn,’high’);

             freqz(b,a,128,1000);

作者:TianSL 

  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Butterworth高通滤波器是一种常用的数字滤波器,用于在数字信号上实现高频滤波。它的优点是在滤除高频信号的同时,保留了信号的幅度响应和相位响应。在C语言中,实现Butterworth高通滤波器可以使用数字滤波器设计工具(如MATLAB、Scilab等)生成所需的滤波器系数,然后将其编写成C语言程序进行实现。 假设我们已经得到了Butterworth高通滤波器的系数:b和a,它们表示滤波器的分子和分母多项式系数。然后,我们可以按照以下步骤在C语言中实现Butterworth高通滤波器: 1. 定义输入,输出和滤波器系数数组 ``` double input[N]; // 输入信号数组 double output[N]; // 输出信号数组 double b[M+1]; // 分子多项式系数数组 double a[M+1]; // 分母多项式系数数组 ``` 2. 初始化滤波器状态变量 ``` double state[M]; // 滤波器状态变量数组 for (int i = 0; i < M; i++) { state[i] = 0.0; } ``` 3. 定义滤波器函数 ``` void butter_highpass_filter(double *input, double *output, int n, double *b, double *a, int m, double *state) { for (int i = 0; i < n; i++) { double y = b[0] * input[i] + state[0]; for (int j = 1; j <= m; j++) { y += b[j] * input[i-j] - a[j] * output[i-j]; } for (int j = m-1; j > 0; j--) { state[j] = b[m-j] * input[i-m+j] - a[m-j] * output[i-m+j] + b[m-j+1] * input[i-m+j+1] - a[m-j+1] * output[i-m+j+1]; } state[0] = b[m] * input[i-m] - a[m] * output[i-m]; output[i] = y; } } ``` 4. 调用滤波器函数 ``` butter_highpass_filter(input, output, N, b, a, M, state); ``` 在调用滤波器函数前,需要确保输入信号数组input和滤波器系数数组b、a的长度和分母多项式系数a[0]的值不为0。调用滤波器函数后,输出信号数组output即为经过Butterworth高通滤波器后的信号。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值