幅值全聚焦成像

clear all;
close all;
%%    
t0=2.50000000000000e-07;
vs0=2538;
% vs0=5336;
A=zeros(200,200);

h = @(x) abs(hilbert(x));

d(1).d=@(x,y) sqrt((x-100)^2+(y-0)^2);
d(2).d=@(x,y) sqrt((x-200)^2+(y-0)^2);
d(3).d=@(x,y) sqrt((x-200)^2+(y-100)^2);
d(4).d=@(x,y) sqrt((x-200)^2+(y-200)^2);
d(5).d=@(x,y) sqrt((x-100)^2+(y-200)^2);
d(6).d=@(x,y) sqrt((x-00)^2+(y-200)^2);
d(7).d=@(x,y) sqrt((x-00)^2+(y-100)^2);
d(8).d=@(x,y) sqrt((x-00)^2+(y-0)^2);
%%
%data1=xlsread('E:\data\700-600\100KHZ\16\16-1.xlsx',1,'A2:I801')/J2:R801
data16_1=xlsread('E:\data\700-600\100KHZ\16\16-1.xlsx',1,'J2:R801');
data0_1=xlsread('E:\data\700-600\100KHZ\0\0.xlsx',1,'J2:R801');
hd_data16_1=h(data16_1-data0_1);
data16_3=xlsread('E:\data\700-600\100KHZ\16\16-1.xlsx',2,'J2:R801');
data0_3=xlsread('E:\data