用matlab设计的IIR滤波器源程序-IIR一阶低通/高通

用matlab设计的IIR滤波器源程序

（1）IIR一阶低通滤波器 P576

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）一阶高通滤波器 P581

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

200多MATLAB经典教程及论文: http://www.elecfans.com/soft/special/matlab/

阅读全文

matlab(228270) matlab(228270)
滤波器(174522) 滤波器(174522)



请按住滑块，拖动到最右边



了解新功能

基于FPGA和MATLAB实现IIR数字滤波器的设计和仿真验证分析

IIR数字滤波器在很多领域中都有着广阔的应用。与FIR数字滤波器相比，IIR数字滤波器可以用较低的阶数获得较高的选择性，而且所用存储单元少，经济效率高。一个N阶IIR数字滤波器的系统函数为：

2020-07-23 17:05:08

2191

IIR滤波器种类和设计

IIR滤波器（Infinite Impulse Response Filter）是一种数字滤波器，具有无限冲激响应特性。

2023-06-20 11:31:36

5822

Matlab辅助设计无限冲激响应滤波器（IIR）算法的方法

了解IIR滤波器的特点，掌握Matlab辅助设计滤波器系数的方法，并实现IIR滤波器滤除高频信号，并在LCD上显示结果。

2023-11-16 17:34:50

683

全文详解IIR滤波器原理与设计方法

IIR滤波器具有很高的滤波效率，在相同幅频响应条件下，所需的滤波器阶数明显比FIR滤波器低。其次，IIR滤波器的设计，可以借鉴模拟滤波器的设计成果。

2024-01-22 09:25:09

526

4-6 无限冲激响应滤波器（IIR）算法（MATLAB辅助，LCD显示）

4-6 无限冲激响应滤波器（IIR）算法（MATLAB辅助，LCD显示）实验目的本次教程是基于创龙教仪DSP教学实验箱：TL6748-PlusTEB完成的本节视频的目的是了解IIR滤波器

2023-06-08 12:04:28

IIR滤波器和FIR滤波器的对比分析介绍

，它具有线性相位、容易设计的优点。这也就说明，IIR滤波器具有相位不线性，不容易设计的缺点。而另一方面，IIR却拥有FIR所不具有的缺点，那就是设计同样参数的滤波器，FIR比IIR需要更多的参数。这也就说明，要增加DSP的计算量。DSP需要更多的计算时间，对DSP的实时性有影响。以下都是低通滤波器的设计。

2019-06-26 06:15:35

IIR滤波器的连接

IIR滤波器的怎么连

2012-05-10 15:53:51

IIR数字滤波器的Matlab和FPGA实现

本帖最后由 eehome 于 2013-1-5 10:01 编辑 IIR数字滤波器的Matlab和FPGA实现

2012-08-20 22:16:49

Matlab设计滤波器

介绍如何利用Matlab设计滤波器，IIR，FIR等等，希望对大家有用

2016-02-25 13:02:12

iir滤波器学习

1.简介在做音频eq的时候，iir滤波器是必不可少的，所以这节主要学习一下iir滤波器，当然，不会偏于理论，而是着重怎么用。iir滤波器主要是用来处理信号的频率的，这是第一个最基本的认识。摆一

2021-08-17 07:08:28

iir滤波器相位校正

iir滤波器后面级联一个全通滤波器进行相位校正，不知道怎样设计一个全通滤波器，以及不会用函数把两个滤波器级联起来，望大神指导、、、、{:4:}

2013-06-09 15:20:23

iir滤波器设计

iir滤波器设计联合matlab仿真

2014-07-10 17:23:30

正在加载...

林超文PCB设计：PADS教程，PADS视频教程	郑振宇老师：Altium Designer教程，Altium Designer视频教程
张飞实战电子视频教程	朱有鹏老师：海思HI3518e教程，HI3518e视频教程
李增老师：信号完整性教程，高速电路仿真教程	华为鸿蒙系统教程，HarmonyOS视频教程
赛盛：EMC设计教程，EMC视频教程	杜洋老师：STM32教程，STM32视频教程
唐佐林：c语言基础教程，c语言基础视频教程	张飞：BUCK电源教程，BUCK电源视频教程
正点原子：FPGA教程，FPGA视频教程	韦东山老师：嵌入式教程，嵌入式视频教程
张先凤老师：C语言基础视频教程	许孝刚老师：Modbus通讯视频教程
王振涛老师：NB-IoT开发视频教程	Mill老师：FPGA教程，Zynq视频教程
C语言视频教程	RK3566芯片资料合集
朱有鹏老师：U-Boot源码分析视频教程	开源硬件专题

搜索历史

用matlab设计的IIR滤波器源程序-IIR一阶低通/高通

评论