椭圆滤波器滤波
椭圆函数滤波器
椭圆函数滤波器
维基百科，自由的百科全书
跳转到：导航, 搜索
椭圆滤波器（Elliptic filter）又称考尔滤波器（Cauer filter），是在通带和阻带等波纹的一种滤波器。
椭圆滤波器相比其他类型的滤波器，在阶数相同的条件下有着最小的通带和阻带波动。
它在通带和阻带的波动相同，这一点区别于在通带和阻带都平坦的巴特沃斯滤波器，以及通带平坦、阻带等波纹或是阻带平坦、通带等波纹的切比雪夫滤波器。
一个低通椭圆滤波器的频率响应的幅度为：
四阶低通椭圆滤波器的频率响应。
其中R n是n阶雅可比椭圆函数（Chebyshev rational functions）。
[编辑]与其他线性滤波器的比较
下图给出了椭圆滤波器与其他常见滤波器的比较，各滤波器的参数一样多。
由图可见，椭圆滤波器比其他滤波器更陡，因此在选择滤波器的时候，椭圆滤波器能够以较低的阶数获得较窄的过渡带宽，但是它在通带和阻带上都有波动。
椭圆滤波器(考尔滤波器)
椭圆滤波器（考尔滤波器）
特点：幅值响应在通带和阻带内都是等波纹的，对于给定的阶数和给定的波纹要求，椭圆滤波器能获得较其它滤波器为窄的过渡带宽，就这点而言，椭圆滤波器是最优的，其振幅平方函数为
式中，RN（Ω，L）为雅可比椭圆函数，L是一个表示波纹性质的参量。
由图可见，在归一化通带内（-1≤Ω≤1），在（0，1）间振荡，而超过ΩL后，在L2和间振荡。
L越大，ΩL也变大。
这一特点使滤波器同时在通带和阻带具有任意衰减量。
下图为典型的椭园滤波器振幅平方函数：
图中ε和A的定义与切比雪夫滤波器相同。
当Ωc、Ωs、ε和A确定后，阶次N的确定方法为：
上面讨论了三种最常用模拟低通滤波器的特性和设计方法，设计时按照指标要求，合理选用。
一般，椭圆滤波器的阶次可最低，切比雪夫次之，巴特沃兹最高，参数的灵敏度则恰恰相反。
巴特沃兹、切比雪夫、椭圆滤波器的滤波性能比较
三种滤波器的性能比较
实际滤波器由于电路实现的限制 ，只能在某些方面(通带特性、阻带 特性、衰减特性、相位特性等)逼近 理想滤波器。常用抗混滤波器有巴 特沃斯、切比雪夫、椭圆；主要特 征如下:
(1) 从幅频特性上看，具有相同阶数 的各类滤波器中，衰减特性依次为椭 圆、切比雪夫、巴特沃斯；通带波纹 依次为巴特沃斯、切比雪夫、椭圆； 巴特沃斯滤波器具有“最平幅度”特性 。
3
>> 0 >> 1 >> 2 >> 3 >> 4 >>
2
>> 0 >> 1 >> 2 >> 3 >> 4 >> Nhomakorabea1
>> 0 >> 1 >> 2 >> 3 >> 4 >>
Let `s START
>> 0 >> 1 >> 2 >> 3 >> 4 >>
巴特沃兹、切比雪 夫、椭圆滤波器的 滤波性能比较
(2) 同一种滤波器如巴特沃斯，阶数 增加，衰减特性改善，相应的实现电 路变得复杂。
(3) 巴特沃斯、切比雪夫、椭圆滤波 器是从滤波器的幅频特性上考虑的， 滤波器的相位特性较差，其中最差的 是椭圆滤波器，切比雪夫滤波器次之 ，巴特沃斯滤波器较好。
在满足相同滤波器幅频响应指 标时： （1）椭圆型的阶数最低，巴 特沃兹型的阶数最高。 （2）就满足滤波幅频响应指 标而言，椭圆型的性价比较高 ，应用比较广泛。
相位逼近情况
巴特沃思和切比雪夫滤 波器在大约3/4的通带上非常 接近线性相位特性，而椭圆 滤波器仅在大约半个通带上 非常接近线性相位特性。
归一化椭圆函数滤波器的正元件值全解及其cad
归一化椭圆函数滤波器的正元件值全解及其
cad
归一化椭圆函数滤波器是一种带有优良频率响应的滤波器。
根据其滤波器结构，其正元件值可分为4组满组件，每组有2个元件，它们分别是抗输入电感（L1，L2，L7，L8），抗源电容（C2，C3，C6，C7），抗源抗输出电容（C1，C4，C5，C8）和抗输出电感（L3，L4，L5，L6）。
传统的椭圆滤波器设计方法是根据不同的阶数，定义不同的椭圆常数（K）值，然后利用大量的数学公式来求出滤波器的正元件的值。
而归一化椭圆函数滤波器利用椭圆函数滤波器的正电容值实现元件值一次性全解，可以有效节省计算时间。
为了开发基于归一化椭圆函数滤波器的CAD，首先需要利用椭圆函数滤波器的正元件值，按照以下步骤计算归一化在特定频率表示的椭圆函数滤波器。
其次，利用椭圆函数滤波器的归一化元件值建立归一化椭圆函数滤波器的模型，并实现归一化椭圆函数滤波器的设计、仿真和电路分析。
最后，根据所需频率参数，通过仿真和计算获得归一化椭圆函数滤波器的最佳设计参数。
通过归一化椭圆函数滤波器的正元件值全解及其cad，可以实现更节省时间、更简便的椭圆滤波器设计。
此外，它还可以确保滤波器的性能，实现优异的频率响应，有效提高滤波器的性能。
椭圆滤波器的实现
摘要在模拟信号的处理中，滤波器起着重要的作用，一个好的滤波电路对一个产品的整个性能有着重要的影响。
然而用传统的方法设计无源的高阶的滤波电路，不仅耗时，计算也比较复杂。
因此本文在总结传统方法的基础上，运用计算机进行辅助设计，不仅省时，而且效率高。
本文对无源低通椭圆滤波器进行了研究。
首先阐述了滤波器设计的理论，在此基础上，了解了常用的滤波器的原理和结构。
接着根据椭圆滤波器的结构和原理，得出了理论上的椭圆滤波器。
然后分别利用滤波器设计软件Filter Solutions和MATLAB软件对椭圆滤波器进行设计和仿真。
在比较两者优点的基础上，对设计出的电路进行了改进。
最后，对仿真后的电路，利用PROTEL 99SE画出电路的原理图，生成PCB板，并对元器件的值进行了修正，以便在实际中选择元器件。
对焊接好的电路，进行了测试，结果表明设计出的椭圆滤波器符合规定的设计要求。
关键词：滤波器，传递函数，椭圆函数AbstractThe filters play the vital role in processing analog signals. A good filter circuit has important influence on the whole performance of a product. However, in traditional methods, designing the passive and higher-order filter circuits wastes a great deal time and the computation is quite complex. Therefore, on the basis of summarizing traditional methods, this paper proposes the novel method which uses computers to aid design. The method makes time saved and efficiency improved.In the paper, the passive low pass elliptic filter is studied. Firstly, the theory about filter design is explained, and usual filters are introduced. On the basis of understanding structure and principle of the elliptic filter, the elliptic filter is theoretically designed.And then separately in the light of filter design software-Filter Solution and MATLAB, elliptic filter is realized and simulated. According to their respective advantage, designed circuit is improved.Finally, Schematic and PCB of designed circuit are charted by using Protel 99. In order to choose proper components, the values of components are amended. And sealed circuit is tested. The result has indicated that proposed elliptic filter satisfies design demands.Keywords: filters; transfer function; elliptic function目录摘要 (I)ABSTRACT ................................................................................................................................... I I 目录.. (III)第一章绪论 (1)1.1滤波器的发展 (1)1.2椭圆滤波器课程设计的现实意义 (2)1.3本课程设计的主要工作 (2)第二章滤波器的理论 (3)2.1滤波器的概念 (3)2.1.1滤波网络及转移函数 (3)2.2滤波器的分类 (4)2.3相位函数和时延函数 (7)2.4有源滤波器 (8)2.4.1有源滤波器的设计 (9)2.5无源滤波器 (10)2.5.1理想低通无源滤波器 (10)2.5.2 巴特沃思和切比雪夫滤波器 (11)2.5.3椭圆滤波器 (12)2.6小结 (13)第三章无源低通椭圆滤波器的实现 (15)3.1无源低通椭圆滤波器的实现方法 (15)3.1.1FILTER SOLUTIONS方法 (15)3.1.2 MATLAB实现椭圆滤波器 (20)3.2灵敏度 (23)3.2.1极点和零点的灵敏度 (24)3.3小结 (24)第四章滤波器电路的仿真与制板 (25)4.1电路制板 (25)4.2仿真 (26)4.3电路分析 (28)4.4小结 (31)第五章总结和展望 (32)致谢.................................................................................................................. 错误!未定义书签。
常见的滤波器函数
附件9-2-1 常见的滤波器函数由于理想滤波器的特性不可能实现，因而在实际滤波器的设计中通常采用某个函数来逼近。
根据逼近函数有很多种，以下介绍根据常用的逼近函数所设计的巴特沃兹滤波器（Butterworth filter ）、切比雪夫滤波器（Chebyshev filter ）和椭圆函数滤波器（elliptic filter ）。
由这些函数所决定的实际滤波器特性各有其突出特点，有的衰减特性在过渡区很陡峭，有的相位特性（即延时特性）较为规律，应用中要根据实际需要来选用。
一、巴特沃兹滤波器巴特沃兹滤波器的特点是通带内的频率响应曲线最大限度平坦，没有起伏，而在阻带则逐渐下降为零。
巴特沃兹滤波器的时域特性也比较好，其脉冲响应具有适当的过冲及振铃。
R p =3dB 的巴特沃兹滤波器幅频特性的数学表达式为：()nn f f H 22c 1lg 101lg 10lg 20Ω+-=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛+-=式中f c 是截止频率，Ω=f /f c 是归一化频率，n 是其阶数。
这个响应在Ω=0处有20lg|H |=0dB ，其后随Ω增大而单调增大，在Ω<1即f <f c 的通带内，曲线增长极其缓慢，比较平稳；在Ω>1即f >f c 的阻带内，曲线增长甚快，比较陡峭。
因为函数Ω2n 在Ω=0处的一阶、二阶直至2n -1阶导数均为0，反映了函数的变化率极小，所以巴特沃兹响应也称为最平坦响应。
阻带曲线增长的速率由n 来决定，n 越大，增长越快。
一阶巴特沃兹滤波器的衰减率为每倍频6分贝。
二阶巴特沃兹滤波器的衰减率为每倍频12分贝，三阶巴特沃兹滤波器的衰减率为每倍频18分贝，如此类推。
图1所示为一阶至五阶巴特沃兹低通滤波器的幅频特性。
f20lg|H |/dB图1 一阶至五阶巴特沃兹低通滤波器二、切比雪夫滤波器在巴特沃兹滤波器中，幅度响应在通带内是最平坦且没有起伏的，在阻带内是单调下降的，然而衰减速度相对较为缓慢。
椭圆函数低通滤波器
GOAL
Goal OptimGoal2 Expr="dB(S(2,1))" SimInstanceName="SP1" Min= Max=-65 Weight= RangeVar[1]="freq" RangeMin[1]=240MHz RangeMax[1]=250MHz
0 -20
dB(S(2,1)) dB(S(1,1))
（七）低通滤波器的设计步骤
根据设计指标在ADS设计向导中建模仿真； ADS设计向导中建模仿真 1、根据设计指标在ADS设计向导中建模仿真； 对电路进行优化； 2、对电路进行优化； 电容换成标称值，对电感进行优化； 3、电容换成标称值，对电感进行优化； 最后对电感用专门软件计算； 4、最后对电感用专门软件计算； 画版图加工； 5、画版图加工； 安装调试； 5、安装调试；
几种常见滤波器的比较
滤波器 类型 LC滤波 器 微带滤 波器 腔体滤 波 适用频率 插损 3GHz以 下 3～ 30GHz 3～ 30GHz 较小 较大 很小 体积 接口 小 小 SMA SMA
较大 波导
（六）低通滤波器的设计指标
1、工作频率：DC～200MHz 工作频率：DC～ 2、插损 ≤1.5dB 纹波≤ 3、纹波≤1dB 带外抑制：@250MHz衰减 衰减≥ 4、带外抑制：@250MHz衰减≥60dB 5、驻波≤1.5 驻波≤ 接头：SMA-506、接头：SMA-50-K
-40 -60 -80 -100 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
freq, GHz
低通滤波器电路版图
（八）低通滤波器的调试步骤 1、根据指标要求对网络分析仪校准； 根据指标要求对网络分析仪校准； 2、对电感进行调试； 对电感进行调试； 3、改变电感线圈间距来微调电感，间距变大，电感变小； 改变电感线圈间距来微调电感，间距变大，电感变小； 4、重点调试滤波器通带内的驻波比。 重点调试滤波器通带内的驻波比。 特别注意： 特别注意： 线圈要绕正确； 1、线圈要绕正确； 2、装配要规范； 装配要规范； 3、主要的问题是电容和电感值偏离仿真值。 主要的问题是电容和电感值偏离仿真值。
椭圆函数LC带通滤波器的应用设计
椭圆函数LC带通滤波器的应用设计
滤波器类型的选择可根据滤波器设计的带宽等指标和具体的应用场合来选择。
相对带宽在20%以下的为窄带滤波器，应选用窄带滤波器的设计方法来设计；相对带宽在40%以上的为宽带滤波器，应选用宽带滤波器的设计方法来设计；而介于两者之间的为中等带宽滤波器。
由上面的指标可以看出本滤波器是窄带带通滤波器。
采用巴特沃斯滤波器来设计可以使通带内具有最大平坦的幅频响应；而切比雪夫滤波器的好处是：带外抑制好，但是带内有一定的波动；本文设计的滤波器要求带外近端抑制良好（可以用切比雪夫滤波器或椭圆函数滤波器来实现，但是从后面的分析看要使用LC滤波器，而用LC滤波器的话，使用切比雪夫形式电路元件的值过于小，很难实现，这个可以用软件仿真来说明），以此可以看出，用椭圆函数滤波器更适合。
微带滤波器通过采用不同的衬底材料可以在很大的频率范围内应用（从几百MHz到几十GHz）；同轴滤波器由于其微小的尺寸，制作精度很难达到；波导滤波器在小信号电平上，它的频率基本是8~100 GHz;陶瓷介质滤波器体积大，形状因子与品质因数较小；LC滤波器适用于本滤波器频段，且较容易制作和调试。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%%%%%%%%%%%%%%%%%%%%%%%%%%%% 参数初始化 %%%%%%%%%%%%%%%%%%%%%%%% clear all; % 关闭以往仿真的数据栈[x, fs, Nbits] = wavread('dajiahao.wav') ; % 把语音信号进行加载入Matlab 仿真软件平台中%sound( x, fs) ; % 对加载的语音信号进行回放n= length( x) ; % 求出语音信号的长度t= 0: 1/fs:( n- 1) /fs; % 得出频点figure(1);plot(t,x);xlabel('时间 (t)');ylabel('幅度（X)');title('语音信号的时域波形.' ) ;x1= fft( x, n) ; % 傅里叶变换Fx1=abs(x1);f= 0: fs/ n: fs* ( n- 1) / n; % 得出频点figure(2)plot( f, Fx1) ;xlabel('频率(f)');ylabel('幅度(x1)');title( '语音信号的频谱.' ) ; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 构造噪声 %%%%%%%%%%%%%%%%%%%%%%%A=5; %噪声的振幅u1=0.01;%噪声振幅的衰减系数u2=0.02;u3=0.03;f1=50;%多频噪声频率f2=200;f3=800;y1=u1*A*sin(2*pi*f1*t)+u1*A*sin(2*pi*f2*t)%两两噪声相加y2=u1*A*sin(2*pi*f3*t)y=y1+y2figure(3) %画出噪声的时域图，取200个点plot(t(1:200),y(1:200));grid on;xlabel('时间 (t)');ylabel('幅度（y)');title('噪声的时域图');Fy=fft(y,n); %对噪声进行傅立叶变换Fy1=abs(Fy);n1=floor(n/2); %对w求采样点数的一半f=(0:n1)*fs/n; %域上的采样点数figure(4) %画出噪声的频谱图plot(f,Fy1(1:n1+1));grid on;xlabel('频率（f）');ylabel(' 幅度（Fy1）');title('噪声的频谱图');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%信号合成%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%z= x'+y; % 噪声信号的叠加w= fft( z, n) ; % 对叠加信号进行频谱变换f= 0: fs/ n: fs* ( n- 1) / n; % 得出频点figure(5);plot(t,z);xlabel('时间 (t)');ylabel('幅度（Z)');title('加噪声后的时域波形.' ) ;figure(6);xlabel('频率 (t)');ylabel('幅度（w)');plot( f, abs(w) ) ; % 绘出频谱图axis( [ 0 8000 0 400] ) ;title( '含噪语音信号的频谱. ') ;pause(10);% sound( z, fs) ;% %%%%%%%%%%%%%%%%%%% 加噪声前后的时域图比较 %%%%%%%%%%%%%%%%%%%% figure(7) %画出加噪声前后的时域比较图subplot(2,1,1);plot(t,x);grid on;xlabel('时间 (t)');ylabel('幅度（x)');title('加噪声前的时域图');subplot(2,1,2);plot(t,z);grid on;xlabel('时间 (t)');ylabel('幅度（z)');title('加噪声后的时域图');% %%%%%%%%%%%%%%%%%%%%%%%%%% 加噪声前后的频域图比较 %%%%%%%%%%%% %%%figure(8) %画出加噪声前后的频域比较图subplot(2,1,1);plot(f,Fx1);grid on;xlabel('频率（f）');ylabel(' 幅度（Fx1）');title('加噪声前的频谱图');subplot(2,1,2);plot(f,abs(w));grid on;xlabel('频率（f）');ylabel(' 幅度（w）');title('加噪声后的频谱图'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%构造椭圆滤波器%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%% 构造800Hz的带阻滤波器%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%f0_stop1=800; %fc=20; % 设定过渡带的频率Rp=1; % 通带区的波纹系数Rs=30; % 阻带区的波纹系数ws_stop1=[(f0_stop1-fc)/(fs/2) (f0_stop1+fc)/(fs/2)]; % 通带的拐角频率wp_stop1=[(f0_stop1-5*fc)/(fs/2) (f0_stop1+5*fc)/(fs/2)]; % 阻带的拐角频率[N_stop1,wn_stop1]=ellipord(wp_stop1,ws_stop1,Rp,Rs); % 求出椭圆滤波器的阶数N及频率参数wc[num_stop1,den_stop1]=ellip(N_stop1,Rp,Rs,wn_stop1,'stop'); % 求出椭圆滤波器带阻数字滤波器的传递函数模型系数[h_stop1,w_stop1]=freqz(num_stop1,den_stop1,fs) % 求出离散系统频率响应的数值figure(9);subplot(2,3,1); % 画出带阻滤波器的幅频特性图plot(w_stop1*fs/(2*pi),20*log10(abs(h_stop1)));xlabel('频率（w_stop1）');ylabel('幅度（h_stop1)');title('带阻滤波器在800Hz处的幅频特性图');grid on;B_stop1=filter(num_stop1,den_stop1,z); %对含噪信号Y1进行带阻滤波% pause(10);% sound(B_stop1,fs); % 读出滤去800Hz噪声后的采样信号subplot(2,3,4); % 画出带阻滤波器滤去800Hz噪声后的时域图plot(t,B_stop1);xlabel('时间 (t)）');ylabel('幅度（B_stop1)');title('用带阻滤波器滤去800Hz噪声后的时域图');grid on;% %%%%%%%%%%%%%% 构造200Hz的带阻滤波器 %%%%%%%%%%%%%%%%%%%%%%%%f0_stop2=200;ws_stop2=[(f0_stop2-fc)/(fs/2) (f0_stop2+fc)/(fs/2)]; % 通带的拐角频率wp_stop2=[(f0_stop2-5*fc)/(fs/2) (f0_stop2+5*fc)/(fs/2)]; % 阻带的拐角频率[N_stop2,wn_stop2]=ellipord(wp_stop2,ws_stop2,Rp,Rs); % 求出椭圆滤波器的阶数N及频率参数wc[num_stop2,den_stop2]=ellip(N_stop2,Rp,Rs,wn_stop2,'stop'); % 求出椭圆滤波器带阻数字滤波器的传递函数模型系数[h_stop2,w_stop2]=freqz(num_stop2,den_stop2,fs) % 求出离散系统频率响应的数值subplot(2,3,2);plot(w_stop2*fs/(2*pi),20*log10(abs(h_stop2)));xlabel('频率（w_stop2）');ylabel('幅度（h_stop2)');title('带阻滤波器在200Hz处的幅频特性图');grid on;B_stop2=filter(num_stop2,den_stop2,B_stop1);%对含噪信号Y1进行带阻滤波% pause(10);% sound(B_stop2,fs); % 读出滤去800Hz和200Hz噪声后的采样信号subplot(2,3,5); % 画出带阻滤波器滤去800Hz和200Hz噪声后的时域图plot(t,B_stop2);xlabel('时间 (t)）');ylabel('幅度（B_stop2)');title('带阻滤波器滤去800Hz和200Hz噪声后的时域图');grid on;%%%%%%%%%%%%%% 构造50Hz的带阻滤波器 %%%%%%%%%%%%%%%%%%%%%%f0_stop3=50ws_stop3=(f0_stop3-fc)/(fs/2) % 阻带的拐角频率wp_stop3=(f0_stop3+5*fc)/(fs/2) % 通带的拐角频率Wws=0.0075;Wwp=0.0375;[N_stop3,wc_stop3]=ellipord(Wwp,Wws,Rp,Rs,'s'); % 求出椭圆滤波器的阶数N及频率参数wc[num_stop3,dem_stop3]=ellip(N_stop3,Rp,Rs,wc_stop3,'high'); % 求出椭圆滤波器高通数字滤波器的传递函数模型系数[h_stop3,w_stop3]=freqz(num_stop3,dem_stop3,fs); % 求出离散系统频率响应的数值subplot(2,3,3); %画出高通滤波器的幅频特性图plot(w_stop3*fs/(2*pi),(abs(h_stop3)));xlabel('频率（w_stop3）');ylabel('幅度（h_stop3)');title('高通滤波器在50Hz处的幅频特性图');grid on;B_stop3=filter(num_stop3,dem_stop3,B_stop2);subplot(2,3,6); % 画出高通滤波器滤去800Hz和200Hz和50Hz噪声后的时域图plot(t,B_stop3);xlabel('时间 (t)）');ylabel('幅度（B_stop3)');title('滤去800Hz和200Hz 和50Hz噪声后的时域图');grid on;% axis([0 10 -60000 60000]);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%用不同的滤波器滤去相应噪声频率后的时域比较图 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%figure(10);subplot(2,2,1);plot(t,z);grid on;xlabel('时间 (t)');ylabel('幅度（z)');title('加噪声后的时域图'); subplot(2,2,2);plot(t,B_stop1);xlabel(' 时间 (t)）');ylabel('幅度（B_stop1)');title('用带阻滤波器滤去800Hz噪声后的时域图');grid on;subplot(2,2,3);plot(t,B_stop2);xlabel(' 时间 (t)）');ylabel('幅度（B_stop2)');title( ' 带阻滤波器滤去800Hz和200Hz噪声后的时域图');grid on; subplot(2,2,4);plot(t,B_stop3);xlabel(' 时间 (t)）');ylabel('幅度（B_stop3)');title( '高通滤波器滤去800Hz和200Hz和50Hz噪声后的时域图');grid on;%%%%%%%%%%%%%%%%%%%%%%%%%% 用不同的滤波器滤去相应噪声频率后的频域比较 %%%%%%%%%%%%figure(11);FY=fft(z,n); %对加噪声后的采样信号1进行傅立叶变换FY1=abs(FY);n1=floor(n/2);f=(0:n1)*fs/n;FY_B_stop1=fft(B_stop1,n); %对滤去800Hz噪声后的采样信号进行傅立叶变换FY1_B_stop1=abs( FY_B_stop1);FY_B_stop2=fft(B_stop2,n); %对滤去800Hz和200Hz噪声后的采样信号进行傅立叶变换FY1_B_stop2=abs( FY_B_stop2);FY_B_stop3=fft(B_stop3,n); %对滤去800Hz和200Hz和50Hz噪声后的采样信号进行傅立叶变换FY1_B_stop3=abs( FY_B_stop3);subplot(2,2,1);plot(f, FY1(1:n1+1));grid on;xlabel('频率（f）');ylabel(' 幅度（FY1）'); title('加噪声后的频谱图');subplot(2,2,2);plot(f, FY1_B_stop1(1:n1+1));grid on;xlabel('频率（f）');ylabel(' 幅度（ FY1_B_stop1）'); title('用带阻滤波器滤去800Hz噪声后的的频谱图');subplot(2,2,3);plot(f, FY1_B_stop2(1:n1+1));grid on;xlabel('频率（f）');ylabel(' 幅度（ FY1_B_stop2）'); title('用带阻滤波器滤去800Hz和200Hz噪声后的的频谱图');subplot(2,2,4);plot(f, FY1_B_stop3(1:n1+1));grid on;xlabel('频率（f）');ylabel(' 幅度（ FY1_B_stop3）'); title('用高通滤波器滤去800Hz和200Hz和50Hz噪声后的的频谱图');pause(5); % 暂停5s%wavplay(B_stop3,fs); % 读出滤去800Hz和200Hz和50Hz噪声后的采样信号%wavwrite(B_stop3,fs,16,'dajiahao椭圆滤波器滤波后的语音信号.wav');% 写出滤波后的语音信号。