matlab - MATLAB 中的光学晶格

Question

我正在编写脚本来绘制以下图片

以下代码可以正常工作并绘制与上面相同的形状，但没有球体。

clear all
PS=zeros(100,100); 
A=2.4; 
B=3.4; 

for i=1:100 

  for j=1:100        
   PS(i,j) = cos((.1*i)*A)*cos((.1*j)*B); 
  end 
end 
surfc(PS)

我的问题是，如何绘制这些球体？

Answer 1

一旦您在晶格上获得了它们的坐标，就可以很容易地绘制原子/球体，也就是说，您需要找到（如维基百科文章中所述）潜在的最小值。我想你可以弄清楚那部分。如果不是我误解了你的问题对不起！

对于演示，我将使用较小的网格以减少原子并用于scatter3绘制球体。然后你可以用精确的最小值替换我的近似值。

这是整个代码：

clear
clc
close all

N = 30;
PS=zeros(N,N); 
A=2.4; 
B=3.4; 

for i=1:N

    for j=1:N 

        PS(i,j) = cos((.1*i)*A)*cos((.1*j)*B); 
    end

end

%// Approximate locations of the atoms. You can calculate this more
%// accurately of course.
xAtom = [1 10 10 19 28 28];
yAtom = [13 2 27 13 2 27];

%// Plot atoms at height of .6.
zAtom = repmat(.6,1,numel(xAtom));

surfc(PS)

hold on

%// Scatter plot. You can customize the parameters.
scatter3(xAtom,yAtom,zAtom,200,'k','filled')

rotate3d on

和输出：