高阶FDTD的加权系数优化方法研究.pdf
《高阶FDTD的加权系数优化方法研究.pdf》由会员分享,可在线阅读,更多相关《高阶FDTD的加权系数优化方法研究.pdf(5页珍藏版)》请在咨信网上搜索。
1、第 卷 第期 年月金 陵 科 技 学 院 学 报J OUR NA LO FJ I N L I NGI N S T I TUT EO FT E CHNO L O G YV o l ,N o M a r,D O I:/j c n k i /n 高阶F D T D的加权系数优化方法研究朱敏(金陵科技学院电子信息工程学院,江苏南京 )摘要:为了提高电磁散射的计算精确度,对高阶时域有限差分法(F D T D)的加权系数进行优化.提出一种基于Y e e氏蛙跳网格的高阶F D T D优化方法,主要通过优化加权系数消除传播方向上沿轴及网格对角线的色散误差.通过雷达截面散射仿真实验验证了提出的高阶F D T D
2、优化方法相较于未优化的高阶F D T D方法具有更高的计算精度,可应用在复杂目标的电磁散射问题上.关键词:高阶F D T D;各向异性;加权系数;色散误差;雷达散射截面中图分类号:TN ;O 文献标识码:A文章编号:X()收稿日期:基金项目:金陵科技学院科研基金孵化项目(j i t f h x m )作者简介:朱敏(),女,安徽蚌埠人,讲师,博士,主要从事计算电磁学、电磁场与微波技术研究.R e s e a r c ho nW e i g h t i n gC o e f f i c i e n t sO p t i m i z a t i o nM e t h o d so fH i g h
3、 e r O r d e rF D T DZ HU M i n(J i n l i n gI n s t i t u t eo fT e c h n o l o g y,N a n j i n g ,C h i n a)A b s t r a c t:I no r d e r t o i m p r o v e t h e c a l c u l a t i o na c c u r a c yo f e l e c t r o m a g n e t i c s c a t t e r i n g,t h ew e i g h t i n gc o e f f i c i e n t so
4、ft h eh i g h e r o r d e rf i n i t e d i f f e r e n c et i m e d o m a i n(F D T D)m e t h o da r eo p t i m i z e d T h eo p t i m i z e dh i g h e r o r d e rF D T D m e t h o db a s e do nY e e sl e a p f r o gg r i di sp r o p o s e d,w h i c hm a i n l ye l i m i n a t e s t h ed i s p e r s
5、 i o ne r r o r a l o n g t h e a x i s a n dg r i dd i a g o n a l i n t h ep r o p a g a t i o nd i r e c t i o nb yo p t i m i z i n gt h ew e i g h t i n gc o e f f i c i e n t s S i m u l a t i o ne x p e r i m e n t so fr a d a rc r o s ss e c t i o nf u r t h e rv e r i f yt h a t t h eo p t
6、i m i z e dh i g h e r o r d e rF D T D m e t h o dp r o p o s e dh a sh i g h e rc a l c u l a t i o na c c u r a c yt h a nt h eu n o p t i m i z e dh i g h e r o r d e rF D T D m e t h o d I tc a nb ea p p l i e dt oe l e c t r o m a g n e t i cs c a t t e r i n gp r o b l e m so f c o m p l e xt
7、 a r g e t s K e yw o r d s:h i g h e r o r d e rF D T D;a n i s o t r o p y;w e i g h t e dp a r a m e t e r;d i s p e r s i o ne r r o r s;r a d a rc r o s ss e c t i o n传统时域有限差分法(f i n i t e d i f f e r e n c et i m e d o m a i n,F D T D)是最先被提出的一种有限差分算法,但受数值稳定性和色散条件的制约,它难以提供较高的计算精度.高阶F D T D在某些特定
8、条件下可降低色散误差和各向异性,但算法本身仍需要完善.为了解决上述问题,学者们相继提出了很多新的高精度差分算法.C o l e提出了一种高精度的非标准差分方法.T a f l o v e和H a g n e s s对高阶F D T D和其他多种算法进行了比较研究和优化.N e h r b a s s等给出了一种在不提高阶数的情况下仍能增加计算精度的优化方法.J u n t u n e n等提出的各向异性F D T D方法在不增加计算量的前提下提高了数值准确度.W a n g和T e i x e i r a提出了一些改进方法,如复杂滤波器对角优化法,以降低色散误差.S a x e n a等提出第
9、期朱敏:高阶F D T D的加权系数优化方法研究了高阶局部一维HO L O D F D T D方法,以提高四步L O D F D T D、五步L O D F D和六步L O D D T D方法的精度,与HO L O D F D T D方法相比在所需方向上提供了更高的精度且花费的计算时间更少.D i n g提出了一种基于各向同性色散有限差分格式的无条件稳定分裂步长S S F D T D方法.此方法对各向同性色散F D T D的加权因子和比例因子进行了改进,可以为单个频率生成较精确的相位速度.与传统的四级分步F D T D方法相比,即使在稳定因子较大的情况下,S S F D T D方法的数值色散误
10、差也显著减小,数值模拟验证了此方法的有效性.W a n g等提出了一种基于跳变交替方向隐式F D T D和传统F D T D计算传输线网络瞬态响应的混合时域有限差分法,混合方法具有较好的精确度.K u a f等给出了一种新的部分隐式三维F D T D方法,此方法的最大色散误差和各向异性误差分别为 和 .这些改进方法在计算精度上有一定的提升,但是缺少对传播方向上色散误差的改进.本文在上述研究的基础上,提出一种基于Y e e氏蛙跳网格的高阶F D T D优化方法,降低数值色散误差,消除在给定频率下各向异性带来的误差,对于宽频带信号也能通过该方法减小色散误差和平均累计的相位差.加权系数的优化本文提出
11、的高阶F D T D优化方法主要是通过优化加权系数来降低数值色散误差.为了简单说明该方法,下面考虑一个在线性各向同性均匀无耗介质中传播的二维T Ez波,高阶F D T D优化方法的M a x w e l l方程在x方向的迭代公式可以改写为如下离散形式:Enx(i/,j,k)Enx(i/,j,k)tHn/z(i/,j/,k)Hn/z(i/,j/,k)yHn/y(i/,j,k/)Hn/z(i/,j,k/)z()tHn/z(i/,j/,k)Hn/z(i/,j/,k)yHn/y(i/,j,k/)Hn/z(i/,j,k/)z()式中,为介电常数,E为电场强度,H为磁场强度,y和z分别为y和z轴上的空间步
12、长,离散均匀即yz,t为时间步长,i、j、k分别为计算区域元胞的标号.为需要优化的加权系数,代表采用的是传统F D T D方法;代表采用的是高阶F D T D方法.下面通过增长因子 来分析本文提出的高阶F D T D优化方法的色散误差.使用傅里叶分析方法 可以得到优化的高阶F D T D增长因子为(rxpxrypy)j(rxpxrypy)e x p j t a n(rxpxrypy)(rxpxrypy)()式中,j,rxcts i n(kxx/)/x,rycts i n(kyy/)/y,c为真空中的光速,kx为沿x轴向的数值相位常数,kxkc o s(),kyks i n(),数值相位kkxk
13、y,为与x轴相关的电磁波的传播角度,px()s i n(kxx/)/,py()s i n(kyy/)/.通过对式()增长因子方程 的计算和整理,可以得到关于加权系数的色散方程为s i n(t/)cts i n(kxx/)x()s i n(kxx/)xs i n(kyy/)y()s i n(kyy/)y()从式()可以看出,不仅可以通过改变加权系数来降低色散误差,还可以通过调节其他的相关参数来优化色散误差.沿轴方向数值色散误差的优化在均匀离散的情况下即xy时,为了消除沿x轴和y轴方向的色散误差,本文对式()进行处理,金陵科技学院学报第 卷这里以x轴方向为例来说明,由式()可得x轴方向的色散方程为
14、s i n(t/)cts i n(kxx/)x()s i n(kxx/)x()求解式()可以得到轴向优化加权系数为s i n(t/)cts i n(kxx/)xs i n(kxx/)xs i n(kxx/)x()式中,在给定的频率下kx/,为给定频率下的波长;网格步长x/N,N 为每波长网格数.这些参数代入式(),就是加权系数关于每波长网格数N 和时间步长的方程,确定其中的两个量就能得到另一个量的解.根据x轴色散方程式(),简化得到以下形式:s i n(s/N)/s s i n(/N)s i n(/N)s i n(/N)()式中,s为C o u r a n t稳定系数.图加权系数随每波长网格数
15、N变化的曲线图为不同s下加权系数随每网格波长数N 变化的曲线,其中Sc,s分别为 Sc、Sc、Sc和 Sc.从图可以看出,当N 固定时随着s的增加而减小,当s固定时随着N 的增加而减小.数值色散误差可以通过数值相速Vn/c来衡量.图给出了加权系数优化前N 时的相对相速在不同s条件下随着传播角度变化的曲线.从图可以看出,色散误差随着s的减小而减小.根据式(),加权系数可以通过相对相速Vn/c来优化,优化后的曲线如图所示.从图可以看出,加权系数优化后的相对相速与传播角度的关系与优化前的一致,优化后的色散误差在相同s下更低.图加权系数优化前N 时相对相速随传播角度变化的曲线图加权系数优化后N 时相对
16、相速随传播角度变化的曲线 对角线上数值色散误差的优化本节采用对角线优化方法(d i a g o n a l l yo p t i m i z e dm e t h o d,D OM)对数值色散的误差进行优化,在二维情况下,将传播角度 代入式(),可得沿 方向上的色散方程为s i n(t/)cts i n(dx/)x()s i n(dx/)x()第期朱敏:高阶F D T D的加权系数优化方法研究式中,kdk/,k/.求解式()得到关于加权系数d的表达式为ds i n(t/)cts i n(dx/)xs i n(dx/)xs i n(dx/)x()式()的加权系数d在不同的C o u r a n
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- FDTD 加权 系数 优化 方法 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。