关于滤波器设计的matlab函数简表

合集下载

Matlab滤波信号处理函数

Matlab滤波信号处理函数2009-12-04 19:32:22| 分类：matlab方法| 标签：|字号大中小订阅1 conv功能：求卷积。

格式：c = conv(a,b)说明：c = conv(a,b)返回向量a、b的卷积c。

举例：a = [1 2 3]b = [4 5 6]c = conv(a,b)c=4 13 28 27 182 impz功能：数字滤波器的冲激响应。

格式：[h,t] = impz(b,a)[h,t] = impz(b,a,n)[h,t] = impz(b,a,n,Fs)impz(b,a)impz(...)说明：[h,t] = impz(b,a)返回系统(b,a)的冲激响应h和相应的时间轴向量t，b、a分别为系统传递函数的分子和分母系数向量。

[h,t] = impz(b,a,n)返回指定的n点冲激响应[h,t] = impz(b,a,n,Fs)指定了冲激响应采样点的频率间隔1/Fs。

Fs 为相对频率，缺省值为1。

impz(b,a)和impz(...)绘制冲激响应的图形。

举例：计算线性系统(b,a)的冲激响应，结果见图1.4.1。

b =[0.2 0.1 0.3 0.1 0.2];a =[1 ?.1 1.55 ?.7 0.3];impz(b,a,50)3 zplane功能：离散系统的零极点图。

格式：zplane(z,p)zplane(b,a)说明：zplane(z,p)和zplane(b,a)绘制系统的零极点图，用“o”表示零点，“x”表示极点。

z、p分别为零点和极点向量，b、a分别为系统传递函数的分子和分母系数向量。

举例：计算线性系统(b,a)的零点和极点，结果见图1.4.2。

b =[0.2 0.1 0.3 0.1 0.2];a =[1.0 -1.1 1.5 -0.7 0.3];zplane(b,a)4 abs功能：求幅值。

格式：y = abs(x)说明：y = abs(x)返回复数向量x的幅值向量y。

使用MATLAB设计FIR滤波器

使⽤MATLAB设计FIR滤波器1. 采⽤fir1函数设计，fir1函数可以设计低通、带通、⾼通、带阻等多种类型的具有严格线性相位特性的FIR滤波器。

语法形式：b = fir1(n, wn)b = fir1(n, wn, ‘ftype’)b = fir1(n, wn, ‘ftype’, window)b = fir1(n, wn, ‘ftype’, window, ‘noscale’)参数的意义及作⽤：b：返回的FIR滤波器单位脉冲响应，脉冲响应为偶对称，长度为n+1；n：滤波器的介数；wn：滤波器的截⽌频率，取值范围为0<wn<1，1对应信号采样频率⼀半。

如果wn是单个数值，且ftype参数为low，则表⽰设计截⽌频率为wn的低通滤波器，如果ftype参数为high，则表⽰设计截⽌频率为wn的⾼通滤波器；如果wn是有两个数组成的向量[wn1wn2]，ftype为stop，则表⽰设计带阻滤波器，ftype为bandpass，则表⽰设计带通滤波器；如果wn是由多个数组成的向量，则根据ftype的值设计多个通带或阻带范围的滤波器，ftype为DC-1，表⽰设计的第⼀个频带为通带，ftype为DC-0，表⽰设计的第⼀个频带为阻带；window：指定使⽤的窗函数，默认为海明窗；noscale：指定是否归⼀化滤波器的幅度。

⽰例：N=41; %滤波器长度fs=2000; %采样频率%各种滤波器的特征频率fc_lpf=200;fc_hpf=200;fp_bandpass=[200 400];fc_stop=[200 400];%以采样频率的⼀半，对频率进⾏归⼀化处理wn_lpf=fc_lpf*2/fs;wn_hpf=fc_hpf*2/fs;wn_bandpass=fp_bandpass*2/fs;wn_stop=fc_stop*2/fs;%采⽤fir1函数设计FIR滤波器b_lpf=fir1(N-1,wn_lpf);b_hpf=fir1(N-1,wn_hpf,'high');b_bandpass=fir1(N-1,wn_bandpass,'bandpass');b_stop=fir1(N-1,wn_stop,'stop');%求滤波器的幅频响应m_lpf=20*log(abs(fft(b_lpf)))/log(10);m_hpf=20*log(abs(fft(b_hpf)))/log(10);m_bandpass=20*log(abs(fft(b_bandpass)))/log(10);m_stop=20*log(abs(fft(b_stop)))/log(10);%设置幅频响应的横坐标单位为Hzx_f=0:(fs/length(m_lpf)):fs/2;%绘制单位脉冲响应%绘制单位脉冲响应subplot(421);stem(b_lpf);xlabel('n');ylabel('h(n)');subplot(423);stem(b_hpf);xlabel('n');ylabel('h(n)');subplot(425);stem(b_bandpass);xlabel('n');ylabel('h(n)');subplot(427);stem(b_stop);xlabel('n');ylabel('h(n)');%绘制幅频响应曲线subplot(422);plot(x_f,m_lpf(1:length(x_f)));xlabel('频率(Hz)','fontsize',8);ylabel('幅度(dB)','fontsize',8);subplot(424);plot(x_f,m_hpf(1:length(x_f)));xlabel('频率(Hz)','fontsize',8);ylabel('幅度(dB)','fontsize',8);subplot(426);plot(x_f,m_bandpass(1:length(x_f)));xlabel('频率(Hz)','fontsize',8);ylabel('幅度(dB)','fontsize',8);subplot(428);plot(x_f,m_stop(1:length(x_f)));xlabel('频率(Hz)','fontsize',8);ylabel('幅度(dB)','fontsize',8);2. 采⽤fir2函数设计，函数算法是：⾸先根据要求的幅频响应向量形式进⾏插值，然后进⾏傅⾥叶变换得到理想滤波器的单位脉冲响应，最后利⽤窗函数对理想滤波器的单位脉冲响应激进型截断处理，由此得到FIR滤波器系数。

matlab自带的滤波器函数

matlab自带的滤波器函数
Matlab自带的滤波器函数可以用于对信号进行滤波处理，常用的函数有：
1. fir1函数：设计一阶低通、高通、带通、带阻滤波器的FIR 数字滤波器，可自定义通带和阻带的截止频率。

2. cheby1函数：设计ChebyshevI型低通、高通、带通、带阻数字滤波器，可自定义通带和阻带的截止频率和最大通带波纹。

3. butter函数：设计Butterworth型低通、高通、带通、带阻数字滤波器，可自定义通带和阻带的截止频率和滤波器阶数。

4. filtfilt函数：对信号进行双向滤波处理，可避免滤波后信号的相位畸变和滞后。

这些函数可以在Matlab的Signal Processing Toolbox中找到，可根据需要选择合适的函数进行滤波处理。

- 1 -。

matlab滤波器设计(源代码)

某合成信号，表达式如下：f=10cos(2pi*30t)+cos(2pi*150t)+5cos(2pi*600t),请设计三个滤波器，分别提取出信号中各频率分量，并分别绘制出通过这三个滤波器后信号的时域波形和频谱这个信号的频率分量分别为30、150和600Hz，因此可分别设计一个低通、带通和高通的滤波器来提取。

