iir低通滤波c语言实现代码的简单介绍

本文目录一览：

• 1、一段matlab低通滤波器程序，求改编成C语言。
• 2、求IIR及FIR数字滤波器的C语言实现。（VC++）
• 3、如何用C语言实现低通滤波器
• 4、【原创】C++实现IIR二阶数字滤波器(一)
• 5、求脉冲响应不变法设计Elliptic型IIR数字低通滤波器

一段matlab低通滤波器程序，求改编成C语言。

这个我刚好做过一个滤波器，事实上对时域信号做FFT，截取一定点数再做逆FFT相当于理想滤波。设计滤波器代码如下：

f1=100;f2=200;%待滤波正弦信号频率

fs=2000;%采样频率

m=(0.3*f1)/(fs/2);%定义过度带宽

M=round(8/m);%定义窗函数的长度

N=M-1;%定义滤波器的阶数

b=fir1(N,f2/fs);%使用fir1函数设计滤波器

%输入的参数分别是滤波器的阶数和截止频率

figure(1)

[h,f]=freqz(b,1,512);%滤波器的幅频特性图

%[H,W]=freqz(B,A,N)当N是一个整数时函数返回N点的频率向量和幅频响应向量

plot(f*fs/(2*pi),20*log10(abs(h)))%参数分别是频率与幅值

xlabel(‘频率/赫兹’);ylabel(‘增益/分贝’);title(‘滤波器的增益响应’);

figure(2)

subplot(211)

t=0:1/fs:0.5;%定义时间范围和步长

s=sin(2*pi*f1*t)+sin(2*pi*f2*t);%滤波前信号

plot(t,s);%滤波前的信号图像

xlabel(‘时间/秒’);ylabel(‘幅度’);title(‘信号滤波前时域图’);

subplot(212)

Fs=fft(s,512);%将信号变换到频域

AFs=abs(Fs);%信号频域图的幅值

f=(0:255)*fs/512;%频率采样

plot(f,AFs(1:256));%滤波前的信号频域图

xlabel(‘频率/赫兹’);ylabel(‘幅度’);title(‘信号滤波前频域图’);

figure(3)

sf=filter(b,1,s);%使用filter函数对信号进行滤波

%参数分别为滤波器系统函数的分子和分母多项式系数向量和待滤波信号输入

subplot(211)

plot(t,sf)%滤波后的信号图像

xlabel(‘时间/秒’);ylabel(‘幅度’);title(‘信号滤波后时域图’);

axis([0.2 0.5 -2 2]);%限定图像坐标范围

subplot(212)

Fsf=fft(sf,512);%滤波后的信号频域图

AFsf=abs(Fsf);%信号频域图的幅值

f=(0:255)*fs/512;%频率采样

plot(f,AFsf(1:256))%滤波后的信号频域图

xlabel(‘频率/赫兹’);ylabel(‘幅度’);title(‘信号滤波后频域图’);

求IIR及FIR数字滤波器的C语言实现。（VC++）

这个问题比较复杂，最近本人也在研究数字滤波，

结合图片说一下

第一个图是fir的流程图，其中Z-1是延迟，是单个采样时间1/fs

n阶的fir滤波器就是选取最近的n+1个样本，然后使他们各自乘以自己的滤波器系数即图中的F（n），[一般其他书的表示是h（n）]

然后相加得到输出的y（n）就是一个输出点

，其中F（n）的得出需要根据采样频率和滤波器的通带和阻带来决定

其中为了改善旁瓣的幅值，一般在采样后给样本或者h（n）加窗，当然可以用“最佳方法”来做

得出h（n）大致方法是先将矩形窗进行DFT，得出h（n），然后对h（n）进行加窗得出h（k），然后将∑h（k）×x（n）=y（n），假如阶数较多可以用傅里叶变换使时域变频域后再将卷积相加，可以利用FFT来改进实时性，提升速度

上面就是fir滤波器的简述

第二个图片上传不了，直接给链接

