自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+
  • 博客(48)
  • 资源 (2)
  • 收藏
  • 关注

原创 matlab中一种直观的构建矩阵的方式

matlab中一种构建矩阵的技巧

2024-04-22 21:21:49 78

原创 群折射率$ n_g$与有效折射率$ n_{eff}$的关系详细推导(假设$ \lambda_0$为真空入射光波长)

群折射率和有效折射率关系公式的详细推导

2023-08-15 15:01:50 1317

原创 matlab设置图像的坐标轴刻度朝外

matlab图像的坐标轴刻度设置朝外

2023-08-01 20:55:04 862

原创 平板光波导中导模的(注意不是泄露模)传播常数β的matlab计算(验证了是对的)

(Slab waveguide)平板光波导的传播常数的计算

2023-07-29 20:58:59 1123

原创 用平面波展开法(PWEM)计算光子晶体条状超原胞的投影能带

(平面波展开法)PWEM算投影能带

2023-05-13 11:04:53 546

转载 Universitatea Babeş-Bolyai网站找到的一些计算电磁学资料(包括PWEM、FDTD、TTM)

一些计算电磁学的代码

2023-05-10 21:27:57 90

原创 美国老自治领大学的Dr. ShangPingGuo分享的平面波展开法代码合集文档

Dr.ShangPingGuo的平面波展开法代码合集文档

2023-05-10 20:58:00 94

原创 二维体光子晶体的平面波展开法代码2

这篇文章接前面《二维体光子晶体的平面波展开法代码》这个是实际网格划分,各标记的含义上篇文章有说明下面是傅里叶系数恢复的相对介电常数分布(虚部),看起来非常大,那是因为imag(EpsMat)趋于0的原因。下面这个图是实空间相对介电常数分布的伪彩图:下面这个图是TE能带图,可以看到和FEM的结果还是有点偏差,但实际上只要加密网格同时遵守前一篇文章所述的nG参数和precisStruct参数的关系就可以获得误差更小的能带图。

2023-05-10 15:35:15 248

原创 matlab实现DFT(取了一些数据计算结果和fft相差很小,但未详细验证)

DFT的matlab实现,和fft进行了少量对比误差非常小,但没详细验证,慎重。

2023-05-06 21:56:02 121

原创 DFT振幅

和离散傅里叶变换有关

2023-05-06 21:35:10 73

原创 matlab根据中心和边长构建正方形

【代码】matlab根据中心和边长构建正方形。

2023-05-04 16:02:18 909

原创 二维体光子晶体的平面波展开法代码

二维平面波展开法

2023-05-04 09:25:48 1037

原创 (PWEM前置知识)解析公式求相对介电常数分布函数的 Fourier 展开系数以及基于此的绘制实空间原胞的相对介电常数分布

(PWEM前置知识)解析公式求相对介电常数分布函数的 Fourier 展开系数以及基于此的绘制实空间原胞的相对介电常数分布

2023-05-03 15:21:32 335

原创 实空间四极矩不变量、体偶极矩不变量的计算

matlab实空间体四极矩计算(参考了github的一个代码)

2023-04-17 09:14:08 296

原创 Matlab绘图为什么设置了相同的字体大小结果出的图字体大小还是不一样

matlab相同字号设置但最终图片字号大小不一致(matlab2018b中已验证过)

2023-02-26 10:06:50 897 1

原创 matlab同一条曲线上plot不同x位置颜色不同

matlab同一条曲线上不同颜色,注意曲线是光滑的,不是scatetr

2022-10-01 11:37:06 2407

原创 寻找子集一维数组在母集一维数组中的位置(注意只适用于特定关系数组,未详细测试)

只适用于特定关系的数组

2022-07-30 09:45:05 219

原创 一维列数组分类加强版(可以对重复元素使用,但不能含有0元素)

一维数组按间隔分类,请仔细阅读注释

2022-07-28 10:14:28 101