以FIR滤波器为例，程序如下：clear;fs=2000;t=(1:1000)/fs;x=10*cos(2*pi*30*t)+cos(2*pi*150*t)+5*cos(2*pi*600*t);L=length(x);N=2^(nextpow2(L));Hw=fft(x,N);figure(1);subplot(2,1,1);plot(t,x);grid on;title('滤波前信号x');xlabel('时间/s');% 原始信号subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw));% 查看信号频谱grid on;title('滤波前信号频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');%% x_1=10*cos(2*pi*30*t)Ap=1;As=60;% 定义通带及阻带衰减dev=[(10^(Ap/20)-1)/(10^(Ap/20)+1),10^(-As/20)];% 计算偏移量mags=[1,0];% 低通fcuts=[60,100];% 边界频率[N,Wn,beta,ftype]=kaiserord(fcuts,mags,dev,fs);% 估算FIR滤波器阶数hh1=fir1(N,Wn,ftype,kaiser(N+1,beta));% FIR滤波器设计x_1=filter(hh1,1,x);% 滤波x_1(1:ceil(N/2))=[];% 群延时N/2，删除无用信号部分L=length(x_1);N=2^(nextpow2(L));Hw_1=fft(x_1,N);figure(2);subplot(2,1,1);plot(t(1:L),x_1);grid on;title('x_1=10*cos(2*pi*30*t)');xlabel('时间/s');subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw_1));% 查看信号频谱grid on;title('滤波后信号x_1频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');%% x_2=cos(2*pi*150*t)Ap=1;As=60;% 定义通带及阻带衰减dev=[10^(-As/20),(10^(Ap/20)-1)/(10^(Ap/20)+1),10^(-As/20)];% 计算偏移量mags=[0,1,0];% 带通fcuts=[80,120,180,220];% 边界频率[N,Wn,beta,ftype]=kaiserord(fcuts,mags,dev,fs);% 估算FIR滤波器阶数hh2=fir1(N,Wn,ftype,kaiser(N+1,beta));% FIR滤波器设计x_2=filter(hh2,1,x);% 滤波x_2(1:ceil(N/2))=[];% 群延时N/2，删除无用信号部分L=length(x_2);N=2^(nextpow2(L));Hw_2=fft(x_2,N);figure(3);subplot(2,1,1);plot(t(1:L),x_2);grid on;title('x_2=cos(2*pi*150*t)');xlabel('时间/s');subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw_2));% 查看信号频谱grid on;title('滤波后信号x_2频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');%% x_3=5*cos(2*pi*600*t)Ap=1;As=60;% 定义通带及阻带衰减dev=[10^(-As/20),(10^(Ap/20)-1)/(10^(Ap/20)+1)];% 计算偏移量mags=[0,1];% 高通fcuts=[500,550];% 边界频率[N,Wn,beta,ftype]=kaiserord(fcuts,mags,dev,fs);% 估算FIR滤波器阶数hh2=fir1(N,Wn,ftype,kaiser(N+1,beta));% FIR滤波器设计x_3=filter(hh2,1,x);% 滤波x_3(1:ceil(N/2))=[];% 群延时N/2，删除无用信号部分L=length(x_3);N=2^(nextpow2(L));Hw_3=fft(x_3,N);figure(4);subplot(2,1,1);plot(t(1:L),x_3);grid on;title('x_3=5*cos(2*pi*600*t)');xlabel('时间/s');subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw_3));% 查看信号频谱grid on;title('滤波后信号x_3频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');。

matlab滤波函数详解

matlab滤波函数详解Matlab作为一种广泛应用于数值计算和数据处理的软件，提供了许多用于信号处理和图像处理的函数。

其中，滤波函数是其中非常重要的一部分，它们在许多应用中都起着关键的作用。

本文将详细介绍Matlab中常见的滤波函数，包括它们的用途、参数设置、使用方法和示例。

一、滤波函数概述滤波函数主要用于对信号进行滤波处理，以消除噪声、突出信号特征或实现其他特定的处理目标。

在Matlab中，常见的滤波函数包括低通、高通、带通、带阻等类型，它们可以根据不同的应用需求选择。

滤波器通常由一组数学函数组成，用于对输入信号进行加权和叠加，以达到滤波的目的。

二、低通滤波函数低通滤波函数用于消除高频噪声，保持低频信号的完整性。

在Matlab中，常用的低通滤波函数包括lfilter和filter等。

lfilter函数适用于线性滤波器，而filter函数适用于任意滤波器设计。

低通滤波函数的参数包括滤波器系数、输入信号和采样率等。

通过调整滤波器系数，可以实现不同的滤波效果。

三、高通滤波函数高通滤波函数用于消除低频噪声，突出高频信号特征。

在Matlab 中，常用的高通滤波函数包括hilbert和highpass等。

hilbert函数适用于频谱分析和高频信号提取，而highpass函数则适用于消除低频噪声。

高通滤波函数的参数包括滤波器系数、采样率和信号类型等。

通过调整滤波器系数，可以实现不同的高通效果。

四、带通滤波函数带通滤波函数用于选择特定频率范围内的信号进行过滤。

在Matlab中，常用的带通滤波函数包括bandpass和butter等。

bandpass函数适用于设计带通滤波器，而butter函数则适用于连续时间滤波器设计。

带通滤波函数的参数包括带外抑制值、带宽和采样率等。

通过调整带宽参数，可以实现不同的带通效果。

五、其他滤波函数除了以上三种常见的滤波函数外，Matlab还提供了其他一些滤波函数，如带阻、限幅、防混叠等类型。

(整理)各类滤波器的MATLAB程序清单.

