plotshade.m

水声模型 很不错的东西

function NewDir = PlotShade(DefaultDir)
% PlotShade()
% plots the shade file
%Original by Mike Porter,
%Modded by Alec Duncan
%Modded by Alec Duncan to plot binary shade files 13/6/02



if nargin < 1
   ChangeDir = 0;
else
   StartDir = pwd;
   DirStat = 1;
   if ~exist(DefaultDir, 'dir')
       DirStat = LibMakeDirectory(DefaultDir);
   end
   if DirStat
       cd(DefaultDir);
       ChangeDir = 1;
   else
       ChangeDir = 0;
   end
end

[FName, Path] = uigetfile('*.shd', 'Select file to plot');

if ChangeDir
   cd(StartDir);
end

Status = 1;
if FName == 0
    if nargin < 1
        NewDir = pwd;
    else
        NewDir = DefaultDir;
    end
    Status = 0;
end

if Status
   NewDir = Path;
   filename = [Path FName];

	% read
	[ pltitl, freq, nsd, nrd, nrr, sd, rd, rr, tlt, Status ] = ReadShadeBin( filename );
end

if Status
	isd = 1;  %index of source to be plotted
   if nrd <= 1
      msgbox(...
         'Propagation code must be run with multiple receiver depths in order for field to be plotted', ...
      'PlotShade:', 'warn');
      return;
   end
   
	zt = rd;
	taker = 1:nrr;
	rt = rr( taker );


    Ans = inputdlg({'Figure number', 'Plot dynamic range (dB)', ...
            'Remove cylindrical spreading (y/n)'}, '', 1, {int2str(gcf), '60', 'n'});
    if isempty(Ans)
       return;
    end
   
    FigNum = str2num(Ans{1});
    DynRange = str2num(Ans{2});
    RemoveCylindrical = strcmpi(Ans{3}, 'y');
   
    figure(FigNum);
    clf;
    if exist('GUI_SetupGraphMenu.m', 'file')
        GUI_SetupGraphMenu('pcolor');
    end
    
    TL = -20*log10(abs(tlt));
    if RemoveCylindrical
	    TL = TL - ones( nrd, 1 ) * 10.0 * log10( rt )';
    end
    
    PltTLMin = floor(min(min(TL))/5)*5;  %Minimum transmission loss for plot - round to nearest 5 dB
    PltTLMax = PltTLMin + DynRange;

	pcolor( rt, zt, TL ); 
    shading interp;
    caxis([PltTLMin PltTLMax]);
    colormap('default');
	Map = colormap;
	colormap(flipud(Map));  %Want small TL = big signal = red

	h = colorbar;
    h2 = get(h, 'ylabel');
    if RemoveCylindrical
        set(h2, 'string', 'Transmission Loss - 10logR (dB)');
    else
        set(h2, 'string', 'Transmission Loss (dB)');
    end
	view( 0, -90 );
	xlabel( 'Range (m)' ); ylabel( 'Depth (m)' );
    title( [deblank( pltitl ) '. f=' num2str(freq) 'Hz, Zs=' num2str(sd(isd)) 'm.'] , ...
      'interpreter', 'none');  
end