原创 检查一维数组中的重复元素以及每个重复元素的重复次数(结果以矩阵形式给出)

matlab函数,检测一列数组中的重复元素和重复次数,结果以矩阵形式给出,第一行是重复元素,第二行是对应元素的重复次数

2022-07-27 21:38:42 507

原创 对输入的一维数组(按顺序列且不重复)按容差tol分类(未详细测试)

主要用于识别简并能带

2022-07-27 18:45:15 118

原创 一点绕另一点旋转angle度的matlab子程序(未详细测试)

某点绕某点旋转angle度的matlab程序

2022-07-25 14:40:57 326

原创 运用小矩阵作为元素构建大矩阵(matlab函数,未详细测试)

如题

2022-07-01 11:09:19 1449

原创 一个有趣的极限(会在求标量格林函数的Ewald级数时且处理它的奇点的时候会遇到)

如题目所示,这就是要表达的意思

2022-06-29 16:23:02 69

原创 一个和二维泊松求和有关的公式(推导Ewald级数中有用,运用了2D泊松求和公式,傅里叶变换的位移性质)

如题所示

2022-06-28 21:24:57 345

原创 matlab复误差函数的高效计算(有官方背书)

matlab复误差函数的高效计算,可以和mathematica的官方函数对比验证下,有权威官方背书,是参考网站的代码改写而来。

2022-06-27 17:40:46 421

原创 赫尔维兹_勒奇超越函数(matlab自编函数)

function [out] = Lerch_fun(z,s,v)%Lerch_fun(exp(1j*2),3,0.1)% Lerch transcendent%表达式:\Phi(z, s, v)=\sum_{n=0}^{\infty} \frac{z^{n}}{(n+v)^{s}}%syms tf=@(t)((t.^(s-1)).*exp(-v*t))./(1-z.*exp(-t));out=(1/gamma(s))*(integral(f,0,Inf));end和Mathematic

2022-03-01 21:29:14 796

原创 带选择程序(未测试)

function [band_sto]=choice_band(vec,val,band_index)%根据band_index选择能带%vec,val是最原始的vec,valval=real(sqrt(diag(val)));%w/c,,real可改,决定画能带的实部还是虚部[~,order]=sort(val);%w/cband_sto=zeros(size(vec,1),length(band_index));for count=1:length(band_index)tp=band.

2022-01-05 17:14:49 184

原创 三角格子BZ离散代码(T-K-M-T)

