matlab 2d pml,二维TM波,PML边界,matlab程序。但是仿真结果有问题,求高手。
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
clc;
clear;
ft=1.5;
fs=20;
NX=100;%X方向网格数
NY=150;%Y方向网格数
NA=20;%吸收层厚度
bxrhx=1:100;
bxirhx=2:150;
byrhy=2:100;
byirhy=1:150;
bxrez=1:100;
byrez=1:150;
hx=zeros(100,151);
hy=zeros(101,150);
ezx=zeros(100,150);
ezy=zeros(100,150);
sigmamy=hx;
sigmamx=hy;
sigmax=ezx;
sigmay=ezy;
sigmamax=0.2;
sigma=sigmamax*(NA:-1:1)/NA;
sigmam=sigmamax*(2*NA-1:-2:1)/2/NA;
sigmamx([1:NA,end:-1:end-NA+1],:)=repmat([sigmam,sigmam]',1,150);
sigmax([1:NA,end:-1:end-NA+1],:)=repmat([sigma,sigma]',1,150);
sigmamy(:,[1:NA,end:-1:end-NA+1])=repmat([sigmam,sigmam],100,1);
sigmay(:,[1:NA,end:-1:end-NA+1])=repmat([sigma,sigma],100,1);
f1hx=exp(-sigmamy(bxrhx,bxirhx));
f1hy=exp(-sigmamx(byrhy,byirhy));
f2hx=zeros(size(hx));
f2hy=zeros(size(hy));
f2hx(sigmamy == 0)=1/ft;
f2hy(sigmamx==0)=1/ft;
f2hx(sigmamy~=0)=(1-exp(-sigmamy(sigmamy~=0)))./sigmamy(sigmamy~=0)/ft;
f2hy(sigmamx~=0)=(1-exp(-sigmamx(sigmamx~=0)))./sigmamx(sigmamx~=0)/ft;
f2hx=f2hx(bxrhx,bxirhx);
f2hy=f2hy(byrhy,byirhy);
f1ezx=exp(-sigmax);
f1ezy=exp(-sigmay);
f2ezx=zeros(size(ezx));
f2ezy=zeros(size(ezy));
f2ezx(sigmax==0)=1/ft;
f2ezy(sigmay==0)=1/ft;
f2ezx(sigmax~=0)=(1-exp(-sigmax(sigmax~=0)))./sigmax(sigmax~=0)/ft;
f2ezy(sigmay~=0)=(1-exp(-sigmay(sigmay~=0)))./sigmay(sigmay~=0)/ft;
centerx=floor((1+NX)/2);
centery=floor((1+NY)/2);
for nn=1:150
ezx(centerx,centery)=ezx(centerx,centery)+sin(2*pi*nn/ft/fs);
ezy(centerx,centery)=ezy(centerx,centery)+sin(2*pi*nn/ft/fs);
hx(bxrhx,bxirhx)=f1hx.*hx(bxrhx,bxirhx)-f2hx.*(ezx(bxrhx,bxirhx)-ezx(bxrhx,bxirhx-1)+ezy(bxrhx,bxirhx)-ezy(bxrhx,bxirhx-1));
hy(byrhy,byirhy)=f1hy.*hx(byrhy,byirhy)+f2hy.*(ezx(byrhy,byirhy)-ezx(byrhy-1,byirhy)+ezy(byrhy,byirhy)-ezy(byrhy-1,byirhy));
ezx=f1ezx.*ezx+f2ezx.*(hy(bxrez+1,:)-hy(bxrez,:));
ezy=f1ezy.*ezy-f2ezy.*(hx(:,byrez+1)-hx(:,byrez));
end
matlab 2d pml,二维TM波,PML边界,matlab程序。但是仿真结果有问题,求高手。相关推荐
- 【图像分割】基于matlab灰狼算法二维oust图像分割【含Matlab源码 2574期】
⛄一. 灰狼算法简介 1 前言 灰狼优化算法(Grey Wolf Optimizer,GWO)由澳大利亚格里菲斯大学学者 Mirjalili 等人于2014年提出来的一种群智能优化算法.该算法受到了灰 ...
- N圆最密堆积、最小外接正方形的matlab求解(二维、三维等圆Packing 问题)
圆形最密堆积.最小外接正方形的matlab求解(二维.三维等圆Packing 问题) 0 前言 1 N个圆的最小外接正方形求解 2 N个球的最小外接立方体求解 惯例声明:本人没有相关的工程应用经验,只 ...
- 【A_star二维路径规划】基于matlab A_star算法无人机二维路径规划(起终点障碍物可设置)【含Matlab源码 1321期】
⛄一.获取代码方式 获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[A_star二维路径规划]基于matlab A_s ...
- Matlab系列之二维图形(上)
Matlab系列之二维图形(上) 简要 绘制基础 plot 文字标注 (1) 添加图形标题 (2)添加坐标轴标注 (3)图例 (4)文本注释 示例 程序 结果 线型.点型和颜色 坐标轴设置 示例 结果 ...
- 利用MATLAB来绘制二维随机变量的联合概率密度图像
本文档通过MATLAB来绘制二维随机变量(X,Y)的联合概率密度函数. 第一种类型,X服从标准正态分布,Y服从均匀分布. [例题]已知随机变量X与Y相互独立,X~N(0,1);Y在区间[0,2]上服从 ...
- matlab contourf设置范围,如何使用matlab contourf绘制二维决策边界
我知道这已经有一段时间了,但我会试一试,以防其他人提出这个问题. 假设我们有一个2D训练集来训练SVM模型,换句话说特征空间是一个2D空间.我们知道,一个核SVM模型导致得分(或决定)的形式的函数: ...
- 偏微分方程的数值解(五): 二维状态空间的偏微分方程的 MATLAB 解法
偏微分方程的数值解系列博文: 偏微分方程的数值解(一):定解问题 & 差分解法 偏微分方程的数值解(二): 一维状态空间的偏微分方程的 MATLAB 解法 偏微分方程的数值解(三): 化工应用 ...
- matlab 如何画二维图形,Matlab 学习 画图篇 一 二维图形
matlab给绘制二维图形提供了很多的函数,把一些绘制二维图形的基本函数做成一张表,如下图所示: 我就按照表的顺序一一记录一些个函数的简略用法. 首先是 1.plot函数 plot函数有很多重载方法, ...
- matlab三维选取二维,基于Matlab绘制二维和三维图形以及其他图形控制函数的使用方法...
Matlab绘图 强大的绘图功能是Matlab的特点之一,Matlab提供了一系列的绘图函数,用户不需要过多的考虑绘图的细节,只需要给出一些基本参数就能得到所需图形,这类函数称为高层绘图函数.此外,M ...
- matlab中饼图绘制程序,MATLAB中绘制二维饼图的函数是(? ? ) (1.5分)
MATLAB中绘制二维饼图的函数是(? ? ) (1.5分) 答:pie 计算机的内存储器比外存储器读写速度 ______ 答:快 新民主主义革命时期,提出了反帝反封建这个任务就等于真正懂得了什么叫反 ...
最新文章
- python科学计算笔记(九)pandas中DataFrame数据操作函数
- hdu2094产生冠军
- DOCKER - 容器抓包
- 企业应用程序中需要捕获的5大Java性能指标
- [react] 为什么说React中的props是只读的?
- 如何用 Python 实现超级玛丽的人物行走和碰撞检测?
- mysql 柱状图统计_js/jquery 进行动态统计 各种柱状图 饼状图 线条图 等
- “单口相声”回归!罗永浩要开发布会了:黑科技!不售票!
- Oracle、 Mysql 、 SQLserver 分页查询
- 未将对象引用设置到对象实例(转)
- [HDU1754]I Hate It线段树裸题
- kodi pvr 不能安装_Kodi看电视直播教程 安装PVR IPTV Simple Client播放m3u8直播源
- 服务器ipsan虚拟化,挂载IPSAN存储+链接克隆模板(修正版)
- hex文件、bin文件、axf文件的区别?
- STP-16-根防护,BPDU防护和BPDU过滤
- 【无标题】打印水仙花(pyth)
- [Games101] Lecture 03-04 Transformation
- 叩响港交所大门,KK集团能否成为“中国版秋叶原”?
- android遥控器适配
- 导出富文本格式word
热门文章
- 《Android群英传:神兵利器》— Android 书籍
- java实现阴历、阳历互相转换
- Android【报错】Description Resource Path Location Type Unparsed aapt error(s)! Cheheck the console for o
- 7-35 混合类型数据格式化输入 (5 分)
- L1-003 个位数统计 (15 分)—团体程序设计天梯赛
- python自动登录百度_Python实现自动登录百度空间的方法
- 基于Cocos2d-x开发guardCarrot--8 《保卫萝卜2》关卡选择场景开发
- YANG-SUITE安装步骤
- 防火墙双机热备配置实例(二)
- 静默安装Oracle11G