各类滤波器的MATLAB程序一、理想低通滤波器IA=imread('lena.bmp');[f1,f2]=freqspace(size(IA),'meshgrid');Hd=ones(size(IA));r=sqrt(f1.^2+f2.^2);Hd(r>0.2)=0;Y=fft2(double(IA));Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA));subplot(2,2,2),imshow(uint8(Ia));figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');二、理想高通滤波器IA=imread('lena.bmp');[f1,f2]=freqspace(size(IA),'meshgrid');Hd=ones(size(IA));r=sqrt(f1.^2+f2.^2);Hd(r<0.2)=0;Y=fft2(double(IA));Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya));figuresubplot(2,2,1),imshow(uint8(IA));subplot(2,2,2),imshow(uint8(Ia));figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');三、B utterworth低通滤波器IA=imread('lena.bmp');[f1,f2]=freqspace(size(IA),'meshgrid');D=0.3;r=f1.^2+f2.^2;n=4;for i=1:size(IA,1)for j=1:size(IA,2)t=r(i,j)/(D*D);Hd(i,j)=1/(t^n+1);endendY=fft2(double(IA));Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya));figuresubplot(2,2,1),imshow(uint8(IA));subplot(2,2,2),imshow(uint8(Ia));figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');四、B utterworth高通滤波器IA=imread('lena.bmp');[f1,f2]=freqspace(size(IA),'meshgrid');D=0.3;r=f1.^2+f2.^2;n=4;for i=1:size(IA,1)for j=1:size(IA,2)t=(D*D)/r(i,j);Hd(i,j)=1/(t^n+1);endendY=fft2(double(IA));Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya));figuresubplot(2,2,1),imshow(uint8(IA));subplot(2,2,2),imshow(uint8(Ia));figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');五、高斯低通滤波器IA=imread('lena.bmp');IB=imread('babarra.bmp');[f1,f2]=freqspace(size(IA),'meshgrid');D=100/size(IA,1);r=f1.^2+f2.^2;Hd=ones(size(IA));for i=1:size(IA,1)for j=1:size(IA,2)t=r(i,j)/(D*D);Hd(i,j)=exp(-t);endendY=fft2(double(IA));Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya));figuresubplot(2,2,1),imshow(uint8(IA));subplot(2,2,2),imshow(uint8(Ia));figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');六、高斯高通滤波器IA=imread('lena.bmp');IB=imread('babarra.bmp');[f1,f2]=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D=0.3;r=f1.^2+f2.^2;for i=1:size(IA,1)for j=1:size(IA,2)t=r(i,j)/(D*D);Hd(i,j)=1-exp(-t);endendY=fft2(double(IA));Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya));figuresubplot(2,2,1),imshow(uint8(IA));subplot(2,2,2),imshow(uint8(Ia));figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');七、梯形低通滤波器IA=imread('lena.bmp');IB=imread('babarra.bmp');[f1,f2]=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D0=0.1;D1=0.4;r=sqrt(f1.^2+f2.^2);Hd=zeros(size(IA));Hd(r<D0)=1;for i=1:size(IA,1)for j=1:size(IA,2)if r(i,j)>=D0 & r(i,j)<=D1Hd(i,j)=(D1-r(i,j))/(D1-D0);endendendY=fft2(double(IA));Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya));figuresubplot(2,2,1),imshow(uint8(IA));subplot(2,2,2),imshow(uint8(Ia));figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');八、梯形高通滤波器IA=imread('lena.bmp');IB=imread('babarra.bmp');[f1,f2]=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D0=0.1;D1=0.4;r=sqrt(f1.^2+f2.^2);Hd=ones(size(IA));Hd(r<D1)=0;for i=1:size(IA,1)for j=1:size(IA,2)if r(i,j)>=D0 & r(i,j)<=D1Hd(i,j)=(D0-r(i,j))/(D0-D1);endendendY=fft2(double(IA));Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya));figuresubplot(2,2,1),imshow(uint8(IA));subplot(2,2,2),imshow(uint8(Ia));figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');九、用其他方法编写的理想低通、理想高通、Butterworth低通、同态滤波程序1、理想低通i1=imread('lena.bmp');i2=imnoise(i1,'salt & pepper',0.1);f=double(i2);k=fft2(f);g=fftshift(k);[N1,N2]=size(g);d0=50;u0=floor(N1/2)+1;v0=floor(N2/2)+1;for i=1:N1for j=1:N2d=sqrt((i-u0)^2+(j-v0)^2);if d<=d0h=1;elseh=0;endy(i,j)=g(i,j)*h;endendy=ifftshift(y);E1=ifft2(y);E2=real(E1);figuresubplot(2,2,1),imshow(uint8(i1));subplot(2,2,2),imshow(uint8(i2));subplot(2,2,3),imshow(uint8(E2));2、理想高通i1=imread('lena.bmp');i2=imnoise(i1,'salt & pepper',0.1); f=double(i2);k=fft2(f);g=fftshift(k);[N1,N2]=size(g);n=2;d0=10;u0=floor(N1/2)+1;v0=floor(N2/2)+1;for i=1:N1for j=1:N2d=sqrt((i-u0)^2+(j-v0)^2);if d<=d0h=0;else h=1;endy(i,j)=g(i,j)*h;endendy=ifftshift(y);E1=ifft2(y);E2=real(E1);figuresubplot(2,2,1),imshow(uint8(i1)); subplot(2,2,2),imshow(uint8(i2)); subplot(2,2,3),imshow(uint8(E2));3、Butterworth低通i1=imread('lena.bmp');i2=imnoise(i1,'salt & pepper',0.1); f=double(i2);k=fft2(f);g=fftshift(k);[N1,N2]=size(g);n=2;d0=50;u0=floor(N1/2)+1;v0=floor(N2/2)+1;for i=1:N1for j=1:N2d=sqrt((i-u0)^2+(j-v0)^2);h=1/(1+(d/d0)^(2*n));y(i,j)=g(i,j)*h;endendy=ifftshift(y);E1=ifft2(y);E2=real(E1);figuresubplot(2,2,1),imshow(uint8(i1)); subplot(2,2,2),imshow(uint8(i2)); subplot(2,2,3),imshow(uint8(E2));4、同态滤波I=rgb2gray(imread('fabric00.bmp')); [M,N]=size(I);T=double(I);L=log(T);F=fft2(L);A=2;B=0.3;for i=1:Mfor j=1:ND(i,j)=((i-M/2)^2+(j-N/2)^2);endendc=1.1;%锐化参数D0=max(M,N);H=(A-B)*(1-exp(c*(-D/(D0^2))))+B;F=F.*H;F=ifft2(F);Y=exp(F);figuresubplot(1,2,1),imshow(I);subplot(1,2,2),imshow(uint8(real(Y)));十、十一、Gabor滤波器。

matlab之经典数字滤波函数介绍

matlab 之经典数字滤波函数介绍南京理工大学仪器科学与技术专业谭彩铭2010-3-121 butter 函数设计一个9阶高通Butterworth 数字滤波器，截止频率为300Hz图1图2下面看一看freqz 函数例如对离散系统传递函数1212(1)(2)(3)()(1)(2)(3)b b z b z H z a a za z----++=++freqz 函数的的主要计算环节是计算()sjw T H e编写下图所示程序验证之图3图4可见h1和h相等图2中出现小于-360度是否表达其他不同的意义？-361度和-1度有什么区别吗？对于正弦波应该是一样的，故理论上说应该没有区别。

butter函数的原理是什么？顾名思义，butter函数的原理是基于Butterworth滤波器。

这里始终要带着这个问题去研究，滤波系数本身有什么特性竟然可使低频的滤掉，高频的通过，其实这里想要寻找的是敏捷控程，理论上这个问题的答案已经很成熟。

2 impinvar函数用冲击响应不变法数字仿真模拟Butterworth滤波器，程序如下。

图53 bilinear 函数用双线性变换法数字仿真模拟Butterworth 滤波器，程序如下。

图64 fir1函数由理想滤波器幅频特性反推滤波系数，得出来的系数数量是无穷多的。

故可采用加窗的方法舍去部分，留下有限的滤波系数数量，使仍能基本达到需要的滤波效果。

图7图85 fir2函数fir2函数的基本原理同fir1函数，它的功能更进一层，可以设计任意形状的频率响应图形。

图9图10、图9中，w 如果是角频率值，将w 转换成频率值时，应该是将w 除以2*pi ，但是程序中为什么除以的是pi 呢？准确地说，图9中的w 并非是角频率值，而是频率值，freqz 函数调用时若没有加入采样率参数，其返回的频率值的范围是0~pi 。

6 fir1函数补充对调用方式b=fir(n,wn)，wn 的取值范围是（0，1），其中1对应于0.5fs （fs 为采样率）。

基于MATLAB设计FIR滤波器

基于MATLAB设计FIR滤波器FIR（Finite Impulse Response）滤波器是一种数字滤波器，它具有有限的冲激响应长度。

