一维FDTD等离子体的Matlab,修正过的一维FDTD等离子体MATLAB代码(公式修正)

修正过的一维FDTD等离子体MATLAB代码

clear;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%系统参数

TimeT=3000;%迭代次数

KE=2000;%网格数目

kc=400;%源的位置

kpstart=600;%等离子体开始位置

kpstop=1000;%等离子体终止位置

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%物理参数

f=1.356e7;%电磁场频率

fpe=1e8;%等离子体频率

c0=3e8;%真空中波速

zdelta=1e-1;%网格大小

dt=zdelta/(2*c0);%时间间隔

u0=2e10;%碰撞频率

wpe=2*pi*fpe;%等离子体圆频率

epsz=1/(4*pi*9*10^9); % 真空介电常数

mu=1/(c0^2*epsz);%磁常数

ex_low_m1=0;

ex_low_m2=0;

ex_high_m1=0;

ex_high_m2=0;

a0=2*u0/dt+(2/dt)^2;

a1=-8/(dt)^2;

a2=-2*u0/dt+(2/dt)^2;

b0=wpe^2+2*u0/dt+(2/dt)^2;

b1=2*wpe^2-8/(dt)^2;

b2=wpe^2-2*u0/dt+(2/dt)^2;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%初始化电磁场

Ex=zeros(1,KE);

Ex_Pre=zeros(1,KE);

Hy=zeros(1,KE);

Hy1=zeros(1,KE);

Hy_Pre=zeros(1,KE);

Dx=zeros(1,KE);

Dx_Pre=zeros(1,KE);

一维FDTD等离子体的Matlab,修正过的一维FDTD等离子体MATLAB代码(公式修正)相关推荐

matlab pdepe函数边界,科学网-使用MATLAB中pdepe函数求解一维偏微分方程-邓浩鑫的博文...
由于自己科研水平较低,记录的各种体会更多的是给自己做个小结,错误之处,欢迎大家指正. 使用MATLAB求解偏微分方程或者方程组,大致有三类方法.第一种是使用MATLAB中的PDE Toolbox,PD ...
fdtd 磁场 matlab,矩形谐振腔电磁场的FDTD分析和Matlab仿真
矩形谐振腔电磁场的FDTD分析和Matlab仿真摘要:目前,电磁场的时域计算方法越来越引人注目.这种方法已经广泛应用到各种电磁问题的分析之中.而将Matlab作为一种仿真工具,用于时域有限差分法,可 ...
MATLAB转移矩阵法仿真,利用转移矩阵和MATLAB求解一维薛定谔方程的一种简捷方法...
引言在半导体量子器件的分析和设计中,需要计算薛定谔方程,转移矩阵方法是常见的求解途径之一[1,2].转移矩阵方法是较为基于抽象的矩阵迭代计算,处理比较琐杂.由于转移矩阵方法的实质仍是求解线性方程组,因 ...
numerov算法matlab,Numerov算法求解一维薛定谔方程研究.doc
Numerov算法求解一维薛定谔方程研究摘要本文主要讲述的是运用MATLAB对一维定态薛定谔方程求解中遇到问题的分析.问题指的是在运用文献[1]所提供的程序,解一维定态薛定谔方程时出现了波函数与 ...
通过有限差分和matlab矩阵运算直接求解一维薛定谔方程,通过有限差分和MATLAB矩阵运算直接求解一维薛定谔方程.doc...
通过有限差分和MATLAB矩阵运算直接求解一维薛定谔方程通过有限差分和MATLAB矩阵运算直接求解一维薛定谔方程文章编号:1672-8785(2010)03-0042-05 通过有限差分和MATL ...
随机游走matlab,从Matlab到Python的一维随机游走
你做了很多不同的事情.在首先,在MATLAB中使用rand,它返回一个介于0和1之间的随机浮点.在python中,使用的是randint,它返回一个随机整数.您正在执行randint(0, 1),这 ...
matlab需要多大运存_提高matlab运行效率
用过Matlab的人都知道,Matlab是一种解释性语言,存在计算速度慢的问题,为了提高程序的运行效率,matlab提供了多种实用工具及编码技巧. 1. 循环矢量化 Matlab是为矢量和矩阵操作而设 ...
matlab的函数要写在哪,matlab函数库在哪
---WORD 格式--可编辑--专业资料--- 在 MATLAB 环境下访问外部函数的共享库文件,必须首先把该库文件加载到内存中.一旦加载成功,就能直接在 MATLAB 中直接请求...... 1 ...
图像分割matlab阈值自动计算,图像分割的阈值算法matlab实现.doc
图像分割的阈值算法matlab实现.doc 图像分割的阈值算法matlab实现[OTSU,1DEntropy,2DEntropy]今天看了几篇论文,实现了一下,没有验证各算法的有效性也没有进行定量比较 ...

一维FDTD等离子体的Matlab,修正过的一维FDTD等离子体MATLAB代码(公式修正)

一维FDTD等离子体的Matlab,修正过的一维FDTD等离子体MATLAB代码(公式修正)相关推荐

最新文章

热门文章