IIR滤波器matlab源程序

最新推荐文章于 2024-07-27 13:30:21 发布
最新推荐文章于 2024-07-27 13:30:21 发布
阅读量540 收藏 2
点赞数
文章标签： matlab

（1）IIR一阶低通滤波器

clear;
fi=1;fs=10;Gc2=0.9;
wc=2*pi*fi/fs;     
omegac=tan(wc/2);
alpha=(sqrt(Gc2)/sqrt(1-Gc2))*omegac;
a=(1-alpha)/(1+alpha);
b=(1-a)/2;
w=0:pi/300:pi;     
Hw2=alpha^2./(alpha^2+(tan(w/2)).^2);
plot(w/pi,Hw2);
grid;
hold on;
（2）一阶高通滤波器 
clear;
fi=1;fs=10;Gc2=0.5;
wc=2*pi*fi/fs;
omegac=tan(wc/2);
alpha=(sqrt(1-Gc2)/(sqrt(Gc2)))*omegac;
a=(1-alpha)/(1+alpha);
b=(1+a)/2;
w=0:pi/300:pi;
Hw2=(tan(w/2).^2)./(alpha^2+(tan(w/2)).^2);
plot(w/pi,Hw2);
grid;
hold on;
（3）Notch 嵌波滤波器 

clear;
Gb2=0.5;
w0=0.35*pi;
deltaw=0.1*pi;
b=1/(1+tan(deltaw/2)*(sqrt(1-Gb2)/sqrt(Gb2)));
B=[1 -2*cos(w0) 1].*b;
A=[1 -2*b*cos(w0) (2*b-1)];
w=0:pi/500:pi;
H=freqz(B,A,w);
plot(w/pi,abs(H));
grid;
（4）Peak 滤波器 

clear;
Ac=3;
Gb2=10^(-Ac/10);
w0=0.35*pi;
deltaw=0.1*pi;
b=1/(1+tan(deltaw/2)*(sqrt(Gb2)/sqrt(1-Gb2)));
B=[1 0 -1].*(1-b);
A=[1 -2*b*cos(w0) (2*b-1)];
w=0:pi/500:pi;
H=freqz(B,A,w);
plot(w/pi,abs(H));
grid;
（5）IIR低通滤波（Butterworth） 

% IIR Lowpass Use Butterworth
% copyright by Etual
clear;
fs=20;fpass=4;fstop=5;
Ap=0.5;As=10;
wp=2*pi*fpass/fs;ws=2*pi*fstop/fs;
omegap=tan(wp/2);omegas=tan(ws/2);
ep=sqrt(10^(Ap/10)-1);
es=sqrt(10^(As/10)-1);
N=ceil(log(es/ep)/log(omegas/omegap));
omega0=omegap/ep^(1/N);
K=floor(N/2);
for i=1:K
    theta(i)=pi*(N-1+2*i)/(2*N);
end
for i=1:K
    G(i)=omega0^2/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a1(i)=2*(omega0^2-1)/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a2(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2);
end
if K<(N/2)
    G0=omega0/(omega0+1);a0=(omega0-1)/(omega0+1);
end
w=0:pi/300:pi;
Hw2=1./(1+(tan(w/2)/omega0).^(2*N));
plot(w/pi,Hw2);
grid;
（6）IIR高通滤波（Butterworth） 

% IIR Hightpass Use Butterworth
% copyright by Etual
clear;
fs=20;fpass=5;fstop=4;
Ap=0.5;As=10;
wp=2*pi*fpass/fs;ws=2*pi*fstop/fs;
omegap=cot(wp/2);omegas=cot(ws/2);
ep=sqrt(10^(Ap/10)-1);
es=sqrt(10^(As/10)-1);
N=ceil(log(es/ep)/log(omegas/omegap));
omega0=omegap/ep^(1/N);
K=floor(N/2);
for i=1:K
    theta(i)=pi*(N-1+2*i)/(2*N);
end
for i=1:K
    G(i)=omega0^2/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a1(i)=-2*(omega0^2-1)/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a2(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2);
end
if K<(N/2)
    G0=omega0/(omega0+1);a0=-(omega0-1)/(omega0+1);
