登录
首页 » Others » 时域有限差分法 (高本庆)

时域有限差分法 (高本庆)

于 2021-05-06 发布
0 427
下载积分: 1 下载次数: 11

代码说明:

时域有限差分学习书籍,有研究FDTD方法的可以学习学习国防科技图书出版基金第二届评审委员会组成人员名誉主任委员怀国模主任委员黄宁副主任娄员殷鹤龄高景德陈芳允曾铎员尤子平朱森无朵英贤按姓氏笔划为序)刘仁何庆芝何因伟何新贵宋家树张汝果范学虹胡万忱柯有安侯迁候正明莫梧生崔尔东秘书长刘琯德时域有限差分法(Fnte- Differerce Time-Domain Method简称FDTD Method)是求解电磁问题的种数值技术,它是在1966年由K.S.Yee第一次提出的。fDTD法直接将有限差分式代替麦克斯韦时域场旋度方程中的微分式,得到关子场分量的有限差分式,用具有相同电参量的空间网格去模拟被研究体,选取合适的场初始值和计算室间的边界条件,可以得到包括时间变量的麦克斯韦方程的四维数值解。通过傅里叶变换可求得三维空间的频域解20多年来FDTD法历经了一个蓬勃发展的过程。最初是用它来求解金属体上的散射问题,用的是笛卡尔坐标系,使空间单元网格呈直角六面体。鉴于当时的计算机容量水平,特别是FDTD技术本身尚有若重要问题未很好解决,使得早期的数值精度不够高,应用范围也不很广,这种局面大约延续到70年代末期。随着FDTD技术的发展,首先需要解决的是有限计算空间的无反射截断问题,早期采用的一种方法是加大边界与散射体间距离,以在边界上构成外向行波,这种方法精度不高、计算空问亦大。直到将波方程的二阶近似用以处理边界上的场值,得到了较好近似的吸收边界条件,才将这个间题的解决向前推进了一大步在直角坐标系中用FDTD技术进行模拟时,光滑曲线形媒质表面将呈锯齿形状这可能产生沿面的表面波,加大了数值色散误差,解决这个问题的有效方法是“共形”技术的提出,这包括:或是使用曲线坐标系使媒质表面与坐标曲线共形,或是在直角坐标系中改变媒质介面上的网格形状,使二者共形,利用共形网格明显提高了计算精度。在类电磁问题中,当媒质结构尺寸比网格尺寸小时(如细线、窄槽或薄介质层等),将使FDTD模拟变得很困难。近来相继出现以麦克斯韦方程的国路积分形式建立相应FDTD算法式,FDTD与其他方法(如积分方程法或矩量法)的混合技术,以及媒质参数竹网格平均技术等,均提供了解决这类特殊问题的途径。FDTD法的特点是很易得到被研究体的近场,但不易一次直接得到远场值。80年代初期提出的利用等效原理将频城近场变换为远场是解决这个问题的好方法,近几年又将此技术发展到时域这二种路径给求解散射间题和天线问题提供了强有力的工具。值得提出的是,早期fTD方法中没有计及媒质的色散特性,即假定被研究媒质的电参数是与频率无关的。实际上自然界中有很多媒质的电参数具有很强的色散特性,近几年已开始注意研究色散媒质中的FDTD算法,为解决这电磁领域内难题铺平了道路。上述几方面问题的进展有力推动了FDTD技术的发展和应用,使它在解决复杂形体结构和多种媒质并存的一类呵题中占有重要的一席之地。今天,它不仅在电磁散射、电磁兼容预测、生物电磁学中得到卓有成效的应用,而且在天线做波技术、光电子学等的应用中意益受到重视。本书内容包括三部分:第一部分论述了FDTD法基本知识和各种FDTD算法,包括各种坐标系和特殊结构媒质的算法,以及FDTD与其他方法相结合的混合法等;第二部分介绍了FDTD法在电磁学各领域中的应用情况,内容涉及电磁散射微波传输线、天线、电磁榘容预測及生物电磁学范畴;第三部分讨论近期发展的色散媒质的FDTD算祛原理及其应用情况。本书在给出各种应用实例时,多将FDID法数据与理论值或实验数据或其他数值解数据(如矩量法等)进行比较。除早期发表的少数例子误差稍大外,大多数均与相比较的数据吻合很好。从目前水平看,在分析…般散射问题中,数值误差约在L%~3%附近(RCS的误差略大些),而在求谐振器本征值问题中,FDTD数值解与理论值误差低于百分之,某些情况能小于千分之一,这个精度是很喜人的。可以无夸张地说,FDTD法与其他数值解从精度上讲是可以媲美的,有的则有胜过。加之FDTD法得到的是时域解,通过傅氏变换可得到频域解。即它具有次时域计算代替频域上逐点计算的潜力。这些均L表示FDTD法具有较明显的优势。本书是作者尝试将已发表在不同场合不同时间的有关资料经整理推演、加工编写成专著。内容取材上亦包括作者及同事们近几年在这方面开展的工作结果。编写本书的目的是希望能起到抛砖引玉的效果。希给初学者提供…个入门途径,给从事这方面工作的同行们提供较系统的参考资料,以便更好地促进FDTD技术的发展和更进一步拓宽它的应用。由于作者才疏学浅,加上时间关系,错误和不足在所难免,上述目的恣难如愿敬请各位专家学者和读者多提供宝贵意见。在本书确立大纲的过程中,得到北京理工大学张德齐教授和楼仁海教授、西安交通大学汪文秉教授的帮助,汪文秉教授对内容安排提出了贵意见。要特别提到的是,在本书编写过程中得到中国科学技术大学旅美陈金元博士的大力支持和帮助陈博士及时寄来最新研究资料(见§7.6以充实本书内容。书稿完稿后,张德齐教授通阅了全稿,并提出宝贵意见国防工业出版社在本书出版过程中给予很多帮助,作者在此对他们表示衷心的感谢作者t993年8月于北京内容简介本书重点介绍」时域有限差分(FDTD)法的基本知识及其在电磁学各个领域内的应用。全书包括二部分(共八章),第一部分论述了FDTD法的基本知识(第…章)和各种FDTD算法,包括各种坐标和特殊结构媒质中FDTD算法及混合算法等(第二章);第二部分介绍FDTD法的爷种应用,内容涉及电磁散射(第三章〉微波传输线和谐振腔(第四章)、天线(第五章)、电磁兼容预测第六章)生物电磁学(第七章)等领域;第三部分介绍近期发展的色散媒质屮FDTD算法原及其应用情祝等(第八章〕。木书可供从事时域计算电磁学理论和应用妍究的人员参考,亦可供有关专业教师、研究生及商年级大学生作选修和参考用书目·录第一章时域有根差分法基本知识1.1支克斯韦方程1)].2FTD基本方程■h吾早4中■暑h音h鲁↓ψ中4山■■骨日◆p4白(生)1.3解的稳定性“……““〔l1)1.4边界条件…(13)1.5激励派的类型和设置……………………………(22)816误差分析………*………(29)§1.7近场一远场变换山中■冒↓d4白■■■e40)氵1.8FDTD数值解步骤…………(45〕参考文献■··■跏■·幽自嚞ψ■■·自·■甲血幽b中■●鲁歌●■甲晉『甲●■目■■46第二章FDTD算法变异48)82.1一般曲线坐标系的FDTD算法4甲●鲁ψψ■“ψ晕P山血曾■■■省口_P■口(48)§2.2正交曲线坐标系的PDTD算法(58}2.3非均匀网格尺寸的FDTD算法…(69)2.4细薄结构媒质的FTD算法……………*………(80)2.5FDTD瞬态积分方程的混合算法:(95§2.6FDTD矩量法时混合算法…………………*………(101)参考文献…………………………………………(105)第三章在电磁散射问题中的应用聊司自■电■■悬■■··ψb■■t最聊10883.1二维散射体■■画■·■唱口■■』烟·■口甲■唱■■即■■■口·吾■昌■■■■即■■■晶』h画p■■■昌〔108§3.2三维敢射体……………(114)3.3RCS计算……〔121)§3.4散射体的时域综合……………………**……(134)参考文献(143)第阿章在徽波传输线和谐振腔中的应用…………………(141徵带和共面传输线咖自b·山山山啬■■山■■ψ·;跏·■■“■士鲁■b■口■■■动●咖【41§4.2徼带不均匀性和徼带元件●日日鲁■日聊聊昏目申聊■与;自语↓ψ■ψ口一·↓c自ψb53)氵4.3波导传贛线元件……………………………………(63)8.4谐振腔的本征值……4…(j72)4.5谐振腔的e值(I名1参考文献…………………………………………(185第五章在天线问题中的应用·……(!8?)§5.1圆柱形单极天线………4■·↓■自·h■■■■■■■■甲4■■■▲187)5.2波导口和喇明线………………………………………(196)§5.3徽带贻片天线200)§5.4天线互朝的计算…ra……(210)参考文献1■幽昏p■■◆(215第六章在电磁兼容预测中的应用…(218§6·]瞬变电驚环境下飞机表面效应………“(216)§6.2EMP对连接有尾气焰导弹体的效应…(22)s63EM场透入导弹导引头的预浏■■■非十(230)6.4飞行体上微带贴片天线的EMP效应……a(233)6.5TEM传导胞腔中不均匀性影响的预测…-4…*4(239参考文献昏■■·聊■●■■■·ψ■山幽■如卓■■ψ卩■■中斷鲁晋■■口自■■导画晋■b山1〔241第七章在生物电磁学中的应用……246§7.1生物组织的电磁特性及人体电磁模型……………"……(246)72平面波照射下人体内的电磁效应号昌◆+卓b4即u(z56)7.3工业加热器对人体的作用…………………………(263)7.4动力电的人体效应…sa…(2687.5高压EMP的生物效.…·…(273)N7.6虾篾电话的人体效应…■■■■■↓鲁号↓個山山……(280)参考文………………………………………(2865第八章色教媒质的FDTD算法日■■d十■T■一一會■■『■■■血288)N8.1离散时域卷积建立FDTD方程……〔288N8.2由德拜方程和频域场到吋域的直接变换建立FDTD方程"sr…:…(295)氵8.3由z变换建立FDTD方程日ψ号■自卓自甲ψ吾■■■日■■■道语晶画『4b即44■(298)4色散媒质的吸收边界条件…………………"………(303)8.5(FD)TD法应用举例……(3058.6表面阻抗概念在fDTD法中的应用……(314)s8.了表面阻抗FDT法应用举例………■·■■b■324)参考文献(召30

下载说明:请别用迅雷下载,失败请重下,重下不扣分!

发表评论

0 个回复

  • FECG 胎儿心电
    胎儿心电信号的预处理以及胎儿心率的提取算法,基于matlab代码实现
    2020-12-03下载
    积分:1
  • sscom官方最新版本
    sscom5.13.1是sscom串口工具的最新版本,解决了usb插拔死机的问题
    2020-12-09下载
    积分:1
  • 图像分割的源代码资源
    图像分割图像分割图像分割图像分割图像分割图像分割图像分割图像分割
    2020-11-29下载
    积分:1
  • 独立子空间
    用matlab 完成独立子空间分析的实现
    2020-12-07下载
    积分:1
  • 基于中颖sh79f166A demo 上手代码
    基于中颖sh79f166A的demo 涉及到定时器 中断 PWM 串口 LCD 按键的一初始化和应用的demo,代码风格规范,合适刚开始接触改单片机的新手使用。
    2020-12-04下载
    积分:1
  • 电磁场与电磁波-(王家礼,朱满座,路宏敏)课后习答案
    电磁场与电磁波-(王家礼,朱满座,路宏敏)课后习题答案
    2020-12-04下载
    积分:1
  • 常用数据校验源代码(CRC8 CRC16 和校验 异或校验)
    文档中包含了CRC8位数据校验 CRC16位数据校验 和校验 异或校验等常用校验方式,使用C++代码编写,很实用。
    2020-12-07下载
    积分:1
  • 视频中人体动作识别
    频中的人体动作识别研究,视频动作识别,视频识别动作视频下载,人体油画视频,人体解剖学视频,德国人体解剖学视频,新鲜人体解剖学视频,真实人体解剖学
    2020-11-06下载
    积分:1
  • 三维装箱的模型与改进遗传算法
    关于三维装箱算法问题, 一些算法理论, 感觉对这方面的应用有一定帮助144效学的实践与认识40着∑(B,*v)≤VB,B,PD,PWy=0或者1v∈{1,2,…,D},y∈{1,2,…,W},z∈{1,2,…,H},j∈{1,2,…,n}(12目标函数是箱子未装填物品的空间最少(亦即空间浪费最少)条件(2)确保子的1个装填空间单元被装填不超过1次即保证物品间不会互相嵌入;(3)式说明上层物品会有支撑,不会悬空(4),(5),(6)式说明物品装箱位置约束;(7),(8},(9)说明物品的摆放问;(11)是箱子的容积约束2這传算法21遗俊法遗传算法(GAs)是建立在达尔文进化论基础上的搜索算法,它从代表问题潜在解的个种群( population)开始,而一个种群则由经过基因(gene)编码 coding)的一定数目的个体individual)组成遗传算法采用了自然进化模型,如选择,交叉变异等计算开始时,一定数目S个个体(父个体1、父个体2……)即种群随机地初始化,并计算每个个体的适应度函数第一代也即初始代产生如果不满足优化准则,开始产生新一代的计算为了产生新一代按照适应度选择个体,父代通过基因重组(交叉)而产生子代所有的子代按一定的概率变异然后重新计算子代的适应度,将子代插入到种群中取代父代构成新的一代循环执行这一过程,直到满足优化准则22算法设计2.21编码方法采用矩阵编码方法,用多维数组(二维矩阵表示染色体结构,数组元素表示染色体基因,编码清晰,易于理解,遗传算子操作方便染色体S=(L,P,Px,Py,T)来表示问题的一个解其中:向量L=(Li,L1,…,Ln)为待装箱物品的一个排列;向量P=(Bn,B1,…,B3n)为对应于排列L的B,一个排列向量Px=(PB,PB,…,PB)为对应于排列L的PB一个排列向量Py=(PB,PB},…,PB为对应于排列L的PBx一个排列矩阵T=(x2=欢面为对应于排列L的装箱物品坐标22适值函数问题的目标是最小化箱子的浪费空间,适应度函数可定义为空间利用率函数(S代表染色体C1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net2期陈德良,等:三维装箱问题的模型与改进遗传算法Fitness(s)(∑B3*v/若∑(B;*v)≤V否则23解的不可行性,罚函数与评估函数由于对染色体作遗传运算时可能获得不可行的子代,惩罚技术是用遗传算法解约束优化问题中最常用的技术,本质上它是通过惩罚不可行解将约束优化问题转化为无约爽问题就本文讨论的问题而言,惩罚项包括:1)物品在装箱时不交叠,即满足约束条件{2},有着∑By≤1g:(S)1,否则2)物品装箱时不能出现悬空即满足约束条件(3),有0若∑B-B+)>0g2(S)=(151,否则3)物品装箱不能超出箱子边界,即满足约束条件(4,(5)和(6),有0若吃+B(Pp++Pwy*吗)≤D1,否则0若+B*(PD*+PWy*m)≤W941,否则(17)95(S)=0,若x+B*h;≤H8)1,否则eat(s=∑9(S)b=1那么,式(14)至(18)任何一个取值为1,都是不可行评估函数eval(S)=Fitness(S)*(5-Genalty (S)24算法步骤)初始化进化代数计数器,随机产生一定数目(大于设定的初始种群规模)的染色体;2)利用式(14)检验初始种群染色体可行性,对不可行解旋转赌轮接受小部分不可行解,与可行解构成初始种群3)对初始种群染色体进行遗传运算;①按照式(14)至(20)计算评估函数:⑩按顺序交叉方法产生子代;④变异算子;4)旋转赌轮选择染色体;)重复3至4)直到完成给定的循环次数;C1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net数学的实践与认识40卷6)确定最好的染色体作为最优解3实验结果我们用C++编程实现了上述算法在配置为CPU24GH/512 Mb ram的微机上,用随机产生的数据进行实验取遗传算法运行参数为:{群体大小进化代数,交叉概率,变异概率}-{100,50,0.85,0.05}用随机产生的数据进行实验,求解20个种类100件物品的装箱问题,得到最好解耗时小于1秒;计算50个种类200件物品的装箱问题,得到最好解耗时小于2秒以下是3类共16件物品的装箱问题.实验数据图2,第!行为箱子尺『;第2至第4行为待装箱物品,每行第1个数据表小序号第2至镌4数据分别为物品尺寸,第5个数据表示物品件数在计算转桌中包含数据依次是:序号,是否装载,物品长,物品宽,物品高,纵向坐标横向坐标,垂向坐标纵向长度,横向长度,垂向长度(图4).从图4可知第12号物品未能装箱,物品装箱的顺序可以从“序号列中得出.绘制的物品装箱示意图见图31421,2,乙2,2,图2实验数据图3装箱示意图文件((格式(Q帮助新 s REPORT耗时:.1 most g sec次数:01615积:7580001016每a0库:92.875989名寸:=280;y=1210;2=300NO: P st Din 1 Din 2 in 3 C xC YPu y Pu 2202002002002001002812B20020012010010020012鲁2020012010020020020020012B100100212021201215024015824815015561111111111115152每000ao00015150240202020055020200201002001215502D0200120100100201205s012020028012015024075000015024815015024075015015024a152002009o002002001002012090020012日未装相物品121501502年图1计算结果o1994-2010ChinaAcademicJournalElectronicPublishingHousealLrightsreservedhttp://www.cnki.net2期陈德良,等:三维装箱问题的模型与改进遗传算法1474结束语装箱问题是一常见而难解的优化问题,利用遗传算法求解时,随机产生的初始解会出现大量的不可行解(装箱物品占用空间出现大量交叠),本文将箱子内部空间划分为一个个立方体单元:算法的第2)步对标准遗传算法做了修改通过剔除大量不可行解提高算法的收敛速度,实验结果表明此算法运算过程及绪果稳定,具有较强的实际应用价值能有效解决复杂的三维装箱何题,今后将继续研究将该方法运用到其它不同的有关装箱问题或组合优化问题中参考文献[1] John J, et al. An improved algctithra for the ncn-guillctine-constraincd cutting-stock problem(JIOperational Resee ch Society, 19 /0,+1: 141-149[2] Coffmau E. G, et al. ver age-case analysis of cutting and packing in two dimensions [J]. Euro. Jof Operatic al Reseaich, 1990, 44: 134-14413) Fabien C, et al. A Two-phase heuristic for the two-dinensional cutting-stock problem [J. Opera-tional Research Society, 1991, 42: 39-744 Martello Silvano, Pisinger, David, and Vigo, Daniele. The Three-Dimensional Bin Packing ProblemJ. Operations Research, 2000 Informs. Vol. 48: 256-267]何大勇,査建中,姜义东遗传算法求解复杂集装箱装载问题方法研究向]软件学报,201,12(9):13801385阿]张德富魏丽军陈青山陈火旺等.三维装箱问题的组合启发式算法软件学报,2007,18(9):20832089A Mixed Integer Programming Model ofThree-Dimensional Bin-Packing Problem and ImprovedGenetic AlgorithmsCHEN De-liang, 2, CHEN Zhi-yaSchool of Traffc &z transportation Engineering, Central South University, Changsha 41076, China)(2. Logistics School, Central South University of Forestry Technology, Changsha 410004, ChinaAbstracts The three-dimensional bin-packing problem is complicated but a high level ofinterest in developing effective way to solve this kinds of NP-hard problem. First a MixedInteger Programming model was worked out in this paper, which resorted to dividing box spaceinto unit cube. Then an improved genetic algorithm was mainly developed. Tests on hundredsof problems show that this algorithm makes the most of volume utilization in minimal timeKeywords: three-dimensional bin-packing problem; space division; mixed integer program-ming model; improved genetic algorithmso01994-2010ChinaAcademicJournalElectronicPublishinghOuse.Allrightsreservedhttp://www.cnki.net
    2020-12-05下载
    积分:1
  • IEICE论文LATEX模板
    IEICE的论文Latex模板,包含mac,unix,win下的英文和日文模板
    2021-05-06下载
    积分:1
  • 696516资源总数
  • 106914会员总数
  • 0今日下载