巴特沃斯数字(精选)低通滤波器
巴特沃斯数字低通滤波器要点说明书
目录1.题目........................................................ .................................. .22.要求........................................................ (2)3.设计原理........................................................ . (2)3.1 数字滤波器基本概念......................................................... (2)3.2 数字滤波器工作原理......................................................... (2)3.3 巴特沃斯滤波器设计原理 (2)3.4脉冲响应不法......................................................... .. (4)3.5实验所用MATLAB函数说明 (5)4.设计思路........................................................ .. (6)5、实验内容........................................................ . (6)5.1实验程序......................................................... . (6)5.2实验结果分析......................................................... . (10)6.心得体会........................................................ . (10)7.参考文献........................................................ . (10)一、题目：巴特沃斯数字低通滤波器二、要求：利用脉冲响应不变法设计巴特沃斯数字低通滤波器，通带截止频率100HZ,采样频率1000HZ，通带最大衰减为0.5HZ，阻带最小衰减为10HZ，画出幅频、相频相应相应曲线。
巴特沃斯低通滤波器分析
巴特沃斯低通滤波器分析LT摘要：本篇论文叙述了数字滤波器的基本原理、IIR数字滤波器的设计方法和IIR数字高通滤波器设计在MATLAB上的实现与IIR数字滤波器在实际中的应用。
无限脉冲响应（IIR）数字滤波器是冲激响应函数h(t)包含无限个抽样值的滤波器，一般是按照预定的模拟滤波器的逼近函数来转换成相应的数字滤波器，现有的逼近函数如巴特沃斯、切比雪夫。
设计IIR数字滤波器在工程上常用的有两种：脉冲响应不变法、双线性变换法。
其设计过程都是由模拟滤波器的系统函数H(s)去变换出相应的数字滤波器的系统函数H(z)。
关键词:数字滤波器 MATLAB 巴特沃斯切比雪夫双线性变换法Abstract:Digital filter is described in this paper basic principles, IIR digital filter design method of IIR digital high-pass filter design in MATLAB realization of IIR digital filter and application in practice. Infinite Impulse response (IIR) digital filter is the impulse response h (t) includes unlimited sampling filter, usually according to the analog filter approximating function to convert into digital filters, such as Butterworth, Chebyshev approximation functions. Design of IIR digital filters there are two commonly used in engineering: impulse response method, the bilinear transformation method. The design process are by the analog filter system function h (s) to transform the digital filter of the system function h (z).Key words: Digital filter MATLAB Butterworth ChebyshevBilinear transformation methed目录目录 (7)1.前言................... 错误!未定义书签。
数字低通巴特沃斯滤波器的设计实验报告
实验报告姓名：学号：实验日期：实验题目：数字低通巴特沃斯滤波器的设计实验目的：掌握IIR数字滤波器的设计方法实验内容：1.设计一个低通巴特沃斯模拟滤波器，绘制滤波器的的幅频响应及零极点图。
指标如下：通带截止频率：WP＝1000HZ, 通带最大衰减：RP=3dB阻带截止频率：Ws＝2000HZ, 阻带最小衰减：Rs=40 dB参考程序butter1.m2. 用冲激响应不变法和双线性变换法将一模拟低通滤波器转换为数字低通滤波器并图释H(S)和H(Z)，采样频率Fs＝1000Hz实验地点：4305机房实验结果：%巴特沃兹滤波器的幅频响应图subplot(1,2,1)；％分两个窗口，幅频图在第一个窗口wp=1000;ws=2000;rp=3;rs=40; ％设置指标[N,wn]=buttord(wp,ws,rp,rs,'s') ％计算巴特沃斯低通滤波器的阶数和3dB截止频率[B,A]=butter(N,wn,'s'); ％代入N和Wn设计巴特沃斯模拟低通滤波器[Z,P,K]=buttap(N); ％计算滤波器的零、极点[h,w]=freqs(B,A,1024); ％计算1024点模拟滤波器频率响应h，和对应的频率点w ％画频率响应幅度图plot(w,20*log10(abs(h)/abs(h(1))))grid;xlabel('频率Hz');ylabel('幅度（dB）');％给x轴和y轴加标注title('巴特沃斯幅频响应') ％给图形加标题axis([0,3000, -40,3]);line([0,2000],[-3,-3]);line([1000,1000],[-40,3]);%绘制巴特沃斯滤波器的极点图subplot(1,2,2) ％在第二个窗口画极点图p=P';q=Z';x=max(abs([p,q]));x=x+0.1;y=x;axis([-x,x,-y,y]);axis('square')plot([-x,x],[0,0]);hold on plot([0,0],[-y,y]);hold on plot(real(p),imag(p),'x') 程序运行结果： N = 7wn = 1.0359e+003 Z = []P = -0.2225 + 0.9749i -0.2225 - 0.9749i -0.6235 + 0.7818i -0.6235 - 0.7818i -0.9010 + 0.4339i -0.9010 - 0.4339i -1.0000 K = 1.0000100020003000-40-35-30-25-20-15-10-50频率Hz幅度（d B ）巴特沃斯幅频响应-2-1012-1.5-1-0.50.511.5b=1:a=[1,1000];w=[0:1000]*2*pi;％模拟频率为2 f ，其中f 取0~1000Hz [h,w]=freqs(b,a,w);％计算模拟滤波器的频率响应 subplot(2,2,1)plot(w/2/pi,abs(h)/abs(h(1)));grid;％画模拟滤波器幅频特性 title('模拟频率响应');xlabel('f(Hz)');ylabel('幅度'); Fs=1000;[bz,az]=impinvar(b,a,Fs); ％冲激响应不变法设计数字滤波器 [bzl,azl]=bilinear(b,a,Fs); ％双线性变换法设计数字滤波器 wz=[0:pi/512:pi]; hz1=freqz(bz,az,wz); hz2=freqz(bzl,azl,wz);subplot(2,2,2);plot(wz/pi,abs(hz1)/hz1(1));grid; ％画出冲激响应不变法滤波器的幅频图， axis([0,1,0,1]) ％数字频率wz 归一化为0~1 title('冲激响应不变法数字频率响应')subplot(2,2,3);plot(wz/pi,abs(hz2)/hz2(1));grid; ％画出双线性变换法滤波器的幅频图， axis([0,1,0,1]) ％数字频率wz 归一化为0~1 title(双线性变换法数字频率响应);500100000.51模拟频率响应f(Hz)幅度00.510.51冲激响应不变法数字频率响应00.510.51双线性变换法数字频率响应结果分析：总结：。
巴特沃斯低通滤波器的设计
巴特沃斯低通滤波器的设计1、巴特沃斯滤波器的介绍巴特沃斯低通滤波器的幅度平方函数定义为2221|()|1NH j C λλ=+其中C 为一常数参数，N 为滤波器阶数，λ为归一化低通截止频率,/p λ=ΩΩ。
式中N 为整数，是滤波器的阶次。
巴特沃斯低通滤波器在通带内具有最大平坦的振幅特性，这就是说，N 阶低通滤波器在0Ω=处幅度平方函数的前2N-1阶导数等于零，在阻带内的逼近是单调变化的。
巴特沃斯低通滤波器的振幅特性如图a 所示。
滤波器的特性完全由其阶数N 决定。
当N 增加时，滤波器的特性曲线变得更陡峭，这时虽然由a 式决定了在p Ω=Ω处的幅度函数总是衰减3dB ，但是它们将在通带的更大范围内接近于1，在阻带内更迅速的接近于零，因而振幅特性更接近于理想的矩形频率特性。
滤波器的振幅特性对参数N 的依赖关系如图a 所示。
设归一化巴特沃斯低通滤波器的归一化频率为λ，归一化传递函数为()H p ，其中p j λ=，则可得：2221()1(1)N NpjH j C pλλ==+- 由于p图a 巴特沃斯低通滤波器的振幅特性221()()()1()a a jsNcH s H s A s j Ω=--=Ω=+Ω所以巴特沃斯滤波器属于全极点滤波器。
2、常用设计巴特沃斯低通滤波器指标p λ：通带截止频率； p α：通带衰减，单位：dB ；s λ：阻带起始频率；s α：阻带衰减，单位：dB 。
说明：（1）衰减在这里以分贝（dB ）为单位；即222110lg10lg 1()NC H j αλλ⎡⎤==+⎣⎦（2）当3dB α=时p C Ω=Ω为通常意义上的截止频率。
（3）在滤波器设计中常选用归一化的频率/C λ=ΩΩ，即1,p sp s ppλλΩΩ===ΩΩ图b 为巴特沃斯低通滤波器指标3、设计巴特沃斯低通滤波器的方法如下：（1）计算归一化频率1p p pλΩ==Ω，ss pλΩ=Ω。
（2） 根据设计要求按照210101pC α=-和lg lg saN λ=其中a =特沃斯滤波器的参数C 和阶次N ；注意当3p dB α=时 C=1。
巴特沃斯低通滤波器传递函数
巴特沃斯低通滤波器传递函数一、引言巴特沃斯滤波器是一种常见的滤波器，它可以用于信号处理、图像处理等领域。
其中，低通滤波器是最基本的一种。
本文将详细介绍巴特沃斯低通滤波器传递函数的计算方法。
二、巴特沃斯低通滤波器1. 巴特沃斯低通滤波器概述巴特沃斯低通滤波器是一种对频率响应有要求的低通滤波器，其传递函数为：H(s) = 1 / (1 + (s/wc)^2n)^0.5其中，s为Laplace变换中的复频率变量，wc为截止频率，n为阶数。
2. 巴特沃斯低通滤波器传递函数推导（1）将传递函数H(s)转化为标准形式：H(s) = 1 / (1 + (s/wc)^2n)^0.5= 1 / [(s/wc)^2n + 1]^0.5= 1 / [(s^2n + wc^2n) / wc^2n]^0.5= wc^n / [(s^2n + wc^2n)^0.5]（2）将复平面上的频率变量s转化为极坐标形式：s = σ + jω= r * e^(jθ)其中，σ为实部，ω为虚部，r为模值，θ为相位角。
（3）将传递函数H(s)中的s用极坐标表示：H(s) = wc^n / [(s^2n + wc^2n)^0.5]= wc^n / [(r^2n * e^(j2nθ) + wc^2n)^0.5]（4）将传递函数H(s)中的分母进行有理化：H(s) = wc^n / [(r^2n * e^(j2nθ) + wc^2n)^0.5] = wc^n * (r^2n * e^(j2nθ) - wc^2n)^-0.5（5）将传递函数H(s)中的极坐标形式转化为直角坐标形式：H(s) = wc^n * cos(nθ) - jwc^n * sin(nθ)----------------------------------(r^2n - wc^2n)^0.5（6）根据频率响应要求，令模值等于1时的频率为截止频率wc，则有：1 = |H(jwc)| = wc^n / (wc^2n - wc^2n)^0.5=> 1 = (wc/wc)^n=> n = 1 / [ln(1/√R)] / [ln(tan(π/4 + fc/fs/2))]其中，R为通带最大衰减，fc为通带截止频率，fs为采样频率。
巴特沃斯低通滤波器的设计精编资料
巴特沃斯低通滤波器的设计巴特沃斯低通滤波器的设计1、巴特沃斯滤波器的介绍巴特沃斯低通滤波器的幅度平方函数定义为2221|()|1NH j C λλ=+其中C 为一常数参数，N 为滤波器阶数，λ为归一化低通截止频率,/p λ=ΩΩ。
式中N 为整数，是滤波器的阶次。
巴特沃斯低通滤波器在通带内具有最大平坦的振幅特性，这就是说，N 阶低通滤波器在0Ω=处幅度平方函数的前2N-1阶导数等于零，在阻带内的逼近是单调变化的。
巴特沃斯低通滤波器的振幅特性如图a 所示。
滤波器的特性完全由其阶数N 决定。
当N 增加时，滤波器的特性曲线变得更陡峭，这时虽然由a 式决定了在p Ω=Ω处的幅度函数总是衰减3dB ，但是它们将在通带的更大范围内接近于1，在阻带内更迅速的接近于零，因而振幅特性更接近于理想的矩形频率特性。
滤波器的振幅特性对参数N 的依赖关系如图a 所示。
设归一化巴特沃斯低通滤波器的归一化频率为λ，归一化传递函数为()H p ，其中p j λ=，则可得：2221()1(1)N Np jH j C pλλ==+-p 图a 巴特沃斯低通滤波器的振幅特性由于221()()()1()a a jsNcH s H s AsjΩ=--=Ω=+Ω所以巴特沃斯滤波器属于全极点滤波器。
2、常用设计巴特沃斯低通滤波器指标pλ：通带截止频率；pα：通带衰减，单位：dB；sλ：阻带起始频率；sα：阻带衰减，单位：dB。
说明：（1）衰减在这里以分贝（dB）为单位；即222110lg10lg1()NCH jαλλ⎡⎤==+⎣⎦（2）当3dBα=时p CΩ=Ω为通常意义上的截止频率。
（3）在滤波器设计中常选用归一化的频率/Cλ=ΩΩ，即1,p sp sp pλλΩΩ===ΩΩ图b 为巴特沃斯低通滤波器指标3、设计巴特沃斯低通滤波器的方法如下：（1）计算归一化频率1p p pλΩ==Ω，ss pλΩ=Ω。
（2） 根据设计要求按照210101pC α=-和lg lg saN λ=其中a =特沃斯滤波器的参数C 和阶次N ；注意当3p dB α=时 C=1。
巴特沃斯低通滤波器
带最小衰减α =30dB，按照以上技术指标设计巴特沃斯低通滤波器。 0.1a s
1a p
1a s
2.4
0242 4.25, 2.4
2.4 10 1 2 f lg 0.0242 lg 0.0242 NN 2 4.25, 55 lgf 2.4 4.25, N N s sp lg 2.4 2.4 2 f p
H（ a s）
N c
(s s
k 0
N 1
k
)
7 j 3
• 例如N=3, 通过下式可以计算出6个极点 5 2 4 j j j j s 3 c 3 s 2 c 3 s 0 c 3 s1 c
s 4 c
j2
s 5 c
要求
f i g u r e ; p l o t ( Q , H a s ) ; a x i s ( [ 0 5]);xlabel('f(kHz)'),ylabel('20lg(abs(H_{a}(j{\Omega})))(dB)');
3 0
- 7 0
• • • • •
L=length(Ha); Yt=Xt(1:L).*Ha; figure;plot(Q,abs(Yt));axis([0 60 0 150]); yt=ifft(Yt); figure;plot(Q,yt);
• 模拟低通滤波器的设计指标 • 构造一个逼近设计指标的传输函数Ha(s) • Butterworth（巴特沃斯）低通逼近
模拟低通滤波器的设计指标及逼近方法(续)
• 模拟低通滤波器的设计指标有αp, Ωp,αs和Ωs。 • Ωp；通带截止频率 • Ωs：阻带截止频率
巴特沃斯数字低通滤波器要点说明
目录1.题目........................................................ .................................. .22.要求........................................................ (2)3.设计原理........................................................ . (2)数字滤波器基本概念......................................................... (2)数字滤波器工作原理......................................................... (2)巴特沃斯滤波器设计原理 (2)脉冲响应不法......................................................... .. (4)实验所用MATLAB函数说明 (5)4.设计思路........................................................ .. (6)5、实验内容........................................................ . (6)实验程序......................................................... . (6)实验结果分析......................................................... . (10)6.心得体会........................................................ . (10)7.参考文献........................................................ . (10)一、题目：巴特沃斯数字低通滤波器二、要求：利用脉冲响应不变法设计巴特沃斯数字低通滤波器，通带截止频率100HZ,采样频率1000HZ，通带最大衰减为，阻带最小衰减为10HZ，画出幅频、相频相应相应曲线。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
目录1.题目...................................................................... (2)2.要求...................................................................... . (2)3.设计原理...................................................................... .. (2)3.1数字滤波器基本概念 (2)3.2数字滤波器工作原理 (2)3.3巴特沃斯滤波器设计原理 (2)3.4脉冲响应不法...................................................................... . (4)3.5实验所用MATLAB函数说明 (5)4.设计思路...................................................................... (6)5、实验内容...................................................................... .. (6)5.1实验程序...................................................................... (6)5.2实验结果分析...................................................................... (10)6.心得体会...................................................................... .. (10)7.参考文献...................................................................... .. (10)一、题目：巴特沃斯数字低通滤波器二、要求：利用脉冲响应不变法设计巴特沃斯数字低通滤波器，通带截止频率100HZ,采样频率1000HZ，通带最大衰减为0.5HZ，阻带最小衰减为10HZ，画出幅频、相频相应相应曲线。
并假设一个信号x(t)=sin(2*pi*f1*t)+sin(2*pi*f2*t),其中f1=50HZ,f2=200HZ。
用此信号验证滤波器设计的正确性。
三、设计原理1、数字滤波器的基本概念所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤波器除某些频率成分的数字器件或程序，因此，数字滤波的概念和模拟滤波相同，只是的形式和实现滤波方法不同。
正因为数字滤波通过数值运算实现滤波，所以数字滤波处理精度高、稳定、体积小、质量轻、灵活、不存在阻抗匹配问题，可以实验模拟滤波器无法实现的特殊滤波功能。
如果要处理的是模拟信号，可通过A\DC 和D\AC,在信号形式上进行匹配转换，同样可以使用数字滤波器对模拟信号进行滤波。
2、数字滤波器的工作原理数字滤波器是一个离散时间系统，输入x(n)是一个时间序列，输出y(n)也是一个时间序列。
如数字滤波器的系统函数为H(Z),其脉冲响应为h(n),则在时间域内存在下列关系y(n)=x(n)h(n)在Z 域内，输入输出存在下列关系Y(Z)=H(Z)X(Z)式中，X(Z),Y(Z)分别为输入x(n)和输出y(n)的Z 变换。
同样在频率域内，输入和输出存在下列关系Y(jw)=X(jw)H(jw)式中，H(jw)为数字滤波器的频率特性，X(jw)和Y(jw)分别为x(n)和y(n)的频谱。
w 为数字角频率，单位rad 。
通常设计H(jw)在某些频段的响应值为1，在某些频段的响应为0.X(jw)和H(jw)的乘积在频率响应为1的那些频段的值仍为X(jw)，即在这些频段的振幅可以无阻碍地通过滤波器，这些频带为通带。
X(jw)和H(jw)的乘积在频段响应为0的那些频段的值不管X(jw)大小如何均为零，即在这些频段里的振幅不能通过滤波器，这些频带称为阻带。
一个合适的数字滤波器系统函数H(Z)可以根据需要输入x(n)的频率特性，经数字滤波器处理后的信号y(n)保留信号x(n)中的有用频率成分，去除无用频率成分。
3、巴特沃斯滤波器设计原理（1）基本性质巴特沃斯滤波器以巴特沃斯函数来近似滤波器的系统函数。
巴特沃斯滤波器是根据幅频特性在通频带内具有最平坦特性定义的滤波器。
巴特沃思滤波器的低通模平方函数表示1()ΩΩ+=Ωc N /22a 11)(j H N=1,2，……（2-6） 下面归纳了巴特沃斯滤波器的主要特征a 对所有的N ，()1aj H 20=Ω=Ω。
b 对所有的N ，()707.0aj 2c =ΩΩH =Ω即()dB 3a lg 20j H c =Ω=ΩΩ c ()Ωj H a 2是Ω的单调下降函数。
d ()Ωj H a 2随着阶次N 的增大而更接近于理想低通滤波器。
如下图2所示，可以看出滤波器的幅频特性随着滤波器阶次N 的增加而变得越来越好，在截止频率Ωc处的函数值始终为1/2的情况下，通带内有更多的频带区的值接近于1；在阻带内更迅速的趋近于零。
图2巴特沃思低通滤波平方幅频特性函数（2）系统函数设巴特沃斯的系统函数为Ha（s），则：（3）设计过程巴特沃思低通滤波技术指标关系式为a p >-20log|Ha(jΩ)|，Ω<ΩPa s <-20log|Ha(jΩ)|，Ω>Ωs其中：Ωp 为通带边界频率，Ωs经过化简整理可得：取满足上式的最小整数N作为滤波器的阶数。
再将N代入可得：或查表求得归一化传输函数H(s),令s/Ωc代替归一化原型滤波器系统函数中的s,即得到实际滤波器传输函数。
4、脉冲响应不变法所谓脉冲响应不变法就是数字滤波器的脉冲响应序列h(n)等于模拟滤波器的响应ha(t)的采样值，即h(n)=ha(t)|t=nT=ha(nT)式中，T为采样周期。
因此数字滤波器的系统函数H(Z)可由下式求得H(z)=Z[h(n)]=Z[ha(nT)]Z[-]表示[-]的内容进行变换，变换的内容请参考相应的数字信号处理材料。
如果已经获得了满足性能指标的模拟滤波器的传递函数Ha(s),求与之对应的数字滤波器的传递函数H(z)的方法是：（1）、求模拟滤波器的单位脉冲响应ha(t)。
式中，L[Ha(s)]表示对Ha(s)的Laplace.逆变换。
Laplace变换内容请参考高等数学的积分变换或信号处理教材。
（2）、求模拟滤波器单位冲激响应ha(t)的采样值，即数字滤波器冲激响应序列h(n)。
（3）、对数字滤波器的冲激h(n)响应进行z变换，得到传递函数H(z)。
由上述方法推论出更直接地由模拟滤波器系统函数Ha(s)求出数字滤波器系统函数H(z)的步骤是：（1）利用部分分式展开将模拟滤波器的传递函数H(z)展开成Ha(s)=Rk\(S-Pk)在MATLAB中这步可通过residue函数实现若调用residue函数的形式为[b,a]=residue(R,P,K)形式。
若为[R,P,K]=residue(a,b)则为上面调用形式的反过程。
(2)将模拟极点Pk变换为数字极点e^pkT即得到数字系统的传递函数H(z)=Rk\(1-e^pkT*z*(-1))式中T为采样间隔。
（3）将上式转换为传递函数形式，可采用[R,P,K]=residue(b,a)。
对于上面的步骤，中已经提供了冲激响应不变法设计数字滤波器的函数，调用格式为[bz,az]=impinvar(b,a[,Fs],Fp)式中，b,a为模拟滤波器分子和分母多项式系数向量；Fs为采样频率（所滤波数据），单位Hz，缺省时为1Hz，为预畸变频率（prewarpedfrequency），是一个“匹配”频率，在该频率上，频率响应在变换前后和模拟频率可精确匹配。
一般设计中不考虑。
bz,az分别为数字滤波器分子和分母多项式系数向量。
5、实验所用MATLAB函数说数。
（1）[N,wc]=buttord(wp,ws,RP,As,’s’)该格式用于计算巴特沃斯模拟滤波器的阶数N和3db截止频率wc。
Wp、ws和wc是实际模拟角频率(rad\s)。
Rp和As为通带最大衰减和最小衰减。
（2）[Z,P,k]=buttap(N)该格式用于计算N阶巴特沃斯归一化模拟低通原型滤波器系统函数的零、极点和增益因子，返回长度为N的列向量Z和P，分别给出N个零点和极点的位置，K表示滤波器增益。
（3）Y=filter(b,a,x)式中b表示系统传递函数的分子多项式的系数矩阵；a表示系统传递函数的分母多项式的系数矩阵；x表示输入序列；filter表示输出序列。
IIR函数实现的直接形式。
(4)[b，a]=butter（N，wc，‘ftype’）计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量b、a。
说明：调用参数N和wc分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值，一般是调用buttord格式（1）计算N和wc。
系数b、a是按照z-1的升幂排列。
(5)[B，A]=butter（N，Ωc，‘ftype’，‘s’）计算巴特沃斯模拟滤波器系统函数的分子、分母多项式系数向量。
说明：调用参数N和Ωc分别为巴特沃斯模拟滤波器的阶数和3dB截止频率(实际角频率)，可调用buttord（2）格式计算N和Ωc。
系数B、A按s的正降幂排列。
tfype为滤波器的类型：◇ftype=high时，高通；Ωc只有1个值。
◇ftype=stop时，带阻；Ωc=[Ωcl,Ωcu]，分别为带阻滤波器的通带3dB下截止频率和上截止频率。
◇ftype缺省时：若Ωc只有1个值，则默认为低通；若Ωc有2个值，则默认为带通；其通带频率区间Ωcl<Ω<Ωcu。
(6)[H,w]=freqz(b,a,N)b和a分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了离散系统频响在0～pi范围内N个频率等分点的值（其中N为正整数），w则包含了范围内N个频率等分点。
调用默认的N时，其值是512。
可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，绘制出系统的频响曲线。
(7)lp2lp函数[bt,at]=lp2lp(b,a,w0)该函数用于实现由低通模拟原型滤波器至低通滤波器的频率变换，可以用传递函数和状态空间进行转换，但无论哪种形式，其输入必须是模拟滤波器原型。