科技写作结课作业(时域有限差分法的Matlab仿真开题报告)
时域有限差分法的Matlab仿真
时域有限差分法的Matlab仿真关键词: Matlab 矩形波导时域有限差分法摘要：介绍了时域有限差分法的基本原理，并利用Matlab仿真，对矩形波导谐振腔中的电磁场作了模拟和分析。
关键词：时域有限差分法；Matlab；矩形波导；谐振腔目前，电磁场的时域计算方法越来越引人注目。
时域有限差分（Finite Difference Time Domain，FDTD）法［1］作为一种主要的电磁场时域计算方法，最早是在1966年由K. S. Yee提出的。
这种方法通过将Maxwell旋度方程转化为有限差分式而直接在时域求解，通过建立时间离散的递进序列，在相互交织的网格空间中交替计算电场和磁场。
经过三十多年的发展，这种方法已经广泛应用到各种电磁问题的分析之中。
Matlab作为一种工程仿真工具得到了广泛应用［2］。
用于时域有限差分法，可以简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多的时间。
下面将采用FDTD法，利用Matlab仿真来分析矩形波导谐振腔的电磁场，说明了将二者结合起来的优越性。
1FDTD法基本原理时域有限差分法的主要思想是把Maxwell方程在空间、时间上离散化，用差分方程代替一阶偏微分方程，求解差分方程组，从而得出各网格单元的场值。
FDTD 空间网格单元上电场和磁场各分量的分布如图1所示。
电场和磁场被交叉放置，电场分量位于网格单元每条棱的中心，磁场分量位于网格单元每个面的中心，每个磁场（电场）分量都有4个电场（磁场）分量环绕。
这样不仅保证了介质分界面上切向场分量的连续性条件得到自然满足，而且还允许旋度方程在空间上进行中心差分运算，同时也满足了法拉第电磁感应定律和安培环路积分定律，也可以很恰当地模拟电磁波的实际传播过程。
1.1Maxwell方程的差分形式旋度方程为：将其标量化，并将问题空间沿3个轴向分成若干网格单元，用Δx，Δy和Δz 分别表示每个网格单元沿3个轴向的长度，用Δt表示时间步长。
时域有限差分法对平面TE波的MATLAB仿真
时域有限差分法对平面TE波的MATLAB仿真姓名：王云璐学号：2011302021（西北工业大学电子信息学院08041103，陕西西安，710072）摘要：本文分析FDTD算法的基本原理及两种典型边界条件的算法特点，通过MATLAB程序对TE波进行了仿真，模拟了高斯制下完全匹配层中磁场分量瞬态分布。
得到了相应的磁场幅值效果图。
最后得出用Matlab语言对FDTD算法编程的几点结论。
关键词：FDTD；MATLAB；PML；平面TE波1引言电磁场的有限差分一般是在时域上进行的，随着计算机技术的发展和应用，近年来，时域计算方法越来越受到重视。
时域有限差分法具有简单、结果直观、网格剖分简单等优点。
近些年FDTD发展的十分迅速，在许多方面都有重要应用，包括天线设计，微波电路设计，电磁兼容分析，电磁散射计算，光子学应用等等。
时域有限差分(FDTD)算法是K.S.Yee于1966年提出的直接对麦克斯韦方程作差分处理，用来解决电磁脉冲在电磁介质中传播和反射问题的算法。
基本思想是：FDTD计算域空间节点采用Yee元胞的方法，同时电场和磁场节点空间与时间上都采用交错抽样；把整个计算域划分成包括散射体的总场区以及只有反射波的散射场区，这两个区域是以连接边界相连接，最外边是采用特殊的吸收边界，同时在这两个边界之间有个输出边界，用于近、远场转换；在连接边界上采用连接边界条件加入入射波，从而使得入射波限制在总场区域；在吸收边界上采用吸收边界条件，尽量消除反射波在吸收边界上的非物理性反射波。
本文主要简述了FDTD算法的基本原理，解的稳定性以及边界条件特点，并用Matlab语言进行对平面TE波进行了编程计算。
2FDTD基本原理时域有限差分法的主要思想是把Maxwell方程在空间、时间上离散化，用差分方程代替一阶偏微分方程，求解差分方程组，从而得出各网格单元的场值。
FDTD空间网格单元上电场和磁场各分量的分布如图1所示。
图1.Yee 氏网格及其电磁场分量分布电场和磁场被交叉放置，电场分量位于网格单元每条棱的中心，磁场分量位于网格单元每个面的中心，每个磁场(电场)分量都有4个电场(磁场)分量环绕。
Matlab仿真实验报告_发布
2、实验目的
通过本次实验，应该掌握： (1) 用傅立叶变换进行信号分析时基本参数的选择。 (2) 经过离散时间傅立叶变换（DTFT）和有限长度离散傅立叶变换（DFT） 后信号频 谱上的区别，前者 DTFT 时间域是离散信号，频率域还是连续的，而 DFT 在两个域中都 是离散的。 (3) 离散傅立叶变换的基本原理、特性，以及经典的快速算法（基 2 时间抽选法） ，体 会快速算法的效率。 (4) 获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨 率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。 (5) 建立 DFT 从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念 的一个典型应用是数字音频压缩中的分析滤波器，例如 DVD AC3 和 MPEG Audio。
~4~
MATLAB 仿真实验报告
N=4000; % N>1/TΔf=1/357*T
n=0:1:N-1; f=fs*n/N;
%DTMF 信号编码 dtmf(1,:)=sin(2*pi*row(1)*T*n)+sin(2*pi*col(1)*T*n); dtmf(2,:)=sin(2*pi*row(1)*T*n)+sin(2*pi*col(2)*T*n); dtmf(3,:)=sin(2*pi*row(1)*T*n)+sin(2*pi*col(3)*T*n); dtmf(4,:)=sin(2*pi*row(2)*T*n)+sin(2*pi*col(1)*T*n); dtmf(5,:)=sin(2*pi*row(2)*T*n)+sin(2*pi*col(2)*T*n); dtmf(6,:)=sin(2*pi*row(2)*T*n)+sin(2*pi*col(3)*T*n); dtmf(7,:)=sin(2*pi*row(3)*T*n)+sin(2*pi*col(1)*T*n); dtmf(8,:)=sin(2*pi*row(3)*T*n)+sin(2*pi*col(2)*T*n); dtmf(9,:)=sin(2*pi*row(3)*T*n)+sin(2*pi*col(3)*T*n); dtmf(10,:)=sin(2*pi*row(4)*T*n)+sin(2*pi*col(2)*T*n); %1 %2 %3 %4 %5 %6 %7 %8 %9 %0
时域有限差分法的Matlab仿真
ｗ ｊ ＝ ｃ
ｎ ｛ 一 ｃ 。 ・ Ｅ 杠 一 《 ￣ ｘ ／ ２ － Ｅ ＋ ． ）
（ ７ ）
式 中系数 Ｃ Ａ， ｃ 曰 ， Ｃ Ｑ 的定义为 ：
Байду номын сангаас
Ｃ Ａ ＝ １ 一 ￣ ａ ｔ １ ／ ： ＋ 警） ； Ｃ Ｂ ＝ （ 等） ／ （ １ ＋ 警） ｃ ￣ － （ １ 一 ＆ ａ ｔ ） ， （ １ ＋ ） ； Ｃ Ｑ ＝ ， （ １ ＋ ）
比 ＋ Ｃ Ｂ 。 川 、 Ｈ Ｉ ＂ ￣ ￣ ／ ２ 川 一 Ｉ ）（ ５ ） Ｅ Ｉ ｉ ｄ ＋ ｌ ＝ Ａ Ｃ 。 ， ： ・ Ｅ 。 ＋ Ｃ Ｂ ｉ ， ｊ ＋ ｌ ／ ２ ＂ （ ｉ ＋ １ ／ ｍ ２ ￣ ／ 一 ｌ ｎ ＊ Ｉ ／ ２ ）（ ６ ）
【 Ｋｅ ｙ ｗｏ ｒ ｄ ｓ ｌ Ｆ Ｄ Ｔ Ｄ； Ｍａ ｔ ｌ ａ ｂ ； Ｅ ｌ ｅ ｃ ｔ ｒ ｉ ｃ ｉ ｆ ｅ ｌ ｄ ｄ ｉ ｓ ｔ ｉ ｒ ｂ ｕ ｔ ｉ ｏ ｎ
０ 引言
时 域有限差 分 （ Ｆ ｉ ｎ ｉ ｔ ｅ Ｄ ｉ ｆ ｅ ｒ ｅ ｎ ｃ ｅ Ｔ ｉ ｍｅ Ｄ ｏ ｍａ ｉ ｎ． Ｆ Ｄ Ｔ Ｄ）法是 Ｋ． Ｓ ． Ｙ ｅ ｅ在 １ ９ ６ ６年 给出的利用有 限差分式把麦克斯 韦（ Ｍａ ｘ ｗ ｅ ｌ １ ） 旋 度方 程 替换为一组差分方程＿ 】 １ 。 并提供所解问题中电磁特性物理意义 的算法 ． 可直接在时域 中求解
【 Ａ ｂ ｓ ｔ ｒ ａ ｃ ｔ 】 Ｔ ｈ ｅ ｂ ａ ｓ ｉ ｃ ｐ ｉ ｒ ｎ ｃ ｉ ｐ ｌ ｅ ｏ ｆ ｉ ｆ ｎ ｉ ｔ ｅ ｄ ｉ ｆ ｆ ｅ ｒ ｅ ｎ ｃ ｅ ｔ ｉ ｍｅ ｄ ｏ ｍａ ｉ ｎ ｉ ｓ ｉ ｎ ｔ ｒ ｏ ｄ ｕ ｃ ｅ ｄ ｉ ｎ ｔ ｈ ｉ ｓ ｐ ａ ｐ ｅ ｒ ． Ｗｉ ｔ ｈ ｔ ｗ ｏ — ｄ ｉ ｍ ｅ ｎ ｓ ｉ ｏ ｎ ａ ｌ ｉ ｆ ｎ ｉ ｔ ｅ ｄ ｉ ｆ ｆ ｅ ｒ ｅ ｎ ｃ ｅ ｔ ｉ ｍ ｅ ｄ ｏ ｍａ ｉ ｎ
MATLAB在离散信号时域变换与运算中的应用与实现(开题报告)
•
会泽百家
至公天下
系统的基本功能
1. 系统由离散信号的变换及离散信号的运算两个子模 块组成； 2. 系统为用户提供交互式的信号选择功能，为用户提 供可视化的离散信号的变换与运算的观察与分析界 面； 3. 通过典型信号将离散信号的变换与运算用生动直观 的二维曲线面图进行可视化表现； 4. 为“信号与系统”课程的教学提供计算机辅助教学 环境。
会泽百家
至公天下
设计意义
以计算机辅助教学的模式进行学习，这样学生能够 更加深刻地了解该课程，从而能够更加好的掌握课程所 涉及到的知识。所以将MATLAB引入到“信号与系统”课 程中的意义的非常巨大的。 • 利用MATLAB高效的数值计算和符号计算功能， 使学生能够从繁琐的数学运算分析中解脱出来 从而把有限的时间用到课程知识的理解与运用 上。 利用MATLAB完备的图形处理功能，实现计算结 果和编程的可视化，可以让信号能够很直观的
实验二 利用MATLAB进行时域分析报告
实验二 利用MATLAB 进行时域分析本实验容包含以下三个部分：基于MATLAB 的线性系统稳定性分析、基于MATLAB 的线性系统动态性能分析、和MATALB 进行控制系统时域分析的一些其它实例。
一、 基于MATLAB 的线性系统稳定性分析线性系统稳定的充要条件是系统的特征根均位于S 平面的左半部分。
系统的零极点模型可以直接被用来判断系统的稳定性。
另外，MATLAB 语言中提供了有关多项式的操作函数，也可以用于系统的分析和计算。
（1）直接求特征多项式的根设p 为特征多项式的系数向量，则MATLAB 函数roots()可以直接求出方程p=0在复数围的解v,该函数的调用格式为：v=roots(p)例3.1 已知系统的特征多项式为： 123235++++x x x x特征方程的解可由下面的MATLAB 命令得出。
>> p=[1,0,3,2,1,1]; v=roots(p) 结果显示：v =0.3202 + 1.7042i 0.3202 - 1.7042i -0.7209 0.0402 + 0.6780i 0.0402 - 0.6780i利用多项式求根函数roots(),可以很方便的求出系统的零点和极点，然后根据零极点分析系统稳定性和其它性能。
（2）由根创建多项式如果已知多项式的因式分解式或特征根，可由MATLAB 函数poly()直接得出特征多项式系数向量，其调用格式为：p=poly(v)如上例中：v=[0.3202+1.7042i;0.3202-1.7042i;-0.7209;0.0402+0.6780i; 0.0402-0.6780i];>> p=poly(v) 结果显示p =1.0000 0.0001 3.00002.0001 0.9998 0.9999由此可见，函数roots()与函数poly()是互为逆运算的。
（3）多项式求值在MATLAB 过函数polyval()可以求得多项式在给定点的值，该函数的调用格式为： polyval(p,v) 对于上例中的p 值，求取多项式在x 点的值，可输入如下命令：>> p=[1,0,3,2,1,1]; x=1polyval(p,x) 结果显示 x = 1 ans = 8(4)部分分式展开 考虑下列传递函数：nn n nn n a s a s a b s b s b den num s N s M +⋅⋅⋅+++⋅⋅⋅++==--110110)()( 式中0a 0≠，但是i a 和j b 中某些量可能为零。
电磁仿真时域有限差分法初步应用实验报告
验
报
刘伶
告
一、实验室名称：物理电子学院仿真实验室 二、实验项目名称：电磁仿真综合实践 I 级 三、实验学时：36 学时 四、实验原理：电磁仿真的时域有限差分法 五、实验目的：
加强学生的计算机综合应用能力、尤其是运用计算机分析和解决专业问 题的能力培养，使学生对独立进行科学研究有初步实践；初步掌握一种纯数 值电磁仿真方法——时域有限差分法；初步学会综合应用一种程序设计语言 进行科学与工程计算；增强科技报告写作能力，学会相关软件使用。
11
收边界后波形就传出。 (2)终端短路时
图 a: 向 Z=L 传输时的时域波形
图 b: 向 Z=0 传输时的时域波形
图 c: 两观测面的电压值，从上到下依次是左端面和右端面
分析： 图 a 和 b 知是电磁波整个空间内随时间变化波形图的两个截图。 由仿真结果可知， 电磁波随时间的流逝在空间上是先沿 Z=L 方向传播，再沿 Z=0 方向传播。由于终端是短路 边界，因此传输进来时值是 1，经端面后反射回去的波幅度与原波形等大，且反向，符合理 论预期。图 c 是两个观测面上的电压随时间的变化情况。由波形可知两端面在电磁波传到 该处时均再现高斯脉冲，到短路边界后脉冲再一次反向返回。
常用的一阶差分格式有三种：
因为上述三种差分格式中， 只有中心差分格式精度是 2 阶， 故可用中心差分代替时间和 空间的差分：
3）Yee 氏网格
为了采用数值差分对电磁场进行计算，在空域进行 Yee 氏网格划分：
由麦克斯伟方程组可以得出电场与磁场之间的 6 组偏微分方程关系式， 用中心差分格式 代替偏导，于是可得下面 6 组差分格式——Yee 差分格式：
n 1/2 Hy
;
Ezn 1 ; Exn 1 ;
matlab模拟的电磁学时域有限差分法
matlab模拟的电磁学时域有限差分法时域有限差分法（FDTD）是一种计算电磁波传播及散射的数值模拟方法。
它是基于麦克斯韦方程组进行仿真的一种方法，而且从计算电磁波传播的实质上来看，FDTD方法是一种求解时域麦克斯韦方程的有限差分方法。
在FDTD方法中，我们将区域空间离散化，并定义电场、磁场等量的格点值。
然后，根据麦克斯韦方程组的时域形式，在各个时刻进行场量的更新。
FDTD方法在实践应用中具有计算时间和空间复杂度低，且适用于复杂的结构和非线性介质等特点，所以在电磁学数值仿真中应用广泛。
我们可以用MATLAB来进行FDTD的电磁学仿真，下面详细介绍MATLAB的使用步骤：1. 建立空间离散化格点在仿真开始前，需要先根据空间大小和仿真目的来建立离散化格点。
对于一个一维的结构，我们可以用以下代码来建立：x = linspace(0,1,N); %建立离散化空间格点Ex = zeros(1,N); %电场，长度为N的全0数组Hy = zeros(1,N); %磁场，长度为N的全0数组其中N为获取离散化格点数量的参数，x为离散化空间格点，Ex和Hy为电场和磁场。
2. 定义电场和磁场边界条件在进行仿真时，需要了解仿真的边界情况并将其定义成特殊的边界条件。
例如，仿真空间内可能存在各种元件、环境等，这些都会对电场和磁场的性质产生影响。
所以，我们需要用特殊边界条件来约束仿真空间内电场和磁场的行为。
在FDTD中，通常采用数值反射边界条件（DNG Boundary）来进行仿真。
例如，在这个边界条件下，在仿真空间内部设置经典的电场边界条件：场强等于零；并在仿真空间外部添加一层基质，该基质的介电常数和磁导率均为负值，并且在该基质中场的强度和方向均反向。
相当于在仿真空间外设置一个虚拟折射界面，能够将场边界反射。
我们设定如下代码：M = 20; % 反射界面层数Ex_low_M1 = 0; %反射界面边界条件Ex_high_M1 = 0; %反射界面边界条件for i = 1:MEx_low_M2(i) = Ex_high_M2(i-1); %反转反射界面内的电场贡献Ex_high_M2(i) = Ex_low_M2(i-1); %反转反射界面内的电场贡献end3. 计算电场的场值FDTD仿真中最核心的内容就是判断时刻要计算的电场场值。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
开题报告
论文题目：基于matlab的时域有限差分法的电磁仿真研究（10分）
学院：电气工程及其自动化学院学号：1103000105姓名：__杨志刚___
一、论文选题的目的和意义（300字以内；15分）
时域有限差分法，因具有多种优点被运用到电磁场理论研究的各个方面，而且其使用成效和应用领域还在迅速扩大和提高，在现代电磁场理论研究中具有很大的重要性和很强的可操作性。
但是同时这种方法也存在一定的缺陷，主要表现在对无边界问题需要吸收边界条件处理，有色散误差，消耗内存大等方面。
本课题在利用时域有限差分法对一些实际的算例进行实验仿真和验证，同时对这种方法在解决实际问题的缺陷进行一定程度的研究和分析。
Matlab作为一种工程仿真工具得到了广泛应用。
用于时域有限差分法，可以简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多的时间。
二、国内外关于该论题的研究现状和发展趋势（500字以内；15分）
时域有限差分方法作为一种典型的全波时域分析方法，因其原理直观、编程简便、实用性强在目前的计算电磁学领域内被人们广泛深入地研究，并取得巨大应用成功的方法。
时域数值技术的一个突出优点是可以给出关于问题空间的丰富的时域信息，而且经过简单的时频变换，即可得到宽带范围的频域信息，相对频域方法显著地节约了计算量。
最近几十年，是电磁场数值计算时域技术蓬勃发展的时期，各具优势和特色的新颖时域算法层出不穷。
但是到目前为止国内关于时域有限差分法中的PML 算法文献较少，其中绝大多数文献集中在综述和应用方面。
而在国际的学报和杂志上对于这方面的文献非常多。
时域有限差分法经过了三十年多年的高速发展之后，仍然还是计算电磁学制高点的研究热潮，而且其应用的范围和成效还在迅速的扩大和提高。
本课题正是利用时域有限差分法的基础理论，利用matlab对一些实际的电磁场问题进行仿真研究。
三、论文的主攻方向、主要内容、研究方法及技术路线（1000字左右；40分）
通过对时域有限差分法理解基础之上，利用matlab仿真软件按照这种方法编程，实现对三种情况下的电磁场情况的仿真研究。
1.时域有限差分法的理论研究：理论研究的主要内容是在麦克斯韦方程组的基础上，利用时域有限差分法的研究思路，全面细致地推导出在直角坐标系中，时域有限差分法二维仿真研究的基本数值计算公式。
并且，根据时域有限差分法的缺点，分析这种方法的数值稳定性条件，吸收边界条件和数值色散分析。
2.时域有限差分法的实现步骤：本课题对于时域有限差分法的研究主要是基于matlab 这一编程软件的仿真研究。
因此，实现仿真研究首先需要搭建matlab 仿真平台，然后从理论上分析时域有限差分法的实现功能和实现步骤，并给出伪代码或者核心思路。
最后对这种方法实现的电磁仿真研究进行评估。
3.基于matlab 的时域有限差分法的验证和应用：选择三种合适的算例，分别是二维方金属柱电磁散射算例、电解槽算例以及等离子体算例，利用理论分析得到的结果在matlab 软件上编程实现仿真进行算例测试，并且对测试的结果进行比较和分析。
时域有限差分法的理论基础
时域有限差分法是在麦克斯韦方程组的微商形式基础上，利用差商代替微商的方法进行离散化处理，将连续型数学模型转化为等价的离散型数学模型，由离散数值构成的离散代数方程组，然后求解出该数学模型的离散数值解。
麦克斯韦方程组
宏观电磁现象的基本规律可以用麦克斯韦方程组表示，如果研究有耗煤质的电磁问题，在求解过程中引入假想的磁荷概念是非常方便的，相应的麦克斯韦方程组式（1）中需引入磁流密度。
已经证明的麦克斯韦方程组修正形式如下所示：ρ
ρ=⋅∇=⋅∇-∂∂-=⨯∇∂∂+=⨯∇D B J t B E t D J H m m ϖϖϖϖϖϖϖ 式中，除了与原来相同参数之外，m J ϖ——磁流密度； m ρ——磁荷密度。
Yee 元胞
Yee 元胞中E ρ，H ρ各分量空间节点与时间步取值的整数和半整数约定如下表所示：
电磁场分量
空间分量取样
时间轴
t 取样 X 坐标 Y 坐标 Z 坐标
E ρ节点 x E
21+i j
k n y E i 21+j k
z E i j 2
1+k H ρ节点 x H i 21+j 2
1+k 21+n y H 2
1+i j 21+k z H 21+i 21+j k
麦克斯韦方程组的时域有限差分法的差分形式
根据（1）中的麦克斯韦方程组以及各项同性介质中的本构关系：H
J E J H B E D m m ρρρρρρρργγμε==== 以E ρ，H ρ为变量，在直角坐标中，展开麦克斯韦第一、二方程，分别为
⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧+∂∂=∂∂-∂∂+∂∂=∂∂-∂∂+∂∂=∂∂-∂∂z γE t z E εy x H x y H y γE t y E εx z H z
x H x γE t x E εz y H y z H ⎪⎪⎪⎪⎩
⎪⎪⎪⎪⎨⎧-∂∂-=∂∂-∂∂-∂∂-=∂∂-∂∂-∂∂-=∂∂-∂∂z H m γt z H μy x E x y E y H m γt y H μx z E z x E x H m γt x H μz y E y z E 令()z,t y,x,f 代表E ρ，H ρ
在直角坐标中的任何一个分量，离散符号取为 ()()()k ,j ,i f
t n ,z k ,y j ,x i f t z y x f n =∆∆∆∆=,,,
直角坐标系中时域有限差分法的二维迭代方法 直角坐标系中时域有限差分法的二维迭代方法，实际上也就是matlab 软件仿真的程序设计的基本思想。
作出matlab 仿真实现的流程图
技术要求：
1．时域有限差分法的程序化设计；
2．二维方金属柱电磁散射算例验证；
3．电解槽算例验证；
4．等离子体算例验证
在对时域有限差分法的理论研究基础之上，利用matlab软件实现程序化的电磁仿真研究。
利用matlab的编程结果分别对二维方金属柱的电磁散射算例，电解槽算例和等离子体算例三种情况进行实验仿真，并对仿真的结果进行比较和分析。
四、论文工作3个月的时间安排（字数不限；10分）
2013年12月20日-2014年02月20日：收集相关文献资料，论文开题。
寒假期间阅读有关计算电磁学和时域有限差分法的相关书籍和资料，对课题主体方法和核心思路有很深的认识。
2014年02月20日-2014年04月12日：基本上完成论文的前两章内容，也就是论文的核心理论，在主要算法的matlab实验上有明确思路和核心伪代码，基本完成后两章的框架。
2014年04月12日-2014年05月22日：完成时域有限差分法的matlab仿真实验验证的全部内容和论文剩余部分的撰写，并对论文进行排版。
2014年05月22日-2014年05月31日：修改论文，完善最后对实验结果优化，宏观调整论文内容，准备结题验收和答辩。
五、论文主要参考文献（字数不限；10分）
[1] 葛德彪,闫玉波.电磁场时域有限差分方法[M].西安:西安电子科技大学出版社,2005.
[2] Yee K S. Numerical solution of initial boundary value problems involving Maxwell equations in isotropic media[J].IEEE Trans, Antennas Propagat, May 1996,AP-14(3): 302-307.
[3] (美)Duane Hanselman,Bruce Littlefield著.精通Matlab 7[M].朱仁峰,译.北京:清华大学出版社, 2006.5.
[4] 杨臻颖.三维散射的FDTD模拟及平面入射波补偿若干问题的研究[J].2006.
指导教师签名：
年月日。