end
w=(0+eps):pi/300:pi;
Hw2=1./(1+(cot(w/2)/omega0).^(2*N));
plot(w/pi,Hw2);
grid;
（7）IIR带通滤波（Butterworth） 

% IIR Bandpass Use Butterworth
% copyright by Etual
clear;
fs=20;fpa=2;fpb=4;fsa=1.5;fsb=4.5;
Ap=0.0877;As=16.9897;
wpa=2*pi*fpa/fs;wpb=2*pi*fpb/fs;wsa=2*pi*fsa/fs;wsb=2*pi*fsb/fs;
c=sin(wpa+wpb)/(sin(wpa)+sin(wpb));
omegap=abs((c-cos(wpb))/sin(wpb));
omegasa=(c-cos(wsa))/sin(wsa);omegasb=(c-cos(wsb))/sin(wsb);
omegas=min(abs(omegasa),abs(omegasb));
ep=sqrt(10^(Ap/10)-1);es=sqrt(10^(As/10)-1);
N=ceil(log(es/ep)/log(omegas/omegap));
omega0=omegap/ep^(1/N);
K=floor(N/2);
for i=1:K
    theta(i)=pi*(N-1+2*i)/(2*N);
end
for i=1:K
    G(i)=omega0^2/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a1(i)=4*c*(omega0*cos(theta(i))-1)/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a2(i)=2*(2*c^2+1-omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a3(i)=-(4*c*(omega0*cos(theta(i))+1))/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a4(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2);
end
if K<(N/2)
    G0=omega0/(1+omega0);a0(1)=-2*c/(1+omega0);a0(2)=(1-omega0)/(1+omega0);
end
w=(0+eps):pi/300:pi;
Hw2=1./(1+((c-cos(w))./(omega0*sin(w))).^(2*N));
plot(w/pi,Hw2);
grid;
（8）IIR带阻滤波（Butterworth） 

% IIR Bandstop Use Butterworth
% copyright by Etual
clear;
fs=20;fpa=1.5;fpb=4.5;fsa=2;fsb=4;
Ap=0.5;As=10;
wpa=2*pi*fpa/fs;wpb=2*pi*fpb/fs;wsa=2*pi*fsa/fs;wsb=2*pi*fsb/fs;
c=sin(wpa+wpb)/(sin(wpa)+sin(wpb));
omegap=abs(sin(wpb)/(c-cos(wpb)));
omegasa=sin(wsa)/(cos(wsa)-c);omegasb=sin(wsb)/(cos(wsb)-c);
omegas=min(abs(omegasa),abs(omegasb));
ep=sqrt(10^(Ap/10)-1);es=sqrt(10^(As/10)-1);
N=ceil(log(es/ep)/log(omegas/omegap));
omega0=omegap/ep^(1/N);
K=floor(N/2);
theta=zeros(1,K);
for i=1:K
    theta(i)=pi*(N-1+2*i)/(2*N);
end
G=zeros(1,K);a1=zeros(1,K);a2=zeros(1,K);
for i=1:K
    G(i)=omega0^2/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a1(i)=2*(omega0^2-1)/(1-2*omega0*cos(theta(i))+omega0^2);
end
for i=1:K
    a2(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2);
end
if K<(N/2)
    G0=omega0/(omega0+1);a0=(omega0-1)/(omega0+1);
end
w=(0+eps):pi/300:pi;
Hw2=1./(1+(sin(w)./(omega0*(c-cos(w)))).^(2*N));
plot(w/pi,Hw2);
grid;
（9）IIR低通滤波（chebyshev 1） 

% IIR Lowpass Use Chebyshev Type 1
% copyright by Etual
clear;
fs=20;fpass=4;fstop=5;
Ap=0.5;As=10;
wp=2*pi*fpass/fs;ws=2*pi*fstop/fs;
omegap=tan(wp/2);omegas=tan(ws/2);
ep=sqrt(10^(Ap/10)-1);
es=sqrt(10^(As/10)-1);
e=es/ep;w=omegas/omegap;
N=ceil(log(e+sqrt(e^2-1))/log(w+sqrt(w^2-1)));
a=log(1/ep+sqrt(1/ep^2+1))/N;
omega0=omegap*sinh(a);
K=floor(N/2);
theta=zeros(1,K);omega=zeros(1,K);
for i=1:K
    theta(i)=pi*(N-1+2*i)/(2*N);
end
for i=1:K
    omega(i)=omegap*sin(theta(i));
end
G=zeros(1,K);a1=zeros(1,K);a2=zeros(1,K);
for i=1:K
    G(i)=(omega0^2+omega(i)^2)/(1-2*omega0*cos(theta(i))+omega0^2+omega(i)^2);
end
for i=1:K
    a1(i)=2*(omega0^2+omega(i)^2-1)/(1-2*omega0*cos(theta(i))+omega0^2+omega(i)^2);
end
for i=1:K
    a2(i)=(1+2*omega0*cos(theta(i))+omega0^2+omega(i)^2)/(1-2*omega0*cos(theta(i))+omega0^2+omega(i)^2);
end
if K<(N/2)
    G0=omega0/(omega0+1);a0=(omega0-1)/(omega0+1);  
else
    H0=sqrt(1/(1+ep^2));
end
f=0:1/300:10;
Hf2=1./(1+ep^2*(cheby(N,tan(pi*f/fs)/omegap)).^2);
plot(f,abs(Hf2));
grid;
（9）IIR低通滤波（chebyshev 1） 

% IIR Lowpass Use Chebyshev Type 2
% copyright by Etual
clear;
fs=20;fpass=4;fstop=5;
Ap=0.5;As=10;
wp=2*pi*fpass/fs;ws=2*pi*fstop/fs;
omegap=tan(wp/2);omegas=tan(ws/2);
ep=sqrt(10^(Ap/10)-1);
es=sqrt(10^(As/10)-1);
e=es/ep;w=omegas/omegap;
N=ceil(log(e+sqrt(e^2-1))/log(w+sqrt(w^2-1)));
a=log(es+sqrt(es^2+1))/N;
omega0=omegas/sinh(a);
K=floor(N/2);
for i=1:K
    theta(i)=pi*(N-1+2*i)/(2*N);
end
for i=1:K
    omega(i)=omegas/sin(theta(i));
end
for i=1:K
    G(i)=(1+omega(i)^-2)/(1-2*omega0^-1*cos(theta(i))+omega0^-2+omega(i)^-2);
end
for i=1:K
    a1(i)=2*(1-omega0^-2+omega(i)^-2)/(1-2*omega0^-1*cos(theta(i))+omega0^-2+omega(i)^-2);
end
for i=1:K
    a2(i)=(1+2*omega0^-1*cos(theta(i))+omega0^-2+omega(i)^-2)/(1-2*omega0^-1*cos(theta(i))+omega0^-2+omega(i)^-2);
end
for i=1:K
    b1(i)=2*(1-omega(i))/(1+omega(i));
end
if K<(N/2)
    G0=omega0/(omega0+1);a0=(omega0-1)/(omega0+1);  
else
    H0=sqrt(1/(1+ep^2));
end
f=(0+eps):1/100:10;
Hf2=(cheby(N,omegas./tan(pi*f/fs))).^2./((cheby(N,omegas./tan(pi*f/fs))).^2+es^2);
plot(f,abs(Hf2));
grid;
（10）chebyshev 中用到的函数 cheby.m 
function CN=cheby(N,x)
if x<=1
    CN=cos(N*acos(x));
else
    CN=cosh(N*log(x+sqrt(x.^2-1)));
end









iteye_3619
关注
matlab iir高通滤波器,用matlab设计的IIR滤波器源程序-IIR一阶低通/高通
weixin_32825955的博客
03-17 921
用matlab设计的IIR滤波器源程序(1)IIR一阶低通滤波器 P576clear;fi=1;fs=10;Gc2=0.9;wc=2*pi*fi/fs;omegac=tan(wc/2);alpha=(sqrt(Gc2)/sqrt(1-Gc2))*omegac;a=(1-alpha)/(1+alpha);b=(1-a)/2;w=0:pi/300:pi;Hw2=alpha^2./(alpha^2+(t...
实验：使用IIR滤波器与非时不变卡尔曼滤波器的IIR卡尔曼滤波器效果对比
xsuwyXyg的博客
04-25 299
实验设计如下：首先，我们将通过数学推导和理论分析，给出IIR滤波器和非时不变卡尔曼滤波器的设计方法和表达式。接着，我们将构造不同类型的信号，包括包络、正弦、噪声等，作为实验输入，通过IIR滤波器和非时不变卡尔曼滤波器进行滤波处理。在本实验中，我们将使用二阶形式的陷波滤波器进行研究。在本实验中，我们将结合IIR滤波器，将非时不变卡尔曼滤波器应用于信号处理任务中，探究其在滤波性能和状态估计准确性方面的表现。综上所述，本实验通过使用IIR滤波器和非时不变卡尔曼滤波器，对滤波器的性能和应用进行了研究。
用MATLAB设计IIR数字滤波器
03-10
用MATLAB设计IIR数字滤波器，用MATLAB设计IIR数字滤波器，用MATLAB设计IIR数字滤波器
iir滤波器matlab程序,IIR滤波器matlab源程序
weixin_29825411的博客
03-20 494
(1)IIR一阶低通滤波器clear;fi=1;fs=10;Gc2=0.9;wc=2*pi*fi/fs;omegac=tan(wc/2);alpha=(sqrt(Gc2)/sqrt(1-Gc2))*omegac;a=(1-alpha)/(1+alpha);b=(1-a)/2;w=0:pi/300:pi;Hw2=alpha^2./(alpha^2+(tan(w/2)).^2);plot(w/pi,H...
MATLAB仿真：数字信号处理IIR数字滤波器设计
最新发布
lengmei1的博客
07-27 367
（2）要求将st中三路调幅信号分离，通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。（2）学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。本实验的数字滤波器的MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。
【matlab】iir滤波器
weixin_43778388的博客
04-04 1万+
文章目录三. 参考资料 三. 参考资料 IIR滤波器设计（调用MATLAB IIR函数来实现）
IIR数字滤波器设计及软件实现 Matlab
12-08
1．实验目的 （1）熟悉用双线性变换法设计IIR数字滤波器的原理与方法； （2）学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。 （3）掌握IIR数字滤波器的MATLAB实现方法。 （3）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。 2．实验原理 设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。基本设计过程是：①先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； ②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。第六章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2和椭圆模拟和数字滤波器。本实验要求读者调用如上函数直接设计IIR数字滤波器。 本实验的数字滤波器的MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。
IIR滤波器matlab源程序.doc
07-05
下面，我们将从给定的MATLAB源程序中生成相关的知识点。 一、IIR滤波器的基本概念 IIR滤波器是一种基于递归公式的信号滤波方法，它可以用来实现低通、高通、带通、带阻等各种频率响应特性的滤波器。IIR滤波器的...
IIR滤波器-巴特沃斯(Butterworth) matlab源程序
01-27
【Matlab】音频信号分析及IIR滤波处理——巴特沃斯(Butterworth) 文中的源程序
MATLAB设计的简单滤波器程序源代码.zip
11-14
这个名为"MATLAB设计的简单滤波器程序源代码.zip"的压缩包包含了一系列用于实现不同滤波器功能的源代码。这些源代码为我们提供了深入理解滤波器设计原理和MATLAB编程技巧的机会。 首先，我们要了解滤波器的基本概念...
IIR 数字低通滤波器(用MATLAB GUI设计.内附源码)
06-27
IIR 数字低通滤波器(用MATLAB GUI设计.内附源码) 使用MATLAB GUI 制作可以更改参数的IIR数字低通滤波器,已生成可独立执行的EXE格式文件.里面有工程文件的原码FIG.用户可以根据自已的要求更改界面,再生成EXE格式，希望对大家有所帮助。
基于MATLAB的IIR滤波器的设计与仿真
07-25
1. 学习双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，利用双线性变换法或脉冲响应不变法设计低通、高通和带通其中一种IIR数字滤波器。 2. 观察双线性变换或脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法或脉冲响应不变法的特点。 3. 了解Butterworth滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。 要求：设计巴特沃斯数字低通滤波器，要求通带边界频fp=2.5kHZ，通带最大衰减Rp=0.5dB；阻带边界频率fs=9kHZ，阻带最小衰减Rs=25dB,采样频率为Fs=30kHZ。
IIR高通滤波器（内含完整的MATLAB代码）
04-11
IIR高通滤波器（内含完整的MATLAB代码）
IIR滤波器 ccs程序，C语言和汇编
05-13
IIR滤波器，结合硬件，可完美实现滤波。如果无相应硬件，可修改选项，改为软件导入。
IIR滤波器设计基础及Matlab设计示例
热门推荐
chenxy_bwave的专栏
06-21 1万+
目录1. 概要2. IIR 和 FIR 滤波器的定义3. IIR滤波器和FIR滤波器的优缺点对比4. 双二次（Biquad）IIR 滤波器5. IIR滤波器设计指标（specifications）5.1 滤波器的频率响应特性5.2 模拟滤波器类型5.3 IIR滤波器结构5.3.1 Direct Form IIR structure5.3.2 Direct Form I IIR structure5.3.3 Direct Form II IIR structure6. 基于Matlab的SOS IIR滤波器设
【毕业设计】78-基于MATLAB的IIR数字滤波器设计与仿真（详细说明书+PPT+原理图+matlab仿真+源代码）
cqtianxingkeji的博客
03-22 1445
数字滤波器在数字信号处理系统中，通过使用离散系统的特性对系统的信号进行处理，从而改变信号的波形。通过滤波的方式使需要通过的信号通过，杂波或者不需要通过的信息抑制。数字滤波器在应用过程中通常使用做筛选的作用。在实际应用中，可以用作工业零件的筛选、音频信号的过滤。通过去除工业零件的粉尘或者过滤信号的杂波从而获取到最优质的零件或者音频信号的目的。在本文研究了在环境MATLAB环境下通过对无限脉冲响应（IIR）数字滤波器进行研究与设计，将其进行应用，对数字滤波器进行模拟，数字转化。
Matlab实现IIR数字滤波器设计
weixin_30662109的博客
12-14 5912
低通滤波器的设计： 1、基本指标 　　低通滤波器顾名思义就是低频率成分通过，高频率成分截止，那么在设计一个低通滤波器时首先要明白想要截止多大的频率。比如想截掉4Hz以上的信号，理想状态下就是将4Hz以上的信号成分全部截止，4Hz以下的信号全部保留。然而事实上几乎不存在这样的滤波器，通常情况下总是在通过频率和截止频率之间存在一个过渡带。通过频率这部分称为通带，允许通过的最大频率为通带截止频率ωp...
9阶iir数字滤波器汇编源程序的编写
06-01
很抱歉，我是一个语言模型，无法直接提供汇编代码。不过，我可以给你一些参考资料和思路： 1. 首先，你需要确定使用哪种架构和指令集。比如，x86架构上的汇编指令集有很多种，包括8086/8088、80286、80386、80486等等，你需要选择其中一种合适的指令集。 2. 其次，你需要了解数字滤波器的原理和设计方法。你可以参考一些经典的滤波器设计书籍，比如《数字信号处理（第四版）》、《数字信号处理与应用》等等。 3. 然后，你需要根据你选择的滤波器类型和参数，计算出其系数。这部分内容比较复杂，需要一定的数学基础。你可以使用MATLAB等工具来辅助计算。 4. 接着，你需要将计算出来的系数转换成汇编代码。这一部分需要一定的汇编语言基础，你需要熟悉汇编指令的格式和用法。 5. 最后，你需要在程序中实现滤波器的功能。这一部分需要将滤波器系数载入滤波器寄存器中，并对输入信号进行滤波处理。你需要注意滤波器的状态保存和恢复等问题。 总之，编写一个高质量的9阶IIR数字滤波器汇编源程序需要具备一定的数学基础、汇编语言基础和数字信号处理基础。如果你对上述知识点不熟悉，建议先学习相关知识再尝试编写。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值