基于MATLAB设计FIR滤波器可以使用signal工具箱中的fir1函数。

fir1函数的语法如下：b = fir1(N, Wn, window)其中，N是滤波器的阶数，Wn是截止频率，window是窗函数。

要设计一个FIR低通滤波器，可以按照以下步骤进行：步骤1：确定滤波器的阶数。

阶数决定了滤波器的截止频率的陡峭程度。

一般情况下，阶数越高，滤波器的陡峭度越高，但计算复杂度也会增加。

步骤2：确定滤波器的截止频率。

截止频率是指在滤波器中将信号的频率限制在一定范围内的频率。

根据应用的需求，可以选择适当的截止频率。

步骤3：选择窗函数。

窗函数是为了在时域上窗口函数中心增加频率衰减因子而使用的函数。

常用的窗函数有Hamming、Hanning等。

窗函数可以用来控制滤波器的幅度响应特性，使得它更平滑。

步骤4：使用fir1函数设计滤波器。

根据以上步骤确定滤波器的阶数、截止频率和窗函数，可以使用fir1函数设计FIR滤波器。

具体代码如下：N=50;%设定阶数Wn=0.5;%设定截止频率window = hanning(N + 1); % 使用Hanning窗函数步骤5：使用filter函数对信号进行滤波。

设计好FIR滤波器后，可以使用filter函数对信号进行滤波。

具体代码如下：filtered_signal = filter(b, 1, input_signal);其中，input_signal是输入信号，filtered_signal是滤波后的信号。

以上，便是基于MATLAB设计FIR滤波器的简要步骤和代码示例。

根据具体需求和信号特性，可以进行相应的调整和优化。

用MATLAB设计滤波器

用MATLAB 设计滤波器1 IIR 滤波器的设计● freqz功能：数字滤波器的频率响应。

格式：[h ，w]=freqz(b ，a ，n)[h ，f]=freqz(b ，a ，n ，Fs)[h ，w]=freqz(b ，a ，n ，'whole')[h ，f]=freqz(b ，a ，n ，'whole'，Fs)h=freqz(b ，a ，w)h=freqz(b ，a ，f ，Fs)freqz(b ，a)说明:freqz 用于计算由矢量"和b 构成的数字滤波器H(z)=A(z)B(z)= n-1--n-1 l)z a(n ....a(2)z l l)z b(n .... b(2)z b(l)++++++++ 的复频响应H(j ω)。

[h ，w]=freqz(b ，a ，n)可得到数字滤波器的n 点的幅频响应，这n 个点均匀地分布在上半单位圆(即0~π)，并将这n 点频率记录在w 中，相应的频率响应记录在h 中。

至于n 值的选择没有太多的限制，只要n>0的整数，但最好能选取2的幂次方，这样就可采用 FFT 算法进行快速计算。

如果缺省，则n=512。

[h ，f]二freqz(b ，a ，n ，Fs)允许指定采样终止频率Fs(以Hz 为单位)，也即在0~Fs/2 频率范围内选取n 个频率点(记录在f 中)，并计算相应的频率响应h 。

[h ，w]=freqz(b ，a ，n ，'whole')表示在0~2π之间均匀选取n 个点计算频率响应。

[h ，f]=freqz(b ，a ，n ，'whole'，Fs)则在O~Fs 之间均匀选取n 个点计算频率响应。

h=freqz(b ，a ，w)计算在矢量w 中指定的频率处的频率响应，但必须注意，指定的频率必须介于0和2π之间。

h=freqz(b ，a ，f ，Fs)计算在矢量f 中指定的频率处的频率响应，但指定频率必须介于 0和Fs 之间。

数字滤波器的设计及其MATLAB实现

