通信原理设计报告new

1、这是我的课程设计的报告,大家可以参考一下,现在已经大四了,程序和报告还保存着,希望对你有所帮助,这基本都是自己做的哦,要自己努力,只能参考哦.蔡泽鹏中南大学目录第一章 概述21.1课程设计的基本原理.21.2课程设计目的.2 1.3课程设计的内容.2 第二章 总体设计及关键技术分析42.1卷积演示实验.42.2采样定理演示.52.3模拟滤波器设计演示.72.4设计切比雪夫I型低通滤波器.102.5双线性变换法设计巴特沃斯低通数字滤波器.112.6用凯塞窗设计高通滤波器.13第三章 程序实现 14 3.1卷积演示实验.143.2采样定理演示.193.3模拟滤波器设计演示.233.4设计切比雪夫I

2、型低通滤波器.263.5双线性变换法设计巴特沃斯低通数字滤波器.273.6用凯塞窗设计高通滤波器.29第四章 结束语.32参考文献33第一章 概述1.1 课程设计的基本原理数字信号处理的直接对象是数字信号，处理的方式是数值运算的方式，它涉及到的内容非常丰富和广泛，它是应用最快、成效最显著的新学科之一，做为信息专业的专业基础课，需要掌握其基本理论和基本的分析方法，通过理论和实践、原理和应用结合让学生掌握知识。而matlab是集数学计算、图形处理和程序语言设计于一体的软件，通过matlab来实现数字信号处理的有关表达，最终实现课程设计的目的。1.2课程设计目的1、通过基于MATLAB的算法仿真实验

3、及分析、基于DSP的算法综合实验等实践活动，进一步领会和深化课堂上学到的有关数字信号处理的基本概念、基本原理以及基本的信号处理操作及滤波器设计方法。 2、掌握线性卷积与圆周卷积软件实现的方法，并验证二者之间的关系；验证奈奎斯特取样定理，加深对时域取样后信号频谱变化的认识；掌握模拟滤波器的频率变换模拟高通、带通、带阻滤波器的设计与对比演示；以及根据窗函数设计滤波器。1.3课程设计内容 1. 设计一卷积演示程序。 （1）可输入任意2待卷积序列x1(n)、x2(n)，长度不做限定。测试数据为： x1(n)=1,1,1,1,0,0,1,1,1,1,0,0，x2(n)=0,1,2,1,0,0,0,1,2

4、,1,0,0； （2）分别动态演示2序列进行线性卷积x1(n)x2(n)和圆周卷积x1(n) x2 (n)的过程；要求分别动态演示翻转、移位、乘积、求和的过程。 （3）圆周卷积默认使用2序列中的最大长度，但卷积前可以指定卷积长度N用以进行混叠分析。 （4）根据实验结果分析2类卷积的关系。 （5）假定时域序列x1(n)、x2(n)的长度不小于10000，序列内容自定义。利用FFT实现快速卷积，验证时域卷积定理，并与直接卷积进行效率对比。2. 演示采样定理，时域采样，频谱周期延拓，同时演示采样频率小于2fc时，产生的混叠效应。 （1）假设频域归一化三角波的频带宽度fc=100Hz，对应的时域信号为

5、： y(t)=fcsinc(fct/4)2/2。 （2）要求输入采样频率fs（根据程序处理需要指定范围）后，在时域演示信号波形、采样脉冲及采样后信号；在频域演示对应的信号频谱、采样脉冲及频域周期拓延。 注：sinc(x)=sin(x)/(x)。 3. 模拟滤波器设计演示从模拟低通滤波器到模拟高通，带通，带阻的幅度特性对比演示。 （1）设计过程详见教材相关内容。 （2）使用巴特沃斯滤波器，其阶数N应该根据实际参数计算（计算公式和方法如教材所述），为方便作图，这里指定阶数为N=5，并假定通带截止频率fp=1，阻带截止频率fs=2。（3）分别用不同颜色曲线绘制通带、过渡带和阻带。要求根据变换关系动态

6、演示低通滤波器和目标滤波器的幅度特性。 4. 设计一切比雪夫I型低通滤波器，各参数要求：fp=100Hz，p=2dB，fs=120Hz，s=60dB。给出所设计滤波器的幅度特性并分析是否满足设计需要。5.使用双线性变换法设计巴特沃斯低通数字滤波器，各设计指标如下：s=0.4，p=0.6，p=0.5 dB，s=50 dB令T=2，给出所设计滤波器的幅度衰减特性及其冲激响应。6. 利用凯塞窗设计高通滤波器，设计指标分别是：s=0.4，p=0.6，p=0.5dB，s=60dB。给出窗函数及所设计滤波器的幅度特性。 第二章 总体设计及关键技术分析2.1 卷积演示实验2.1.1基本的原理1、线性卷积：线

