基于遗传算法的复杂无源滤波器参数设计

　　第5步,单次过程结束判断

　　当F

　　当F≤Er且gen≥MAXGEN时,gen=gen+1,返回第3步.

　　当F≥Er时且gen≥MAXGEN时,返回第2步.

　　第6步,全过程结束判断

　　dcnt=dcnt-1;当dcnt>0时,返回第2步;否则,停机.

　　3 数值实验例

　　图2为一带通无源滤波器电路结构,通频带要求在950～1050 rad/s之间.为此,每隔5 rad/s作一次采样,采样点的幅度大于0.85;设定低频截止频率为800 rad/s,幅度小于1e-5;高频截止频率为1300rad/s,幅度小于1e-5.建立如下优化模型:

　　s.t. XL=[0,0,...,0]

　　其中

　　X=[X1,X2,...,X19]=[L,C3,C4,L6,L7,C7,L9,L10,C12,L13,C13,L15,L16,C18,L19,C19,L21,L22,C24];

　　R1=10×Min(+[1.0e-5-L1]-0.0,0),对应W1=800rad/s

　　Rj=1.0×Min(+[Lj-0.85]-0.05,0);j=2,3,...,22,对应Wj=(5×j+940)rad/s

　　R23=10×Min(+[1.0e-5-L23]-0.0,0),对应W23=1300rad/s

　　在NEC4800/210Ⅱ工作站完成上述算法.算法的参数设置为Er=1e-6,α=0.1,Pc=Pc(0)=0.6,Pm(0)=0.1,POP_SIZE=40,MAXGEN=2000,程序语言为UNIX-C.dcnt取10,得到10组设计值,皆能使频率特性满足要求.平均世代数为1508代,平均时间为9.8min.其中一组结果为:

　　X*=[0.079212.65140.075213.01570.10583.3131

　　0.179313.33860.07260.13343.3260.172215.12180.06330.08761.92880.333310.31710.0900].

　　按此参数设计后,滤波器频率特性较好地达到了预期要求,如图3所示.

　　本文简要分析了无源滤波器参数设计存在困难的主要原因.对无源滤波器的常用电路结构,提出了计算频率特性的简易迭代法,并将求解满足指标要求的参数设计值的问题转化为优化模型的求解过程,使的原来难以描述和解决的设计问题变得明确和简单.

　　在遗传算法方面,我们并不照搬前人的方法[3],例如,把最小目标函数的求解要求转化为进化的驱动力而不是刻意求得最优解,不仅减少了计算时间,还可提供设计者多种可选择的方案.在遗传算法的变异阶段,采用的是以现有的基因为中心的左右变异方法,有别于常用的在变量范围内变异的方法.因本法中心不断灵活飘移,求解过程出现两个极端,即快速找到合适解或走入死区——永远找不到解.我们通过设定最大世代数消除了后者的可能性,总体上提高了效率.

　　本法有效地克服了无源滤波器参数设计的困难,具有普遍性的价值和意义.其基本思想及原理亦可在电气、电子、自动化等较复杂系统的参数设计领域得到广泛的应用.