matlab在一阶无源带阻滤波器的应用.docx

书目1设计相关理论311串联谐振电路3电路模型分析3电路谐振条件4电路频率响应4电路品质因数512无源滋波器8低通滤波器8高通注波器9常通途波器9常阻逐波器92设计内容102.1原理分析112.2数学建模112.3卷数推导122.4可行性分析132. 5程序设计14程序框图14程序代码15函数说明及运用说用163设计结果213.1程序运行结果和图表213. 2结果分析233. 3设计结论234心得体会24参考文献26摘要本次课程设计是探讨一阶带阻漉波器，以UO为响应，求域率响应函数,画出其幅频特性和相频特性曲线.首先，通过对波波器背景理论的学习与理解，对设计题目进行分析并做出数学模型，然后借助Mauab软件，画出选目所要求的幅频、相频图，并通过调整和设计参数画出/比较合理的图形.在设计中，我深刻理解到J'波波器的工作原理，而且发觉志向和实际中的曲线有不小差距，与此同时，也娴熟MaUab的相关操作和应用，对信号理论和R1.C串联谐振有更新的相识。最终，本文对此次课程设计中的收获和体会进行了总结.关健词：滤波微幅频相频Matlab串联谐振一阶无源带阻滤波器的分析和设计1设计相关理论1.1 串联谐振电路电路频率响应的最全要的特征是其在幅度特性上所呈现的峰值点（或尖峰点、谐振峰值点谐振的概念应用于科学和工程的多个领域之中。任何有更共鞭极点对的系统都会产生谐振，这是振荡产生的根源。谐振峰值的现象用在通信网络中可以进行频率识别。在至少有一个电容和一个电感的任何电路中都可能产生储能由一种形式到另一种形式转换的谐振振荡。谐振是R1.C电路中的一种状态，该电路中电容和电感的电抗大小是相等的，结果呈现出纯电阻的阻抗性质。串联或并联谐振电路的传递函数有很高的频率选择性，所以在设计和制作滤波器过程中是很有用的，其他很多应用中包括收音机的选台和电视机的选频道等。1.1.1 电路模型分析考虑如图所示的串联R1.C频率电路.其输入阻抗是:Z=H()=R+j1.+7IjC1.1.2 电路谐振条件若传递函数的虚部为零,则产生谐振，即:Im(Z)=1.=OC满意上述条件的值称为谐振频率，因此谐振的条件是:o°l=1o=rads又因0=2f0,所以:在谐振条件下，有：1、阻抗是纯电阻Z=R.换宫之，1.C的串联组合相当于短路，整个电压都加在电阻R上.2、电压VS和电流是同相的,所以功率因数为13、传递函数H(3)=Z(3),其运值最小.4、电感上电压和电容上的电压可能比原电压高得多.可由下列关系证明：VmV,=-01.=QVmIVCl="TTTF=QVmc式中Q是品质因数.1.1.3 电路频率响应RI.C电路的电流的频率响应是:=ll=可+31.-嬴其图形如图所示。该图说明白频率轴为对数坐标时的I大小的对称特性,1.1.4 电路品质因数R1.C电路所消耗的平均功率是:P()=I2R谐振频率时,电路消耗最大功率是:/、v'>2P(0)=rr令3=5,32时，消耗的功率是最大功率的一半，即:P(J=P(2)=2R4R则称Sl.32为半功率(点)频率.半功率频率可以通过令Z=MR求得.即:解s,得至小32=t+fe）+由此可以得到半功率点频率与谐振频率的关系是:0='2即谐振频率是半功率领率的几何平均值。一般来讲，频率响应并不对称于谐振频率，所以巴，32也不是对称于°的。但是，认为半功率点频率对称于谐振频率是一个比较合理的近似。图中谐振曲线的峰值取决于电阻R,而该曲线的宽度取决于其频带宽度B.带宽定义为两个半功率点频率之差：B=2-j这种频带宽度的定义是几种常用的定义之一。严格地说，上式所定义的带宽称为半功率带宽，是半功率点频率之间的谐振曲线的频带宽度。谐振曲线的“锐度”用品质因数Q这个量来度量。电路谐振时，电路中的电抗能量在电感和电容之间来回振荡。因质因数将电路储存的很大（峰值）能量与每振荡个周期所消耗的能量间的关系连系起来，定义为：C电路存储的峰值能量Q=2”，振荡一个周期电路所消耗的能曼品质因数也是电路中储能的性能与耗能性能之间关系的一个度址。在串联R1.C电路中,储能的峰值是1.22,一个周期的耗能是：（-R）2f,所以；2t,2211f1.Q=211.=2R（l0或（OO1.1Q=T=砌品质因数是无量纲的值，带宽B和品质因数Q之间的关系由将式代入,并利用式关系得到:R0BTQ或者B=02CR谐振电路的品质因数是其谐振频率与带宽之比.注：上述公式R(>qCR只适用于R1.C串联电路,电路的Q值越高，其带宽越小如图所示，Q值越高，电路的频率选择性越强，而其频带也越窄，R1.C电路的选频特性是电路对每个频率的响应，并解除其它频率的一种实力，假如要选用或要排斥的频带很窄品侦因数必需很高，反之若频带较宽，则谐振电路的品质因数应当低些。品质因数也是电路电路性（或谐振“锐度”）的i个度量。通常谐振电路总是工作在其谐振频率处或其邻近处，当电路的品痂因数等于或大于10时,称为高Q电路。在高Q电路（Q10）的全部实际应用中，其半功率频率都可认为是对称于谐振频率的且可挖地按下式表示：BB1*0-y2g>0+高Q电路常常在通信网络中运用.综上所述，谐振电路可以用以下五个相关的参数来表征：两个半功率点的频率<山和谐振频率3。，带宽B和品质因数Q.1.2 无源滤波器戏波器的概念从起先就始终是电子工程进展中的一个主要组成部分，没有电子泄波器的参加，某些技术成果将是不行能的.戏波器是一个被细心设计的电路，它只允许一些特定频率的信号通过,而阻挡或衰减其他频率的信号。作为一个减率选择装置，灌波器可以用来将信号的顼谱限制在某个指定的频带宽度范国内，在无线电接收机和电视机中，从空间许很多多广播信号中选出所须耍的信号频道所用的电路就是泄波器.假如战波器电路的组成元件只有无源元件R.1.和C,则称为无源泄波器：假如除J'无源元件R.1.和C外，还有有源器件（例如晶体管、运算放大器等），则称为有源造波器.本文探讨无源浊波器.1.2.1 低通泣波器低通泄波器是只通过从直流到截止频率3。的泄波器.一个RC电路取其电容器上的电压作为输出，就构成一个典型的低通灌波器.一个R1.电路以电阻两端作为输出也构成了一个低通港波港。1.2.2 高通滤波的高通泄波器是通过全部高于其截止领率3。的戏波器“一个RC电路以电阻两端作为输出，就构成r一个高通港波淞.，若Rl电路以电感两端作为输出也构成一个虑通/波器。1.2.3 带通滤波器带通泄波器是通过频带<1<<2;中全部缴率的滤波器，R1.C串联谐振电路以电阻两端为输出，就是-个带通灌波罂。12.4带阻淖波器个波波罂阻挡两个指定值(O)IflIO2)之间的频带通过，则称之为带阳滤波器、带止滤波罂或陷波罂。当个串段R1.C谐振电路取其1.C串联的组合为输出时；即构成了带阳滤波器。Cr如图所示，其传递函数是：Uoj(1.-lC)HG,D=.=R+j(31.3C)由式可见:H(O)=I,H()=l图所示为H()的幅频特性，其中心频率为:t0°=衣同祥，4波器的半功率频率、带宽和品质因数，仍旧可以用上述串联谐振电路的公式来计算,例如:R21.+S=+O+R0B=-=1.Q这里的。称为抑制预率，而对应的带宽(B=2-3)称为抑制带宽II带阻戏波器是抑制或消退在频率<3<32内全部频率的泄波器.留意：对于相同R1.和C值的带通泄波器的传递函数与带阳泄波器的传递函数相加得到的结果是在任何频率下都为1.这是因为这两个电路其中一个电路的特性是另一个的反特性，当然，一般状况下，上述结果是不成立的。无源戏波器的最大增版是1，耍想得到大于的增益，应当用有源滤波器.2设计内容本次课程设计的内容是探讨一阶带阻波波器。如图2.1,以UO为响应，求频率响应函数，画出其幅频响应(幅频特性)H(3)和相频的响应(相频特性)0()o图2.12.1原理分析如图2.1.一个申联R1.C谐振电路取其1.C串联的组合作为输出，构成J'带阻滤波器.其传递函数，即频率响应函数是：UOj(1.-lC)H()=-=R+j(1.-itoc)中心频率为：°=K带宽B和品质因数Q之间的关系：R0B=1.Q2 .2数学建模频率响应函数：Uoj(1.-lC)H3)=R+j(1.-lC)令x=)1.-±,则：jxjRx+X2HR=西丁南寿×iRx=+iR2+x2R2+X2于是,有:x2(R2÷X2)(R2+X2)2故我低频响应(幅频特性)为:3 =昌*”相频响应(相频特性):,R8()=tan1-1 R=tan'1-2.3参数推导经过反复尝试和推算，分析得出当R=In,1.=WmH,C=lF时得出的曲线是比较合理的。由:求得半功率点频率:1=951.249219725039rad/s2=1051.24921972504radso电路的谐振频率：3。=表=100°radst,带宽：B=（i）2-<<）=100rads.品质因数：W0Q=R=Uh2.4 可行性分析由于Q=10,所以此电路属于高Q电路，通常用于通信电路。Q表征一个储能器件（如电感线圈、电容等）、谐振电路中所储能量同每周期损耗能量之比的一种质量指标：电抗元件的Q值等于它的电抗与其等效串联电阻的比值；元件的Q值愈大，用该元件组成的电路或网络的选择性愈佳。Q值越大，曲线越尖锐，通频带越窄，电路的选择性越好。在恒压源供电时，电路的品质因数、选择性与通频带只确定于电路本身的参数，与信号源无关。提高品质因数可以选用带有磁心的电感线圈。电感线圈中带有磁心时，可使线圈圈数及我电阻大大削减，有利于Q值的提裔。另外，还可以采纳依据工作频率选择绕制线圈的导线、选用优质竹架，削减介质损耗、合理选择屏蔽罩的尺寸等等一系列措施来提高品质因数，提海电路性能。对于电路所用到的器材，电阻、电感、电容，都是常见值，比较好找。而且电路结构简洁，简洁搭建。但是由于这个波波器不含放大电路，放大倍数最大就是1.不能对有用的信号加以放大处理，这是一个缺陷。它只能把谐振频率旁边的频段阻碍掉，留下有用的频率。总体来讲，作为一个一阶无源带阻灌波器，这次设计还是可行的。可以做出成品，并且电路拥有相应的功能，而且还比较稳定。但从波波器的角度看，也存在不足，比如：能放大、功能过少、好用价值低等等。若仅仅作为设计还是可以的，假如要用于工业生产和运用，本设计还得加以改进和调试，方可更经济更有效的运用.2.5 程序设计Mablab是矩阵试验室(Matrix1.alxxatory)的简称，是美国MalhWOrkS公司出品的商业数学软件，用丁尊法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括Mablab和Simulink两大部分.MabIab是由美国nuthworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境“它将数值分析、矩阵计克、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于运用的视窗环境中，为科学探讨、工程设计以及必需进行有效数值计算的众多科学领域供应了一种全面的解决方案，并在很大程度上摆脱，传统非交互式程序设计语言(如C、FOEan)的媪辑模式，代表了当今国际科学计兑软件的先进水平。2.5.1 程序框图2.5.2 程序代码»figure,x,y=fplot("cdtFi*10-2)*(l(ix*10-4)(l÷i*x*10'2)÷,l,i*x*104l1,0,2000,1¢-4);size()fplot<*abs"i*10*(-2*<l<i*10(-4)")(l*i*x*10(-2+<l<i*10*(-4)'S,2000,le7)场福姨特性曲线fisure,x,y=fplot(,angle(i*x*10-(-2H<l(i*x*10'(-4.)(Hi*x10'(-2>+<l<i*xM0(-4)',0,2000,le-4).size(x)fpl"('an8h"i*x*10>yXl(i*x*l<(7)"(Ni*x*l(-23"(i*x*10l-4>)>)',0,2000,k7)、画相涮性曲线ans=19201sns=5691»X=SoIVe('(0.01*-l(*ll(7)2)=T)%利用半功率频率求解明代入R1.C值并化简，解方程1051.2492197250392863848606074161-951.24921972503928638486060741613951.24921972503928638486060741613-1051.2492197250392863848606074161»*l=95l.24921972503928638486060741613¥2=1051.2492197250392863848606074161B=v2-wl*求出半功至巨角颜至和阳帚宽Vl=951.24921972503910S1.2492197250499,9999999999999»B=IO0«带亮妁为10。VO=Isqrt(10(-2)*(l(T)力求出谐方领军Q=W/B%求出昆后因效B=100W=100O»x=O:l:2000,yl=l.<t(x<wl)+l.*(x>w2);y2=abs(i.*x.*10(-2)*(l.(i.*x.*10H).(l+i.*x*10(-2)*(l.(i.+x+10(*4);plot(x,yl,x,y2)遍出理想医和实除图2.5.3iablab曲数说明及运用说明(1) fpk>t函数函数功能：在指定的范围内绘制函数图像运用方法：fpot在指定的范困内绘制函数图像，函数必需是y=f(X)的形式，其中X是个指定范闱limits的向量，y是和X有相同大小的向盘并包含在点X处的值。假如对个给定的x值.函数返回多于一个值，则y是每列包含f(X)的每个重量的矩阵。fplot<fun.limits)在指定的范围limits内画出函数名为fun的图像。其中IimitS是个指定X轴范闱的向AUxniinXmax或者是X轴和y轴范围的向员xminXmaXyminymax.fun可能为：1、M文件函数名字。2、可能传递给eval函数的带变能X的字符串,比如，Sin(x)'.'diric(x,10),'sin<x).COS(X)'.3、M文件函数或匿名函数的函数句柄。对于向员X的每个元素，函数f(X)必需返回个行向量。比如，假如f(x)返回向SUfl(X)，(x>,(X).输入参员为X=IXl;x2),则函数f(X)返回矩阵：fl(xl)f2(xl)B(xl)fl(x2)f2(x2)(x2)例如：fplot(fun.limits.1.ineSpec)用指定的线型1.ineSPeCI面出函数fun。fplot(fun.limits,tol)用相对误差值为col画出函数fun。fplot(fun.limits,tol.1.ineSpec)用指定的相对误差值为tol和指定的线型1.ineSPeC画出函数fun。fplot(fun.Iims.)允许可选参数tol.n和1.ineSpec以随意组合方式输入。fplot<axes_handle.)用指定句柄axes.handle代替当前坐标轴句柄来画图。X.Y=fplot(fun>limits.)返回懂坐标与做坐标的值赋给X和丫，此时fplot不给出图形，若想画山，可用吩咐plot(X.Y)。应用举例：画1-2,2区间的双曲正切函数：fnch=(3)tanh;fplot(fnch,(-2,2)(2) PIot函数函数功能：plot是绘制二维图形的最基本函数，它是针对向量或矩阵的列来绘制曲线的。运用方法：运用PlOt函数之前，必需首先定义好曲线上每一点的X及y坐标，常用格式为：1. plot(x)：当X为向室时,以X元素的值为纵坐标，X的序号为横坐标值绘制曲线。当X为实矩阵时，则以其序号为横坐标,按列绘制每列元素值相对于其序号的曲线,当X为mxn矩阵时，就由n条曲线。2. plot(x.y)：以X元素为横坐标值，y元素为纵坐标值绘制曲线。3、plot(X,yl.X,y2.)以公共的X元素为横坐标值，以yl,y2,元素为纵坐标值绘制多条曲线。(3) ezplot函数函数功能：在Mablab的吩咐窗口中键入helpezplot吩咐或者docezplot即可获得本函数的帮助信息。ezplot即easytousefunciionplotter(它是个易用的元函数绘图函数,特殊是在绘制含有符号变应的函数的图像时，ezpbt要比plot更便利。因为plot绘制图形时要指定自变室的范闹，而ezplot无需数据打算，干脆绘出图形。运用方法：ezplot的调用格式:1、ezplot(fun)2、ezplot(fun.min.max)3、ezplot(fun2)4、ezplot(fun2.xmin.xmax.ymin.ymax)5.ezplot(fun2.min.nax)6、ezplot(funx.funy)7-ezplot(funx.funy,tnin.tmax)8、ezplot(.figure-landle)9. ezplot(axesIandle.)10»h=ezplot<.)(4) solve函数函数功能：主要是用来求解线性方程组的解析解或者精确解。对于得出的结果是符号变量，可以通过VPaO得出陋意位数的数值解，运用方法：solve函数的语法定义主要有以下四种：solve(eq)solve(eq,var)solve(eql.eq2eqn)g=solve<eql.eq2.»eqn.vail,var2.,varn)eq代表方程，Var代表的是变量。(5) abs函数函数功能：求发数实部与虚部的平方和的算术平方根。运用方法：格式：abs(x)例如：x=l+jsyabs(x)：»>=1.4142(6) angle函数：函数功能：求复数的相角。运用方法：格式:angle<x)3设计结果3.1程序运行结果和图表图3.11帕顿特性曲线图3.12相频特性曲线图3.13帝阳漉波器的志向和实际频率晌应3.2结果分析从幅频特性曲线上能够看出，当角频率从Orad/s向100orad/s渐渐增大的过程当中，起初，曲线有很长一段时间趋于平缓,下降幅度很小,H(3)的幅值接近于I。当角猱率渐渐接近F半功率点频率100Orad/s的时候，曲线起先快速下降，始终到谐振频率(,>o=10O0rads时，H(3)幅值变为0。接者，曲线在o=100rads以右的旁边区间内又急剧上升，直至出(3)|的幅值接近丁1,然后曲线又犹如起先阶段那样变得平缓，上升幅度也明显减小。明显视察出在角频率从Orad/S起先到2O00rads的范围内，曲线大致关F3。=100orad/s对称。通过幅频特性，可以发觉这个电路对100Orad/s的角频率有很好的阻挡作用，带宽B=100rad/s很率，品质因数Q=F=Io很高.电路的阻挡实力很强.再看相频特性曲线，在从Orad/s达到谐振频率“>o=ooorads以前，曲线先从零起先缓慢下降，趋势平缓，然后在快达到中心频率旁边又急剧下滑到谐振点处。(3)降为最低，接着0(3)突然跳变到正最大值，在/。=100orad/s以右的旁边区间内又快速下降再几乎以与刚起先阶段相同的趋势缓慢下降，。(3)及终趋近零。这反映出电路当角菽率从OraWS到谐振频率o=100Orad人期间呈现出容性，在谐振顼率他)=10rad.'s以后呈现出感性，在谐振点o=100OraWs的时刻，相角为零，该电路中电容和电感的电抗大小是相等的,结果显示出纯电阻特性，这就正好与幅频特性曲线对应起来.再在后面的志向图与实际图中，我们也可以看出在阻带带宽中，幅频曲线快速下降，达到了阻挡对应频率的特性。另外，得出电路的中心频率<。0=100OraWs,半功率频率3=951.249219725039rads,2=1051.24921972504rad/s3.3设计结论本次美干一阶滤波器的课程设计做出了个网挡频率特性比较好Il图像比较合理的无源带阻波波器，对对应频率有良好的网挡效果。它有着较高的Q值与较窄的带宽B,实现r对相应谐振频率的阻碍作用.但也存在不足，就是这种泄波器只能阻碍肯定的频率，却不能对有用的频率进行放大处理。因此，无源/波器主要有三个方面的限制：第一，增益不能大于1,无源器件不能增加网络的能量：其次，可能须要笨重和较费的电感元件：第：,在低丁音频范围工作时，戏波器性能不好.可见，与有源滤波器相比，这种一阶无源滤波器还是有很大的不足的。不过，这种泄波器电路简洁，原理易懂,简洁搭建和集成,所需元器件较少，在高频工作时，无源滤波器是很有用的。除了滤波器，在设计中借助了软件Mablab.Mablab作为线性系统的种分析和仿真工具，是理工科高校生应当驾驭的技术工具，它作为种编程语言和可视化工具，可解决工程、科学计算和数学学科中的很多问题。通过应用，我发觉这个软件对代码编写的格式要求比较严格，略微不对就得不到正确的结果。所以，在输入代码的过程中应当特殊当心留意，要段段的、一块块地编写、运行，以免最终整个程序运行不出来,这样可以避开因某些小地方出错导致最终花费大量时间重笈忙乱检查而影响效率，不简洁改正。最终，就是课程设计应当不断地修改和学试，再加上不断地学习和实践，如此反发才能谢试出正确的、合理的图像,得到比较完满的结果，才能让我们感觉到自己真正的学到学问。4心得体会经过为期周的暑假基础强化训练，我最终完成了我的课程设计。此课程设计是对我们专业学问的次考察，特殊是对我们电路原理和数字电了技术基础相关理论驾驭程度的考验。在收获学问的同时也收获了成熟。在课程设计的过程中，我通过查找大量资料，清教同学，再加上自己不褥的努力，不仅培育了自身独立思索、动手操作的实力，在和同学探讨问题时的沟通沟通实力上也有了提高。更重要的是，在此期间，我学会了很多学习的方法，而这乂怡恰是在今后的学习工作生活中最珍贵最好用的东西，让我受益匪浅。要面对社会的挑战，只有不断的学习、实践，再学习、再实践。我选择的是第21题一一关于阶带阳谑波落的设计，题目要求求出频率响应函数.并画出电路的幅频特性和相领特性曲线。确定好目标后便起先查找各种相关资料。我通过网络、书籍、课本、笔记找到了很多有用的资料，为设计电路打下了坚实的基础。因为要设计好个电路必需了解跟它相关的各方面学问。在这个过程中，我学会了Mablab的运用。它真是个强大的软件，让数学计算画图变成编程！经过查阅它的运用资料和观看相关运用视频，我驾驭了不少Mablab的吩咐，熟识了它的殷操作，而11能画出二维函数图像。在不断的摸索中，我渐渐学会在不了解个新事物新东西的时候，如何上网查阅资料去学习和运用它的方法，为我供应更多更好的思路和思维方式，拓宽了我的视野。经过以上充分的打算，我已积累了不少资料，对所选的题目也也许仃些了解，然后就是在这个基附上，综合已有的资料来更透彻的分析题目。经过多次调整和尝试画图，并加以计算,我最终找到了比较合适的参数。在此之上，我用MabIab画出了比较合理的图形，求出了方程组的解。接下来，我乂求出表征此电路的几个重要相关参数，以此来评价我设计出的泄波器的质量。招意到在编写源代码的时候最好是编写个个的小模块分别进行调试，这样可以避开因某些小地方出错导致最终花物大量时间重纪忙乱检查而影响效率。渐渐地，就这样一步步地做下去之后，我突然发觉要做出整个程序来并不难，但是肯定要细心和耐性地每完成一小块儿就检查并修正，如此不断的修改调试，再修改再调试，就做出了正确完整的程序，虽然这个滤波罂设计的还不是很完备。当把整个任务完成后，我细细思索，才发觉，我们课堂上所学的理论学问有多么重要,而敏搬运用所学学问也同样重要，感到困雄的缘由是思想没有完全打开，自己对学问的驾驭程度不够，娴熟运用的实力还相当欠缺。我知道像这样的课程设计常常有，所以每次有这样的机会时，肯定要好好利用，争取每次都要更上层楼，平常也要多思索，石新的想法时可以记录下来，留意身边的电器元件，报长思索，这才是我们课改的真正目的。总之，不管起先你会的多少，你都会觉得困难比较多，没有什么头绪，真是万事开头难.不知道如何入手。我还得出个结论：学问必需通过应用才能实现其价值！有些东西以为学会了，但真正到用的时候才发觉是两回事，所以我认为只仃到真正会用的时候才是真的学会了。参考文献1于歆杰.电路基础.第三版.清华高校出版社，20082康光华.电子电工基础.高等基化出版社,20063赵光宙.信号分析与处理B.其次版.机械工业出版社，20064孙祥,徐流美，吴清.MAT1.AB7.0基础教程.第一版.清华高校出版社，20055王正盛.MAT1.AB与科学计算.第一版.国防工业出版社，2011