7、性时不变系统（Linear Time-Invariant System, or L. T. I系统）输入、输出间的关系为：当系统输入序列为，系统的单位脉冲响应为，输出序列为，则系统输出为：或 上式称为离散卷积或线性卷积。2、圆周卷积D F T设两个有限长序列和，均为点长 如果则 N上式称为循环卷积或圆周卷积注：为序列的周期化序列；为的主值序列。编程计算时，可表示如下：3、两个有限长序列的线性卷积序列为点长，序列为点长，为这两个序列的线性卷积，则为且线性卷积的最大长，也就是说当和时。2.1.2设计思想首先建立一个基本的框架，制作一个菜单，其中包括主程序菜单和子程序的菜单，子程序菜单可以选择回到主

8、程序菜单选择功能。菜单的框架完成后，实现可以任意输入两个序列，然后分别制作动态演示序列的线性卷积的程序、动态演示序列的圆周卷积、以及验证时域卷机定理以及比较运行速率的程序。结合上面建立的框架完成菜单选择以及功能的调用，让整个设计完美。2.1.3 设计流程图：开始输入x1（n）、x2（n）序列调用程序主菜单、选择子菜单动态演示2序列的线性卷积动态演示序列的圆周卷积验证时域卷机定理、比较速率退出程序结束2.14 关键技术分析本个设计主要要实现动态的演示，为实现动态演示，有很多不同的方法，而我采用的是for循环加上pause，在for循环中的一次执行中，实现一次绘图，表示某一时刻的状态，用pause

9、暂停等待for循环中下一次执行，实现动态演示。这是主要的一个设计，另外一个重要的技巧是如何将循环卷积表示出来，因为循环卷积的结果是周期的，其循环卷积的方法和线性卷积不同，并且循环卷积中要处理几种情况：如当x1序列和x2序列之间的最大长度比卷积N长度大时是一种画图方法，比它小时是另外一种方法，然而困难的是当循环卷积长度比序列最大长度小时还要分析x1和x2序列各自长度与N的关系而做不同的情况分析，具体的解决方法见程序实现。2.2采样定理演示2.2.1基本的原理奈奎斯特取样定理指出：为了使实信号取样后能够不失真还原，取样频率必须大于信号最高频率的两倍。若为有限带宽的连续信号，其频谱为，以T为取样间隔

10、对理想取样，得到理想取样信号。的频谱为： 也就是说，一个连续信号经过理想取样后，它的频谱将沿着频率轴，从开始，每个一个取样频率重复出现一次，即频谱产生周期延拓。2.2.2设计思想首先通过时域来表示出其波形的变化，从而观察采样前后的波形做比较，第二步通过频域检验其是否混叠，观察采样前后的频域，再通过改变不同的采样频率来分析采样定理。以此设计思想编写好相关程序之后，建立整体框架，实现完成一次采样后可以重新输入采样频率，以便进行分析。2.2.3设计流程图：开始输入采样频率fs调用程序并且绘制出时域信号、采样信号以及采样之后的信号的图形继续执行程序，并且在figure2中绘画出时俞信号、采样信号以及采

11、样后的信号的频谱图，如果混叠则画出混叠部分是否继续？1/0结束是/1否/02.2.4 关键技术分析本次的重点在于如何实现时域原信号的频谱，就是说如何画出三脚波，由于在实验中通过fft不能很好的表示出其频谱，即傅立叶变换，如果按照傅立叶变换求出其频谱也很难实现，所以在实验中采用画波形，得到其原信号的波形之后的关键就在于如何动态地表示出其采样过程以及混叠的过程。动态过程同样使用for循环以及hold on和pause，通过判断fc和fs的关系看是否混叠，混叠则绘制混叠部分，通过数学计算，求出其关系式，在程序中直接应用就可以实现绘制混叠部分的说明，具体见第三章的程序实现部分。2.3模拟滤波器设计演示

