【无标题】

clear all;
close all;

%%
%data1=xlsread(‘E:\data\700-600\100KHZ\16\16-1.xlsx’,1,‘A2:I801’)/J2:R801

data1=xlsread(‘E:\data\600-600\30-1.xlsx’,1,‘A2:F9000’);
data2=xlsread(‘E:\data\600-600\30-1.xlsx’,2,‘A2:F9000’);
data3=xlsread(‘E:\data\600-600\30-1.xlsx’,3,‘A2:F9000’);
data4=xlsread(‘E:\data\600-600\30-1.xlsx’,4,‘A2:F9000’);
data5=xlsread(‘E:\data\600-600\30-1.xlsx’,5,‘A2:F9000’);
data6=xlsread(‘E:\data\600-600\30-1.xlsx’,6,‘A2:F9000’);
hdata1=h(data1);
hdata2=h(data2);
hdata3=h(data3);
hdata4=h(data4);
hdata5=h(data5);
hdata6=h(data6);

plot(data1(:,3));
%%
t0=data1(5,1)-data1(4,1);
vs0=5324;
A=zeros(300,300);

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

d(1).d=@(x,y) sqrt((x-00)^2+(y-0)2);
d(2).d=@(x,y) sqrt((x-150)^2+(y-0)2);
d(3).d=@(x,y) sqrt((x-300)^2+(y-00)2);
d(4).d=@(x,y) sqrt((x-000)^2+(y-300)2);
d(5).d=@(x,y) sqrt((x-150)^2+(y-300)2);
d(6).d=@(x,y) sqrt((x-300)^2+(y-300)2);
%%
%r c行列

for r=1:1:300
for c=1:1:300
for i=[1,2,3,4,5,6]
for j=1:6
if i~=j
distance=d(i).d(r,c)+d(j).d(r,c);
row=round(distance1e-3/(vs0t0));
if 2000<row && row<6300
switch(i)
case 1
A(r,c)=A(r,c)+hdata1(row+100,j+1);
case 2
A(r,c)=A(r,c)+hdata2(row+100,j+1);
case 3
A(r,c)=A(r,c)+hdata3(row+100,j+1);
case 4
A(r,c)=A(r,c)+hdata4(row+100,j+1);
case 5
A(r,c)=A(r,c)+hdata5(row+100,j+1);
case 6
A(r,c)=A(r,c)+hdata6(row+100,j+1);
end
end
end
end
end
end
end

% 画图
figure(2);

% clims = [14.3e-9 14.5e-9];
% imagesc(A,clims)
imagesc(A)
colormap(‘jet’)
colorbar