用离散小波包求解介质柱的散射.pdfVIP

用离散小波包求解介质柱的散射 魏兴昌 粱昌洪 (西安电子科技大学微波电信工程系 西安710071) 【摘要】 矩量法求解电大尺寸介质柱敬射时生成大型线性方程组，本文用离散小 波包方法对其快速求解。比较了该方法与离散小波变换方法的时间复杂度，分析和计 算结果表明，在电大尺寸介质柱散射分析中，用离散小波包可以最大限度稀疏化系数 矩阵，有效减少其求解时闯。 【关键词】 离散小波色 矩量法 电磁散射 矩量法(MOM)将电场积分方程(EFlE)变为线性方程组，是求解电磁散射的一个有 效方法，而对电大尺寸敬射体。MOM产生了大尺寸、稠密的系数矩阵，其求解十分费 时。长期以来，人们采用了各种高频近似方法解决这个问题，但由于大多数散射体既 表现出宏观的电大尺寸，又具有细节上的复杂结构．高频近似法显得理论模型粗糙， 计算精度低。近几年来，离散小波变换(DWT)的发展提供了一种把稠密矩阵快速稀疏 化的方法”】，把它与解方程组的共轭梯度法(CG)【21相结合，可以快速求解大型方程 的系数矩阵稀疏化后用CG法求解，比较了该DWP方法与D、ⅣT方法的时间复杂度． 计算了一个二维介质柱散身j，分析和计算结果表明，该方法可以有效地减少计算时间。 一、离教小波变换和离散小波包 根据多分辨分析理论．离散小波变换可以逐步地把信号不同频段的高频分量从信 号中提取出来，而离散小波包可以避一步对提取出的高频分量进行分解，因而能用更 少的项数展开信号。下面简介对Ⅳ维向量x进行DWT和DWP分解的过程Ⅱ】： 设小波滤波器系数为h。，g}=(一lrh。．其中．1≤k≤P时，hk≠0，其余k， 魄=0。定义，矿2x，，矩阵，V。，G∥其元索分别为： 14。U，，J。n(，一2帅l。d。 m G。(i，J)=gfH抽m ”7 利用MaUat算法，对Z进行～次分橱后的向量为 阡|器1Ix] ㈦ 其中，向量j，d分别代表向量x的低、高频分量。51,}于DWT，在下一步变换中只对s 替换s．如此选代下去，最终得到x的DWT分 一融 喇铲 陋啦 ≯T● 对于DWP，在下一步变换中有三种选择：l、 一42 = 替 ，．．．．，、。．L 耐甜 P．．．．．，，．L‰瓯 此n ，●●●●●，●j p 替龇，、同时对t删i姚即NNNISS!ds=f铅m 换s、d，如此迭代下去。可见，DWP的迭代过程对应于。个、：叉树，如煳1所不，}冬{ 准则为变换后向量的代价函数小于变换前向量的代价函数，也称为最好莲选取准则。 代价函数的定义有好几种Ⅲ，一般定义代价函数越小，对应的向量越容易稀疏化。 综上所述，Ⅳ维向量Ⅳ的DwT、DWP分解向量可甑用分解矩阵71表示为TX，由 小波滤波器的性质可以得到r矩阵的性质： l、r是正交、稀疏矩阵； 2、我们选取矩阵最大元素的口倍作为门限值，当矩阵元素值小于门限值时，认为 它是零，对于稠密矩阵c，在门限值的作用下，Cr‘，7C将是稀疏矩阵(对于DWP， 我们先选取～个检验向量Y，生成对于y是最好基的r矩阵．再与C相乘)。在此定义 矩阵的压缩率=矩阵中非零元素个数，矩阵中总共元素个数。 二、MOM系数矩阵的稀琉化 用MOM求解介质柱的散射，可以得到如下的线性方程组 CE=g (3) 其中．C是N×N的系数矩阵，￡是待求的介质柱内总电场向量，g是入射电场向量。 的，我们得到坐标交换下与(3)等价的方程组 C’E、=g、