12、2.3.1基本的原理巴特沃斯滤波器：特点：具有通带内最大平坦的振幅特性，且随f单调 ，其幅度平方函数具有如下形式： 式中，N为整数，称为滤波器的阶数，N越大，通带和阻带的近似性越好，过渡带也越陡。如图。 图 巴特沃兹filter 振幅平方函数 1) 在通带，分母/c1，随着N增加，/c1，A(2)快速下降。 3) =c时， ，幅度衰减 ，相当于3db衰减点。 2.3.2 设计思想 依据低通滤波器和高通滤波器、带通滤波器和带阻滤波器的关系绘制动态图形对比其幅度特性。首先从所给的参数求出所需要的绘制低通滤波器的条件参数，然后通过各滤波器之间的关系表达式画图。由于他们之间的关系是通过横做标的关系来表

13、达的，即通过转换横坐标就可以动态的画图。例如低通与高通之间，低通归一化频率w和高通归一化频率W之间的关系是w=1/W，所以，低通中w点对应于高通中1/W点的值是相同的，由此画图。2.3.3 设计流程图 开始输入数字选择菜单，进行滤波器的转换动态表示1、低通到高通的变换2、低通到带通的变换3、低通到带阻的变换4、退出程序结束2.3.4 关键技术设计本实验关键的设计在于如何实现低通于其他的转换，在本设计中是通过改变横坐标的关系来设计的，在具体的设计中，采用到的设计如下：1 模拟高通滤波器设计： 确定高通滤波器的技术指标：通带下限频率p，阻带上限频率s，通带最大衰减p，阻带最小衰减s。 确定相应低通

14、滤波器的设计指标：按1/，将高通滤波器的边界频率转换成低通滤波器的边界频率，各相设计指标为：低通滤波器通带截止频率p1/p；低通滤波器阻带截止频率s1/s；通带最大衰减仍为p，阻带最小衰减仍为s 设计归一化低通滤波器G(p)。 求模拟高通的H(s)。将G(p)按1/转换成归一化高通H(q)，为去归一化，将q=s/c代入H(q)中，得H(s)G(p)|p=c/s2 模拟带通滤波器设计思路： 确定模拟带通滤波器的设计指标，即：带通上限频率u，带通下限频率l；下阻带上限频率s1，上阻带下限频率s2 通带中心频率02=lu，通带宽度B=u-l与以上边界频率对应的归一化边界频率如下：还需要确定的技术指标

15、有：通带最大衰减p，阻带最小衰减s。 确定归一化低通技术要求； 设计归一化低通G(p)； 将G(p)转换成带通H(s)。3 模拟带通阻滤波器设计思路： 确定模拟带通滤波器的设计指标，即：下通带上限频率l，上通带下限频率u 阻带下限频率s1，阻带上限频率s2 阻带中心频率02=lu，通带宽度B=u-l 它们相应的归一化边界频率为： 归一化模拟低通技术要求； 设计归一化模拟低通G(p)； 将G(p)转换成带阻滤波器H(s)2.4 设计切比雪夫I型低通滤波器2.4.1 基本的原理误差值在规定的频段上等波纹变化。 巴特沃兹滤波器在通带内幅度特性是单调下降的，如果阶次一定，则在靠近截止 处，幅度下降很多

16、，或者说，为了使通带内的衰减足够小，需要的阶次N很高，为了克服这一缺点，采用切比雪夫多项式来逼近所希望的 。切比雪夫I型低通滤波器的 在通带范围内是等幅起伏的，所以在同样的通常内衰减要求下，其阶数较巴特沃兹滤波器要小。 2.4.2 设计思想可以通过直接调用有关切比雪夫I型低通滤波器的库函数，matlab本身就包含有这些函数供用户使用，另一种方法是按照公式一步一步地往下做，但是由于所学知识有限以及时间问题在本设计中采用了第一种方法，调用函数。然后通过定点绘制说明线来表达其特性。2.4.3 设计流程图开始输入有关参数fp、fs、ap、as调用库函数，求出N以及过度频率wc，并且求出等差频率所对应的

17、幅度值绘制幅度特性图结束选择继续进行不继续2.4.4 关键技术设计 由于采用了调用库函数的方法，所以关键的部分就不在于求出其转换的表达式了，重点在于如何更好的表达切比雪夫I型低通滤波器的特性，所以在本设计中田加了说明线，例如在其幅度的增益图象中，说明其通带的部分、过渡带以及阻带部分。2.5 双线性变换法设计巴特沃斯低通数字滤波器2.5.1 基本的原理1.确定技术指标ap，as，Wp，Ws2.求滤波器阶数N其中 3求归一化极点再将代入式，得到归一化传输函数。4将去归一化，得到实际的，即2.5.2 设计思想本设计中，由于我是选择先做这个双线性，所以采用了利用公式一步一步地做出来。首先利用双线性变换

