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

最新推荐文章于 2024-07-14 10:01:20 发布
最新推荐文章于 2024-07-14 10:01:20 发布
阅读量1.1k 收藏
点赞数
文章标签： matlab实现一阶低通滤波器

用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

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

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

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

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

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

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:///soft/special/matlab/

NetNinja
关注
一阶低通滤波器 matlab,(1条消息)RC一阶低通滤波器
weixin_35642839的博客
03-16 1493
由于工作原因一直在接触滤波器，但是对其详细概念和原理一知半解，所以一直想好好学习一下。最近终于抽出时间搜遍各种资料，但总觉得知识点太散，不太方便系统的理解，所以自己整理了一下，希望有所帮助。第一次编辑发博客，实在不易，word编辑好了，想直接拷贝上来还不行，公式无法复制，只好一个个截图，实在尴尬！是不是我没掌握技巧。文章参照和汇集多位大神之作，如有侵权还请给予谅解，谢谢！@长弓的坚持@Jason...
matlab iir高通滤波器,用matlab设计的IIR滤波器源程序-IIR一阶低通/高通
weixin_32825955的博客
03-17 932
用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...
matlab一阶及二阶低通、高通、陷波滤波器等滤波实例
罗伯特祥的博客
03-28 956
【代码】matlab一阶及二阶低通、高通、陷波滤波器等滤波实例。
数字 一阶低通滤波器 matlab 幅频响应代码
12-15
数字 一阶低通滤波器 matlab 幅频响应代码，详情见博客：数字 一阶低通滤波器 详细分析 冰三点水
一阶数字低通滤波器设计matlab
qczhao_10的博客
12-12 2638
一阶数字低通滤波器设计matlab
基于MATLAB的Simulink搭建一阶低通滤波器
m0_73931287的博客
10-25 4616
基于MATLAB的Simulink搭建一阶低通滤波器
基于Matlab中firpm函数设计FIR滤波器实现一阶RC低通滤波器
05-04
本m文件基于matlab使用firpm函数计算出滤波器系数，从而完成了使用FIR滤波器拟合一阶RC电路的功能。这种思想可以扩展到其他的传输函数上。
UANNHXT.rar_matlab 一阶低通_滤波器_通讯编程
07-15
设计低通滤波器的matlab源程序 IIR一阶低通滤波器,高通滤波器,Notch 嵌波滤波器
IIR.rar_IIR数字滤波器_iir滤波_巴特沃斯低通_带通滤波器_高通滤波 c
09-24
在压缩包中的"IIR"文件可能包含了这些实现的源代码、头文件或者示例程序，通过阅读和分析这些代码，你可以深入理解IIR滤波器的理论和实践应用。对于初学者，这是一个很好的学习资源，不仅可以了解滤波器的设计，还...
matlab中低通滤波器,设计数字低通滤波器(用matlab实现)
weixin_33462804的博客
03-16 3901
DSP设计滤波器报告姓名：张胜男 班级：07级电信(1)班 学号：一低通滤波器的设计(一)实验目的：掌握IIR数字低通滤波器的设计方法。(二)实验原理：1、滤波器的分类滤波器分两大类：经典滤波器和现代滤波器。经典滤波器是假定输入信号中的有用成分和希望取出的成分各自占有不同的频带。这样，当通过一个线性系统(即滤波器)后可讲欲去除的成分有效的去除。现代滤波器理论研究的主要内容是从含有噪声的数据记录(又...
低通滤波器 MATLAB程序
11-17
低通滤波器的MATLAB程序 ，值得借鉴。附带快速傅里叶变换的子程序
MATLAB设计的简单滤波器程序源代码-综合文档
05-21
这个压缩包包含的源代码可能提供了一种实现滤波器设计的方法。接下来，我们将详细探讨MATLAB中滤波器设计的基本概念、步骤以及可能涉及的源代码结构。 1. **滤波器类型**： - **低通滤波器**：允许低频信号通过，...
【学习总结】一阶低通滤波与matlab仿真
林梦不远的博客
04-25 1271
1、所谓低通滤波器就是允许低频信号通过，不允许高频信号通过。滤波是将信号中特定波段频率滤除的操作，一阶低通滤波也就是一阶惯性滤波，是使用软件实现普通硬件RC低通滤波的功能。2、一阶低通滤波法采用本次采样值与上次滤波输出值加权，得到有效滤波值，使输出对输入有反馈作用。，其中k为滤波系数，在（0, 1）上，X(n)为本拍采样值，Y(n - 1)是上一拍滤波输出值，Y(n)是本次滤波输出值。
一阶低通滤波(IIR)
最新发布
weixin_45052358的博客
07-14 629
递归性质：由于 Y(n) 依赖于 Y(n−1)，这个滤波器具有递归性质，即当前输出不仅与当前输入有关，还与过去的输出有关。当 a 接近 0 时，滤波器的截止频率较低，更多的高频成分被滤除；当 a 接近 1 时，截止频率较高，允许更多的高频成分通过。平滑作用：公式中的 (1−a)⋅Y(n−1) 部分表示上一次输出对当前输出的贡献。当 a 较小时，上一次输出的权重较大，这使得输出信号更加平滑，因为高频的突变会被抑制。这个公式中的 a 是一个介于 0 和 1 之间的系数，它决定了滤波器的截止频率和滤波效果。
matplotlib绘制一阶低通滤波器频率响应
2301_81614815的博客
12-08 1799
H(jω)=uoui=1jωCR+1jωC=11+jωRCH(j\omega )=\frac{u_o}{u_i}=\cfrac{\frac{1}{j\omega C}}{R+\frac{1}{j\omega C}}=\cfrac{1}{1+j\omega RC}H(jω)=ui​uo​​=R+jωC1​jωC1​​=1+jωRC1​通常截止频率时输出功率为传导频率的一半，即输出电压为输入电压的1/2\sqrt{1/2}1/2​ 当ω0=2πf0=1RC\omega_0=2\pi f_0=\frac{1}{R
Matlab 一阶低通滤波/一阶IIR
笔记
07-18 2936
LPF 一阶低通滤波 同样拿那个炸机气压计日志分析，可以更好对比 新建变量，把高度存进去BaroAlt.mat 写获取高度值函数 写LPF函数 function output = LPF(x) persistent lastx Init if isempty(Init) lastx = x; Init = 1; end alpha = 0.9; output = alpha*lastx + (1 - alpha)*x; lastx = output; 写个调用脚本RunLPF.m
一阶数字低通滤波器-C语言/matlab实现
yuyan7045的博客
10-18 1万+
1.一阶硬件滤波器 原理： 计算： 实验： 2.一阶数字滤波器 原理： 计算： 实验： 应用及注意事项： 3.一阶数字滤波器代码实现 代码： 实验数据：
用MATLAB设计低通滤波器
热门推荐
菜的不行的IT民工
05-28 3万+
滤波器有很多种设计方式。最简单的，就是用已有的buttord计算出最符合条件的巴特沃思滤波器的阶数n和截止频率，再用butter计算n阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量b、a。用freqz函数画出滤波器的幅频、相频图。用filter滤波。例：　　规定：wp为通带截止频率；ws为阻带截止频率；ap为通带最大衰减（dB）；as为阻带最大衰减(dB)；wc为3dB截止频率;fn为采样...
MATLAB实现低通滤波器(附完整代码)
corn1949的博客
02-08 9706
MATLAB实现低通滤波器(附完整代码)
MATLAB实现IIR数字滤波器的源码解析
本资源提供了IIR滤波器的MATLAB源程序，这些程序可以用来设计低通、高通、带通和带阻滤波器等多种类型的IIR滤波器。通过这些程序，用户能够直观地了解IIR滤波器的设计原理和实现方法。源代码的编写将涉及MATLAB的...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值