设计低通数字滤波器，要求在通带内频率低于0.2pirad时，允许幅度误差在1dB以内，在频率0.3pi rad～pi rad之间的阻带衰减大于15dB，用脉冲响应不变法设计数字滤波器，T=1: 切比雪夫I型模拟滤波器的设计子程序：function [b,a]=afd_chb1(Omegap,Omegar,Ar)if Omegap<=0error('通带边缘必须大于0')endif(Dt<=0)|(Ar<0)error('通带波动或阻带衰减必须大于0');endep=sqrt(10^(Dt/10)-1);A=10^(Ar/20);OmegaC=Omegap;OmegaR=Omegar/Omegap;g=sqrt(A*A-1)/ep;N=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)));fprintf('\n***切比雪夫I型模拟低通滤波器阶数=%2.0f\n',N);[b,a]=u_chblap(N,Dt,OmegaC);设计非归一化切比雪夫I型模拟低通滤波器原型程序：function [b,a]=u_chblap(N,Dt,OmegaC)[z,p,k]=cheb1ap(N,Dt);a=real(poly(p));aNn=a(N+1);p=p*OmegaC;a=real(poly(p));aNu=a(N+1);k=k*aNu/aNn;b0=k;B=real(poly(z));b=k*B;直接形式转换成级联形式子程序：function [C,B,A]=sdir2cas(b,a)Na=length(a)-1;Nb=length(b)-1;b0=b(1);b=b/b0;a0=a(1);a=a/a0;C=b0/a0;p=cplxpair(roots(a));K=floor(Na/2);if K*2==NaA=zeros(K,3);for n=1:2:NaArow=p(n:1:n+1,:);Arow=poly(Arow);A((fix(n+1)/2),:)=real(Arow);elseif Na==1A=[0 real(poly(p))];elseA=zeros(K+1,3);for n=1:2:2*KArow=p(n:1:n+1,:);Arow=poly(Arow);A((fix(n+1)/2),:)=real(Arow);endA(K+1,:)=[0 real(poly(p(Na)))];endz=cplxpair(roots(b));K=floor(Nb/2);if Nb==0B=[0 0 poly(z)];elseif K*2==NbB=zeros(K,3);for n=1:2:NbBrow=z(n:1:n+1,:);Brow=poly(Brow);B((fix(n+1)/2),:)=real(Brow);endelseif Nb==1B=[0 real(poly(z))];elseB=zeros(K+1,3);for n=1:2:2*KBrow=z(n:1:n+1,:);Brow=poly(Brow);B((fix(n+1)/2),:)=real(Brow);endB(K+1,:)=[0 real(poly(z(Nb)))];End计算系统函数的幅度响应和相位响应子程序：function [db,mag,pha,w]=freqs_m(b,a,wmax)w1=0:500;w=w1*wmax/500;h=freqs(b,a,w);mag=abs(h);db=20*log10((mag+eps)/max(mag));pha=angle(h);脉冲响应不变法程序：function [b,a]=imp_invr(c,d,T)[R,p,k]=residue(c,d);p=exp(p*T);[b,a]=residuez(R,p,k);b=real(b).*T;数字滤波器响应子程序：function [db,mag,pha,grd,w]=freqz_m(b,a);[H,w]=freqz(b,a,1000,'whole');H=(H(1:501))';w=(w(1:501))';mag=abs(H);db=20*log10((mag+eps)/max(mag));pha=angle(H);grd=grpdelay(b,a,w);直接转换成并联型子程序：function [C,B,A]=dir2par(b,a)M=length(b);N=length(a);[r1,p1,C]=residuez(b,a);p=cplxpair(p1,10000000*eps);x=cplxcomp(p1,p);r=r1(x);K=floor(N/2);B=zeros(K,2);A=zeros(K,3);if K*2==Nfor i=1:2:N-2br=r(i:1:i+1,:);ar=p(i:1:i+1,:);[br,ar]=residuez(br,ar,[]);B((fix(i+1)/2),:)real(br');A((fix(i+1)/2),:)real(ar');end[br,ar]=residuez(r(N-1),p(N-1),[]);B(K,:)=[real(br') 0];A(K,:)=[real(ar') 0];elsefor i=1:2:N-1br=r(i:1:i+1,:);ar=p(i:1:i+1,:);[br,ar]=residuez(br,ar,[]);B((fix(i+1)/2),:)real(br);A((fix(i+1)/2),:)real(ar);endEnd比较两个含同样标量元素但（可能）有不同下标的复数对及其相位留数向量子程序：function I=cplxcomp(p1,p2)I=[];for i=1:length(p2)for j=1:length(p1)if(abs(p1(j)-p2(i))<0.0001)I=[I,j];endendendI=I';双线性变换巴特沃斯低通滤波器设计：巴特沃思模拟滤波器的设计子程序：function [b,a]=afd_butt(wp,ws,Rp,rs)if wp<=0error('通带边缘必须大于0');endif ws<=wperror('阻带边缘必须大于通带边缘');endif(Rp<=0)|(Rs<0)error('通带波动或阻带衰减必须大于0');endN=ceil((log10((10^(Rp/10)-1)/(10^(Rs/10)-1)))/(2*log10(wp/ws))); fprintf('\n***Butterworth Filter Order=%2.0f\n',N);OmegaC=wp/((10^(Rp/10)-1)^(1/(2*N)));[b,a]=u_buttap(N,OmegaC)设计非归一化巴特沃思模拟低通滤波器原型子程序：function [b,a]=u_buttap(N,OmegaC)[z,p,k]=buttap(N);p=p*OmegaC;k=k*OmegaC^N;B=real(poly(z));b0=k;b=k*B;a=real(poly(p));直接型到级联型形式的转换：function [b0,B,A]=dir2cas(b,a)b0=b(1);b=b/b0;a0=a(1);a=a/a0;b0=b0/a0;M=length(b);N=length(a);if N>Mb=[b,zeros(1,N-M)];a=[a,zeros(1,M-N)];elseNM=0;endk=floor(N/2);B=zeros(k,3);A=zeros(k,3);if k*2==Nb=[b,0];a=[a,0];endbroots=cplxpair(roots(b));aroots=cplxpair(roots(a));for i=1:2:2*kbr=broots(i:1:i+1,:);br=real(polt(br));B((fix(i+1)/2),:)=br;ar=aroots(i:1:i+1,:);ar=real(polt(ar));A((fix(i+1)/2),:)=ar;Endfunction [db,mag,pha,grd,w]=freqz_m(b,a)[h,w]=freqz(b,a,1000,'whole');h=(h(1:501))';w=(w(1:501))';mag=abs(h);db=20*log10((mag+eps)/max(mag));pha=angle(h);grd=grdelay(b,a,w);设计一个巴特沃思高通滤波器，要求通带截止频率为0.6pi，通带内衰减不大于1dB，阻带·起始频率为0.4pi，阻带内衰减不小于15dB，T=1：>> wp=0.6*pi;ws=0.4*pi;>> Rp=1;Rs=15;T=1;>> [N,wn]=buttord(wp/pi,ws/pi,Rp,Rs) 计算巴特沃思滤波器阶数和截止频率N =4wn =>> [b,a]=butter(N,wn,'high'); 频率变换法计算巴特沃思高通滤波器>> [C,B,A]=dir2cas(b,a)C =0.0751B =1.0000 -2.0000 1.00001.0000 -2.0000 1.0000A =1.0000 0.1562 0.44881.0000 0.1124 0.0425>> [db,mag,pha,grd,w]=freqz_m(b,a);>> subplot(2,1,1);plot(w/pi,mag);>> subplot(2,1,2);plot(w/pi,db);椭圆带通滤波器的设计--ellip函数的应用：>> ws=[0.3*pi 0.75*pi]; 数字阻带边缘频率>> wp=[0.4*pi 0.6*pi]; 数字通带边缘频率>> Rp=1;Rs=40;>> Ripple=10^(-Rp/20); 通带波动>> Attn=10^(-Rs/20); 阻带衰减>> [N,wn]=ellipord(wp/pi,ws/pi,Rp,Rs) 计算椭圆滤波器参数N =4wn =0.4000 0.6000>> [b,a]=ellip(N,Rp,Rs,wn); 数字椭圆滤波器的设计>> [b0,B,A]=dir2cas(b,a) 级联形式实现b0 =0.0197B =1.0000 1.5066 1.00001.0000 0.9268 1.00001.0000 -0.9268 1.00001.0000 -1.5066 1.0000A =1.0000 0.5963 0.93991.0000 0.2774 0.79291.0000 -0.2774 0.79291.0000 -0.5963 0.9399>> figure(1);>> [db,mag,pha,grd,w]=freqz_m(b,a);>> subplot(2,2,1);plot(w/pi,mag);>> grid on;>> subplot(2,2,3);plot(w/pi,db);grid on;>> subplot(2,2,2);plot(w/pi,pha/pi);grid on;>> subplot(2,2,4);plot(w/pi,grd);设计一个巴特沃思带阻滤波器，要求通带上下截止频率为0.8pi、0.2pi，通带内衰减不大于1dB，阻带上起始频率为0.7pi、0.4pi，阻带内衰减不小于30dB：>> wp=[0.2*pi 0.8*pi];>> ws=[0.4*pi 0.7*pi];>> Rp=1;Rs=30;>> [N,wn]=buttord(wp/pi,ws/pi,Rp,Rs);>> [b,a]=butter(N,wn,'stop');>> [C,B,A]=dir2cas(b,a)C =0.0394B =1.0000 0.3559 0.99941.0000 0.3547 1.00401.0000 0.3522 0.99541.0000 0.3499 1.00461.0000 0.3475 0.99601.0000 0.3463 1.0006A =1.0000 1.3568 0.79281.0000 1.0330 0.46331.0000 0.6180 0.17751.0000 -0.2493 0.11131.0000 -0.6617 0.37551.0000 -0.9782 0.7446>> [db,mag,pha,grd,w]=freqz_m(b,a); >> subplot(2,1,1);plot(w/pi,mag);>> subplot(2,1,2);plot(w/pi);数字低通---数字带阻：function [bz,az]=zmapping(bZ,aZ,Nz,Dz) bzord=(length(bZ)-1)*(length(Nz)-1); azord=(length(aZ)-1)*(length(Dz)-1);bz=zeros(1,bzord+1);for k=0:bzordpln=[1];for i=0:k-1pln=conv(pln,Nz);endpld=[1];for i=0:bzord-k-1pld=conv(pld,Dz);endbz=bz+bZ(k+1)*conv(pln,pld); endfor k=0:azordpln=[1];for i=0:k-1pln=conv(pln,Nz);endpld=[1];for i=0:azord-k-1pld=conv(pld,Dz);endaz=az+aZ(k+1)*conv(pln,pld); endall=az(1);az=az/az1;bz=bz/az1;线性相位FIR滤波器的幅度特性：function pzkplot(num,den)hold on;axis('square');x=-1:0.01:1;y=(1-x.^2).^0.5;y1=-(1-x.^2).^0.5;plot(x,y,'b',x,y1,'b');num1=length(num);den1=length(den);if(num1>1)z=roots(num);elsez=0;endif(den1>1)p=roots(den);elsep=0;endif(num>1&den1>1)r_max_z=max(abs(real(z)));i_max_z=max(abs(imag(z)));a_max_z=max(r_max_z,i_max_z);r_max_p=max(abs(real(p)));i_max_p=max(abs(imag(p)));a_max_p=max(r_max_p,i_max_p);a_max=max(a_max_z,a_max_p);elseif (num1>1)r_max_z=max(abs(real(z)));i_max_z=max(abs(imag(z)));a_max=max(r_max_z,i_max_z);elser_max_p=max(abs(real(p)));i_max_p=max(abs(imag(p)));a_max=max(r_max_p,i_max_p);endaxis([-a_max a_max -a_max a_max]);plot([-a_max a_max],[0 0],'b');plot([0 0],[-a_max a_max],'b');plot([-a_max a_max],[a_max a_max],'b');plot([a_max a_max],[-a_max a_max],'b');Lz=length(z);for i=1:Lz;plot(real(z(i)),imag(z(i)),'bo');endLp=length(p);for j=1:Lpplot(real(p(j)),imag(p(j)),'bx');endtitle('The zeros-pole plot');xlabel('虚部');ylabel('实部');function [Hr,w,a,L]=Hr_Type1(h)M=length(h);L=(M-1)/2;a=[h(L+1) 2*h(L:-1:1)];n=[0:1:L];w=[0:1:500]'*pi/500;Hr=cos(w*n)*a';设计I型线性相位FIR滤波器：>> h=[-4 1 -1 -2 5 6 5 -2 -1 1 -4];>> M=length(h);n=0:M-1;>> [Hr,w,a,L]=Hr_Type1(h);>> amax=max(a)+1;>> amin=min(a)-1;>> subplot(2,2,1);stem(n,h);>> axis([-1 2*L+1 amin amax]);text(2*L+1.5,amin,'n'); >> xlabel('n');ylabel('h(n)');title('脉冲响应');>> subplot(2,2,3);stem(0:L,a);>> axis([-1 2*L+1 amin amax]);>> xlabel('n');ylabel('a(n)');title('a(n) 系数');>> subplot(2,2,2);plot(w/pi,Hr);>> grid on;text(1.05,-20,'频率pi');>> xlabel('频率');ylabel('Hr');title('I 型振幅响应');>> subplot(2,2,4);pzkplot(h,1);>> title('零极点分布');function [hr,w,b,L]=Hr_Type2(h)M=length(h);L=M/2;b=2*h(L:-1:1);n=[1:1:L];n=n-0.5;w=[0:1:500]'*pi/500;hr=cos(w*n)*b';II型线性相位FIR滤波器：>> h=[-4 1 -1 -2 5 6 5 -2 -1 1 -4];>> M=length(h);n=0:M-1;>> [Hr,w,b,L]=Hr_Type2(h);Warning: Integer operands are required for colon operator when used as index. > In Hr_Type2 at 2>> bmax=max(b)+1;bmin=min(b)-1;>> subplot(2,2,1);stem(n,h);axis([-1 2*L+1 bmin bmax]);text(2*L+1.5,bmin,'n');xlabel('n');ylabel('h(n)');title('脉冲响应');>> subplot(2,2,3);stem(1:L,b);axis([-1 2*L+1 bmin bmax]);xlabel('n');ylabel('b(n)');title('b(n) 系数');>> subplot(2,2,2);plot(w/pi,Hr);grid on;text(1.05,-20,'频率pi');xlabel('频率');ylabel('Hr');title('II 型振幅响应');>> subplot(2,2,4);pzkplot(h,1);title('零极点分布');function [hr,w,c,L]=Hr_Type3(h)M=length(h);L=(M-1)/2;b=2*h(L+1:-1:1);n=[1:1:L];w=[0:1:500]'*pi/500;hr=cos(w*n)*c';用MA TLAB编程绘制各种窗函数的形状。

