%初始化三角模式中的三个点
rng(0);

delta = 5*sqrt(3);
delta = 1;
%随机选择一个方向探索
phi_0 = (rand()*2-1)*pi;
delta_phi_0 = [2*pi/3 -2*pi/3];
n = 300;
x = zeros(n,2);
xinit = [0 30];
x(1,:) = delta*[cos(phi_0) sin(phi_0)]+xinit;
a = phi_0 + delta_phi_0;
x(2:3,:) = delta*[cos(a)' sin(a)']+xinit;
rpos = [-0.5,0.5];
env = WaveEnv();
yfun = @(angles,k)env.getIntensityByYP(rpos(1),rpos(2),angles(1),angles(2),0/100);

yaw = atand(rpos(2)/rpos(1));
yaw = 135;
ap = 0:10:180;
ints = zeros(length(ap),1);
for i = 1:length(ap)
    ints(i) = yfun([yaw,ap(i)],0);
end
plot(ap,ints)