In [2]:

    

run_num_all = [79,87,94,103,115,124,131];
ping_num{1} = 4:9:877;    % run 79
ping_num{2} = 2:2:1000;   % run 87
ping_num{3} = 3:3:899;    % run 94
ping_num{4} = 3:5:911;    % run 103
ping_num{5} = 2:2:926;    % run 115
ping_num{6} = 2:2:973;    % run 124
ping_num{7} = 2:2:1000;   % run 131


    

In [3]:

    

run_num_all


    

    




run_num_all =

    79    87    94   103   115   124   131

In [4]:

    

ping_num


    

    




ping_num =

  1×7 cell array

  Columns 1 through 4

    [1×98 double]    [1×500 double]    [1×299 double]    [1×182 double]

  Columns 5 through 7

    [1×463 double]    [1×486 double]    [1×500 double]

In [5]:

    

iR=1;
A(iR) = summarize_run_fcn(run_num_all(iR),ping_num{iR},1);


    

In [6]:

    

A


    

    




A = 

  struct with fields:

    energy_in_bnd: [98×1 double]
            si_W1: [98×1 double]
            si_W2: [98×1 double]
           max_W1: [98×1 double]
           max_W2: [98×1 double]
        ping_time: [98×1 double]

In [7]:

    

iR=2;
A(iR) = summarize_run_fcn(run_num_all(iR),ping_num{iR},1);


    

In [16]:

    

clear A


    

In [17]:

    

A(1) = summarize_run_fcn(run_num_all(1),ping_num{1},0);
for iR=2:length(run_num_all)
    A(iR) = summarize_run_fcn(run_num_all(iR),ping_num{iR},0);
end


    

In [18]:

    

A


    

    




A = 

  1×7 struct array with fields:

    energy_in_bnd
    si_W1
    si_W2
    max_W1
    max_W2
    ping_time

In [29]:

    

figure
for iR=1:length(run_num_all)
    subplot(311)
    plot(A(iR).ping_time,10*log10(A(iR).max_W1)+5*(iR-1))
    hold on
    ylim([60 120])
    subplot(312)
    plot(A(iR).ping_time,10*log10(A(iR).max_W2)+5*(iR-1))
    hold on
    ylim([60 120])
    subplot(313)
    plot(A(iR).ping_time,10*log10(A(iR).energy_in_bnd)+5*(iR-1))
    hold on
    ylim([160 220])
end


    

In [39]:

    

figure
for iR=2:length(run_num_all)
    subplot(311)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).max_W1),1))
    hold on
    ylim([70 85]); grid on;
    subplot(312)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).max_W2),1))
    hold on
    ylim([60 85]); grid on;
    subplot(313)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).energy_in_bnd),1))
    hold on
    ylim([165 195]); grid on;
end


    

In [42]:

    

figure
for iR=2:length(run_num_all)
    subplot(311)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).max_W1),5))
    hold on
    ylim([70 85]); grid on;
    xlim([18 32])
    subplot(312)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).max_W2),5))
    hold on
    ylim([60 85]); grid on;
    xlim([18 32])
    subplot(313)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).energy_in_bnd),5))
    hold on
    ylim([165 195]); grid on;
    xlim([18 32])
end


    

In [55]:

    

figure
for iR=[2,3,5,6,7]
    subplot(311)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).max_W1),5)-max(smooth(10*log10(A(iR).max_W1),5)))
    hold on
    ylim([-10 0]); grid on;
    xlim([18 32])
    legend(num2str(run_num_all([2,3,5,6,7])'),'location','eastoutside')
    subplot(312)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).max_W2),5)-max(smooth(10*log10(A(iR).max_W2),5)))
    hold on
    ylim([-10 0]); grid on;
    xlim([18 32])
    legend(num2str(run_num_all([2,3,5,6,7])'),'location','eastoutside')
    subplot(313)
    plot(A(iR).ping_time,smooth(10*log10(A(iR).energy_in_bnd),5)-max(smooth(10*log10(A(iR).energy_in_bnd),5)))
    hold on
    ylim([-15 0]); grid on;
    xlim([18 32])
    legend(num2str(run_num_all([2,3,5,6,7])'),'location','eastoutside')
end
subplot(311)
title('Max echo level in W1')
subplot(312)
title('Max echo level in W2')
subplot(313)
title('Total energy')


    

    




Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)
Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)
Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)
Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)
Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)
Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)
Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)
Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)
Warning: Ignoring extra legend entries.
> In legend>set_children_and_strings (line 658)
  In legend>make_legend (line 334)
  In legend (line 282)






    

In [ ]: