野草算法优化提取广义切比雪夫滤波器耦合矩阵
合集下载
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
收稿日期: 2 0 1 3- 0 5- 2 5 作者简介:田雨波( 1 9 7 1 —) , 男, 教授, 博士, 研究方向为群集智能技术和计算智能技术研究及其电磁应用. E m a i l : t i a n y u b o @j u s t . e d u . c n
能优化方法得到任意拓扑结构的耦合矩阵. 文献 [ 3 ] 给出单纯形法和梯度法相结合的方案来提取 4 ] 采用粒子群算法优化提取广义 耦合矩阵, 文献[ C h e b y s h e v 线性相位滤波器耦合矩阵, 以上两种方 法均为直接优化滤波器的传输函数之法, 因而对初 始值的选取比较敏感. 野草算法( i n v a s i v ew e e d s o p t i m i z a t i o n , I WO ) 是 一种模拟杂草繁殖过程的新型的随机: I n v a s i v e w e e d s o p t i m i z a t i o n( I WO )a l g o r i t h mi s i n t r o d u c e dt o e x t r a c t t h e c o u p l i n g m a t r i x o f t h e g e n e r a l i z e dC h e b y s h e vf i l t e r .F i r s t l y ,t h eI WOi s u s e dt oc o n t r o l t h er o t a t i o na n g l ei nt h ep r o c e s s o f c o u p l i n gm a t r i x s i m i l a r t r a n s f o r m a t i o n .B y o p t i m i z i n g t h e o b j e c t i v e f u n c t i o n ,z e r o s u p p r e s s i o no f c o u p l i n g m a t r i x e l e m e n t s a t a r b i t r a r yp o s i t i o n s i s r e a l i z e d .T h e n ,s e v e r a l t o p o l o g i e so f f i l t e r a r eg i v e na n dt r a n s m i s s i o nz e r oi ss e t t oa n yf i n i t e p o i n t . T h ee x t r a c t i o no f c o u p l i n g m a t r i x d e m o n s t r a t e s t h e f e a s i b i l i t y o f t h e p r o p o s e dm e t h o d . F i n a l l y , b a s e do nt h e e x t r a c t e dc o u p l i n gm a t r i x , m i c r o w a v ef i l t e r i sf a b r i c a t e d ,a n dt h et e s t i n gr e s u l t sa r ei na g r e e m e n t w i t ht h ed e s i g n e do n e s . K e yw o r d s :i n v a s i v ew e e d s o p t i m i z a t i o n ;c o u p l i n gm a t r i x ;f i l t e r 通信系统对微波滤 随着现代通信技术的发展, 波器的性能指标及体积大小提出了越来越严格的 广义切比雪夫滤波器通过引入有限频率的传 要求. 输零点, 在不增加滤波器阶数的同时, 能有效提高 滤波器的性能, 因此被广泛应用于各种微波通信系 统中. 文献[ 1- 2 ] 中分别提出了 N阶及 N+ 2阶耦 合矩阵的综合理论, 但是提取出来的耦合矩阵对应 的拓扑结构很难直接用于工程实际中, 需要根据实 际的拓扑结构对其进行相似变换. 由此可见, 如何 提取出可行的耦合矩阵是广义切比雪夫滤波器设 计的重点及难点. 传统的矩阵相似变换方法只能适 用于特定的拓扑结构, 很难推广使用, 故而采用智
T 置矩阵 R 来实现的: T M1 =R ·M0·R
( 3 )
式中: M0为初始矩阵, M1为相似变换后的矩阵. 6 阶 旋转矩阵 R ( 支点为[ 3 , 5 ] , 旋转角度为 θ )如式( 3 ) 定义, 其中 c o s ( ) , s i n ( ) . θ θ r =c r =s
1 0 0 R = 0 0 0
对称零点结构归一化幅频响应fig8normalizedamplitudefrequencyresponsesymmetriczerosstructure从以上4例可以看出对任意的拓扑结构使用野草算法来优化提取耦合矩阵经过迭代都能得到良好的收敛效果实现了较快的收敛速度和较好的收敛精度并且优化结果不依赖于算法初值的选取预先设定的传输零点的位置与所提取出的耦合矩阵的传输函数的零点位置基本吻合得到耦合矩阵对应的广义切比雪夫滤波器在物理上是可实现的
n
图1 野草算法流程 F i g . 1 F l o w c h a r t o f I WO
2 广义切比雪夫滤波器耦合矩阵的相 似变换
文献[ 1 ] 详细介绍了综合法提取广义切比雪夫 滤波器耦合矩阵的方法. 综合方法得到的广义切比 雪夫滤波器耦合矩阵一般来说无法直接使用, 因为 耦合矩阵的元素都是非零元素. 由于耦合矩阵的特 征值决定了滤波器的传输和反射特性, 只要保证耦 合矩阵的特征值不变, 滤波器的频率响应特性就不 会改变. 根据矩阵理论, 对矩阵进行相似变换, 即在 不改变矩阵特征值的前提下, 对矩阵的元素作消零 处理. 一个 N阶耦合矩阵 M0 的相似变换, 是通过在 M0 的左边和右边分别乘以 N阶旋转矩阵 R及其转
第2 7卷第 6期 V o l 2 7 N o 6 江苏科技大学学报( 自然科学版) D e c . 2 0 1 3 2 0 1 3年 1 2月 J o u r n a l o f J i a n g s uU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y( N a t u r a l S c i e n c eE d i t i o n ) d o i : 1 0 . 3 9 6 9 / j . i s s n . 1 6 7 3- 4 8 0 7 . 2 0 1 3 . 0 6 . 0 1 3
5 ] 优化算法 [ . I WO算法自 2 0 0 6年提出以来, 得到
了广泛的应用, 并显现出强大的优势和潜力, 在电
6 - 1 0 ] 磁学领域内亦获得了应用 [ . 文中采用野草算
优化旋转矩阵的角 法, 在电路综合理论的基础上,
5 7 0
江苏科技大学学报( 自然科学版)
2 0 1 3年
度, 以实现广义切比雪夫滤波器耦合矩阵的消元. 最 后, 依据提取的耦合矩阵设计微波滤波器, 并给出相 应的测试结果.
T i a n gY u b o ,M aY u a n
( S c h o o l o f E l e c t r o n i c a n dI n f o r m a t i o n , J i a n g s uU n i v e r s i t yo f S c i e n c ea n dT e c h n o l o g y , Z h e n j i a n gJ i a n g s u 2 1 2 0 0 3 , C h i n a )
; i 为按适度值由低到高的每个父代杂 例因子, 取为 1 草所对应的序号)的正态分布进行扩散. 按以上扩 散方式, 适应度好的父代杂草在其较近的位置产生 子代, 而适应度较差的父代杂草在其较远的位置产 生子代, 这样可以进一步提高算法的寻优能力. 4 )竞争排斥. 通过数代的繁殖后, 后代的数目 将超过环境资源的可承受能力. 通过最大种群数目 P 在算法迭代过程中, 种群 m a x 对种群数量进行控制. 中的所有杂草和其后代按适度值从大到小依次排
野草算法优化提取广义切比雪夫滤波器耦合矩阵
田雨波,马 远
( 江苏科技大学 电子信息学院, 江苏 镇江 2 1 2 0 0 3 ) 摘 要:介绍了野草算法优化提取广义切比雪夫滤波器耦合矩阵的方法. 首先用野草算法控制耦合矩阵相似变换时的旋 转角度, 通过对目标函数的优化归零, 实现对耦合矩阵任意位置元素的消零. 接着给出广义切比雪夫滤波器的几种拓扑结 构, 设置传输零点在任意有限远处, 并通过对其耦合矩阵的提取, 计算耦合矩阵对应于归一化原型的幅频响应, 验证了这种 方法的可行性. 最后以优化提取的耦合矩阵为基础设计了微带滤波器, 测试结果与设计指标相吻合. 关键词:野草算法;耦合矩阵;滤波器 中图分类号: T N 8 文献标志码:A 文章编号: 1 6 7 3- 4 8 0 7 ( 2 0 1 3 ) 0 6- 0 5 6 9- 0 6
E x t r a c t i o no f c o u p l i n gma t r i xo f g e n e r a l i z e dC h e b y s h e vf i l t e r b a s e do ni n v a s i v ew e e d s o p t i mi z a t i o n
取适度值最优的前 P 其余个体将被去 序. m a x个个体, 除. 根据上述野草算法的主要步骤使用 M A T L A B 编程实现, 具体的算法流程如图 1 .
1 野草算法基本理论及改进
I WO是受野草启发而提出的、 基于种群的数值 优化计算方法, 其主要执行过程有如下 4 个步骤: 1 ) 初始化种群. 对一些参数进行初始化设置: 问 题的维数 D , 随机初始化 G个杂草植株, 最大迭代次 数i t e r , 最大种群数目 P , 最大和最小可生成种 m a x m a x 和S , 非线性指数 n , 区间步长初始值 σ 子数 S m a x m i n i n i t 和最终值 σ 及初始搜索空间 X . f i n a l 2 ) 生长繁殖. 繁殖过程按照自然界中的繁殖法 则, 给予不可行的个体较小的生存和繁殖的机会, 根 据适度值及种群所有个体的最小和最大适度值来确 定种群中的成员能够播散的种子数, 计算公式如下: N u m ( i )= ( s e e d e e d ) f i t n e s s ( i )-W o r s t F i t n e s s ) ( m a x -s m i n + ( B e s t F i t n e s s -W o r s t F i t n e s s ) s e e d m i n ( 1 ) s e e d e e d 式 中: m a x 和 s m i n 为 最 大 和 最 小 种 子 数, B e s t F i t n e s s 和W o r s t F i t n e s s 为族群种子的最好和最差 适度值,f i t n e s s ( i )为族群中第 i 个种子的适度值. 3 )空间扩散. 产生的种子个体按照正态分布在 选取的 其父代杂草个体附近的 D维空间进行扩散, 正态分布均值为 0 , 标准差为 σ . 每轮迭代对应不 c u r 同的标准差, 随着迭代的进行, 标准差从 σ i n i t开始至 σ σ f i n a l结束. c u r计算公式如下: σ c u r= ( i t e r t e r ) m a x -i )+σ ( σ σ i n i t- f i n a l f i n a l n ( i t e r ) m a x ( 2 ) 为了提高算法的性能, 对上述空间扩散的过程 进行适当的改进. 对父代杂草按其适度值进行排序 ( 假定父代杂草的个数为 n f ) , 每个父代杂草按均值 为0 , 标准差为 σ =σ β c u r+ i -1 ( 其中 β 为比 σu r n f c
能优化方法得到任意拓扑结构的耦合矩阵. 文献 [ 3 ] 给出单纯形法和梯度法相结合的方案来提取 4 ] 采用粒子群算法优化提取广义 耦合矩阵, 文献[ C h e b y s h e v 线性相位滤波器耦合矩阵, 以上两种方 法均为直接优化滤波器的传输函数之法, 因而对初 始值的选取比较敏感. 野草算法( i n v a s i v ew e e d s o p t i m i z a t i o n , I WO ) 是 一种模拟杂草繁殖过程的新型的随机: I n v a s i v e w e e d s o p t i m i z a t i o n( I WO )a l g o r i t h mi s i n t r o d u c e dt o e x t r a c t t h e c o u p l i n g m a t r i x o f t h e g e n e r a l i z e dC h e b y s h e vf i l t e r .F i r s t l y ,t h eI WOi s u s e dt oc o n t r o l t h er o t a t i o na n g l ei nt h ep r o c e s s o f c o u p l i n gm a t r i x s i m i l a r t r a n s f o r m a t i o n .B y o p t i m i z i n g t h e o b j e c t i v e f u n c t i o n ,z e r o s u p p r e s s i o no f c o u p l i n g m a t r i x e l e m e n t s a t a r b i t r a r yp o s i t i o n s i s r e a l i z e d .T h e n ,s e v e r a l t o p o l o g i e so f f i l t e r a r eg i v e na n dt r a n s m i s s i o nz e r oi ss e t t oa n yf i n i t e p o i n t . T h ee x t r a c t i o no f c o u p l i n g m a t r i x d e m o n s t r a t e s t h e f e a s i b i l i t y o f t h e p r o p o s e dm e t h o d . F i n a l l y , b a s e do nt h e e x t r a c t e dc o u p l i n gm a t r i x , m i c r o w a v ef i l t e r i sf a b r i c a t e d ,a n dt h et e s t i n gr e s u l t sa r ei na g r e e m e n t w i t ht h ed e s i g n e do n e s . K e yw o r d s :i n v a s i v ew e e d s o p t i m i z a t i o n ;c o u p l i n gm a t r i x ;f i l t e r 通信系统对微波滤 随着现代通信技术的发展, 波器的性能指标及体积大小提出了越来越严格的 广义切比雪夫滤波器通过引入有限频率的传 要求. 输零点, 在不增加滤波器阶数的同时, 能有效提高 滤波器的性能, 因此被广泛应用于各种微波通信系 统中. 文献[ 1- 2 ] 中分别提出了 N阶及 N+ 2阶耦 合矩阵的综合理论, 但是提取出来的耦合矩阵对应 的拓扑结构很难直接用于工程实际中, 需要根据实 际的拓扑结构对其进行相似变换. 由此可见, 如何 提取出可行的耦合矩阵是广义切比雪夫滤波器设 计的重点及难点. 传统的矩阵相似变换方法只能适 用于特定的拓扑结构, 很难推广使用, 故而采用智
T 置矩阵 R 来实现的: T M1 =R ·M0·R
( 3 )
式中: M0为初始矩阵, M1为相似变换后的矩阵. 6 阶 旋转矩阵 R ( 支点为[ 3 , 5 ] , 旋转角度为 θ )如式( 3 ) 定义, 其中 c o s ( ) , s i n ( ) . θ θ r =c r =s
1 0 0 R = 0 0 0
对称零点结构归一化幅频响应fig8normalizedamplitudefrequencyresponsesymmetriczerosstructure从以上4例可以看出对任意的拓扑结构使用野草算法来优化提取耦合矩阵经过迭代都能得到良好的收敛效果实现了较快的收敛速度和较好的收敛精度并且优化结果不依赖于算法初值的选取预先设定的传输零点的位置与所提取出的耦合矩阵的传输函数的零点位置基本吻合得到耦合矩阵对应的广义切比雪夫滤波器在物理上是可实现的
n
图1 野草算法流程 F i g . 1 F l o w c h a r t o f I WO
2 广义切比雪夫滤波器耦合矩阵的相 似变换
文献[ 1 ] 详细介绍了综合法提取广义切比雪夫 滤波器耦合矩阵的方法. 综合方法得到的广义切比 雪夫滤波器耦合矩阵一般来说无法直接使用, 因为 耦合矩阵的元素都是非零元素. 由于耦合矩阵的特 征值决定了滤波器的传输和反射特性, 只要保证耦 合矩阵的特征值不变, 滤波器的频率响应特性就不 会改变. 根据矩阵理论, 对矩阵进行相似变换, 即在 不改变矩阵特征值的前提下, 对矩阵的元素作消零 处理. 一个 N阶耦合矩阵 M0 的相似变换, 是通过在 M0 的左边和右边分别乘以 N阶旋转矩阵 R及其转
第2 7卷第 6期 V o l 2 7 N o 6 江苏科技大学学报( 自然科学版) D e c . 2 0 1 3 2 0 1 3年 1 2月 J o u r n a l o f J i a n g s uU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y( N a t u r a l S c i e n c eE d i t i o n ) d o i : 1 0 . 3 9 6 9 / j . i s s n . 1 6 7 3- 4 8 0 7 . 2 0 1 3 . 0 6 . 0 1 3
5 ] 优化算法 [ . I WO算法自 2 0 0 6年提出以来, 得到
了广泛的应用, 并显现出强大的优势和潜力, 在电
6 - 1 0 ] 磁学领域内亦获得了应用 [ . 文中采用野草算
优化旋转矩阵的角 法, 在电路综合理论的基础上,
5 7 0
江苏科技大学学报( 自然科学版)
2 0 1 3年
度, 以实现广义切比雪夫滤波器耦合矩阵的消元. 最 后, 依据提取的耦合矩阵设计微波滤波器, 并给出相 应的测试结果.
T i a n gY u b o ,M aY u a n
( S c h o o l o f E l e c t r o n i c a n dI n f o r m a t i o n , J i a n g s uU n i v e r s i t yo f S c i e n c ea n dT e c h n o l o g y , Z h e n j i a n gJ i a n g s u 2 1 2 0 0 3 , C h i n a )
; i 为按适度值由低到高的每个父代杂 例因子, 取为 1 草所对应的序号)的正态分布进行扩散. 按以上扩 散方式, 适应度好的父代杂草在其较近的位置产生 子代, 而适应度较差的父代杂草在其较远的位置产 生子代, 这样可以进一步提高算法的寻优能力. 4 )竞争排斥. 通过数代的繁殖后, 后代的数目 将超过环境资源的可承受能力. 通过最大种群数目 P 在算法迭代过程中, 种群 m a x 对种群数量进行控制. 中的所有杂草和其后代按适度值从大到小依次排
野草算法优化提取广义切比雪夫滤波器耦合矩阵
田雨波,马 远
( 江苏科技大学 电子信息学院, 江苏 镇江 2 1 2 0 0 3 ) 摘 要:介绍了野草算法优化提取广义切比雪夫滤波器耦合矩阵的方法. 首先用野草算法控制耦合矩阵相似变换时的旋 转角度, 通过对目标函数的优化归零, 实现对耦合矩阵任意位置元素的消零. 接着给出广义切比雪夫滤波器的几种拓扑结 构, 设置传输零点在任意有限远处, 并通过对其耦合矩阵的提取, 计算耦合矩阵对应于归一化原型的幅频响应, 验证了这种 方法的可行性. 最后以优化提取的耦合矩阵为基础设计了微带滤波器, 测试结果与设计指标相吻合. 关键词:野草算法;耦合矩阵;滤波器 中图分类号: T N 8 文献标志码:A 文章编号: 1 6 7 3- 4 8 0 7 ( 2 0 1 3 ) 0 6- 0 5 6 9- 0 6
E x t r a c t i o no f c o u p l i n gma t r i xo f g e n e r a l i z e dC h e b y s h e vf i l t e r b a s e do ni n v a s i v ew e e d s o p t i mi z a t i o n
取适度值最优的前 P 其余个体将被去 序. m a x个个体, 除. 根据上述野草算法的主要步骤使用 M A T L A B 编程实现, 具体的算法流程如图 1 .
1 野草算法基本理论及改进
I WO是受野草启发而提出的、 基于种群的数值 优化计算方法, 其主要执行过程有如下 4 个步骤: 1 ) 初始化种群. 对一些参数进行初始化设置: 问 题的维数 D , 随机初始化 G个杂草植株, 最大迭代次 数i t e r , 最大种群数目 P , 最大和最小可生成种 m a x m a x 和S , 非线性指数 n , 区间步长初始值 σ 子数 S m a x m i n i n i t 和最终值 σ 及初始搜索空间 X . f i n a l 2 ) 生长繁殖. 繁殖过程按照自然界中的繁殖法 则, 给予不可行的个体较小的生存和繁殖的机会, 根 据适度值及种群所有个体的最小和最大适度值来确 定种群中的成员能够播散的种子数, 计算公式如下: N u m ( i )= ( s e e d e e d ) f i t n e s s ( i )-W o r s t F i t n e s s ) ( m a x -s m i n + ( B e s t F i t n e s s -W o r s t F i t n e s s ) s e e d m i n ( 1 ) s e e d e e d 式 中: m a x 和 s m i n 为 最 大 和 最 小 种 子 数, B e s t F i t n e s s 和W o r s t F i t n e s s 为族群种子的最好和最差 适度值,f i t n e s s ( i )为族群中第 i 个种子的适度值. 3 )空间扩散. 产生的种子个体按照正态分布在 选取的 其父代杂草个体附近的 D维空间进行扩散, 正态分布均值为 0 , 标准差为 σ . 每轮迭代对应不 c u r 同的标准差, 随着迭代的进行, 标准差从 σ i n i t开始至 σ σ f i n a l结束. c u r计算公式如下: σ c u r= ( i t e r t e r ) m a x -i )+σ ( σ σ i n i t- f i n a l f i n a l n ( i t e r ) m a x ( 2 ) 为了提高算法的性能, 对上述空间扩散的过程 进行适当的改进. 对父代杂草按其适度值进行排序 ( 假定父代杂草的个数为 n f ) , 每个父代杂草按均值 为0 , 标准差为 σ =σ β c u r+ i -1 ( 其中 β 为比 σu r n f c