科技写作结课作业(时域有限差分法的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仿真中最核心的内容就是判断时刻要计算的电场场值。
实验二用MATLAB实现线性系统的时域分析
实验二用MATLAB实现线性系统的时域分析线性系统是一种重要的数学模型，用于描述许多自然和工程系统的行为。
在实际应用中，对线性系统进行时域分析是非常重要的，以了解系统的稳定性、响应和性能特性。
MATLAB是一种功能强大的数学软件，被广泛用于线性系统的建模和分析。
首先，我们将介绍线性系统的时域分析的基本概念和方法。
然后，我们将学习如何使用MATLAB进行线性系统的时域分析，并通过具体的例子来演示。
时域分析是研究系统在时间上的响应，主要包括系统的因果性、稳定性、阶数、零极点分布、阻尼特性和幅频特性等。
其中，系统的因果性表示系统的输出只依赖于输入的过去和现在，与未来的输入无关；系统的稳定性表示系统的输出有界，不会无限增长或发散；系统的阶数表示系统差分方程的最高阶导数的次数。
在MATLAB中，线性系统可以用传输函数、状态空间或差分方程的形式表示。
传输函数是输入输出之间的比例关系，常用于分析系统的频率特性；状态空间是通过一组状态变量和状态方程描述系统的，可以用于分析系统的稳定性和阻尼特性；差分方程是通过相邻时刻的输入和输出之间的关系来描述系统的，可以用于分析系统的因果性和稳定性。
下面，我们以传输函数为例，介绍如何在MATLAB中进行线性系统的时域分析。
首先，我们需要定义传输函数。
MATLAB提供了tf函数来定义传输函数，其语法为：G = tf(num, den)，其中num是传输函数的分子多项式的系数，den是传输函数的分母多项式的系数。
接下来，我们可以使用MATLAB中提供的各种函数和命令来进行时域分析。
例如，可以使用step函数来绘制系统的阶跃响应曲线，语法为：step(G)；可以使用impulse函数来绘制系统的冲激响应曲线，语法为：impulse(G)；可以使用initial函数来绘制系统的零状态响应曲线，语法为：initial(G, x0)，其中x0是系统的初始状态。
此外，还可以使用MATLAB中的函数和命令来计算系统的阶数、零极点分布、频率响应等。
时域有限差分法对平面TE波的MATLAB仿真概论
时域有限差分法对平面TE波的MATLAB仿真摘要时域有限差分法是由有限差分法发展出来的数值计算方法。
自1966年Yee 在其论文中首次提出时域有限差分以来，时域有限差分法在电磁研究领域得到了广泛的应用。
主要有分析辐射条线、微波器件和导行波结构的研究、散射和雷达截面计算、分析周期结构、电子封装和电磁兼容的分析、核电磁脉冲的传播和散射以及在地面的反射及对电缆传输线的干扰、微光学元器件中光的传播和衍射特性等等。
由于电磁场是以场的形态存在的物质，具有独特的研究方法，采取重叠的研究方法是其重要的特点，即只有理论分析、测量、计算机模拟的结果相互佐证，才可以认为是获得了正确可信的结论。
时域有限差分法就是实现直接对电磁工程问题进行计算机模拟的基本方法。
在近年的研究电磁问题中，许多学者对时域脉冲源的传播和响应进行了大量的研究，主要是描述物体在瞬态电磁源作用下的理论。
另外，对于物体的电特性，理论上具有几乎所有的频率成分，但实际上，只有有限的频带内的频率成分在区主要作用。
文中主要谈到了关于高斯制下完全匹配层的差分公式的问题，通过MATLAB 程序对TE波进行了仿真，模拟了高斯制下完全匹配层中磁场分量瞬态分布。
得到了相应的磁场幅值效果图。
关键词：时域有限差分完全匹配层MATLAB 磁场幅值效果图目录摘要 (1)目录 (2)第一章绪论 (3)1.1 课题背景与意义 (3)1.2 时域有限差分法的发展与应用 (3)2.1 Maxwell方程和Yee氏算法 (6)2.2 FDTD的基本差分方程 (8)2.3 时域有限差分法相关技术 (10)2.3.1 数值稳定性问题 (10)2.3.2 数值色散 (10)2.3.3 离散网格的确定 (11)2.4 吸收边界条件 (11)2.4.1 一阶和二阶近似吸收边界条件 (12)2.4.2 二维棱边及角顶点的处理 (15)2.4.3 完全匹配层 (17)2.5 FDTD计算所需时间步的估计 (21)第三章MATLAB的仿真的程序及模拟 (23)3.1 MATLAB程序及相应说明 (23)3.2 出图及结果 (26)3.2.1程序部分 (26)3.2.2 所出的效果图 (27)第四章结论 (29)参考文献 (30)第一章绪论1.1 课题背景与意义20世纪60年代以来，随着计算机技术的发展，一些电磁场的数值计算方法逐步发展起来，并得到广泛应用，其中主要有：属于频域技术的有限元法（FEM）、矩量法(MM)和单矩法等；属于时域技术方面的时域有限差分法(FDTD)、传输线矩阵法(TLM)和时域积分方程法等。
有限差分法解决平面弹性问题的matlab建模仿真
有限差分法解决平面弹性问题的matlab建模仿真用有限差分法解决弹性平面问题的建模与仿真的报告建模与仿真的目的 :解决弹性平面问题——偏心拉伸应力问题，并与实验数值理论解相比较。
建模与仿真的基本原理：有限差分法建模与仿真的方法途径：经过理论推导出所求变量所满足的关系，通过Matlab 软件进行数学建模，将物理语言和数学语言转化成计算机语言，进行计算。
应用软件进行绘图，将抽象的数学表示法转化成直观的形象具体的图像，进行展示，使人们对该建模所研究的问题有更清晰而感性的认识。
报告论文结构：（一）基本原理介绍（差分法）（二）偏心拉伸问题的理论推导（三）建模与仿真的具体方法（应用Matlab 软件）（四）结果展示（五）建模程序特点介绍（六）与实验结果、理论结果的对比（七）前进展望与程序改进（八）任务分工*************************************************************** ****************一、基本原理——有限差分法的介绍有限差分法使用差分方程代替微分方程，从而把基本方程微分方程限组和边界条件求解化为代数方程组的求解。
在将偏微分方程和相应的边界条件转换成有限差分方程之前，我们先将函数的导数改用差分表示，从而导出常用的差分公式。
二、关于偏心拉伸应力问题的理论推导对于边界外一行的结点，我们适当的选择一个基点A ，使得φA =0，(δφ/δx)A =0，(δφ/δy)A =0 ,这样我们可以通过推导得到边界点满足的方程：(δφ/δy)B =∫X B A ds (δφ/δx)B =?∫Y B A ds φB =∫(y B ?y)X B A ds+∫(x ?x B )Y BA ds 而对于边界外一圈的虚结点，有公式：φ13=φ9+2h x (δφ/δx)A φ14=φ10+2h y (δφ/δx)B因而我们可以用内点和边界点表示出边界外一圈虚点，由此平面上的任意点的应力函数都可以用平面上其他十二个点的应力函数表示，列出这些方程，由外及内，依次解出各点的应力函数。
(完整版)利用MATLAB进行时域分析
自动控制原理与系统课程实验报告实验题目：利用MATLAB进行时域分析班级：机电1131班姓名：刘润学号：38号一、实验目的及内容时域分析法是一种直接在时间域中对系统进行分析的方法，具有直观、准确的优点，并且可以提供系统时间响应的全部信息。
在此实验中，主要介绍时域法进行系统分析，包括一阶系统、二阶系统以及高阶系统,以及系统的性能指标。
通过实验，能够快速掌握、并利用MATLAB及控制系统箱对各种复杂控制系统进行时域分析。
二、实验设备三、实验原理典型的二阶系统在不同的阻尼比的情况下，它们的阶跃响应输出特性的差异是很大的。
若阻尼比过小，则系统的振荡加剧，超调量大幅度增加；若阻尼比过大，则系统的响应过慢，又大大增加了调整时间，下面通过此实验课题分析输出响应变化规律：已知二阶振荡环节的传递函数为：G（s）=ωn*ωn/（s*s+2*ζ*ωn*s+ωn*ωn），其中ωn=0.4，ζ从0变化到2，求此系统的单位阶跃响应曲线，并分析当ζ发生变化时，二阶系统的响应有什么样的变化规律。
四、实验步骤编出程序如下图：五、实验结果画出图表如下图：六、结果分析（1）当ξ=0（无阻尼）（零阻尼）时：无阻尼时的阶跃响应为等幅振荡曲线。
如图ξ=0曲线。
（2）当0<ξ<1（欠阻尼）时：对应不同的ξ，可画出一系列阻尼振荡曲线，且ξ越小，振荡的最大振幅愈大。
如图ξ=0.4曲线。
(3)当ξ=1（临界阻尼）时：临界阻尼时的阶跃响应为单调上升曲线。
如图ξ=1曲线。
（4）当ξ>1（过阻尼）时：过阻尼时的阶跃响应也为单调上升曲线。
不过其上升的斜率较临界阻尼更慢。
如图ξ=1.6曲线七、教师评语。
有限差分法求解偏微分方程MATLAB
南京理工大学课程考核论文课程名称：高等数值分析论文题目：有限差分法求解偏微分方程*名：**学号： 1成绩：有限差分法求解偏微分方程一、主要内容1.有限差分法求解偏微分方程，偏微分方程如一般形式的一维抛物线型方程：22(,)()u uf x t t xαα∂∂-=∂∂其中为常数具体求解的偏微分方程如下：22001(,0)sin()(0,)(1,)00u u x t x u x x u t u t t π⎧∂∂-=≤≤⎪∂∂⎪⎪⎪=⎨⎪⎪==≥⎪⎪⎩2.推导五种差分格式、截断误差并分析其稳定性；3.编写MATLAB 程序实现五种差分格式对偏微分方程的求解及误差分析；4.结论及完成本次实验报告的感想。
二、推导几种差分格式的过程：有限差分法（finite-difference methods ）是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。
有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。
推导差分方程的过程中需要用到的泰勒展开公式如下：()2100000000()()()()()()()......()(())1!2!!n n n f x f x f x f x f x x x x x x x o x x n +'''=+-+-++-+- (2-1)求解区域的网格划分步长参数如下：11k k k kt t x x h τ++-=⎧⎨-=⎩ (2-2) 2.1 古典显格式2.1.1 古典显格式的推导由泰勒展开公式将(,)u x t 对时间展开得 2,(,)(,)()()(())i i k i k k k uu x t u x t t t o t t t∂=+-+-∂ (2-3) 当1k t t +=时有21,112,(,)(,)()()(())(,)()()i k i k i k k k k k i k i k uu x t u x t t t o t t tuu x t o tττ+++∂=+-+-∂∂=+⋅+∂ (2-4)得到对时间的一阶偏导数1,(,)(,)()=()i k i k i k u x t u x t uo t ττ+-∂+∂ (2-5) 由泰勒展开公式将(,)u x t 对位置展开得223,,21(,)(,)()()()()(())2!k i k i k i i k i i u uu x t u x t x x x x o x x x x∂∂=+-+-+-∂∂ (2-6)当11i i x x x x +-==和时，代入式(2-6)得2231,1,1122231,1,1121(,)(,)()()()()(())2!1(,)(,)()()()()(())2!i k i k i k i i i k i i i i i k i k i k i i i k i i i iu uu x t u x t x x x x o x x x x u u u x t u x t x x x x o x x x x ++++----⎧∂∂=+-+-+-⎪⎪∂∂⎨∂∂⎪=+-+-+-⎪∂∂⎩(2-7) 因为1k k x x h +-=，代入上式得2231,,22231,,21(,)(,)()()()2!1(,)(,)()()()2!i k i k i k i k i k i k i k i ku uu x t u x t h h o h x xu u u x t u x t h h o h x x +-⎧∂∂=+⋅+⋅+⎪⎪∂∂⎨∂∂⎪=-⋅+⋅+⎪∂∂⎩ (2-8) 得到对位置的二阶偏导数2211,22(,)2(,)(,)()()i k i k i k i k u x t u x t u x t u o h x h+--+∂=+∂ (2-9) 将式(2-5)、(2-9)代入一般形式的抛物线型偏微分方程得21112(,)(,)(,)2(,)(,)(,)()i k i k i k i k i k i k u x t u x t u x t u x t u x t f x t o h h αττ++---+⎡⎤-=++⎢⎥⎣⎦(2-10)为了方便我们可以将式(2-10)写成11122k kk k k k i i i i i i u u u u u f h ατ++-⎡⎤--+-=⎢⎥⎣⎦(2-11) ()11122k k k k k k i i i i i i u u uu u f h τατ++----+= (2-12)最后得到古典显格式的差分格式为()111(12)k k k k k i i i i i u ra u r u u f ατ++-=-+++ (2-13)2r h τ=其中，古典显格式的差分格式的截断误差是2()o h τ+。
有限差分法及matlab实现
矩阵（数组）是计算机中重要的数据结构，为了方便用矩阵去存储数据，我们网格去划分空间，从而不 仅使方程化为简单的有限差分形式，而且这样的数据类型在 计算机中易于储存和运算。那么 h=k=l=1，并且令 f(x,y,z)=u(x,y,z) 就有
U i1, j , k U i−1, j , k U i , j 1, kU i , j −1, kU i , j , k1 U i , j , k−1
上述六式相加
2
3
U i1, j , k U i−1, j , k U i , j 1, kU i , j −1, kU i , j , k1 U i , j , k−1 2 2 2 ∂ U ∂ U ∂ U ⋯ = 6U i , j , k ∂ x2 ∂ y 2 ∂ z 2 2 代入 ∇ U = ρ U i1, j , k U i−1, j , k U i , j 1, kU i , j −1, kU i , j , k1 U i , j , k−1 = 6U i , j , k∇ 2 U ⋯=6U i , j ,k ρ ⋯
此即拉普拉斯方程的有限差分形式。
----- (2)式
这里，我们通过有限差分的方法，把偏微分方程在三阶精度下简化为形式易于计算的代数方程。从而使 之易于在计算机上实现。
注：有时我们需要解二维静电场，则方程退化为：
U i1, j U i−1, j U i , j 1 U i , j −1 =4U i , j
改为
v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v2(i-1,j)+v2(i,j-1))/4 解释执行后，k=222;
时域有限差分和时域有限元电磁数值计算的研究的开题报告
时域有限差分和时域有限元电磁数值计算的研究的开题报告一、课题背景和研究意义随着计算机技术的不断发展和应用领域的不断扩大，电磁数值计算技术在电磁场求解中得到了广泛应用。
其中，时域有限差分（FDTD）和时域有限元（FEM）方法是常用的电磁数值计算方法，具有计算精度高、适用范围广等优点，广泛应用于电磁场计算、天线设计、散射特性计算、电磁兼容性分析等领域。
本文将重点研究FDTD和FEM方法在电磁数值计算中的应用，探讨两种方法的原理和计算特点，分析它们的优缺点和适用范围，并对它们在电磁场计算、天线设计等方面的应用进行研究和分析。
该研究对于提高电磁数值计算的准确性和效率，推动电磁场计算技术的发展具有重要意义。
二、研究内容和研究方法1. 理论分析：对FDTD和FEM方法的理论基础进行深入分析，探讨它们的数学原理、计算特点以及在电磁场计算中的应用。
2. 编程实现：通过编程实现FDTD和FEM方法，将电磁场的边值问题转化为差分方程或有限元方程，在计算机上进行求解，并得到电磁场的分布情况。
3. 数值模拟与分析：通过数值模拟和实验验证，对比分析FDTD和FEM方法的精度和计算效率，并对它们在电磁场计算、天线设计等方面的应用进行研究和分析。
三、预期研究成果1. 对FDTD和FEM方法的原理和计算特点进行深入分析，分析它们的优缺点和适用范围。
2. 通过数值模拟和实验验证，对比分析FDTD和FEM方法的精度和计算效率，并探讨在电磁场计算、天线设计等方面的应用。
3. 提出针对FDTD和FEM方法的优化算法，进一步提高电磁场计算的精度和效率。
四、研究计划和进度安排1. 第1-2个月：查阅相关文献，深入理解FDTD和FEM方法的原理和计算特点。
2. 第3-4个月：编写FDTD和FEM的数值求解算法，通过计算机程序对其进行验证。
3. 第5-7个月：通过数值模拟和实验验证，对比分析FDTD和FEM方法的精度和计算效率，并对它们在电磁场计算、天线设计等方面的应用进行研究和分析。
电磁场实验一_有限差分法的matlab实现
电磁场与电磁波实验报告实验项目：_______有限差分法__ ____ 班级：_____ __12电子2 ____ __ 实验日期：__2014年12月23日姓名：___ _ __陈奋裕 __ __ 学号：___ ___1215106003 _____ 组员姓名：___ _ __ __ __ 组员学号：___ ___ _____ 指导教师：_ ____张海 ______一、实验目的及要求1、学习有限差分法的原理与计算步骤；2、学习用有限差分法解静电场中简单的二维静电场边值问题；3、学习用Matlab 语言描述电磁场与电磁波中内容，用matlab 求解问题并用图形表示出了，学习matlab 语言在电磁波与电磁场中的编程思路。
二、实验内容理论学习：学习静电场中边值问题的数值法中的优先差分法的求解知识； 实践学习：学习用matlab 语言编写有限差分法计算二维静电场边值问题；三、实验仪器或软件电脑（WIN7）、Matlab7.11四、实验原理 基本思想是把连续的定解区域用有限个离散点构成的网格来代替， 这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似， 积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组 ， 解此方程组就可以得到原问题在离散点上的近似解。
然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。
简单迭代法： 这一方法的求解过程是，先对场域内的节点赋予迭代初值(0),i j ϕ，这里上标(0)表示0次（初始）近似值。
然后按Laplace 方程(k 1)(k)(k)(k)(k),1,,11,,11[]4i j i j i j i j i j ϕϕϕϕϕ+--++=+++(i,j=1,2,…) 进行反复迭代（k=0，1,2,…）。
若当第N 次迭代以后，所有的内节点的相邻两次迭代值之间的最大误差不超过允许范围，即(N)(N-1),,max|-|<W i j i j ϕϕ这里的W 是预设的允许误差，此时即可终止迭代，并将第N 次迭代结果作为内节点上电位的最终数值解。
- 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.
指导教师签名：
年月日。