matlab中关于数字滤波器的函数介绍

matlab中关于数字滤波器的函数介绍MATLAB下的数字信号处理实现示例一信号、系统和系统响应1、理想采样信号序列（1）首先产生信号x(n),0<=n<=50n=0:50; %定义序列的长度是50A=444.128; %设置信号有关的参数a=50*sqrt(2.0)*pi;T=0.001; %采样率w0=50*sqrt(2.0)*pi;x=A*exp(-a*n*T).*sin(w0*n*T); %pi是MATLAB定义的π，信号乘可采用“.*”close all %清除已经绘制的x(n)图形subplot(3,1,1);stem(x); %绘制x(n)的图形title(‘理想采样信号序列’);（2）绘制信号x(n)的幅度谱和相位谱k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n’*k);magX=abs(X); %绘制x(n)的幅度谱subplot(3,1,2);stem(m agX);title(‘理想采样信号序列的幅度谱’);angX=angle(X); %绘制x(n)的相位谱subplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱’)（3）改变参数为：1,0734.2,4.0,10==Ω==TAαn=0:50; %定义序列的长度是50A=1; %设置信号有关的参数a=0.4;T=1; %采样率w0=2.0734;x=A*exp(-a*n*T).*sin(w0*n*T); %pi是MATLAB定义的π，信号乘可采用“.*”close all %清除已经绘制的x(n)图形subplot(3,1,1);stem(x); %绘制x(n)的图形title(‘理想采样信号序列’);k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n’*k);magX=abs(X); %绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’);angX=angle(X); %绘制x(n)的相位谱subplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱’)2、单位脉冲序列在MatLab中，这一函数可以用zeros函数实现：n=1:50; %定义序列的长度是50x=zeros(1,50); %注意：MATLAB中数组下标从1开始x(1)=1;close all;subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);k=-25:25;X=x*(exp(-j*pi/12.5)).^(n’*k);magX=abs(X); %绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);angX=angle(X); %绘制x(n)的相位谱subplot(3,1,3);stem(ang X) ; title (‘单位冲击信号的相位谱’)3、矩形序列n=1:50x=sign(sign(10-n)+1);close all;subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);k=-25:25;X=x*(exp(-j*pi/25)).^(n’*k);magX=abs(X); %绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);angX=angle(X); %绘制x(n)的相位谱su bplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱’)4、特定冲击串：)3()2(5.2)1(5.2)()(?+?+?+=nnnnnxδδδδn=1:50; %定义序列的长度是50x=zeros(1,50); %注意：MATLAB中数组下标从1开始x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1;close all;subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);k=-25:25;X=x*(exp(-j*pi/12.5)).^(n’*k);magX=abs(X); %绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);angX=angle(X); %绘制x(n)的相位谱subplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱’)5、卷积计算：Σ+∞?∞=?=?=mmnhmxnhnxny)()()()()(在MATLAB中。

滤波器设计matlab相关函数

fir1,调用方式： B=fir1(N,Wn); B=fir1(N,Wn,’high’); B=fir1(N,Wn,’stop’); N为滤波器的阶数，因此滤波器的长度为N+1，
Wn为通带截止频率，其值在0-1之间，1对应采样频率的一半；b是设计好的滤波器的系数h(n)。
与FIRDF相关的matlab函数（三）
多项式的系数向量；Wn是通带截止频率，范围在0-1之间， 1对应采样频率的一半。第四个用来设计模拟滤波器。
IIRDF设计实例
设计一个LPDF，数字采样频率为1000Hz，给定技术指标通带截止频率fp=100Hz,阻带截止频率fs=300Hz，通带最大衰减为3db，阻带最小衰减为20db。
Wp=0.2*pi; Ws=0.4*pi; fs=1000; rp=3;rs=20; [n,Wn]=buttord(Wp/pi,Ws/pi,rp,rs);%确定滤波器阶数 [bz,az]=butter(n,Wp/pi);%利用butter函数设计滤波器 [h,w]=freqz(bz,az,128,fs);%计算幅频响应 Plot(w,abs(h)); grid
向量；B,A分别是转换后的H(s)的分子和分母多项式的系数向量。
与IIRDF相关的matlab函数（四）
Bilinear：实现双线性变换，并由模拟滤波器H(s)得到数字滤波器H(z)。调用格式：
[Bz,Az]=bilinear(B,A,Fs) B和A分别是H(s)的分子和分母多项式的系
数向量；Bz和Az分别是H(z)的分子和分母多项式系数向量；Fs是采样频率。
Design Filter部分
Filter Order(滤波器阶数)选项，定义滤波器的阶数，包括Specify Order(指定阶数)和 Minimum Order(最小阶数)。在Specify Order中填入所要设计的滤波器的阶数(N阶滤波器，Specify Order=N-1)。如果选择 Minimum Order，则MATLAB根据所选择的滤波器类型自动使用最小阶数。

matlab滤波器设计命令

matlab滤波器设计命令Matlab滤波器设计命令滤波器是数字信号处理中常用的工具，用于去除信号中的噪声、频率干扰或其他不需要的成分。

Matlab提供了一系列有用的滤波器设计命令，使用户能够轻松设计并应用各种类型的滤波器。

在本文中，我们将详细介绍Matlab中常用的滤波器设计命令，包括滤波器设计函数、滤波器类型和设计过程。

I. Matlab中常用的滤波器设计函数在Matlab中，有几种函数可用于设计滤波器，其中最常用的函数是`designfilt`函数和`fir1`函数。

1. designfilt函数`designfilt`函数是Matlab中最灵活和功能强大的滤波器设计函数之一，可用于设计各种类型的IIR和FIR滤波器。

它的基本语法如下：`filt = designfilt(FilterType, 'PropertyName', PropertyValue, ...)`其中，`FilterType`代表滤波器类型，包括低通滤波器(Lowpass)、高通滤波器(Highpass)、带通滤波器(Bandpass)、带阻滤波器(Bandstop)等。

`PropertyName`和`PropertyValue`是可选的参数，用于设置滤波器的各种属性，如阶数(Order)、截止频率(CutoffFrequency)、通带和阻带的最大衰减(MaximumAttenuation)等。

下面是一个使用`designfilt`函数设计低通滤波器的例子：Fs = 1000; 采样频率Fpass = 20; 通带截止频率Fstop = 30; 阻带截止频率designfilt('lowpassiir', 'FilterOrder', 4, 'PassbandFrequency', Fpass, 'StopbandFrequency', Fstop, 'SampleRate', Fs)该命令将设计一个4阶的低通IIR滤波器，其通带截止频率为20Hz，阻带截止频率为30Hz，采样频率为1000Hz。

MATLAB FIR滤波器

Matlab中，函数fir1（）和fir2（）利用加窗傅里叶级数法设计FIR滤波器。

函数fir1（）用来设计传统的LP（低通）、HP（高通）、BP（带通）、BS（带阻）和多频带FIR滤波器；而函数fir2（）用来设计具有任意幅度响应的的FIR滤波器。

一、函数fir1的各种形式如下：b= fir1(N,Wn)b= fir1(N,Wn,'ftype')b= fir1(N,Wn,window)b= fir1(N,Wn,'ftype',window)b= fir1(...,'normalization')其中，1、N为滤波器节点个数；2、Wn（0<Wn<1）为归一化截止频率；3、ftype为滤波器类型（默认（缺省时）是low）：•'high' for a highpass filter with cutoff frequency Wn.•'stop' for a bandstop filter, if Wn = [w1 w2]. The stopband frequency range is specified by this interval.•'DC-1' to make the first band of a multiband filter a passband.•'DC-0' to make the first band of a multiband filter a stopband.4、window为所加窗的类型（默认的是hamming窗），e.g hamming、chebwin、blackman、hanning、kaiser等。

5、normalization为是否将滤波器的幅度进行归一化：•'scale' (default): Normalize the filter so that the magnitude response of the filter at the center frequency of the passband is 0 dB.•'noscale': Do not normalize the filter.fir1（）举例：设计一个低通滤波器二、函数fir2的各种形式如下：b = fir2(n,f,m)b = fir2(n,f,m,window)b = fir2(n,f,m,npt)b = fir2(n,f,m,npt,window)b = fir2(n,f,m,npt,lap)b = fir2(n,f,m,npt,lap,window)其中，向量f是指定频率点的幅度响应样本，与m定义的幅度响应样本对应；f和m具有相同的长度，并且f的第一个和最后一个分量分别是0和1；可以对f中的频点进行复制，从而跳变地逼近幅度响应指标。

matlab实现滤波器

虽然matlab提供了很多产生低通滤波器的函数，而且也提供了将低通转换为高通、带通等滤波器的方法函数，以及数字化的函数。

但是为了简化设计及设计者方便考虑，matlab还提供了更为简便的产生各种滤波器的方法。

1 besself功能：贝塞尔（Bessel）模拟滤波器设计。

格式：[b,a] = besself(n,Wn)[b,a] = besself(n,Wn,'ftype')[z,p,k] = besself(...)[A,B,C,D] = besself(...)说明：besself函数可以设计模拟低通、高通、带通和带阻贝塞尔（Bessel）滤波器。

[b,a] = besself(n,Wn)返回截止频率为Wn（单位为弧度/秒）的n阶贝塞尔模拟低通滤波器，b、a分别为滤波器传递函数的分子和分母系数向量（降幂排列）。

当Wn为二元向量，即Wn=[W1 W2] (W1<W2)时，[b,a] = besself(n,Wn)返回一个2n阶模拟带通滤波器，其通带为W1<ω< W2。

[b,a] = besself(n,Wn,'ftype')用于设计高通和带阻滤波器，即·ftype=high时，返回截止频率为Wn的高通滤波器；·ftype=stop时,Wn=[W1 W2] (W1<W2)，返回阻带为W1<ω< W2的模拟带阻滤波器。

[z,p,k] = besself(...)得到滤波器的零-极点增益模型，z、p、k分别为零点向量、极点向量和增益系数。

[A,B,C,D] = besself(...)得到滤波器的状态空间模型。

2 butter功能：巴特沃思（Butterworth）模拟／数字滤波器设计。

格式：[b,a] = butter(n,Wn)[b,a] = butter(n,Wn,'ftype')[b,a] = butter(n,Wn,'s')[b,a] = butter(n,Wn,'ftype','s')[z,p,k] = butter(...)[A,B,C,D] = butter(...)说明：butter函数可以设计模拟或数字的低通、高通、带通和带阻Butterworth 滤波器。

Matlab中产生信号滤波器的函数

Matlab中产⽣信号滤波器的函数1 filter功能：直接II型滤波器。

格式：y = filter(b,a,x)[y,zf] = filter(b,a,x)[...] = filter(b,a,x,zi)[...] = filter(b,a,x,zi,dim)说明：y = filter(b,a,x)计算输⼊信号x的滤波输出y，b、a分别为滤波器传递函数的分⼦和分母系数向量（按降幂排列）。

[y,zf] = filter(b,a,x)返回最终的状态向量zf。

[...] = filter(b,a,x,zi)指定了滤波器的初始条件zi。

举例：计算低通滤波器的冲激响应x = [1 zeros(1,100)];[b,a] = butter(12,400/1000);y = filter(b,a,x);stem(y)2 filtic功能：为滤波器的直接II型实现选择初始条件。

格式：z = filtic(b,a,y,x)z = filtic(b,a,y)说明：z = filtic(b,a,y,x)为直接II型滤波器(b,a)选择初始条件z，b、a分别为滤波器传递函数的分⼦和分母系数向量（按降幂排列）。

x和y分别为以前的输⼊和输出。

z = filtic(b,a,y)中假定x为0。

3 fftfilt功能：基于FFT和重叠相加法的FIR滤波。

格式：y = fftfilt(b,x)y = fftfilt(b,x,n)说明：y = fftfilt(b,x)返回信号x经滤波器b滤波后的结果y。

b为滤波器的系数向量。

函数所采⽤的FFT长度为nfft = 2^nextpow2(n)，数据块的长度为nfft –length(b)+1，其中nextpow2(n)函数返回⼤于n的最⼩2的幂。

4 filtfilt功能：零相位数字滤波。

格式：y = filtfilt(b,a,x)说明：y = filtfilt(b,a,x)对输⼊信号x做前向和反向处理，实现零相位数字滤波。

关于滤波器设计的matlab函数简表

高斯窗
chebwin
切比雪夫窗
tukeywin
Tukey窗
巴特沃斯滤波器设计
cheby1
切比雪夫Ⅰ型滤波器设计
cheby2
切比雪夫Ⅱ型滤波器设计
ellip
椭圆滤波器设计
besself
贝塞尔滤波器设计
模拟滤波器频带变换
lp2bp
低通→带通
lp2bs
低通→带阻
lp2hp
低通→高通
lp2lp
低通→低通
滤波器离散化（由模拟滤波器得到数字滤波器）
bilinear
脉冲响应不变法
关于滤波器设计、实现的
（求幅频、相频响应）
abs
求模值
angle
求相角
freqs
模拟滤波器的频率响应
freqz
数字滤波器的频率响应
grpdelay
群延迟
impz
脉冲响应（离散的）
zplane
画出零极点图
fvtool
滤波器可视化工具
滤波器实现（求输入信号通过滤波器的响应）
关于滤波器设计实现的matlab函数分类函数名功能说明滤波器分析求幅频相频响应absangle求相角freqs模拟滤波器的频率响应freqz数字滤波器的频率响应grpdelay群延迟impz脉冲响应离散的zplane画出零极点图fvtool滤波器可视化工具滤波器实现求输入信号通过滤波器的响应convconv2卷积二维卷积filter求信号通过滤波器的响应iir滤波器阶数估算buttord巴特沃斯滤波器阶数估算cheb1ord切比雪夫型滤波器阶数估算cheb2ord切比雪夫型滤波器阶数估算ellopord椭圆滤波器阶数估算iir数字滤波器设计求系统函数hzbuttercheby1cheby2ellip模拟低通滤波器原型归buttap模拟低通巴特沃斯滤波器原型cheb1ap模拟低通切比雪夫型滤波器原型cheb2ap模拟低通切比雪夫型滤波器原型ellipap模拟低通椭圆滤波器原型besselap模拟低通贝塞尔滤波器原型模拟低通滤波器设计求系统函数hsbutter巴特沃斯滤波器设计cheby1切比雪夫型滤波器设计cheby2切比雪夫型滤波器设计ellip椭圆滤波器设计besself贝塞尔滤波器设计模拟滤波器频带变换lp2bp滤波器离散化由模拟滤波器得到数字滤波器bilinear脉冲响应不变法impinvar双线性变换法fir滤波器设计fir1基于窗函数的fir滤波器设计fir2基于窗函数的任意响应fir滤波器设窗函数boxcar矩形窗rectwin矩形窗bartlett三角窗triang三角窗hanningblackman布莱克曼窗kaisergausswin高斯窗chebwintukeywintukey

滤波器matlab

>> n=31;%定义滤波器阶数32fs=12.8*10^3;fc1=49;fc2=51;w1=2*pi*fc1/fs;w2=2*pi*fc2/fs;%参数转换，将模拟滤波器的技术指标转换为数字滤波器的技术指标window=hanning(n+1);%使用hanning窗函数q=fir1(n,[w1/pi w2/pi],hanning(n+1));%滤波器时域函数,使用标准响应的加窗设计函数fir1w=linspace(0,pi,512);h1=freqz(q,1,512);%进行512个点的傅里叶变换figure(2);plot(w/pi,20*log10(abs(h1)));title('滤波器频谱图');xlabel('频率');ylabel('幅度');grid ;设计FIR低通滤波器，系统频率为50MHz，通带截止频率Fpass为1MHz，阻带截止频率Fstop为4MHz，通带最大衰减Apass为1dB，阻带最小衰减Astop 为30dB。

程序和必要的程序注释谢谢最佳答案只要用一个公式就行。

library IEEE;use IEEE.STD_LOGIC_1164.ALL;use IEEE.STD_LOGIC_ARITH.ALL;use IEEE.STD_LOGIC_UNSIGNED.ALL;entity fir isPort (clk: in std_logic;reset: in std_logic;inpx: in std_logic_vector(11 downto 0);outy: out std_logic_vector(11 downto 0));end fir;architecture beh of fir issignal x0,x1,x2,x3: std_logic_vector(11 downto 0);constant c0:integer :=-1234*32768/1000;constant c1:integer :=2345*32768/10000;constant c2:integer :=5*32768;constant c3:integer :=-3*32768/10000;signal p0,p1,p2,p3:integer;signal sum: integer;beginsample_delay_line:process(clk)beginif rising_edge(clk) thenif reset='1' thenx3 <=(others=>'0');x2 <=(others=>'0');x1 <=(others=>'0');x0 <=(others=>'0');elsex3 <=x2;x2 <=x1;x1 <=x0;x0 <=inpx;end if;end if;end process;p0 <= conv_integer(x0)*c0;p1 <= conv_integer(x1)*c1;p2 <= conv_integer(x2)*c2;p3 <= conv_integer(x3)*c3;sum <=p0+p1+p2+p3;outy <=conv_std_logic_vector(sum/32768,12);end beh;4.1 数字滤波器简介数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。