《基于某matlab地低通滤波器》由会员分享，可在线阅读，更多相关《基于某matlab地低通滤波器（10页珍藏版）》请在金锄头文库上搜索。
1、word基于matlab的低通滤波器摘要:调用MATLAB信号处理工具箱中滤波通过观察滤波器输入输出信号的时域波形与其频谱，建立数字滤波的概念。应用最广泛的是双线性变换法。根本设计过程是：先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； 设计过渡模拟滤波器；将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。关键词:滤波器，matlab，c语言，声音Abstract: call MATLAB signal processing toolbox filtering through the observation
2、 filter of input and output signals time domain waveform and spectrum, establish the concept of digital filter. One of the most widely applied is double linear transformation method. The basic design process is: (1) to a given digital filters index converted into analog filter transition index; (2)
3、the design transition simulation filter; (3) transition simulation filter system function will be converted into digital filter system function. MATLAB signal processing toolbox digital filter function design of IIR is bilinear transformation method.Keywords: filter, matlab, the c language, the voic
4、e一.任务：用matlab软件设计IIR模拟、数字以与各种窗函数的FIR低通滤波器二.设计目的：1了解matlab软件的用途以与用法；2了解用冲激响应不变法设计模拟低通滤波器；3了解用脉冲响应不变法设计的巴特沃思数字低通滤波器；4了解基于汉明窗函数的FIR低通滤波器的设计。三.设计内容：3.1用冲激响应不变法设计模拟低通滤波器3.1.1 设计内容：要求按照设计指标设计无限冲激响应IIR巴特沃什模拟低通滤波器。3.1.2 设计原理：低通滤波器的技术要求用图形表示如下：通带通带通带通带 过渡带 阻带如上图表示了一个频域设计、一维低通滤波器的技术要求图。ws和wp 分别称为通带截止频率和阻带截止频率
5、。通带频率X围为0wwp，阻带频率X围wsw。从wp到ws称为过渡带，用w表示，在过渡带里，幅频特性单调下降。在通带和阻带内的衰减一般用数 dB表示，通带内允许最大衰减是s，阻带内允许最小衰减是p。定义分别为：Matlab信号处理工具箱中提供了设计巴特沃思模拟滤波器的函数buttord、buttap和butter，格式如下：用于计算巴特沃思模拟低通滤波器的阶N和3dB截止频率Wc (即本书中的符号)。其中，Wp 和Ws分别是滤波器的通带截止频率和阻止截止频率，单位为rad/s；Rp和Rs分别是通带最大衰减系数和阻带最小衰减系数，单位为dB。用于计算N阶巴特沃思归一化=1模拟低通滤波器系统函数的
6、零、极点和增益因子，返回长度为N的向量z和p分别给出N个零点和极点，G是滤波器增益。得到的滤波器系统函数形式如下：其中，和分别是向量z和p的第k个元素。如果要从零、极点得到系统函数的分子和分母多项式系数向量B和A，可以调用结构转换函数：结构转换后系统函数的形式为：其中，M是向量B的长度，N是向量A的长度，分别是向量B和A的第k个元素。用于计算巴特沃什模拟滤波器系统函数中分子和分母多项式系数向量B和A，其中N和分别是滤波器的阶和3dB截止频率，返回向量B和A中的元素和分别是上面的表示式中的分母和分子系数。ftype缺省时，设计低通滤波器，S缺省时，设计数字滤波器。3.1.3 设计程序如下：1巴特
7、沃思模拟低通滤波器设计指标如下：wp=0.4*pi, ws=0.45*pi, Rp=5,Rs=10;实验所用到的程序如下wp=0.4*pi;ws=0.45*pi;Rp=5;Rs=10;N,wc=buttord(wp,ws,Rp,Rs,s);B,A=butter(N,wc,s);w=0:0.001:6;H,W=freqs(B,A,w);H=20*log10(abs(H);plot(w,H),grid on;xlabel(频率/(rad/s);ylabel(幅度/dB)经matlab 仿真以后得到的滤波器如如下图所示：图1.1 巴特沃思模拟低通滤波器运行后在matlab控制窗口上显示如下结果： N
8、=7 2为了方便观察故设计了不同阶次N的滤波器，即Butterworth模拟原型低通滤波器，其程序代码如下：n=0:0.01:2;for ii=1:4switch iicase 1,N=2;case 2,N=5;case 3,N=10;case 4,N=20;endz,p,k=buttap(N); b,a=zp2tf(z,p,k); H,w=freqs(b,a,n); magH2=(abs(H).2; hold on;plot(w,magH2); endxlabel(w/wc);ylabel(|H(jw)|2);title(Butterworth模拟原型低通滤波器);text(1.5,0.18
9、,n=2) text(1.3,0.08,n=5)text(1.16,0.08,n=10)text(0.93,0.98,n=20)grid on;仿真后结果图形如图：图1.2 Butterworth模拟原型低通滤波器3利用冲激响应不变法，把系统函数为的模拟滤波器变换成等价的数字滤波器，采样间隔T=1S。先经过计算得到：，如此有如下式子Matlab实验参考程序如下:B=1,1;A=1,5,6;T=1;Fs=1/T;Bz,Az=impinvar(B,A,Fs); %用冲激响应不变法将模拟滤波器变换成数字滤波器运行结果：Bz =1.0000 -0.2209Az = 1.0000 -0.1851 0.0
10、067得到的数据和先前计算得到的数据一样，故利用冲激响应不变化法把模拟滤波器设计成数字滤波器正确。3.1.4 实验结果与分析1步骤1 中产生的低通模拟滤波器的系统函数，经计算为：经过仿真后的结果，与上式比照，所得结果进分析一样。此处为N=7的低通模拟滤波器2在Butterworth模拟原型低通滤波器中，可以看到，滤波器的幅频特性随着滤波器阶次N的增加而变得越来越好。3步骤3中，利用冲激响应不变法，把系统函数的模拟滤波器变换成等价的数字滤波器，在matlab仿真软件中验证了冲激响应不变法的正确性。3.2用脉冲响应不变法设计的巴特沃思数字低通滤波器3.2.1设计内容：设计巴特沃思数字低通滤波器和椭
11、圆数字低通滤波器，要求通带边界频率fp=2.1kHZ，通带最大衰减Rp=0.5dB；阻带边界频率fs=8kHZ，阻带最小衰减Rs=30dB,采样频率为Fs=20kHZ。3.2.2设计步骤：用MATLAB进展典型的数字滤波器的设计，一般步骤如下:1将设计指标归一化处理。如果采用双线性变换法，还需进展预畸变。2根据归一化频率，确定最小阶数N和频率参数Wn。可供选用的阶数择函数有:buttord，cheblord，cheb2ord，ellipord等。3运用最小阶数N设计模拟低通滤波器原型。模拟低通滤波器的创建函数有:buttap,cheblap, cheb2ap,ellipap和besselap，
12、这些函数输出的是零极点式形式，还要用zp2tf函数转换成分子分母多项式形式。如果想根据最小阶数直接设计模拟低通滤波器原型，可用butter,chebyl,cheby2,ellip,bessel等函数，只是注意要将函数中的Wn设为1。4根据第2步的频率参数Wn，模拟低通滤波原型转换模拟低通、高通、带通、带阻滤波器，可用函数分别是：lp21p,lp2hp,lp2bp,lp2bs。5运用脉冲响应不变法或双线性变法把模滤波器转数字滤波器，调用的函数是impinvar和bilinear。脉冲响应不变法适用于采样频率大于4倍截止频率的锐截止低通带通滤波器，而双线性变换法适合于相位特性要求不高的各型滤波器。
13、6根据输出的分子分母系数，用tf函数生成Hz的表达式，再用freqz函数验证设计结果。3.2.3 设计程序如下：fp=2100;fs=8000;Fs=20000;Rp=0.5;Rs=30;T=1/Fs; %设计指标W1p=fp/Fs*2;W1s=fs/Fs*2;%求归一化频率N,Wn=buttord(W1p,W1s,Rp,Rs,s);%确定butterworth的最小介数N和频率参数Wnz,p,k=buttap(N); %设计模拟低通原型的零极点增益参数bp,ap=zp2tf(z,p,k); %将零极点增益转换成分子分母参数bs,as=lp2lp(bp,ap,Wn*pi*Fs); %将低通原型
14、转换为模拟低通bz,az=impinvar(bs,as,Fs); %用脉冲响应不变法进展模数变换sys=tf(bz,az,T); %给出传输函数H(Z)H,W=freqz(bz,az,512,Fs); %生成频率响应参数subplot(2,1,1);plot(W,20*log10(abs(H); %绘制幅频响应grid on; %加坐标网格xlabel(频率/Hz); ylabel(振幅/dB);subplot(2,1,2); plot(W,abs(H); grid on; xlabel(频率/Hz); ylabel(振幅/H);运行后的波形如下：图2.1 巴特沃思和椭圆数字低通滤波器N =4bz = 0.0000 0.0999 0.1914 0.0252az= 1.0000 -1.4336 1.0984 -0.4115 0.06273.3基于汉明窗函数的FIR低通滤波器的设计3.3.1设计具体低通滤波器要求：设低通滤波器通带截止频率为，阻带截
