- 博客(49)
- 资源 (2)
- 收藏
- 关注
原创 给定正多边形外切圆半径和边数,给出正多边形各顶点的xy坐标的matlab函数
给定正多边形外切圆半径和边数,给出正多边形各顶点的xy坐标的matlab函数
2024-06-15 20:45:31 128
原创 群折射率$ n_g$与有效折射率$ n_{eff}$的关系详细推导(假设$ \lambda_0$为真空入射光波长)
群折射率和有效折射率关系公式的详细推导
2023-08-15 15:01:50 2235
原创 二维体光子晶体的平面波展开法代码2
这篇文章接前面《二维体光子晶体的平面波展开法代码》这个是实际网格划分,各标记的含义上篇文章有说明下面是傅里叶系数恢复的相对介电常数分布(虚部),看起来非常大,那是因为imag(EpsMat)趋于0的原因。下面这个图是实空间相对介电常数分布的伪彩图:下面这个图是TE能带图,可以看到和FEM的结果还是有点偏差,但实际上只要加密网格同时遵守前一篇文章所述的nG参数和precisStruct参数的关系就可以获得误差更小的能带图。
2023-05-10 15:35:15 322
原创 matlab实现DFT(取了一些数据计算结果和fft相差很小,但未详细验证)
DFT的matlab实现,和fft进行了少量对比误差非常小,但没详细验证,慎重。
2023-05-06 21:56:02 147
原创 (PWEM前置知识)解析公式求相对介电常数分布函数的 Fourier 展开系数以及基于此的绘制实空间原胞的相对介电常数分布
(PWEM前置知识)解析公式求相对介电常数分布函数的 Fourier 展开系数以及基于此的绘制实空间原胞的相对介电常数分布
2023-05-03 15:21:32 434
原创 Matlab绘图为什么设置了相同的字体大小结果出的图字体大小还是不一样
matlab相同字号设置但最终图片字号大小不一致(matlab2018b中已验证过)
2023-02-26 10:06:50 1027 1
原创 检查一维数组中的重复元素以及每个重复元素的重复次数(结果以矩阵形式给出)
matlab函数,检测一列数组中的重复元素和重复次数,结果以矩阵形式给出,第一行是重复元素,第二行是对应元素的重复次数
2022-07-27 21:38:42 561
原创 matlab复误差函数的高效计算(有官方背书)
matlab复误差函数的高效计算,可以和mathematica的官方函数对比验证下,有权威官方背书,是参考网站的代码改写而来。
2022-06-27 17:40:46 517
原创 赫尔维兹_勒奇超越函数(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 877
原创 带选择程序(未测试)
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 210
原创 三角格子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 264
原创 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 6195
原创 matlab angle函数
angle采用的算法:angle(z) = imag(log(z)) = atan2(imag(z), real(z)).
2021-12-29 11:19:01 2310
原创 一个基本成功的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 645
原创 一段判断点是否在多边形内的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 541
原创 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 828
原创 画格点程序
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 135
原创 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 2125
原创 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 273
原创 一个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 1852
智能优化算法—免疫算法笔记.md
2020-09-14
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人