function [k0]=BZ_dispersed()global a0delta=0.01;%可调参数T=[0,0];M=[0,2*pi/(sqrt(3)*a0)];K=[2*pi/(3*a0),2*pi/(sqrt(3)*a0)];TK=K-T;tk=TK(1)+TK(2)*1i;list_tk=(T(1)+T(2)*1i)+[0:delta:1]*tk;KM=M-K;km=KM(1)+KM(2)*1i;list_km=(K(1)+K(2)*1i)+[0:delta:

2022-01-01 17:31:14 233

原创 matlab画圆函数

function [] = circle(x,y,r)%画圆函数% circle(0,0,4) ;grid on% x y是中心,r是半径rectangle('Position',[x-r,y-r,2*r,2*r],'Curvature',[1,1])axis equal%为了修饰曲线的颜色,宽度,圈盘填充颜色等,可以设置其他参数等,例如%'edgecolor','b',其中edgecolor表示边框颜色,后面的b是颜色参数值;%facecolor','r',其中facecolor表示内部填

2021-12-31 11:24:04 6056

原创 matlab angle函数

angle采用的算法:angle(z) = imag(log(z)) = atan2(imag(z), real(z)).

2021-12-29 11:19:01 2181

原创 匿名函数支持输入向量,匿名函数可以组成向量,

################

2021-12-23 19:32:16 267

原创 一个基本成功的TM模式FDFD代码

%Main_FDFD(0)%TM mode%version 2.0function []=Main_FDFD(para0)clearvars -except para0clc;close all%para0为1则画场图%c_const=299792458;%单位mglobal R a0R=1e-3;%圆柱半径 ,单位mPx=5e-3;%超原胞长度Py=5e-3;%超原胞宽度a0=Px;%晶格常数x_center=0.5*Px;%圆柱中心y_center=0.5*Py;%圆柱中心

2021-12-23 16:55:20 605

原创 一段判断点是否在多边形内的matlab代码

function [out]=shape_f(xx,yy,xp,yp)%xx、yy是输入的要判断的点,%xp 、yp是多边形的顶点(顺时针或逆时针,首尾相连)%xp、yp是向量%如果在多边形内或边上,out为1,否则为0out=double(inpolygon(xx,yy,xp,yp));end...

2021-12-23 16:52:38 488

原创 comsol提取的csv数据画能带图

clear;clc;c_const=299792458;%单位mX=csvread('F:\table_desk\COMSOL.VS.紧束缚\test_band1.csv',1);Y=csvread('F:\table_desk\COMSOL.VS.紧束缚\test_band2.csv',1);figure(2)scatter((X),abs(real(Y)).*(5e-3)/(c_const),'r.')axis tight

2021-12-23 16:48:39 715

原创 画格点程序

function []=point_figure()a1=[10,0];a2=[0,10];site=[];for m=-6:6 for n=-6:6 vec=m*a1+n*a2; site=[site;sq(vec(1),vec(2),4)]; endendscatter(site(:,1),site(:,2),'bo','filled') axis equal endfunction [ou

2021-11-01 17:12:07 93

原创 matlab确定特定一行在矩阵中的位置

A=[122 51 41]B=[122 49 41; 122 51 41; 121 49 42; 122 49 42][~,ind]=ismember(A,B,'rows')

2021-10-24 09:57:33 2041

原创 一个clear的高级用法

clearvars -except stoM

2021-10-22 16:06:12 94

原创 Floquet-傅里叶方法里面一个积分的验证小程序

%Floquet-傅里叶-方法里面积分的一个验证小程序function []=main_jifen()clear;clcfor jj=-100:1:100 disp(['jj=',num2str(jj),';',newline]) mm=jifen(jj) disp([newline,'--------------'])endendfunction [out1]=jifen(delta)syms vac w t k T%delta=3;m=2*vac*c.

2021-09-16 21:19:46 227

原创 一个matlab小函数,寻找矩阵最小n个值的位置(没有详细测试)

function [out]=find_M(M,n)%这个函数是找到矩阵M中最小的n个点的位置[u,~]=size(M);sto_v=[];for ii=1:u sto_v=[sto_v,M(ii,:)];endsort_v=unique(sort(sto_v));if length(sort_v)<n sto2=sort_v(:); disp('n超过了矩阵不重复元素数')elsesto2=sort_v(1:n);endclear iisto_si

2021-09-03 10:51:25 1742

原创 正倒空间基矢转换

matlab代码:%正空间、倒空间转换clear;clc;close all%%%%%%%%%%%%%%a=sym(1,'f');%晶格常数%%%%%%%%%%%%%% a1=[sqrt(3)/2,1.5,0].*a;%正空间基矢 a2=[-sqrt(3)/2,1.5,0].*a;%三维%a1=[a,0,0];a2=[0,a,0];a3=[0,0,1];%单位向量A=abs(a3*cross(a1,a2)');s=2*pi./A;b1(1,1)=a2(1,2);b1(1,2)=-a

2021-08-02 16:37:38 1227 3

量子力学贝里相位笔记

量子力学贝里相位笔记

2020-11-19

智能优化算法—免疫算法笔记.md

很详细的免疫算法实现步骤,用于智能优化,与原始版本免疫算法不同,但不是代码,要的可以下载,智能优化算法可以实现无导数优化,实现问题的最优化,但上传的算法是单目标优化算法,要的就下载吧。

2020-09-14

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除