fid=fopen('xishu20m.dat','r');
data=fscanf(fid,'%f');
line=size(data,1)/101;
data=(reshape(data,101,line));
C=20*log10(data);

R=(9500:10:10500)/1000;
Z=15:0.25:25;
figure(1)
pcolor(R,Z,transpose(C))
colormap(gray);
xlabel('距聲源水平距離(km)');ylabel('海深(m)');title('信道歸一化互相關衰落');
shading interp;
colorbar


% figure(2)
% R=(9500:10:10500)/1000;
% Z=45:0.25:55
% surfc(R,Z,transpose(C))
%  xlabel('距聲源水平距離(km)');ylabel('海深(m)');zlabel('信道歸一化互相關衰落');
% colorbar
%