聚焦多波束成像matlab

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:）