双通带切比雪夫带通滤波器的设计matlab
- 格式:doc
- 大小:311.00 KB
- 文档页数:26
目录目录 (1)Abstract (3)1 绪论 (4)2 IIR数字滤波器设计的原理与方法 (5)2.1 IIR数字滤波器设计的原理 (5)2.2 IIR 数字滤波器设计的基本方法 (7)3 IIR带通滤波器的MATLAB 设计 (9)3.1 IIR带通滤波器的设计流程 (10)3.2 IIR带通滤波器的设计步骤 (11)心得与体会 (22)参考文献 (23)摘要数字滤波是数字信号分析中最重要的组成部分之一,数字滤波与模拟滤波相比,具有精度和稳定性高、系统函数容易改变、灵活性高、不存在阻抗匹配问题、便于大规模集成、可实现多维滤波等优点。
数字滤波器的作用是利用离散时间系统的特性对输入信号波形(或频谱)进行加工处理,或者说利用数字方法按预定的要求对信号进行变换。
从广义讲,数字滤波是由计算机程序来实现的,是具有某种算法的数字处理过程。
MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。
它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。
关键词: IIR 数字滤波器,MATLAB,仿真2AbstractDigital filter digital signal analysis is the most important part of digital filter and simulation filter, it is characterized by high precision and stability, system function to change, flexibility, high impedance matching problem does not exist, easy for large scale integrated, may realize the multidimensional filtering, etc. The role of digital filter is using discrete time the characteristics of the system of the input signal waveform (or spectrum) process, or using a digital method according to the requirements of the scheduled to signal transform. Broadly speaking, digital filter is by the computer program to realize, is has some kind of algorithm digital processing process.MATLAB is released by the American mathworks company mainly face of scientific calculation, visualization and interactive program design of the high-tech computing environment.It will numerical analysis, calculation , scientific data visualization and nonlinear dynamic system and simulation, and many other strong function integration in an easy to use Windows environment, for scientific research, engineering design, and to effectively the numerical calculation many fields of science provides a comprehensive solution, and to a large degree from the traditional the interactive programming language (such as C, Fortran) edit mode, which represents the current international scientific computing software advanced level.Keywords: IIR digital filters, MATLAB, the simulation31 绪论在现代通信系统中,由于信号中经常混有各种复杂成分,所以很多信号的处理和分析都是基于滤波器而进行的。
matlab 切比雪夫带通滤波器实现
本文介绍如何使用matlab实现切比雪夫带通滤波器。
切比雪夫
滤波器是一种数字滤波器,可在给定的频率范围内阻止不需要的频率分量。
切比雪夫滤波器的特点在于它对幅频响应的最大偏差是可控的,因此被广泛地应用于信号处理、图像处理、通信等领域。
要在 matlab 中实现切比雪夫带通滤波器,需要先确定以下参数:通带频率范围、阻带频率范围、通带最大衰减度、阻带最小衰减度。
然后,使用 matlab 中提供的 cheb1ap 函数来计算切比雪夫滤波器
的传递函数。
具体步骤如下:
1. 确定通带频率范围、阻带频率范围、通带最大衰减度、阻带
最小衰减度,将这些参数赋值给对应的变量。
2. 使用 cheb1ap 函数计算切比雪夫滤波器的传递函数。
cheb1ap 函数的基本调用格式为 [n,wn]=cheb1ap(Wp,Ws,Rp,Rs),其中 Wp 和 Ws 分别是通带和阻带的归一化频率,Rp 和 Rs 分别是通
带最大衰减度和阻带最小衰减度。
函数返回的 n 和 wn 分别表示滤
波器的阶数和角频率。
3. 将传递函数转换为离散时间域上升通带滤波器的差分方程,
使用 tf2zp 函数将差分方程转换为零极点形式。
4. 使用 zp2sos 函数将零极点形式转换为二阶序列滤波器表示。
5. 使用 sosfilt 函数对信号进行滤波处理。
6. 将滤波结果可视化,比较滤波前后的信号,检查滤波效果。
使用 matlab 实现切比雪夫带通滤波器需要一定的数学基础和编程经验,但是掌握了这种滤波器的应用方法可以为信号处理和通信方面的工作提供很大的便利。
课 程 设 计20011 年 7月 1日 设计题目 学 号专业班级 指导教师 学生姓名 张腾达 吴晔 陈丽娟 杨蕾通信电子电路课程设计 ——数字滤波器的设计 张静 20080302 光信息08-3班实验组员 张静 胡磊 艾永春 赵亚龙王宏道 胡进娟 马丽婷设计要求:某系统接收端接收到的信号为y=cos(2π*60t)+1.2cos(2π*140t)+2sin(2π*220t)+1.5sin(2π*300t)(A) 发现此信号夹杂了一个正弦噪声noise=1.5sin(2π*300t),请设计一个低通滤波器将此噪声滤除,从而恢复原信号。
(B) 发现此信号夹杂了一个正弦噪声noise= cos(2π*60t)+1.5sin(2π*300t) ,请设计一个带通滤波器将此噪声滤除,从而恢复原信号。
(C) 发现此信号夹杂了一个正弦噪声noise=1.2cos(2π*140t)+2sin(2π*220t),请设计一个带阻滤波器将此噪声滤除,从而恢复原信号。
(D) 发现此信号夹杂了一个正弦噪声noise= cos(2π*60t),请设计一个高通滤波器将此噪声滤除,从而恢复原信号。
要求:(1)请写出具体的MATLAB程序,并详细解释每条程序(2)画出滤波前后信号的频谱图(3)画出所设计滤波器的幅频和相频特性图,并写出具体参数参数计算:根据题目要求,开始选取Wp=2*60π,Ws=2*140π。
后来经老师指点,为了将阻带里的信号更好的滤除,通带里的信号更好的保持,达到较好的滤波效果,通带截止频率选取:Wp=2*70π>2*60π,阻带截止频率选取:Ws=2*120π<2*140π,输入信号为:y=cos(2π*60t)+1.2cos(2π*140t)+2sin(2π*220t)+1.5sin(2π*300t) 可知信号最高频率为2*300*π/(2π)=300Hz。
由奈奎斯特抽样定理得,fs>=2*300=600(Hz),这里为了得到更好的抽样效果,同时简化计算,选取fs=1000Hz。
MATLAB带通滤波器1. 简介带通滤波器是一种数字信号处理中常用的滤波器。
它可以选择特定的频率范围内的信号并传递,同时抑制其他频率范围的信号。
在MATLAB中,可以使用信号处理工具箱中的函数来设计和实现带通滤波器。
本文档将介绍如何使用MATLAB设计和使用带通滤波器,包括滤波器的设计方法和常见的应用场景。
2. 带通滤波器的设计带通滤波器的设计过程可以分为以下几个步骤:2.1 滤波器类型选择MATLAB中提供了多种带通滤波器类型的设计方法,包括巴特沃斯滤波器、切比雪夫滤波器和椭圆滤波器等。
根据需求选择合适的滤波器类型。
2.2 滤波器规格确定确定滤波器的通带范围、阻带范围和过渡带宽等规格参数。
2.3 滤波器设计根据滤波器类型和规格参数,使用相应的MATLAB函数进行滤波器设计。
常用的函数包括butter、cheby1和ellip等。
2.4 滤波器特性分析设计完成的滤波器可以通过频率响应、相位响应和零极点分布等特性进行分析。
MATLAB提供了函数来绘制和分析滤波器的特性曲线。
3. MATLAB中的带通滤波器函数MATLAB提供了多个函数用于设计和实现带通滤波器,下面介绍其中几个常用的函数:3.1 butter函数butter函数可用于设计巴特沃斯滤波器。
它的语法为:[b, a] = butter(n, Wn, 'bandpass')其中,n表示滤波器的阶数,Wn为通带范围,可以是一个长度为2的向量表示最低频率和最高频率的范围。
'bandpass'表示带通滤波器。
3.2 cheby1函数cheby1函数可用于设计切比雪夫滤波器。
它的语法为:[b, a] = cheby1(n, Rp, Wn, 'bandpass')其中,n表示滤波器的阶数,Rp为通带中允许的最大衰减量,Wn为通带范围,可以是一个长度为2的向量表示最低频率和最高频率的范围。
'bandpass'表示带通滤波器。
matlab 切比雪夫带通滤波器实现带通滤波器是一种常见的数字信号处理技术,它可用于提取目标频率段内的信号,滤除无用频率段内的干扰信号。
在数字信号处理领域中,存在着多种带通滤波器的实现方式,而切比雪夫带通滤波器是其中一种方法。
切比雪夫带通滤波器与其他滤波器不同之处在于,它能够提供更为陡峭的通带和阻带边缘。
这是由于它声称的Chebyshev响应函数,其中的“Chebyshev”指的是切比雪夫多项式。
切比雪夫滤波器将极陡峭的边缘带给了滤波器,以便按照严格的信号频率进行过滤。
切比雪夫滤波器的设计,需要知道带通滤波器的通带上限fc,下限fc’,以及有选择性地保留其特定频率范围内的信号。
为了从信号中切除这些频率外的信息,切比雪夫设计要求用户指定这些频率范围的阻带跨越带宽bw,以及在通带中允许的最大衰减或输出变弱的增益tol_db。
为了设计一个切比雪夫带通滤波器,我们需要首先输入所需的频率截止值和带宽,然后计算其非标准化参数,最终,选择适当的滤波器阶数进行计算,生成相应的滤波器系数。
Matlab中可以使用signal库中的函数设计切比雪夫带通滤波器,具体步骤如下:1.输入所需的通带上限fc、下限fc’,以及阻带跨越带宽bw和允许的最大衰减或输出变弱的增益tol_db。
2.计算非标准化参数epsilon,并将其应用于滤波器阶数N的计算。
3.使用design()函数和signal.fir1()函数找到设计参数N和alpha。
4.通过signal.fir1()函数生成切比雪夫带通滤波器系数。
例如,以下代码可实现一个带有通带上限50Hz和下限20Hz,阻带跨越带宽10Hz和最大衰减40dB的切比雪夫带通滤波器:fc1 = 20;fc2 = 50;fs = 200;bw = 10;tol_db = 40;%计算非标准化参数epsilondelta_omega_p = 2*pi*bw/fs;delta_omega_s1 = 2*pi*(fc1-bw/2)/fs;delta_omega_s2 = 2*pi*(fc2+bw/2)/fs;epsilon = sqrt(10^(tol_db/10)-1);A =log10((1/epsilon)+sqrt(1+(1/epsilon)^2))/log10(1/delta_omega_ p);%求滤波器阶数和alpha[N, alpha] = cheb1ord(delta_omega_p, delta_omega_s2,tol_db, tol_db);%通过signal.fir1()函数生成切比雪夫带通滤波器系数h = fir1(N, [delta_omega_s1 delta_omega_s2]/(2*pi),'bandpass', chebwin(N+1, alpha));使用这个滤波器可以滤除一些干扰信号,保留20Hz到50Hz之间的信号,得到目标带通信号。
matlab 切比雪夫带通滤波器实现
本文将介绍如何使用MATLAB实现切比雪夫带通滤波器。
切比雪夫带通滤波器是一种数字滤波器,用于将一定范围内的频率信号从输入信号中滤除。
使用 MATLAB 对信号进行数字滤波可以得到更加精确的结果,并且可以方便地进行可视化分析。
步骤:
1. 定义需要滤波的信号。
可以使用 MATLAB 自带的信号生成函数,如 sin、cos 等,也可以导入自己的信号数据。
2. 定义切比雪夫带通滤波器的参数。
主要包括通带频率、阻带频率、通带最大衰减、阻带最小衰减等。
3. 使用 MATLAB 自带的函数 butter 设计数字滤波器。
其中,第一个参数为滤波器的阶数,第二个参数为通带或阻带的频率,第三个参数为滤波器类型,可以选择带通、带阻、低通或高通等。
4. 使用设计好的滤波器对输入信号进行滤波。
可以使用MATLAB 自带的函数 filter 实现。
5. 对滤波后的信号进行可视化展示。
可以使用 MATLAB 自带的绘图函数 plot 进行绘制。
6. 分析滤波效果。
通过观察滤波后的信号波形和频谱图,可以评估滤波器的效果是否符合要求。
总之,使用 MATLAB 实现切比雪夫带通滤波器可以帮助我们更加精确地处理信号,并且可以通过可视化手段来分析滤波效果。
- 1 -。
开始↓读入数字滤波器技术指标↓将指标转换成归一化模拟低通滤波器的指标↓设计归一化的模拟低通滤波器阶数N 和截止频率↓模拟域频率变换,将H(P)变换成模拟带通滤波器H(s)↓用双线性变换法将H(s)转换成数字带通滤波器H(z)↓输入信号后显示相关结果求相应的幅频响应与相频响应↓50100150-202tx 1(t )x1的波形50100150-202tx 2(t )x2的波形50100150-202t x (t )输入信号x 的波形10203040-0.01-0.00500.0050.01ty滤波器输出y 的波形clc;clear all ;结束%数字滤波器的技术指标Rp = 1; % 通带最大衰减Rs = 40;% 阻带最小衰减OmegaS1_1=350; % 通带截止频率OmegaS1_2=550;% 通带截止频率OmegaP1_1=400; % 阻带截止频率OmegaP1_2=500;% 阻带截止频率Fp=2000; % 抽样频率Wp1=2*pi*OmegaP1_1/Fp; % 模数频率变换Wp2=2*pi*OmegaP1_2/Fp;Ws1=2*pi*OmegaS1_1/Fp;Ws2=2*pi*OmegaS1_2/Fp;OmegaP1=2*Fp*tan(Wp1/2); % 非线性变换OmegaP2=2*Fp*tan(Wp2/2); % 非线性变换OmegaS1=2*Fp*tan(Ws1/2); % 非线性变换OmegaS2=2*Fp*tan(Ws2/2); % 非线性变换OmegaP0=sqrt(OmegaP1*OmegaP2);% 等效中心频率Bw=OmegaP2-OmegaP1; % 带通滤波器的通带宽度Eta_P0=OmegaP0/Bw; % 归一化处理Eta_P1=OmegaP1/Bw; % 归一化处理Eta_P2=OmegaP2/Bw; % 归一化处理Eta_S1=OmegaS1/Bw; % 归一化处理Eta_S2=OmegaS2/Bw; % 归一化处理Lemta_P_EquivalentLowPass=Eta_P2/(Eta_P2^2-Eta_P0^2); % 转换成低通参数Lemta_S1_EquivalentLowPass=-Eta_S1/(Eta_S1^2-Eta_P0^2); % 转换成低通参数Lemta_S2_EquivalentLowPass=Eta_S2/(Eta_S2^2-Eta_P0^2); % 转换成低通参数Lemta_S_EquivalentLowPass=min(Lemta_S1_EquivalentLowPass,Lemta_S2 _EquivalentLowPass); % 取最小值% E求滤波器阶数[N, Wn]=cheb2ord(Lemta_P_EquivalentLowPass, Lemta_S_EquivalentLowPass, Rp, Rs,'s');% 滤波器设计[num1,den1]=cheby2(N,Rs,Wn,'s');[num2,den2]=lp2bp(num1,den1,OmegaP0,Bw);[num,den]=bilinear(num2,den2,Fp);[Z,P,K]=cheb1ap(N,Rp);w=linspace(1,1000,100)*2*pi;[M1,N1]=zp2tf(Z,P,K); %将零极点形式转换为传输函数形式[M,N]=lp2bp(M1,N1,OmegaP0,Bw); %对低通滤波器进行频率变换转换为带通滤波器% 计算增益响应w = 0:pi/255:pi;h = freqz(num,den,w);g = 20*log10(abs(h));%绘制切比雪夫带通滤波器幅频特性figure;plot(w/pi,g);gridaxis([0 1 -60 5]);xlabel('\频率/\pi'); ylabel('增益/dB'); title('切比雪夫II型带通滤波器幅频响应');%Plot the poles and zeros[z,p,k]=tf2zp(num,den);figure;zplane(z,p); %绘制传输函数零极点title('?传输函数的零极点')f1=450;f2=600;t=0:0.0001:1x1=sin(2*pi*f1*t);x2=sin(2*pi*f2*t);x=x1+x2;figure;subplot(2,2,1)%绘制x1的波形plot(x1);grid on;axis([0,50*pi,-3,3]);xlabel('t');ylabel('x1(t)');title('x1的波形');subplot(2,2,2)%绘制x2的波形plot(x2);grid on;axis([0,50*pi,-3,3]);xlabel('t');ylabel('x2(t)');title('x2的波形');subplot(2,2,3)%绘制输入x的波形plot(x);grid on;axis([0,50*pi,-3,3]);xlabel('t');ylabel('x(t)');title('输入信号x的波形')%X=fft(x);y=filter(num,den,x);%数字滤波器输出subplot(2,2,4);plot(real(y));grid on;axis([0,15*pi,-0.01,0.01]);xlabel('t');ylabel('y');title('滤波器输出y的波形');附录:PPpppp5. 用双线性变换法设计IIR数字带通滤波器例21-3采用双线性变换法设计一个切比雪夫Ⅰ型数字带通滤波器,要求:通带wp1=0.3p,wp2=0.7p,Rp=1 dB;阻带ws1=0.2p,ws2=0.8p,As=20 dB解程序如下:wp1=0.4*pi;wp2=0.5*pi;ws1=0.35*pi;ws2=0.55*pi;Rp=1;As=40;T=0.0005;Fs=1/T;Omgp1=(2/T)*tan(wp1/2);Omgp2=(2/T)*tan(wp2/2);Omgp=[Omgp1,Omgp2];Omgs1=(2/T)*tan(ws1/2);Omgs2=(2/T)*tan(ws2/2);Omgs=[Omgs1,Omgs2];bw=Omgp2-Omgp1;w0=sqrt(Omgp1*Omgp2);bw=Omgs2-Omgs1;w0=sqrt(Omgs1*Omgs2); %[ZK(]模拟滤波器阻带带宽和中心频率[n,Omgn]=cheb2ord(Omgp,Omgs,Rp,As,'s') %计算阶数n和截止频率[z0,p0,k0]=cheb2ap(n,As); %设计归一化的模拟原型滤波器[n,Omgn]=cheb1ord(Omgp,Omgs,Rp,As,'s')[z0,p0,k0]=cheb1ap(n,Rp);ba1=k0*real(poly(z0));aa1=real(poly(p0));[ba,aa]=lp2bp(ba1,aa1,w0,bw);[bd,ad]=bilinear(ba,aa,Fs)[H,w]=freqz(bd,ad);dbH=20*log10((abs(H)+eps)/max(abs(H)));subplot(2,2,1),plot(w/2/pi*Fs,abs(H),'k');ylabel('|H|');title('幅度响应');axis([0,Fs/2,0,1.1]);set(gca,'XTickMode','manual','XTick',[0,fs,fp,Fs/2]);set(gca,'YTickMode','manual','YTick',[0,Attn,ripple,1]);gridsubplot(2,2,2),plot(w/2/pi*Fs,angle(H)/pi*180,'k');ylabel('\phi');title('相位响应');axis([0,Fs/2,-180,180]);set(gca,'XTickMode','manual','XTick',[0,fs,fp,Fs/2]);set(gca,'YTickMode','manual','YTick',[-180,0,180]);gridsubplot(2,2,3),plot(w/2/pi*Fs,dbH);title('幅度响应( dB)');axis([0,Fs/2,-40,5]);ylabel('dB');xlabel('频率(\pi)');set(gca,'XTickMode','manual','XTick',[0,fs,fp,Fs/2]);set(gca,'YTickMode','manual','YTick',[-50,-20,-1,0]);gridsubplot(2,2,4),zplane(bd,ad);axis([-1.1,1.1,-1.1,1.1]);title('零极图');程序运行结果如下:n = 3Omgn =1.0e+003 * 1.0191 3.9252bd =0.0736 0.0000 -0.2208 0.0000 0.2208 -0.0000 -0.0736ad =1.0000 0.0000 0.9761 0.0000 0.8568 0.0000 0.2919 采用双线性变换法设计一个切比雪夫Ⅱ型数字带通滤波器,其它条件不变,则需要修改下面几句程序:bw=Omgs2-Omgs1;w0=sqrt(Omgs1*Omgs2); %[ZK(]模拟滤波器阻带带宽和中心频率[n,Omgn]=cheb2ord(Omgp,Omgs,Rp,As,'s') %计算阶数n和截止频率[z0,p0,k0]=cheb2ap(n,As); %设计归一化的模拟原型滤波器采用阻带截止频率来计算W0和BW,是因为切比雪夫Ⅱ型模拟低通原型是以阻带衰减As为主要设计指标的。
课程设计题目Matlab课程设计——利用MATLAB 结合双线性变换法设计一个数字切比雪夫带阻IIR滤波器学院信息工程学院专业通信工程班级通信0803姓名王欢指导教师魏洪涛学号:01208103408142011 年1月月14日日武汉理工大学《Matlab课程设计》报告课程设计任务书学生姓名:王欢专业班级:通信0803指导教师:魏洪涛工作单位:信息工程学院题目: Matlab课程设计——利用MATLAB仿真软件系统结合双线性变换法设计一个数字切比雪夫带阻IIR滤波器。
初始条件:Matlab基础知识、计算机要求完成的主要任务:1.方案的理论设计2.方案的安装、调试3.设计报告的撰写时间安排:指导教师签名: 2010 年 1月日系主任(或责任教师)签名: 2010 年 1月日目录摘要 (3)ABSTRACT (4)1 数字滤波器 (5)1.1概述 (5)1.2特点 (6)1.3分类 (6)1.4设计原理 (6)2 双线性变换法 (8)2.1简介 (8)2.2对比 (11)2.2.1优点 (11)2.2.2缺点 (12)3 切比雪夫滤波器 (13)3.1概述 (13)3.2切比雪夫滤波器的种类 (13)3.2.1 I型切比雪夫滤波器 (13)3.2.2 II型切比雪夫滤波器 (13)3.3特点 (13)4 用MATLAB实现切比雪夫IIR带阻滤波器 (14)4.1程序流程图 (14)4.2MATLAB程序代码 (15)4.3仿真结果 (16)5 学习小结 (17)6 参考文献 (17)摘要随着信息时代和数字世界的到来,数字信号处理已成为当今一门极其重要的学科和技术领域。
目前数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。
在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器(DF,Digital Filter)。
数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。
切比雪夫II 型带通IIR 数字滤波器设计1.设计思路(1) 数字—模拟指标转换。
利用双线性变换的频率预畸变公式2tan 2ωT =Ω,把所要求的数字滤波器)(z H 数字频率指标转换为相应的模拟滤波器)(s H 的模拟频率指标。
(2) 低通模拟指标转换。
将模拟滤波器)(s H 的频率指标归一化原型低通滤波器)(p H LP 的频率指标。
(3) 模拟滤波器的设计,得到低通滤波器的归一化传输函数)(p H LP 。
(4) 模拟频率变换,将模拟低通滤波器归一化传输函数)(p H LP 转换成所需要的模拟滤波器传输函数)(z H 。
(5) 模拟—数字滤波器变换。
利用双线性变换得到所要求的数字滤波器传输函数11112)()(--+-==z z T s s H Z H 。
2.设计要求及方案设计一带通切比雪夫II 型IIR 滤波器,要求如下:通带上下边沿频率分别为300 Hz 和400 Hz ,通带最大衰减dB p3=α,阻带上下边沿频率分别为200Hz 和500 Hz ,阻带最小衰减 dB s 18=α,采样频率Hz f s 2000=。
2.1带通滤波器的边缘频率为Hz f p 3001=,Hz f p 4002= Hz f s 2001=,Hz f s 5002=给定的系统取样频率为Hz f s 2000=,相应的数字频率为ππω3.0211=⨯=sp p f f ππω4.0222=⨯=sp p f fππω2.02=⨯=s sl sl f f ππω5.0222=⨯=ss s f f2.2采用双线性变换,以获得相应的模拟带通滤波器的边缘频率s rad f p s p /1.20382tan211==Ωω s rad f p s p /2.29062tan222==Ωωs rad f s s s /7.12992tan211==Ωω s rad f s s s /0.40002tan222==Ωωs rad B p p /1.86812=Ω-Ω=2.3归一化低通滤波器的技术指标1=p λ 11.31212≈Ω-ΩΩ-Ω=p p s s s λ通带最大衰减dB p3=α阻带最小衰减 dB s 18=αss s B s p w p 1001200002202+=Ω+=λ 用MTALAB 算法设计归一化切比雪夫II 型低通模拟滤波器>> [N2,wp2]=cheb2ord(wp,ws,Rp,As,'s'); %Chebyshev Ⅱ型滤波器参数计算(模拟域); >> [B2,A2]=cheby2(N2,Rp,wp2,'s'); %计算H(s)的系数B 和A ;>> [Z,P,K]=cheby2(N2,As,wp2,'s'); %构造Chebyshev Ⅱ型滤波器(零极点模型) >> subplot(2,2,1);>> [Z,P,K]=cheby2(N2,As,wp2,'s'); %构造Chebyshev Ⅱ型滤波器(零极点模型) >> [H,W]=zp2tf(Z,P,K); %将零极点模型转化成传递函数的模型 >> figure(1);>> [P,Q]=freqs(H,W); %在Figure1上显示滤波器的幅频响应及相频响应>> figure(2); >> subplot(2,2,1);>> fk=0:12000/512:12000;wk=2*pi*fk; %在Figure2上显示幅频特性曲线>> Hk=freqs(B2,A2,wk);>> plot(fk/1000,20*log10(abs(Hk)));grid on >> xlabel('Frequency(Hz)') >> ylabel('Magnitude Response')2.4将)(p H LP 转化为带通滤波器的系统sB s p LP w pp H s H 202)()(Ω+==λ2.5用双线性变换法将)(s H 转换成数字滤波器)(z H ,即1112)()(-+-==z z s s H Z H()()()()()()()()()()()()()()()⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-+++-+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-+++-∙⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-+++-+⎪⎪⎪⎪⎪⎭⎫⎝⎛-+++-+-+++-=--------------------6359.011219001148720.0112190011411949.111219001143331.01121900114)5389.01121900114(442.21111212111211112121112111121zz z z z z z z z z z zz z z z z z zz3.利用MATLAB 一步编写切比雪夫II 型带通数字滤波器Matlab 总程序如下:>> W1=300;W2=400;rp=3;rs=18;Fs=2000; >> WP=[200,300];WS=[50,450];>> [N,Wn]=cheb2ord(WP/(Fs/2),WS/(Fs/2),rp,rs);>> [P,Q]=cheby2(N,rp,Wn,'bandpass'); %创建Chebyshev 带通滤波器 >> ylabel('幅度'); >> figure(1);>> freqz(P,Q); %显示产生滤波器的幅频及相频曲线 >> [H,W]=freqz(P,Q); >> figure(2);>> plot(W*Fs/(2*pi),abs(H));grid; >> xlabel('频率/Hz');>> ylabel('幅度')>> title('数字滤波器幅频响应|H(ejOmega)| ');仿真出的幅频特性曲线如下图1.1所示:图1.1:幅频特性曲线相频特性及幅度特性曲线如下图1.2所示:图1.2:相频特性及幅度特性曲线。
中北大学课程设计说明书学生姓名:学号:学院:信息与通信工程学院专业:电子信息工程题目:信息处理信息实践:信号的产生与带通滤波器的设计指导教师:职称:2016 年1 月4 日中北大学课程设计任务书2015/2016 学年第一学期学院:信息与通信工程学院专业:电子信息工程学生姓名:学生姓名:学号:学生姓名:学号:课程设计题目:信息处理信息实践:信号的产生与带通滤波器的设计起迄日期:2016年1月4日~2016年1月22日课程设计地点:学院楼201实验室指导教师:王玉下达任务书日期: 2016 年1月4 日课程设计任务书课程设计任务书目录一基于USB总线A/D卡的报告 (1)B总线介绍 (1)B接口电路设计 (1)B接口的数据采集系统的设计实现 (1)4.A/D转换电路 (2)5.A/D采集卡的应用 (2)二设计内容 (3)三设计原理 (3)四MATLAB软件介绍 (4)五设计过程 (5)1.信号的产生 (5)2.对连续输入信号进行采样,进行频谱分析 (7)3.切比雪夫带通滤波器 (8)4.对信号进行滤波 (10)六结果分析 (12)七心得体会及总结 (13)参考文献 (13)附录:总程序代码 (15)一基于USB总线A/D卡的报告B总线介绍通用串行总线标准USB是1995年微软、康柏、DEC、IBM等公司为解决传统总线不足的问题而推出的一种新型通信标准。
该标准具有高带宽、安装方便、易扩展诸多优点,逐渐成为现代数据传输的主要手段。
基于USB数据采集卡充拥有USB总线的上述优点,有效优化解决传统数据采集卡的安装麻烦、扩展不便等弊端。
USB总线为通用串行总线,USB接口位于PS/2接口和串并口之间,允许外设在开机状态下热插拔,最多可串接下来127个外设,传输速率可达480Mb/S,P它可以向低压设备提供5伏电源,同时可以减少PC机I/O接口数量。
USB是基于通用连接技术,实现外设的简单快速连接,达到方便用户、降低成本、扩展PC连接外设范围的目的。
图1-1 LSI 系统若y(n)、x(n)的傅立叶变换存在,则输入输出的频域关系是:Y()()*()j j j e X e H e ωωω= (B.7)当输入信号x(n)通过滤波器系统h(n)后,其输入y(n)中不再具有|ω|<c ω的频率域成分,仅使|ω|>c ω的信号成分通过。
因此,滤波器的形状不同,其滤波后的信号也不同样。
若滤波器的输入、输出都是离散时间信号,那么该滤波器的单位冲激响应h(n)也必然是离散的,这种滤波器称为为数字滤波器(DF, Digital Filter )。
当用硬件实现一个DF 时,所需的元件是延迟器、乘法器和加法器:而运用MATLAB 软件时,它仅需线性卷积程序便可以实现。
众所周知,模拟滤波器(AF, Analog Filer)只能用硬件来实现,其元件是电阻R ,电感L ,电容C 及运算放大器等。
因此DF 的实现要比AF 容易得多,且更容易获得较抱负的滤波性能。
滤波器的种类滤波器的种类很多,分类方法也不同,可以从功能上分,也可以从实现方法上分,或从设计方法上分等。
但总的来说,滤波器可分为两大类,即经典滤波器和现代滤波器。
经典滤波器是假定输入信号x(n)中的有用成分和无用成分(如噪声)各自占有不同的频带,当x(n)通过滤波器后可将无用成分有效滤去。
假如信号中的有用成分和无用成分的频带互相重叠,那么经典滤波器将无法滤除信号中的无用成分。
现代滤波器理论研究的重要内容是从具有噪声的数据记录(又称为时间时间系列)中估计出信号的某些特性或信号自身。
一旦信号被估计出,那么估计出的信号与原信号相比较高的信噪比。
现代滤波器把信号和噪声都视为随机信号,运用它们的记录特性(如自相关函数、功率谱等)推导出一套最佳的估值算法,然后用硬件或软件予以实现。
现在滤波器理论源于维纳在20世纪40年代及其以后的工作,因此维纳滤波器便是这一类滤波器的典型代表。
此外,尚有卡尔曼滤波器、线性预测滤波器、自适应滤波器等。
matlab 切比雪夫带通滤波器实现切比雪夫带通滤波器是一种数字滤波器,用于去除信号中的噪声和杂波。
该滤波器可以滤波信号的特定频率段内的噪声,并在此频率段外保留信号的有效信息。
本文将介绍matlab中如何实现切比雪夫带通滤波器。
切比雪夫带通滤波器是一种设计采用无限脉冲响应(IIR)数字滤波器的方法。
它是根据Chebyshev多项式来设计滤波器的,其设计目标是让通带中的波形误差最小化,同时让阻带的响应最大化。
切比雪夫带通滤波器的特征是在通带中有强烈的波纹条件,而在阻带中则有一定的振荡。
二、matlab实现在matlab中,实现切比雪夫带通滤波器需要用到cheby1函数。
cheby1函数的语法如下:[b,a] = cheby1(n,Rp,Wn)其中n是滤波器的阶数,Rp是通带最大衰减量,Wn是通带截止频率。
该函数的输出是滤波器的分子和分母系数。
以下是一个实现切比雪夫带通滤波器的例子:% 生成信号t = 0:0.001:1; % 采样时间f1 = 10; % 信号频率f2 = 50; % 基波频率x = sin(2*pi*f1*t) + 0.5*sin(2*pi*f2*t);% 设计滤波器fs = 1000; % 采样频率fcuts = [20,30]; % 通带截止频率Rp = 1; % 通带最大衰减量[n, Wn] = cheb1ord(fcuts/(fs/2), [25,30]/(fs/2), Rp); % 计算阶数和截止频率[b, a] = cheby1(n, Rp, Wn); % 设计滤波器% 滤波信号y = filter(b, a, x);% 绘制滤波器figurefreqz(b, a)% 绘制原始信号和滤波后的信号figureplot(t, x)hold onplot(t, y)legend('原始信号', '滤波后信号')在上面的例子中,我们首先生成了一个带有两个频率分量的信号,然后使用cheby1函数设计了一个通带在20Hz到30Hz之间的切比雪夫带通滤波器。
数字切比雪夫滤波器的设计及matlab仿真数字切比雪夫滤波器是一种常见的数字滤波器,它是由切比雪夫多项式设计而成的。
切比雪夫多项式是一类代数多项式,在理论计算和实际应用中具有广泛的应用。
数字切比雪夫滤波器的设计和matlab仿真是数字信号处理领域中重要的研究内容。
数字切比雪夫滤波器的设计基本步骤如下:首先确定滤波器类型,例如低通滤波器、高通滤波器、带通滤波器、带阻滤波器等。
然后,通过切比雪夫多项式计算出滤波器的传递函数H(s),转换成离散传递函数H(z),并确定滤波器的通带截止频率、阻带截止频率和通带最大衰减量。
最后在matlab中实现数字切比雪夫滤波器的设计和仿真。
在数字切比雪夫滤波器的设计中,需要注意以下几点:1.滤波器类型的选择要符合实际应用中的需求和要求,例如在音频处理中常用的是低通滤波器,用于去除高频噪声和杂音。
2.切比雪夫多项式的阶数越高,滤波器的设计越精细,但会使设计过程变得更加困难和计算量更大,需要权衡设计精细程度和计算成本。
3.通带截止频率和阻带截止频率的设置应该根据实际信号的频率分布情况来确定,保证滤波器能够有效去除噪声信号。
4.通带最大衰减量的选择应该综合考虑信号处理的精度和计算性能,保证在去除噪声的同时尽可能保留信号的有效信息。
数字切比雪夫滤波器的matlab仿真是设计过程中必不可少的环节之一。
在matlab中,可以使用自带的函数cheby1、cheby2、chebwin 等来进行数字切比雪夫滤波器的设计和仿真。
同时,matlab还提供了各种可视化工具,如fdatool、filter designer等,方便用户进行滤波器性能分析和优化。
在数字信号处理应用中,数字切比雪夫滤波器被广泛应用于图像处理、音频处理、信号解调等领域。
其设计和优化方法对于提高数字信号处理的效率和精度具有重要的意义。
因此,深入研究数字切比雪夫滤波器的设计和优化方法,不仅对于学术研究有所帮助,也对于工程实践具有重要的指导意义。
目录摘要 (1)ABSTRACT (2)1 MATLAB简介 (3)1.1 MATLAB的概况 (3)1.2 MATLAB产生的历史背景 (3)1.3 MATLAB的语言特点 (4)1.4 MATLAB界面介绍 (5)1.4.1工作环境窗口 (5)1.4.2当前路径窗口 (6)2. 数字滤波器 (7)2.1数字滤波器的概念 (7)2.2数字滤波器的分类 (7)2.3数字滤波器的设计要求 (9)3. IIR数字滤波器的设计 (9)3.1设计原理 (9)3.2设计步骤 (10)4.切比雪夫IIR数字滤波器 (11)5.双线性变换法设计IIR数字滤波器 (12)5.1设计原理 (12)5.2双线性变换的优缺点 (14)6. 数字切比雪夫高通IIR滤波器设计 (17)6.1设计流程图 (17)6.2程序设计及运行结果 (18)6.2.1设计过程 (18)6.2.2运行结果 (20)7小结 (22)8参考文献 (23)摘要在现代通信系统中,由于信号中经常混有各种复杂成分,所以很多信号分析都是基于滤波器而进行的,而数字滤波器是通过数值运算实现滤波,具有处理精度高、稳定、灵活、不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊滤波功能。
数字滤波器根据其冲激响应函数的时域特性,可分为两种,即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。
实现IIR滤波器的阶次较低,所用的存储单元较少,效率高,精度高,而且能够保留一些模拟滤波器的优良特性,因此应用很广。
Matlab软件以矩阵运算为基础,把计算、可视化及程序设计有机融合到交互式工作环境中,并且为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。
尤其是Matlab中的信号处理工具箱使各个领域的研究人员可以直观方便地进行科学研究与工程应用。
本文首先介绍了数字滤波器的概念,分类以及设计要求。
接着又逐步介绍了IIR数字滤波器,切比雪夫IIR 数字滤波器,以及用双线性变换法设计IIR数字滤波器。
目录1 课题综述 (1)1.1 课题来源 (1)1.2预期目标 (1)1.3 面对的问题 (1)1.4 需要解决的关键技术 (1)2 系统分析 (2)2.1 涉及的基础知识 (2)2.2 解决的基本思路 (2)2.3 总体方案 (2)2.4 功能模块框图 (2)3 详细设计 (3)3.1 巴特沃斯低通滤波特性(MATLAB) (3)3.2 巴特沃斯高通滤波特性(MATLAB) (4)3.3 切比雪夫低通滤波特性(MATLAB) (4)3.4 切比雪夫高通滤波特性(MATLAB) (4)4 程序调试 (5)4.1 巴特沃斯低通滤波特性 (5)4.2 巴特沃斯高通滤波特性 (6)4.3 切比雪夫低通滤波特性 (8)4.4 切比雪夫高通滤波特性 (9)5 运行与测试 (10)5.1 选择音频文件(WAV) (10)5.2 滤波后音频特点 (10)6 全文代码设计 (10)总结 (14)致谢 (15)参考文献 (16)1 课题综述1.1 课题来源随着数字集成电路,设备和系统技术的快速进步,通过数字方法进行信号处理已变得越来越有吸引力。
大规模一般用途的计算机和特殊用途硬盘的高效性,已使得实时滤波既实用又经济。
目前主要有两类滤波器,模拟滤波器和数字滤波器,它们在物理组成和工作方式上完全不同,而模拟滤波器的技术发展已相当成熟,所以研究的重点基本上放在了数字滤波器上。
滤波器的功能是用来移除信号中不需要的部分,比如随机噪声;或取出信号中的有用部分,如位于某段频率范围内的成分。
目前随着计算机技术和数字信号处理器芯片的发展,使我们更为便利地识别和提取各种各样的信号。
因此研究不同数字滤波器的设计原理和稳定性分析对于满足军事、航空、民营等等各个领域的信号处理要求具有十分重要的意义。
1.2 预期目标能够理解并掌握无限脉冲响应数字滤波器(IIR)的机理,分析IIR滤波器的结构特性,观察IIR滤波器的频域特性。
能够通过MATLAB实现巴特沃斯,切比雪夫的高通低通滤波器的幅频响应,相频响应,以及零极点的图像的勾画。
matlab 切比雪夫滤波器设计ftool [matlab 切比雪夫滤波器设计ftool]引言:在信号处理领域中,滤波器是一种常用的工具,可以对信号进行去噪、平滑等处理。
而切比雪夫滤波器是一种具有优秀频率响应特性的滤波器,被广泛运用于信号处理和通信系统中。
本文将介绍如何使用MATLAB中的FTOOL工具进行切比雪夫滤波器的设计。
第一部分:MATLAB FTOOL工具简介FTOOL是MATLAB中的一种可视化滤波器设计工具。
它提供了一种直观且简便的方式来设计和分析滤波器,用户可以通过简单的拖拽和点击操作来完成滤波器的设计。
第二部分:切比雪夫滤波器原理与特性切比雪夫滤波器是一种频率选择性滤波器,它可以实现在频率域内对信号进行有选择地衰减。
切比雪夫滤波器的特点是在通带内具有较大的纹波,但可以通过调整滤波器阶数和纹波限制来进行折衷。
第三部分:切比雪夫滤波器设计过程在MATLAB中,使用FTOOL工具进行切比雪夫滤波器的设计非常简便。
下面将介绍具体的设计步骤:Step 1:打开FTOOL工具在MATLAB命令窗口中输入"ftool"命令,即可打开FTOOL工具。
Step 2:选择滤波器类型在FTOOL工具界面中,选择滤波器类型为"Chebyshev"。
切比雪夫滤波器有两种类型,一种是第一类切比雪夫滤波器(Type I),一种是第二类切比雪夫滤波器(Type II)。
选择相应的类型后,单击"Next"按钮。
Step 3:选择滤波器参数在切比雪夫滤波器设计中,有两个重要的参数需要设置,分别是滤波器阶数和纹波限制。
滤波器阶数决定了滤波器的复杂度和性能,而纹波限制则影响了在通带内的纹波程度。
根据实际需求设置相应的参数,并单击"Next"按钮。
Step 4:选择滤波器响应类型切比雪夫滤波器可以设计为低通、高通、带通或带阻滤波器。
在这一步中,根据自己的需求选择合适的滤波器响应类型,并单击"Next"按钮。
中北大学课程设计说明书学生姓名:学号:学院:信息与通信工程学院专业:电子信息工程题目:信息处理信息实践:信号的产生与带通滤波器的设计指导教师:职称:2016 年 1 月 4 日中北大学课程设计任务书2015/2016 学年第一学期学院:信息与通信工程学院专业:电子信息工程学生姓名:学生姓名:学号:学生姓名:学号:课程设计题目:信息处理信息实践:信号的产生与带通滤波器的设计起迄日期:2016年1月4日~2016年1月22日课程设计地点:学院楼201实验室指导教师:王玉下达任务书日期: 2016 年1月 4 日目录一基于USB总线A/D卡的报告 (1)B总线介绍 (1)B接口电路设计 (1)B接口的数据采集系统的设计实现 (1)4.A/D转换电路 (2)5.A/D采集卡的应用 (2)二设计内容 (3)三设计原理 (3)四 MATLAB软件介绍 (4)五设计过程 (5)1.信号的产生 (5)2.对连续输入信号进行采样,进行频谱分析 (7)3.切比雪夫带通滤波器 (8)4.对信号进行滤波 (10)六结果分析 (12)七心得体会及总结 (13)参考文献 (13)附录:总程序代码 (15)一基于USB总线A/D卡的报告B总线介绍通用串行总线标准USB是1995年微软、康柏、DEC、IBM等公司为解决传统总线不足的问题而推出的一种新型通信标准。
该标准具有高带宽、安装方便、易扩展诸多优点,逐渐成为现代数据传输的主要手段。
基于USB数据采集卡充拥有USB总线的上述优点,有效优化解决传统数据采集卡的安装麻烦、扩展不便等弊端。
USB总线为通用串行总线,USB接口位于PS/2接口和串并口之间,允许外设在开机状态下热插拔,最多可串接下来127个外设,传输速率可达480Mb/S,P它可以向低压设备提供5伏电源,同时可以减少PC机I/O接口数量。
USB是基于通用连接技术,实现外设的简单快速连接,达到方便用户、降低成本、扩展PC连接外设范围的目的。
B接口电路设计图1 电路设计如图1所示,R3是上拉电阻器,它可使USB口的D+端上拉到DS2490S的VB端,表示USB主机系统是高速设备,同时这个上拉电阻器告诉主机有USB设备插入。
该上拉电阻器的设置对适配器的影响很大,它的负载值和1-Wire网络的总长决定1-Wire总线电压上升到5 V的速度。
经过实验测试选择R3的阻值为27Ω±lO%。
R1、R2为USB数据线保护电阻器。
L、L2具有禁止高频干扰并且减弱EMI辐射的功能。
LF33CV为3.3 V电压稳压器,与周围元件C1、C2组成强上拉部分,给EEPROM或温度传感器等器件提供额外的电源。
B接口的数据采集系统的设计实现数据采集系统使用采集卡进行数据采集,然后经过A/D转换器供计算机加工处理。
基于USB接口的数据采集与频谱分析系统本系统结构由硬件部分和软件部分组成,硬件部分主要有计算机、I/O接口设备。
计算机作为硬件平台的核心可采用台式机,系统采用的I/O设备为A/D数据采集卡,该采集卡是一种基于USB总线数据采集产品,可与带USB接口的各种台式计算机、笔记本电脑、工控机连接构成高性能的数据采集测量系统.整个系统主要由4部分组成:USB接口芯片及外围电路、控制电路、数据缓冲电路和A/D转换电路。
USB接口芯片选择了Cypress公司的EZ-USB 2131Q,该芯片内嵌8051控制器,因此整个系统以EZ-USB控制器为核心,由EZ-USB经控制电路实现对A/D转换电路和数据缓冲电路的控制,模拟信号转换后的数据送入数据缓冲器,当数据缓冲器存满之后,通知EZ-USB控制器,由主机取出数据。
整个系统框图如图2所示。
图2系统框图4.A/D转换电路 A/D转换芯片采用了MAXIM公司的MAX122,芯片是12b的高速的A/D转换器。
在完全转换模式下,他的转换时间可以达到2.6μs,采样率为333kS/s。
MAX122有5种工作模式,在数据采集系统中,采用了模式2即连续转换模式。
在这种模式下,每次转换需要13~14个时钟脉冲节拍,转换可以不间断地进行,但是需要提供开始转换使能信号,并且要保证使能信号和时钟信号同步,读信号和片选始终处于有效状态。
数据输出使能信号一直有效,在转换结束时产生新的数据。
5.A/D采集卡的应用数据采集就是把来自各种传感器的信号数据实时地、准确地测量或汇集起来,用计算机进行实时处理或记录存储,实时完成测试和控制功能。
数据采集系统结构通过微机的标准接口连接各种功能模块、仪器仪表和传感器,组成测量系统。
数据采集监测技术已成为日益重要的检测技术,广泛应用于需要监控湿度、温度、压力等工农业中。
数据采集过程是工业控制等系统中的重要过程,通常主要采用一些功能独立的单片机系统来采集,是测控系统不可缺少的一部分,因为数据采集的性能特点直接影响到整个系统的稳定性。
在生产过程中这一系统可对现场的工艺参数进行监视、采集以及记录,从而提高产品质量,提供信息和手段,降低生产成本。
而在科学研究中应用该系统可以获得大量动态然后研究瞬间物理过程,也成为获取科学奥秘重要手段之一。
总之,无论在哪个领域,数据采集系统应用的越及时,工作效率就愈高,获得经济利益就愈大,这是毋庸置疑的。
二设计内容产生三种不同频率的正弦信号并叠加,通过FFT变换得到频域图;并设计在200-300Hz的带通滤波器,并进行频谱分析;再对合成信号进行滤波,对滤波前后信号进行对比。
三设计原理1.利用MATLAB软件产生三个连续信号,并对三个连续信号进行合成。
2.对信号进行抽样,进行频谱分析。
(1)时域采样(奈奎斯特采样)定理:为了避免产生混叠现象,能从抽样信号无失真地恢复出原信号,抽样频率必须大于或等于信号频谱最高频率的两倍。
本设计中信号最高频率是500Hz,抽样频率采用2000Hz。
(2)频谱分析:频谱分析是指对信号进行频域谱的分析,观察其频域的各个分量的功率大小,其理论基础是傅立叶变换,现在一般采用数字的方法,也就是将时域信号数字化后做FFT,可以得到频域的波形。
3.带通滤波器滤波的设计现代生活中,为了滤除谐波干扰,获得所需要的高精度的模拟信号,经常要用到滤波器对信号进行滤波。
典型的模拟滤波器有巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器和椭圆(Ellipse)滤波器等。
其中,巴特沃斯滤波器又叫最平坦响应滤波器,顾名思义,它的响应最为平坦,通带内没有波纹,其频率响应在通带和阻带中都是单调的,且在靠近零频处最平坦,而在趋向阻带时衰减单调增大,巴特沃斯响应能够最大化滤波器的通带平坦度。
该响应非常平坦,非常接近DC信号,然后慢慢衰减至截止频率点为-3dB,最终逼近-20ndB/decade的衰减率,其中n为滤波器的阶数。
切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快,但频率响应的幅频特性不如后者平坦。
切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。
本次课程设计将使用切比雪夫1型带通滤波器对信号进行滤波。
滤波器的结构框图如下图3所示:图3 滤波器的结构框图四 MATLAB软件介绍MATLAB软件是由美国Mathworks公司推出的用于数值计算和图形处理的科学计算系统环境。
在MATLAB环境下,用户可以集成地进行程序设计、数值计算、图形绘制、输入输出、文件管理等各项操作。
MATLAB提供了一个人机交互的数学系统环境,该系统的基本数据结构是矩阵,在生成矩陈对象时,不要求作明确的维数说明。
与利用C语言或FORTRAN语言作数值计算的程序设计相比,利用MATLAB可以节省大量的编程时间。
MATLAB因其提供计算、图形、文字处理的统一环境而深受欢迎。
MTALAB系统由五个主要部分组成,下面分别加以介绍。
(1)MATALB语言是高层次的矩阵/数组语言。
具有条件控制、函数调用、数据结构、输入输出、面向对象等程序语言特性。
利用它既可以进行小规模编程,完成算法设计和算法实验的基本任务,也可以进行大规模编程,开发复杂的应用程序。
(2)MATLAB工作环境这是对MATLAB提供给用户使用的管理功能的总称。
包括管理工作空间中的变量据输入输出的方式和方法,以及开发、调试、管理M 文件的各种工具。
(3)MATLAB图形系统的基础,包括完成2D和3D数据图示、图像处理、动画生成、图形显示等功能的高层MATLAB命令,也包括用户对图形图像等对象进行特性控制的低层MATLAB命令,以及开发GUI应用程序的各种工具。
(4)MATLAB数学函数库这是对MATLAB使用的各种数学算法的总称。
包括各种初等函数的算法,也包括矩阵运算、矩阵分析等高层次数学算法。
(5)MATLAB应用程序接口(API)这是MATLAB为用户提供的一个函数库,使得用户能够在MATLAB环境中使用C程序或Fortran程序,包括从MATLAB中调用于程序(动态链接),读写MAT文件的功能。
可以看出MATLAB是一个功能十分强大的系统,是集数值计算、图形管理、程序开发为一体的环境。
五设计过程1.信号的产生首先产生一个连续输入信号,包含低频(f=30Hz),中频(f=250Hz),高频(f=500Hz)分量,并实现三个信号的合成。
(1)程序代码f1=30;f2=250;f3=500; %产生三种频率t=(1:100)/2000;x1=sin(2*pi*t*f1);figure(1);subplot(2,1,1);plot(x1); %绘制x1(t)的图形xlabel('t');ylabel('x1(t)');title('低频信号');grid;x2=sin(2*pi*t*f2);subplot(2,1,2);plot(x2); %绘制x2(t)的图形xlabel('t');ylabel('x2(t)');title('中频信号');grid;x3=sin(2*pi*t*f3);figure(2);subplot(2,1,1);plot(x3); %绘制x3(t)的图形xlabel('t');ylabel('x3(t)');title('高频信号');grid;x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3);%实现三种信号的叠加subplot(2,1,2);plot(x); %绘制x(t)的图形xlabel('t');ylabel('x(t)');title('合成信号');grid;(2)程序运行结果如图4、图5所示图4 低频中、频信号波形图图5 高频、合成信号波形图(3)结果分析图4、5所示内容为按照设计要求作出的三种频率分量以及合成信号。