;amp;z=0tn=baiduimagedetailword=%D2%BB%BD%D7iir%C2%CB%B2%A8%C6%F7in=12708cl=2cm=1sc=0lm=-1pn=0rn=1di=2607528304ln=1054fr=

图中的Z-1是延时，iir滤波器也叫无限冲击响应滤波器，是有反馈的，

图中的是一阶的，相对fir滤波器来说，iir滤波器可以用较低的阶数来获得较好的滤波特效。但是其相位特性较差。

鉴于实用性，还是建议楼主去图书馆借书看，百度不可能得到确实的方案，

楼主可以去借“数字信号处理”的书，国外的中译本就有详细介绍fir和iir以及fft还有其他变换，国内的dsp大都几乎是dsp用户手册的中译本，对上述问题都是很简陋地带过，不予置评。

本人推荐一本书在上面的dsp专栏有下载，40多M，叫DSP算法、应用和设计，本人有这本实体书，写的较好

如何用C语言实现低通滤波器

float middle_filter(float middle_value [] , intcount)

{

    float sample_value, data;

    int i, j;

    for (i=1; i for(j=count-1; j=i,–j){

        if(middle_value[j-1]=middle_value[j]{

            data=middle_value[j-1];

            middle_value[j-1]=middle_value[j]

            middle_value[j]=data;

        }

    }

    sample_value=middle_value(count-1)/2];

    return(sample_value);

}

【原创】C++实现IIR二阶数字滤波器(一)

如图Fig1所示，是IIR二阶数字滤波器的数学计算公式

转换到离散域，计算公式如下

滤波器主要有以下几种：高通/低通/带通。下分别说明此三种滤波器的系统的求取方法。

通常，对一个滤波器的要求，我们主要给出以下技术规格：中心频率frequency，采样频率sampleRate，增益dBgain，品质因数Q。

为计算方便，先定义以下几个值：

高通滤波器系数的计算：

采样原始信号为一个正弦信号和一个直流分量的相加。直流分量的幅值为1，正弦信号周期为1s，幅值也为1。

初始化采样信号

首先定义数字滤波器的结构体

初始化滤波器，这里我们选择二阶带通滤波器做测试

其中，Init_Filter函数定义如下

DSP滤波器函数

使用滤波器函数滤波原始采样信号

m_signal采样信号如下，由于包含直流分量1，所以它的幅值由正弦信号的（-1,1）沿Y轴向上偏移1位变成（0,2）

滤波后的信号m_signal_filter如下所示

通过上述实验验证了我们给出的二阶带通滤波器算法的有效性，但是我们滤波后的信号其实并不是真正的正弦信号，而是一个无限逼近正弦的信号，这是由于以下几个原因造成的

1）采样时间：采样速率越高，信号失真越小

2）计算精度：计算机处理信号时在存储浮点数会出现一定误差

3）算法精度：算法本身就有精度误差

改进：

1）读者们可以自行提高采样速率，去验证采样速率对信号误差的影响，评判标准可以使用均方根误差

2）读者们还可以通过更改中心频率，品质因子，选择不同的滤波器类型测试对信号的影响

求脉冲响应不变法设计Elliptic型IIR数字低通滤波器

1-1一.试用MATLAB设计一巴特沃斯低通数字滤波器，要求通带截至频率Wp=30HZ，阻带截至频率为Ws=35HZ，通带衰减不大于0.5DB，阻带衰减不小于40DB，抽样频Fs=100HZ。

代码为：

fp = 30;

fs = 35;

Fs = 100;

wp = 2*pi*fp/Fs;

ws = 2*pi*fs/Fs;

wp = tan(wp/2);

ws = tan(ws/2); % 通带最大衰减为0.5dB，阻带最小衰减为40dB

[N, wn] = buttord(wp, ws, 0.5, 40, ‘s’); % 模拟低通滤波器极零点

[z, p, k] = buttap(N); % 由极零点获得转移函数参数

[b, a] = zp2tf(z, p, k); % 由原型滤波器获得实际低通滤波器

[B, A] = lp2lp(b, a, wp);

[bz, az] = bilinear(B, A, .5);

还有三句。。。给分后给你