jack_BOB_LUO的博客

matlab中高精度计算函数vpa

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

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

如题

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

如题所示

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

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

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

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')

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

矩形介质柱的kkk空间相对介电常数计算：α(G)={αAf+αB(1−f),G=0(αA−αB)P(G),G≠0\alpha(G)=\left\{\begin{array}{l}\alpha_{\mathrm{A}} f+\alpha_{\mathrm{B}}(1-f), G=0 \\\left(\alpha_{\mathrm{A}}-\alpha_{\mathrm{B}}\right) P(G), \quad G \neq 0\end{array}\right.α(G)={αA​f+αB​(1

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%下面的子程序可自己编

把平面波展开法当做一个响应器，即有输入和输出，先谈它的输入，输入有：1.介质柱的相对介电常数2.背景的相对介电常数3.填充率（介质柱横截面积比上正格子中二维原胞面积，三角格子光子晶体的原胞是一个菱形，正方格子光子晶体的原胞是一个矩形）4.圆形介质柱的半径R(但是介质柱可能不是圆形甚至是不规则的形状，这个问题留在后面讨论)，这个半径一般表示为晶格常数a的倍数.5.还有一个参数和计算过程中运用的平面波的数目有关，我们假定它是n，平面波展开法中平面波的数目越多计算结果就会越精确，但我不确定这是否是绝对的