toeplitzindextest.m

通过计算晶体结构来确定其光谱分布

%ToeplitzIndexTest.m
clear;
close all;
x1=2;
y1=x1;
x=(-x1:x1);
y=(-y1:y1);
t=size(x);
Nx=t(2);
t=size(y);
Ny=t(2);
Index=zeros(Ny,Nx,2);
ir=0;
for iy=1:Ny
    for ix=1:Nx
        ir=ir+1;
        il=0;
        for iy1=1:Ny
            for ix1=1:Nx
                il=il+1;
                Index(ir,il,1)=iy-iy1;
                Index(ir,il,2)=ix-ix1;
            end
        end
    end
end
hN=0;
hN=hN+1;
figure(hN);
surfl(Index(:,:,1));
title('Index 1');
hN=hN+1;
figure(hN);
surfl(Index(:,:,2));
title('Index 2');
hN=hN+1;
figure(hN);
image(Index(:,:,1));
title('Index 1');
hN=hN+1;
figure(hN);
image(Index(:,:,2));
title('Index 2');