18、的公式w=2*tan(0.5*W)/T,其中w表示模拟角频率，而W表示数字角频率，再通过s=j*和双线性变换中s和z的关系，可以得到s=(2*(1-exp(-j.*W)./(T*(1+exp(-j.*W)，其中的s是数字滤波器的s，s=jw得到结果。这样利用原理可以求出它的几个极点s，带进公式，通过选取一定范围内的频率就可以求得对应的幅度值，结果就出来了。2.5.3 设计流程图开始按照公式求出N以及各极点根据N和极点已经使用s=(2(1-exp(-jW)/(T(1+exp(-jW)带入，进行频率采样绘制图形结束2.5.4 关键技术设计 由于这个程序要按照公式一步一步地执行，所以关键在于如何实现

19、极点与频率W的转换，后来看到了s=j*W，（这个s是数字滤波器的）这就是关键所在，从参考资料上可以知道，s和z-1的关系，而z-1=e-sT,所以便可以将W和模拟滤波器中的极点s的关系求出来了，这就是本实验的关键所在。有了这个关系之后使用for循环便可以将幅度函数和W的关系求出来了。2.6用凯塞窗设计高通滤波器2.6.1 基本的原理凯塞窗函数的时域形式可表示为 其中， 是第1类变形零阶贝塞尔函数， 是窗函数的形状参数，由下式确定： 其中， 为凯塞窗函数的主瓣值和旁瓣值之间的差值(dB)。改变的取值，可以对主瓣宽度和旁瓣衰减进行自由选择。的值越大，窗函数频谱的旁瓣值就越小，而其主瓣宽度就越宽。2

20、.6.2 设计思想 按照公式首先要求出a和，然后求出w（k）,这是前面必须要求出来的，然后通过相关的转换求得结果，其思想与2.5.2基本相同。2.6.3 关键技术设计 这个设计的关键在于1、如何求得w（k）；2、如何求得实际的冲击响应。如何求得w（k）在设计思想已经大概讨论过了，如何求得冲击响应，可以通过hn=fir1(N,wc/pi,high,kaiser(N+1,b)来求得第三章 程序实现3.1卷积原理演示3.1.1 程序的实现本设计中线性卷积的主要实现主要由如下片段实现：p=length(x1);q=length(x2);n=p+q-1;a=0:q-1;y2(a+1)=x2(q-a);f

21、or n=1:p+q-1k=-q+n:1:-1+n;subplot(3,1,2)stem(k,y2)title(x2(n-m);axis(-16,16,0,24);这一部分是实现翻转并移位，在设计中是将序列x2进行翻转和移位-y=conv(x1,x2);t=1:1:nh(t)=y(t);subplot(3,1,3)t=0:n-1;stem(t,h);title(线性卷积y(n)axis(-16,16,0,24);pause(1)end以上整个部分就是实现线性卷积的过程-对于循环卷积，要求我们进行判断并根据情况做不同的分析：p=length(x1);q=length(x2);k=max(p,q)

22、;%pN if pq&pN n=0:1:N-1;x11(n+1)=x1(n+1)x22(n+1)=x2(n+1);else disp(错误，x1的长度要比x2短)end由于原程序比较长，在执行分情况讨论做循环卷积时程序很长所以不显示出来。而第三部分，验证时域卷积定理和效率对比，通过观察通过fft实现出来的线性卷积和直接进行线性卷积运行的时间长短来进行效率的对比，程序长度可以通过rand来实现足够大的长度的随机数来进行验证。3.1.2 结果及分析1、开始运行程序，会进入主菜单，按照提示进行选择：请输入x1:1,1,1,1,0,0,1,1,1,1,0,0请输入x2:0,1,2,1,0,0,0,1,

23、2,1,0,01、动态演示2序列的线性卷积2、动态演示2序列的圆周卷积3、验证时域卷机定理4、退出请选择菜单项目:1 选择1之后可以见到其动态的过程，其最后的结果如下图1： 图1 线性卷积的结果2、运行完成之后可以看到子菜单 1、重新演示2、返回主采单修改2序列的值3、返回主菜单请选择，输入：31、动态演示2序列的线性卷积2、动态演示2序列的圆周卷积3、验证时域卷机定理4、退出请选择菜单项目:2请输入执行循环卷积的数N:10 其最后的结果如下图2所表示： 图2 循环卷积的结果3、选择返回主菜单并修改x1和x2进行第三步实验，如下运行：请输入x1:rand(1,10000)请输入x2:rand(

24、1,10000)1、动态演示2序列的线性卷积2、动态演示2序列的圆周卷积3、验证时域卷机定理4、退出请选择菜单项目:3 其结果如下图所示： 由于是采用随机10000个进行设计，所以得到的图在粘贴时有点不稳定，并且在每个图都有说明标题，所以图下面就不标面图号了。通过分析，第一个图出来比第二个图慢，说明通过fft求线性卷积会比直接卷积快，通过后面2个图可以验证时域卷积定理。3.2 采样定理演示3.2.1 程序的实现程序中主要包括时域部分的分析和频域部分的分析，在时域部分主要将时间定位在t=-0.05:0.001:0.05;将采样时间T=1/fs的应用，就可以实现时域的分析n=0.1/T+1h=on

25、es(1,n);stem(t1,h,)title(采样序列h(t) %这部分是实现采样序列的时域画图而频谱的分析则首先画出原序列的频谱x=0 1 0;t=-0.1 0 0.1;plot(t,x)；title(时域信号的频谱Y(jw)以下的这一部分就是实现频域采样的程序部分wc=0.1;ws=(fs/fc)*wc;N=fix(0.5/ws)for i=0:1:Nsubplot(3,1,2)stem(-ws*i,1,);axis(-0.5,0.5,0,1);hold on;stem(ws*i,1,);axis(-0.5,0.5,0,1)subplot(3,1,3)t1=-ws*i-wc,-ws*i

26、,-ws*i+wc;t2=ws*i-wc,ws*i,ws*i+wcx=0,1,0;plot(t1,x,t2,x);hold on;axis(-0.5,0.5,0,1);pause(1)end3.2.2 结果及分析1、运行程序，程序进如采样频率输入，如下 请输入采样频率fs:150其运行结果如下图3和图4 图3 采样fc=150HZ时的时域部分 图4 采样fc=150HZ时的频域部分 重新选择输入采样频率1、选择1改变采样频率fs，继续2、按其他输入退出程序请选择：1请输入采样频率fs:300 其结果运行如图5和图6 图5 采样fc=300HZ时的时域部分 图6 采样fc=300HZ时的频域部分

27、分析：从图3到图6对比可以知道，当采样频率大于等于2倍的模拟频率时才能保证采样不混碟，从而正常的恢复到原信号。3.3模拟滤波器设计演示3.3.1 程序的实现 在本设计中对每个转换都使用3个情况，即通带、过渡带和阻带，将它们分别讨论分别绘制图形，如，低通到高通m=1时为通带、m=2时过渡带、m=3时阻带部分，现将它们的部分程序列出来： if m=1for i=0:1:10 w1=0:0.1:i*0.1; if i=0 w2=3.2; elseif i=1 w2(2)=3.1 elseif i=2 w2(3)=3.0 elseif i=3 w2(4)=2.8 else w2(i+1)=1/w1(i

28、+1); endG=1./(j.*w1).5+3.236.*(j.*w1).4+5.2361.*(j.*w1).3+5.2361.*(j.*w1).2+3.2361.*(j.*w1)+1);G1=abs(G);subplot(2,1,1)plot(w1,G1);hold onaxis(0,3.5,0,1.2)subplot(2,1,2)plot(w2,G1); hold onaxis(0,3.5,0,1.2)pause(0.4)endend这一部分就是实现完整的低通到高通中的通带部分，为了使其图形更接近资料中的图形，所以在开始的时候故意取了一些固定的值。-for i=10:1:20 w1=1:

29、0.1:0.1*i; w2(i-10+1)=1/w1(i-10+1);这部分就是实现m=2，过度带的主要函数。-for i=20:1:30 w1=2:0.1:0.1*i; w2(i-20+1)=1/w1(i-20+1);这部分就是实现阻带部分的主要函数。3.3.2 结果以及分析运行程序可以看到主菜单：*主菜单*1、观察低通与高通的频域特性2、观察低通与带通的频域特性3、观察低通与阻带的频域特性0、退出程序请选择功能：1其运行结果如图7所示： 图7 低通到高通的转变结果图7第一图表示低通、第二个表示高通，蓝色的线表示通带部分，红色表示过渡带，绿色表示的是阻带部分。 图8 低通到带通的图形转变图8

30、第一图表示低通、第二个表示带通，蓝色的线表示阻带部分，红色表示过渡带，绿色表示的是通带部分。 图9 低通到带阻的图形转变图9中第一个图表示的是低通滤波器的幅度特性，第个图是带阻的幅度特性，其中绿色表示的是通带部分，红色表示是过度代部分，蓝色的是阻带部分，其实阻带在负半轴还有一个跟它一样的波形，只是没有表示出来。3.4 切比雪夫I型低通滤波器的幅频特性设计3.4.1 程序的实现程序中主要的部分就是求Cn，程序如下：if abs(w/wp)=1cn=cos(N*acos(w/wp);elsecn=N*log(w/wp+sqrt(w/wp).2-1);cn=(exp(cn)+exp(-cn)/2;e

31、nd所用到的库函数有N,Wc=cheb1ord(wp,ws,ap,as,s);b,a=cheby1(N,ap,Wc,s); w=linspace(1,400,100)*2*pi;H=freqs(b,a,w);这一部分是实现求出阶数N以及截止频率Wc，b,a=cheby1(N,ap,Wc,s)是求得系统函数H(z)的分子和分母的多项式的系数，H=freqs(b,a,w)用于取得数字滤波器的频率所对应的幅度。3.4.2 结果及分析运行程序后可以看到菜蛋提示，按照要求输入请输入有关的参数fp的值100fs的值120ap的值2as的值50其结果运行如下图10： 图10 切比雪夫I型低通滤波器的幅频特性

32、 图10中第一个图是带增益即分贝db的形式求其特性图，图2是幅度即频谱的绝对值的特性。3.5 双线性变换法设计巴特沃斯低通数字滤波器3.5.1 程序的实现这个程序的实现主要是求N和数字系统函数，其各自的主要程序部分如下：首先求得阶数NUs=2*tan(0.5*Ws)/T;Up=2*tan(0.5*Wp)/T;Ksp=Us/Up;ksp=sqrt(10(0.1*ap)-1)/(10(0.1*as)-1);N=-log10(ksp)/log10(Ksp);N=fix(N)+1;-这一部分是实现将数字滤波器的频率W与模拟滤波器中的系统函数相互转换得到数字滤波器的系统函数Uc=Up*(10(0.1*a

33、p)-1)(-1/(2*N);W=linspace(0,1,100)*pi;s=(2*(1-exp(-j.*W)./(T*(1+exp(-j.*W);Ha=UcN;for i=0:1:N-1;sk(i+1)=Uc*exp(j*pi*(0.5+(2*i+1)./(2*N);Ha=Ha./(s-sk(i+1);endH=abs(Ha);3.5.2 结果以及分析 图11 双线性变换法实现的数字低通滤波器的幅度特性图12 双线性变换法实现的数字低通滤波器的冲击响应从图11可以看出，这个通过双线性变换法求得的低通滤波器是符合要求的数字滤波器，图12是它的冲击响应，由于冲击响应是虚数的，所以分别做了在平面

34、上和绝对值这2个图示。3.6利用凯塞窗设计高通滤波器3.6.1 程序的实现这个设计可以采用库函数法也可以采用公式法，本设计采用利用公式的做法，其部分重要的程序如下：以下这部分是求函数中幅度函数w(n)for n=0:63b=a*sqrt(1-(2*n/63-1)2);x=1;for k=1:20s=1;for i=1:k s=s*i;endx=(b/2)(k)/s)(2)+x;endy=1;for k=1:20s=1;for i=1:ks=s*i;endy=(a/2)(k)/s)(2)+y;end w(n+1)=x/y;end-这一部分是设计实际单位取样响应函数wc=(wp+ws)/2;N=-20*log10(min(1-10(-ap/20),10(-as/20);N=N+rem(N,2);hn=fir1(N,wc/pi,high,kaiser(N+1,b);3.6.2 结果以及分析 图13 凯塞窗窗函数以及其高通滤波器的幅度特性 图14 凯塞窗下的高通滤波器的理想和实际的冲击响应 分析：图13体现了凯塞窗的窗口函数的频域特性以及其对应的高通滤波器的幅度特性，图15展现了其理想和实际中冲击响应