clc
clear?all
close?all
jay=sqrt（-1）;
M=10;
N=10;
f0=200*1000;
fs=20*f0;
T=0.05;
t=0:1/fs:T;
c=1500;
lamda=c/f0;
d=lamda/2;
st=cos（2*pi*f0*t）;
sst=st*st‘/length（st）;
%%?方陣
k=1
p?=?zeros（3M*N）;
for?iii?=?1?:?M
????for?jjj?=?1:?N
????????p（:k）?=?[-（M-1）/2*d+（iii-1）*d-（N-1）/2*d+（jjj-1）*d0]‘;
????????k?=?k+1;
????end
end??????????%陣元位置
figure（1）;
plot3（p（1:）p（2:）p（3:）‘ko‘）;
hold?on;
xlabel（‘\it?x‘）;
ylabel（‘\it?y‘）;
zlabel（‘\it?z‘）;

%?%%?圓陣
%?p=zeros（3M）;
%?R=20;
%?k=1;
%?for?kk=1:M
%?????p（:k）=[R*cos（kk*2*pi/M）R*sin（kk*2*pi/M）0]‘
%?????k=k+1;
%?end
%?figure（1）
%?plot3（p（1:）p（2:）p（3:）‘ko‘）;
%?hold?on;
%?xlabel（‘\it?x‘）;
%?ylabel（‘\it?y‘）;
%?zlabel（‘\it?z‘）;

%設置成像區域
I=100;
J=100;
di=0.1;
Rf=3;
A=zeros（3I*J）;
l=1;
for?i?=?1?:?I
????for?j?=?1:J
????????A（:l）?=?[（-（I-1）/2）*di+（i-1）*di（-（J-1）/2）*di+（j-1）*di-Rf]‘;
????????l=?l+1;
????end
end?????????%成像單元位置
figure（2）;
plot3（A（1:）