光子晶体matlab,Matlab光子晶体程序讨论

clear all

clc;

a=1.0; %晶格常数。

R=0.3*a;

f=pi*R^2/a^2;

Nc=12; %背景介质。

Nh=1.0; %填充介质介电常数。

omiga_pcs=[];

W2=[];

W1=[];

omg=[];

%求解omiga_2D结束

N2d=sqrt((Nc^2*(1+f)+Nh^2*(1-f))*Nh^2)./(Nc^2*(1-f)+Nh^2*(1+f));

%设置一些常数。

n2=1.0;

n3=1.0;

m=0;

h=0.5;

c=1;

mu0=4*pi*1.0e-7;

eps0=8.85*1e-12;

c=1/sqrt(mu0*eps0);

xx=load('F:\郭丽丽\新程序\shuai2.txt');

omiga_pcs=0.1:0.001:0.835;

num=size(omiga_pcs)

dert=0.001;

Fid=fopen('shuai3.txt','w');

for ii=1:736

for jj=1:1000

omiga_pcs(ii)=omiga_k(ii,2)-jj*dert;

Npcs(ii)=N2d*omiga_k(ii,2).\omiga_pcs(ii);

W1(ii)=(sqrt(N2d^2-Npcs(ii)^2)).*omiga_pcs(ii).*h/c ;

W2(ii)=m*pi+atan(sqrt((Npcs(ii)^2-n2^2).\(N2d^2-Npcs(ii)^2)))...

+atan(sqrt((Npcs(ii)^2-n3^2).\(N2d^2-Npcs(ii)^2)));

dert2(ii)=abs(abs(W1(ii))-abs(W2(ii)));

if dert2(ii)<=1.0e-7

omg(ii)=omiga_pcs(ii);

else

omg(ii)=0;

end

pq=[ii,omiga_pcs(ii)];

fprintf(Fid,'%6d %6.4f\r\n',pq);

plot(omiga_k(:,1),omiga_k(ii,2),'.r');

hold on;

plot(omiga_k(:,1),omiga_pcs(ii),'g+');

axis([0 0.5 0.1 0.5]);

title('光子晶体薄板能带计算图(TE模)');

xlabel('wave vector(2\pi/a)');

ylabel('Normalized frequency(c/a)');

set(gca,'Xtick',1:31);

text(0,-0.2,'\Gamma');

text(10,-0.2,'X');

text(20,-0.2,'M');

text(30,-0.2,'\Gamma');

% set(gca,'XtickLable',{'\Gamma','X','M','\Gamma'});

% toc;

光子晶体matlab,Matlab光子晶体程序讨论相关推荐

matlab ask fsk调制解调,Matlab的ASK-FSK--PSK程序仿真
<Matlab的ASK-FSK--PSK程序仿真>由会员分享,可在线阅读,更多相关<Matlab的ASK-FSK--PSK程序仿真(12页珍藏版)>请在人人文库网上搜索. 1. ...
牛顿法的matlab实现例题,【MATLAB|MATLAB牛顿插值法例题与程序Word版】
『易坊知识库摘要_MATLAB|MATLAB牛顿插值法例题与程序Word版』x12345678y22.523.324.421.7025.228.524.825.4二.数学原理假设有n+1个不同的节点及 ...
【 MATLAB 】DFT性质讨论（一）线性、循环反转、共轭与实序列的对称性的MATLAB实现
上篇博文通过在理论上讨论了DFT的三个性质:[ MATLAB ]DFT性质讨论(一)线性.循环反转与共轭分别讨论: 一.线性给出一个例子,给出x1和x2,x3 = 0.3*x1+0.8*x2; 之 ...
Matlab之M程序与M函数
M程序点击MATLAB指令窗工具条上的NewFile图标 ,就可打开如图所示的MATLAB文件编辑调试器MATLABEditor/Debugger.其窗口名为untitled,用户即可在空白窗口中编 ...
MATLAB调用C++程序
MATLAB调用C++程序_仰起脸笑得像满月-CSDN博客_matlab调用c++程序 Matlab与C/C++混合编程 Matlab调用C函数 - NZACH - 博客园
matlab debug出现k,MATLAB下的程序调试
文章主要内容摘自<MATLAB 7.0从入门到精通>,求是科技编著. 程序错误一般分为两种:语法错误和逻辑错误.对于语法错误通常MATLAB会报错,并指出错误所在位置方便用户纠正.对于逻辑 ...
matlab调用Java程序时出现 Java.lang.OutOfMemoryErrot: GC overhead limit exceeded
matlab调用Java程序时出现 java.lang.OutOfMemoryError: GC overhead limit exceeded JDK1.6.0_37和JDK_1.7.0_60版本, ...
matlab程序改为m文件名,在MATLAB中，程序文件的扩展名为.m，所以程序文件也称为M文件...
在MATLAB中,程序文件的扩展名为.m,所以程序文件也称为M文件答:√ 磷酸果糖激酶-2催化6-磷酸果糖生成的产物是答:2,6-二磷酸果糖 <人类行为的经济学分析>的作者是( ). ...
matlab安装无效距离过远,求助matlab的远程序
求助matlab的远程序 function varargout = a1(varargin) % A1 M-file for a1.fig % A1, by itself, creates a new ...
matlab实用小程序,matlab100常用小程序
matlab编程实例100例_计算机软件及应用_IT/计算机_专业资料.1-32 是:图形应用篇 33-66 是:界面设计篇 67-84 是:图形处理篇 85-100 是:数值分析篇实例 1.... ...

光子晶体matlab,Matlab光子晶体程序讨论

光子晶体matlab,Matlab光子晶体程序讨论相关推荐

最新文章

热门文章