jack_BOB_LUO的博客

matlab中高精度计算函数vpa

【代码】对colormap的那个N行三列矩阵进行线性插值的代码(三维scatter3中的确可以看到是正确插好值了的)

【代码】紧束缚开链驻波动画。

给定正多边形外切圆半径和边数，给出正多边形各顶点的xy坐标的matlab函数

matlab中一种构建矩阵的技巧

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

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

二维平面波展开法

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

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

只适用于特定关系的数组

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

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

主要用于识别简并能带

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

如题

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

如题所示

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

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

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.

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表示内部填

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

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

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

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

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

clearvars -except stoM

%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.

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

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

function []=main_surper_cell()n=10;a = 1; %%%%%%%%%%%%%%%定义原胞基矢%%%%%%%%%%%%%% a1 = a* [ sqrt(3)/2 1/2 ]; a2 = a* [-sqrt(3)/2 1/2 ]; a3 = a* [ 0 0 1]; %%%%%%%%%%%%%%%定义三个格点的位置%%%%%%%%%%%%%% u1 = a*[sqrt(3)/6 0 ]; u2 = a*[5*sqrt(3)/12 a/4 ]; u3 =

function [k0]=BZ_dispersed()global aa=1;%临时参数delta=0.01;%可调参数T=[0,0];M=(2*pi/a)*[1/(4+2*sqrt(2)),1/(4+2*sqrt(2))];X=(2*pi/a)*[1/(4+2*sqrt(2)),0];MT=T-M;mt=MT(1)+MT(2)*1i;list_mt=(M(1)+M(2)*1i)+[0:delta:1]*mt;TX=X-T;tx=TX(1)+TX(2)*1i;list_tx

*

function []=main()a=5;%晶格常数d1=a*[1,0];%正空间基矢d2=a*[0.5,-sqrt(3)/2];para=20;sto=[];for u=0:para for v=0:para D=u*d1+v*d2; [p]=triangle_fun(D,0.5*a); sto=[sto;p]; endendscatter(sto(:,1),sto(:,2),'o')end%下面的子程序可自己编

dot函数在对复数向量操作的时候要注意的地方