数字信号处理实验 matlab版 线性相位FIR数字滤波器
数字信号处理MATLAB实验1
实验一熟悉MATLAB环境一、实验目的(1)熟悉MATLAB的主要操作命令。
(2)学会简单的矩阵输入和数据读写。
(3)掌握简单的绘图命令。
(4)用MATLAB编程并学会创建函数。
(5)观察离散系统的频率响应。
二、实验内容认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。
在熟悉了MATLAB基本命令的基础上，完成以下实验。
上机实验内容：(1)数组的加、减、乘、除和乘方运算。
输入A=[1234]，B=[3456]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。
(2)用MATLAB实现以下序列。
a）x(n)=0.8n0≤n≤15b）x(n)=e(0.2+3j)n0≤n≤15c）x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π)0≤n≤15(n)=x(n+16),绘出四个d）将c)中的x(n)扩展为以16为周期的函数x16周期。
(n)=x(n+10),绘出四个e）将c)中的x(n)扩展为以10为周期的函数x10周期。
(3)x(n)=[1,-1,3,5]，产生并绘出下列序列的样本。
a）x 1(n)=2x(n+2)-x(n-1)-2x(n)b）∑=-=51k 2)k n (nx (n) x (4)绘出下列时间函数的图形，对x轴、y轴以及图形上方均须加上适当的标注。
a)x(t)=sin(2πt)0≤t≤10sb)x(t)=cos(100πt)sin(πt)0≤t≤4s(5)编写函数stepshift(n0,n1,n2)实现u(n-n0),n1<n0<n2,绘出该函数的图形，起点为n1，终点为n2。
(6)给定一因果系统)0.9z 0.67z -1)/(1z 2(1H(z)-2-1-1+++=求出并绘制H(z)的幅频响应与相频响应。
(7)计算序列｛8-2-123｝和序列｛23-1-3｝的离散卷积，并作图表示卷积结果。
Matlab设计FIR数字滤波器
FIR数字滤波器专业：学号：XX：一课题目的：1学会使用Matlab的各项功能。
2学会把自己在课堂上学习的知识运用到实践当中。
3了解利用Matlab设计FIR数字滤波器的基本方法。
4在课程设计的过程中掌握程序编译及软件设计的基本方法。
5提高自己对于新知识的学习能力及进行实际操作的能力。
二课题要求：在信号处理过程中所处理的信号往往混有噪音，从接受到的信号中消除或减弱噪音是信号处理过程中十分重要的问题。
根据有用信号和噪音的不同特性，提取有用信号的过程称为滤波，实现滤波功能的系统称为滤波器。
而数字滤波器又是滤波器中运用极为广泛的一种滤波器。
数值滤波技术是数字信号处理的一个重要组成部分，滤波器的设计是信号处理的核心问题之一。
FIR数字滤波器在保证幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性。
要求通过网络及各种资料解决实际问题设计一个符合要求的FIR数字滤波器。
三课题内容：数字滤波器和模拟滤波器有着相同的滤波概念,根据其频率响应特性可分为低通、高通、带通、带阻等类型。
与模拟滤波器相比,数字滤波器除了具有数字信号处理固有优点外,还有滤波精度高、稳定性好、灵活性强等优点。
在数字信号处理中,由于信号中经常混有各种复杂成分，所以很多信号分析都是基于滤波器而进行的，FIR数字滤波器在数字信号处理中发挥着重要作用，采用Matlab软件对FIR数字滤波器进行仿真设计，简化了设计中繁琐的计算。
设计中采用窗函数法，频率采样法和优化设计方法，通过调用Matlab函数设计FIR数字滤波器。
绘制出滤波器的特性图。
利用所设计的滤波器对多个频带叠加的正弦信号进行处理，对比滤波前后的信号时域和频域图，验证滤波器的效果。
最后录制一段语音信号，并对录制的信号进行采样和加噪，绘制出采样后语音信号的时域波形和频谱图，然后用所设计的滤波器对加噪后的信号进行滤波，绘制出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化。
实验3 用MATLAB窗函数法设计FIR滤波器
实验10 用MATLAB 窗函数法设计FIR 滤波器一、实验目的㈠、学习用MA TLAB 语言窗函数法编写简单的FIR 数字滤波器设计程序。
㈡、实现设计的FIR 数字滤波器，对信号进行实时处理。
二、实验原理㈠、运用窗函数法设计FIR 数字滤波器与IIR 滤波器相比，FIR 滤波器在保证幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性。
设FIR 滤波器单位脉冲响应)(n h 长度为N ，其系统函数)(z H 为∑-=-=10)()(N n n zn h z H)(z H 是1-z 的)1(-N 次多项式，它在z 平面上有)1(-N 个零点，原点0=z 是)1(-N 阶重极点。
因此，)(z H 永远是稳定的。
稳定和线性相位特性是FIR 滤波器突出的优点。
FIR 滤波器的设计任务是选择有限长度的)(n h ，使传输函数)(ωj e H 满足技术要求。
主要设计方法有窗函数法、频率采样法和切比雪夫等波纹逼近法。
本实验主要介绍用窗函数法设计FIR 数字滤波器。
图7-10-1 例1 带通FIR 滤波器特性㈡、 用MATLAB 语言设计FIR 数字滤波器例1：设计一个24阶FIR 带通滤波器，通带为0.35<ω<0.65。
其程序如下b=fir1(48,[0.35 0.65]);freqz(b,1,512)可得到如图7-10-1 所示的带通FIR滤波器特性。
由程序可知，该滤波器采用了缺省的Hamming窗。
例2：设计一个34阶的高通FIR滤波器，截止频率为0.48,并使用具有30dB波纹的Chebyshev窗。
其程序如下Window=chebwin(35,30);b=fir1(34,0.48,'high',Window);freqz(b,1,512)可得到如图7-10-2 所示的高通FIR滤波器特性。
图7-10-2 例2 高通FIR滤波器特性例3：设计一个30阶的低通FIR滤波器，使之与期望频率特性相近，其程序如下 f=[0 0.6 0.6 1];m=[1 1 0 0];b=fir2(30,f,m);[h,w]=freqz(b,1,128);plot(f,m,w/pi,abs(h))结果如图7-10-3所示。
FIR数字滤波器的设计与matlab实现
W
(e
j
)
{0.5WR
()
0.25[WR
(
2 )
N 1
WR
(
2
)]}e N 1
j(
N 1) 2
W ()e ja
W
()
0.5WR
()
0.25[WR
(
2 )
N 1
WR
(
2 )]
N 1
(6-2-12b)
W () 是三项矩形窗的幅度响应 WR() 的移位加权和，
n0
即
H (z) z(N1)H (z1)
则有
H (z) 1 [H (z) z(N1)H (z1)] 2
1
N 1
h(n)[ z n
z (N 1) z n ]
2 n0
z ( N 1) / 2
N 1
h(n)
1
(n N 1)
[z 2
(n N 1)
❖ 滤波器设计任务的中心就是求得系统函数。 数字滤波器的系统函数最主要的特征有三个： 幅度平方响应、相位响应和群延迟。
❖ IIR滤波器可以用较少的阶数获得较好的幅度响应， 但由于其结构存在反馈，可能造成系统的不稳定， 其优异幅度响应一般是以相位的非线性为代价的， 非线性相位会引起频率色散。
❖ FIR系统的最主要特性之一就是可以构成具有线性 相位特性的滤波器。所谓线性相位特性是指滤波 器对不同频率的正弦波所产生的相移和正弦波的 频率成直线关系。因此，在滤波器通带内的信号 通过滤波器后，除了由相频特性的斜率决定的延 迟外，可以不失真地保留通带以内的全部信号。
实验四 参考 基于MATLAB的FIR数字滤波器设计
实验4 基于MATLAB 的FIR 数字滤波器设计一、 实验目的（1） 加深对数字滤波器的常用指标理解。
（2） 学习数字滤波器的设计方法。
二、 实验原理 低通滤波器1()1()P P P SS H H δδδπ⎧-≤Ω≤+Ω≤Ω⎪⎨Ω≤Ω≤Ω≤⎪⎩低通滤波器的常用指标：1、 通带边缘频率P Ω，2、 阻带边缘频率SΩ ，3、 通带起伏P δ，4、通带峰值起伏])[1(log2010dB p pδα--=，5、阻带起伏sδ，最小阻带衰减])[(log 2010dB s S δα-=。
三、 要求 （1）在MA TLAB 中，熟悉函数fir1、kaiserord 、remezord 、remez 的使用；B = fir1(n,Wn,'high','noscale')设计滤波器；[n,Wn,beta,ftype] = kaiserord(f,a,dev)估计滤波器阶数；[n,fo,ao,w] = remezord (f,a,dev ,fs)计算等波纹滤波器阶数n 和加权函数w(ω)； B=remez(n,f,a)进行等波纹滤波器的设计（2）阅读附录中的实例，学习FIR 滤波器的设计方法及其在MA TLAB 中的实现；（3）给出FIR 数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。
四、 实验内容利用MA TLAB 编程，分别用窗函数法和等波纹滤波器法设计两种FIR 数字滤波器，指标要求如下：通带边缘频率：π45.01=ΩP ，π65.02=ΩP ，通带峰值起伏：][1dB p≤α。
阻带边缘频率：π3.01=ΩS ，π75.02=ΩS ，最小阻带衰减：][40dB S ≥α。
附录：例1 用凯塞窗设计一FIR 低通滤波器，通带边界频率π3.0=Ωp ，阻带边界频率π5.0=Ωs ，阻带衰减 不小于50dB 。
S Ω-PΩ-P SΩPassband StopbandTransitionbandFig 1 Typical magnitudespecification for a digital LPF解首先由过渡带宽和阻带衰减来决定凯塞窗的N和π2.0=Ω-Ω=∆Ωps，，上图给出了以上设计的频率特性，(a) 为N=30直接截取的频率特性(b)为凯塞窗设计的频率特性。
数字信号处理实验报告-FIR滤波器的设计与实现
数字信号处理实验报告-FIR滤波器的设计与实现在数字信号处理中，滤波技术被广泛应用于时域处理和频率域处理中，其作用是将设计信号减弱或抑制被一些不需要的信号。
根据滤波器的非线性抑制特性，基于FIR（Finite Impulse Response）滤波器的优点是稳定，易设计，可以得到较强的抑制滤波效果。
本实验分别通过MATLAB编程设计、实现、仿真以及分析了一阶低通滤波器和平坦通带滤波器。
实验步骤：第一步：设计一阶低通滤波器，通过此滤波器对波型进行滤波处理，分析其对各种频率成分的抑制效果。
为此，采用零极点线性相关算法设计滤波器，根据低通滤波器的特性，设计的低通滤波器的阶次为n=10，截止频率为0.2π，可以使设计的滤波器被称为一阶低通滤波器。
第二步：设计平坦通带滤波器。
仿真证明，采用兩個FIR濾波器組合而成的阻礙-提升系統可以實現自定義的總三值響應的設計，得到了自定義的總三值響應函數。
实验结果：1、通过MATLAB编程，设计完成了一阶低通滤波器，并通过实验仿真得到了一阶低通滤波器的频率响应曲线，证明了设计的滤波器具有良好的低通性能，截止频率为0.2π。
在该频率以下，可以有效抑制波形上的噪声。
2、设计完成平坦通带滤波器，同样分析其频率响应曲线。
从实验结果可以看出，此滤波器在此频率段内的通带性能良好，通带范围内的信号透过滤波器后，损耗较小，滞后较小，可以满足各种实际要求。
结论：本实验经过实验操作，设计的一阶低通滤波器和平坦通带滤波器具有良好的滤波特性，均已达到预期的设计目标，证明了利用非线性抑制特性实现FIR滤波处理具有较强的抑制滤波效果。
本实验既有助于深入理解FIR滤波器的设计原理，也为其他应用系统的设计和开发提供了指导，进而提高信号的处理水平和质量。
基于matlab的fir数字滤波器的设计
一、引言数字滤波器是数字信号处理中至关重要的组成部分，它能够对数字信号进行滤波处理，去除噪音和干扰，提取信号中的有效信息。
其中，fir数字滤波器作为一种常见的数字滤波器类型，具有稳定性强、相位响应线性等特点，在数字信号处理领域得到了广泛的应用。
本文将基于matlab软件，探讨fir数字滤波器的设计原理、方法和实现过程，以期能够全面、系统地了解fir数字滤波器的设计流程。
二、fir数字滤波器的基本原理fir数字滤波器是一种有限长冲激响应（finite impulse response, FIR）的数字滤波器，其基本原理是利用线性相位特性的滤波器来实现对数字信号的筛选和处理。
fir数字滤波器的表达式为：$$y(n) = \sum_{k=0}^{M}h(k)x(n-k)$$其中，y(n)为输出信号，x(n)为输入信号，h(k)为滤波器的系数，M为滤波器的长度。
fir数字滤波器的频率响应特性由其系数h(k)决定，通过设计合适的系数，可以实现对不同频率成分的滤波效果。
三、fir数字滤波器的设计方法fir数字滤波器的设计方法主要包括窗函数法、频率抽样法、最小最大法等。
在matlab中，可以通过信号处理工具箱提供的fir1函数和firls函数等来实现fir数字滤波器的设计。
下面将分别介绍这两种设计方法的基本原理及实现步骤。
1. 窗函数法窗函数法是fir数字滤波器设计中最为常见的方法之一，其基本原理是通过对理想滤波器的频率响应进行窗函数加权来满足设计要求。
在matlab中，可以使用fir1函数实现fir数字滤波器的设计，其调用格式为：h = fir1(N, Wn, type)其中，N为滤波器的阶数，Wn为滤波器的截止频率，type为窗函数的类型。
通过调用fir1函数，可以灵活地设计出满足特定要求的fir数字滤波器。
2. 频率抽样法频率抽样法是fir数字滤波器设计中的另一种重要方法，其基本原理是在频域上对理想滤波器的频率响应进行抽样，并拟合出一个最优的滤波器。
数字信号处理 实验八 FIR数字滤波器的MATLAB实现
实验八 FIR 数字滤波器的MATLAB 实现一、 实验目的：学习用窗函数法设计FIR 数字滤波器的原理及其设计步骤；学习编写数字滤波器的设计程序的方法，并能进行正确编程；根据给定的滤波器指标，给出设计步骤。
二、 实验仪器：电脑一台，MA TLAB6.5或更高级版本软件一套。
三、 实验内容：（一） 实验原理及实例分析 1. FIR 低通数字滤波器的设计原理如果系统的冲激响应)(n h d 为已知，则系统的输入/输出关系为：)()()(n h n x n y d ⋅=对于低通滤波器，只要设计出低通滤波器的冲激响应函数，就可以由上式得到系统的输出了。
假设所希望的数字滤波器的频率响应为)(jwd e H ，它是频域的周期函数，周期为2π，那么它与)(jwd eH 相对应的傅立叶系数为dw eeH n h jnwjwd d )(21)(⎰-=πππ以)(n h d 为冲激响应的数字滤波器将具有频域响)(jwd e H 。
但是将)(n h d 作为滤波器脉冲响应有两个问题：（1） 它是无限长的，与FIP 滤波器脉冲响应有限长这一前提不一致 （2） 它是非因果的，0,0)(<≠n n h d 对此，要采取以下的措施，（1） 将)(n h d 截短 （2）将其往右平移，由此得到)(2n h 的实际频域响应jnwN n jwd en heH ∑-==12)()(，与理想频域响应)(jwd eH 相近，但不完全一致。
理论证明上述现象是对)(n h d 进行简单截短处理的必然结果，一般称为吉布斯现象，为尽可能的减少吉布斯现象，应对)(n h d 进行加窗截取，即以)()()(n W n h n h N d ⋅=作为FIR 滤波器的系数。
常用的窗函数有矩形窗、海明窗和布莱克曼窗等。
2. 用窗函数法设计FIR 滤波器MA TLAB 设计FIR 滤波器有多种方法和对应的函数，见表8－1。
表8－1 MA TLAB 设计FIR 滤波器的方法和函数方法描述函数窗方法 使用窗函数和逆傅立叶变换实现 fir1,fir2,kaiserord 等 多带方法 包含子带频率域firls ，remez 等 最小二乘法 使用最小二乘法将整个频率域上的错误几率压缩到最小fircls ，fircls1等 任意响应法 使用任意响应，包括非线性相位以及复滤波器cremez 等 余弦法使用三角函数的低通响应firrcos 等窗函数方法不仅在数字滤波器的设计中占有重要的地位，同时可以用于功率谱的估计，从根本上讲，使用窗函数的目的就是消除由无限序列的截短而引起的Gibbs 现象所带来的影响。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验23 线性相位FIR数字滤波器(完美格式版，本人自己完成，所有语句正确，不排除极个别错误，特别适用于山大，勿用冰点等工具下载，否则下载之后的word格式会让很多部分格式错误，谢谢)XXXX学号姓名处XXXX一、实验目的1 加深对线性相位FIR数字滤波器特性的理解。
2 掌握线性相位滤波器符幅特性和零极点分布的研究方法。
3 了解用MATLAB研究线性相位滤波器特性时程序编写的思路和方法。
二、实验内容1 线性相位FIR滤波器的特性2 第一类线性相位滤波器(类型Ⅰ)3 第二类线性相位滤波器(类型Ⅱ)4 第三类线性相位滤波器(类型Ⅲ)5 第四类线性相位滤波器(类型Ⅳ)6 线性相位FIR数字滤波器零点分布特点三、实验环境MATLAB7.0四、实验原理1.线性相位FIR滤波器的特性与IIR滤波器相比，FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性。
设FIR滤波器单位脉冲响应h(n)长度为N，其系统函数为∑-=-=1Nnnzh(n)H(z)当滤波器的系数N满足一定的对称条件时，就可以获得线性相位。
线性相位FIR滤波器共分为四种类型，分别为：(1)类型Ⅰ，系数对称，即h(n)＝h(N－1－n)，N为奇数。
(2)类型Ⅱ，系数对称，即h(n)＝h(N－1－n)，N为偶数。
(3)类型Ⅲ，系数反对称，即h(n)＝－h(N－1－n)，N为奇数。
(4)类型Ⅳ，系数反对称，即h(n)＝－h(N－1－n)，N为偶数。
对于上述四类线性相位FIR滤波器，参考文献［1］中提供了一段通用程序，对考虑正负号的幅度频率特性(简称符幅特性)进行求解，程序名为amplres.m，程序如下：function[A,w,type,tao]=amplres(h)N=length(h);tao=(N-1)/2;L=floor((N-1)/2);n=1:L+1;w=[0:500]*2*pi/500;if all(abs(h(n)-h(N-n+1))<1e-10)A=2*h(n)*cos(((N+1)/2-n)'*w)-mod(N,2)*h(L+1); type=2-mod(N,2);elseif all(abs(h(n)+h(N-n+1))<1e-10)&(h(L+1)*mod(N,2)==0) A=2*h(n)*sin(((N+1)/2-n)'*w); type=4-mod(N,2);elseerror('错误:这不是线性相位滤波器!') end另外，FIR 滤波器的H(z)是z －1的(N －1)次多项式，它在z 平面上有(N －1)个零点，原点z ＝0是(N －1)阶重极点。
由于线性相位FIR 滤波器的极点都在原点处，因此不存在稳定性的问题。
但有必要研究零点对滤波器特性的影响。
下面分别研究和分析上述问题。
2.第一类线性相位滤波器(类型Ⅰ)例23-1 已知FIR 线性相位系统h ＝［3，－1，－5，4，6，4，－5，－1，3］，要求描绘系统的冲激响应和符幅特性。
解 程序如下：h=[3,-1,-5,4,6,4,-5,-1,3]; M=length(h);n=0:M-1;[A,w,type,tao]=amplres(h);type subplot(2,1,1),stem(n,h); ylabel('h(n)');xlabel('n'); subplot(2,1,2),plot(w/pi,A); ylabel('A');xlabel('＼pi');MATLAB 命令窗显示：type = 1 由图23-1可见，这是一个第一类线性相位滤波器。
滤波器的系数N 为奇数(该题N ＝9)，且h(n)＝h(N －1－n)，幅度特性关于w ＝p 对称，在w ＝0和w ＝p 处可以取任何值；可以用于实现低通、高通、带通、带阻等各种滤波特性。
02468-10010h (n )n0.511.52-2002040A＼pi图23-1 类型Ⅰ滤波器冲激响应和符幅特性3.第二类线性相位滤波器(类型Ⅱ)例23-2 已知FIR 线性相位系统h ＝［3， －1，－5，4，4，－5，－1，3］，要求描绘系统的冲激响应和符幅特性。
解 程序如下：h=[3,-1,-5,4,4,-5,-1,3];M=length(h);n=0:M-1;[A,w,type,tao]=amplres(h);type subplot(2,1,1),stem(n,h); ylabel('h(n)');xlabel('n'); subplot(2,1,2),plot(w/pi,A); ylabel('A');xlabel('＼pi');MATLAB 命令窗显示： type = 2 由图23-2可见，这是一个第二类线性相位滤波器。
滤波器的系数N 为偶数(该题N ＝8)，且h(n)＝h(N －1－n)，幅度特性关于w ＝p 反对称，在w ＝0处可以取任何值，在w ＝p 处必定为0；不能用于实现高通和带阻滤波器。
01234567-55h (n )n00.51 1.52-5050A＼pi图23-2 类型Ⅱ滤波器冲激响应和符幅特性4.第三类线性相位滤波器(类型Ⅲ)例23-3 已知FIR 线性相位系统h ＝［3，－1，－5，4，0，－4，5，1，－3］，要求描绘系统的冲激响应和符幅特性。
解 程序如下：h=[3,-1,-5,4,0,-4,5,1,-3]; M=length(h);n=0:M-1;[A,w,type,tao]=amplres(h);type subplot(2,1,1),stem(n,h); ylabel('h(n)');xlabel('n'); subplot(2,1,2),plot(w/pi,A); ylabel('A');xlabel('＼pi');12345678-505h (n)n00.51 1.52-50050A＼pi图23-3 类型Ⅲ滤波器冲激响应和符幅特性MATLAB 命令窗显示： type = 3 由图23-3可见，这是一个第三类线性相位滤波器。
滤波器的系数N 为奇数(该题N ＝9)，且h(n)＝－h(N －1－n)，幅度特性关于w ＝p 反对称，在w ＝0和w ＝p 处都必定为0；只能实现带通滤波器。
5.第四类线性相位滤波器(类型Ⅳ)例23-4 已知FIR 线性相位系统h ＝［3，－1，－5，4，－4，5，1，－3］，要求描绘系统的冲激响应和符幅特性。
解 程序如下：h=[3,-1,-5,4,-4,5,1,-3]; M=length(h);n=0:M-1;[A,w,type,tao]=amplres(h);type subplot(2,1,1),stem(n,h); ylabel('h(n)');xlabel('n'); subplot(2,1,2),plot(w/pi,A); ylabel('A');xlabel('＼pi');MATLAB 命令窗显示： type = 4 由图23-4可见，这是一个第四类线性相位滤波器。
滤波器的系数N 为偶数(该题N ＝8)，且h(n)＝－h(N －1－n)，幅度特性关于w ＝p 对称，在w ＝0必定为0，w ＝p 处可以取任何值；不能用于实现低通、带阻滤波器。
02468-505h (n )n00.51 1.52-2020A＼pi图23-4 类型Ⅳ滤波器冲激响应和符幅特性6.线性相位FIR 数字滤波器零点分布特点线性相位FIR 滤波器的零点分为三种情况： (1)在zi ＝1或zi ＝－1位的零点单独出现；(2)实数零点或在单位圆上的复数零点按z1与1z 1成对出现；(3)非单位圆上的复数零点，则按z1、z*1、1z 1、*1z 1⎪⎪⎭⎫⎝⎛ 四个点同时出现。
例23-5 求例23-1线性相位系统的零极点分布图，并证明实数零点和复数零点成组出现的特点。
解 程序如下：h=[3,-1,-5,4,6,4,-5,-1,3]; M=length(h);n=0:M-1; rz=roots(h) for i=1:8 r(i)=1/rz(i); end r'zplane(h,1)程序执行结果如下和图23-5所示。
可以看出，序号为1、2、6、7的四个复数为一组；在实轴上的3、8两个实数为一组；处于单位圆上的4、5为一组。
rz = 1.1674 + 0.7868i 1.1674 - 0.7868i -1.3732 -0.5391 + 0.8423i -0.5391 - 0.8423i -0.7283 0.5891 + 0.3970i 0.5891 - 0.3970i ans = 0.5891 + 0.3970i 0.5891 - 0.3970i -0.7283 -0.5391 + 0.8423i -0.5391 - 0.8423i -1.3732 1.1674 + 0.7868i 1.1674 - 0.7868i rz = 1.1674 + 0.7868i 1.1674 - 0.7868i -1.3732 -0.5391 + 0.8423i -0.5391 - 0.8423i -0.7283 0.5891 + 0.3970i 0.5891 - 0.3970i ans = 0.5891 + 0.3970i 0.5891 - 0.3970i -0.7283 -0.5391 + 0.8423i -0.5391 - 0.8423i -1.3732 1.1674 + 0.7868i 1.1674 - 0.7868i-1.5-1-0.50.51-1-0.50.518Real PartI m a g i n a r y P a r t图23-5 第一类线性相位滤波器零极点分布图例23-6 求例23-2、例23-3、例23-4的线性相位系统的零极点分布图。
解 程序如例23-5，只需输入不同的h 值，结果如图23-6(a)、(b)、(c)所示。
-1-0.500.51-1-0.500.517Real PartI m a g i n a r y P a r t-1-0.500.51-1-0.500.518Real PartI m a g i n a r y P a r t（a ） (b)-1.5-1-0.500.51-1-0.500.517Real PartI m a g i n a r y P a r t(c)图23-6 线性相位滤波器零极点分布图五、实验过程1 由序列h0＝［3，－1，2，－3］为基础，构成四种类型的线性相位FIR 滤波器，即 ①h1＝［3，－1，2，－3，5，－3，2，－1，3］； ②h2＝［3，－1，2，－3，－3，2，－1，3］； ③h3＝［3，－1，2，－3，0，3，－2，1，－3］； ④h4＝［3，－1，2，－3，3，－2，1，－3］分别求它们的冲激响应和符幅特性，并在同一张图纸上描绘出来，进行比较。