Untitled

clc;clf;clear all;close all;
beam =3;
if(beam ==1)
    [A,B,C]=xlsread('Moments.xlsx',1);
    ht= A(1,(1:44))';
    dates=datetime(strrep(string(((B(4:613))')),'''',''));
    fundasnrs=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',2);
    myextractedSNR=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
    
    [A,B,C]=xlsread('Moments.xlsx',3);
    DopplerWidths=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',4);
    myextractedDWidth=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
    
    [A,B,C]=xlsread('Moments.xlsx',5);
    Powers=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',6);
    myextractedPWR=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
    
    [A,B,C]=xlsread('Moments.xlsx',7);
    NoiseLevels=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',8);
    myextractedNoise=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
elseif(beam ==2)
    [A,B,C]=xlsread('Moments.xlsx',9);
    ht= A(1,(1:44))';
    dates=datetime(strrep(string(((B(4:613))')),'''',''));
    fundasnrs=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',10);
    myextractedSNR=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
   
    
    [A,B,C]=xlsread('Moments.xlsx',11);
    DopplerWidths=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',12);
    myextractedDWidth=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
    
    [A,B,C]=xlsread('Moments.xlsx',13);
    Powers=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',14);
    myextractedPWR=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
    
    [A,B,C]=xlsread('Moments.xlsx',15);
    NoiseLevels=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',16);
    myextractedNoise=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
    
elseif(beam ==3)
    [A,B,C]=xlsread('Moments.xlsx',17);
    ht= A(1,(1:44))';
    dates=datetime(strrep(string(((B(4:613))')),'''',''));
    fundasnrs=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',18);
    myextractedSNR=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
       
    [A,B,C]=xlsread('Moments.xlsx',19);
    DopplerWidths=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',20);
    myextractedDWidth=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
    
    [A,B,C]=xlsread('Moments.xlsx',21);
    Powers=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',22);
    myextractedPWR=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
    
    [A,B,C]=xlsread('Moments.xlsx',23);
    NoiseLevels=A((2:611),(1:44))';
    
    [A,B,C]=xlsread('Moments.xlsx',24);
    myextractedNoise=double(strrep(string(((C((4:613),(2:45)))')),'''',''));
end
xlsp = (linspace(dates(1),dates(numel(dates)),4));
ylsp = (linspace(ht(1),ht(numel(ht)),4));
f12=figure(1); set(f12, 'PaperType', 'a4'); set(f12, 'PaperOrientation', 'portrait');set(f12, 'Resize', 'on');set(f12, 'Position',[0 0 (2480*(95/100)) (3508*(80/100))]);
subplot(5,2,3);
pcolor(dates,ht,fundasnrs);
set (gca, 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');
set(gca, 'clim', [-10 25]);
set(gca, 'XTick', xlsp);
set(gca, 'YTick', ylsp);
set(gca,'TickDir','out'); % The only other option is 'in'
set(gca, 'XTickLabels', [{'05:50' '09:43' '13:35' '17:28'}]);
ytickformat('%.1f')
if(beam ==1)
    title('(a) Fundamental Signal to Noise Ratio (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
elseif(beam ==2)
    title('(a) Fundamental Signal to Noise Ratio (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
elseif(beam ==3)
    title('(a) Fundamental Signal to Noise Ratio (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
end
y = ylabel('Range (km)', 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal'); y.FontWeight='Normal';
shading interp;
colormap('jet');
cBarTicks = round(linspace(-10, 25,5));
colorbar('Ticks',cBarTicks,"TickDirection","out",'FontSize',12,'FontName','Times New Roman', 'FontWeight', 'Normal');box off;
set(gcf,'renderer','opengl') ;
subplot(5,2,5);
pcolor(dates,ht,DopplerWidths);
set (gca, 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');
set(gca, 'clim', [0 1.5]);
set(gca, 'XTick', xlsp);
set(gca, 'YTick', ylsp);
set(gca,'TickDir','out'); % The only other option is 'in'
set(gca, 'XTickLabels', [{'05:50' '09:43' '13:35' '17:28'}]);
ytickformat('%.1f')
if(beam ==1)
    title('(b) Doppler Width (Hz)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
elseif(beam ==2)
    title('(b) Doppler Width (Hz)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
elseif(beam ==3)
    title('(b) Doppler Width (Hz)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
end
y = ylabel('Range (km)', 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal'); y.FontWeight='Normal';
shading interp;
colormap('jet');
h=colorbar;
t=get(h,'Limits');
T=linspace(t(1),t(2),5)
T = 1×5
         0    0.3750    0.7500    1.1250    1.5000
set(h,'Ticks',T);
TL=arrayfun(@(x) sprintf('%.2f',x),T,'un',0);
set(h,'TickLabels',TL)
%colorbar;
box off;
subplot(5,2,7);
pcolor(dates,ht,Powers);
set (gca, 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');
set(gca, 'clim', [75 130]);
set(gca, 'XTick', xlsp);
set(gca, 'YTick', ylsp);
set(gca,'TickDir','out'); % The only other option is 'in'
set(gca, 'XTickLabels', [{'05:50' '09:43' '13:35' '17:28'}]);
ytickformat('%.1f')
if(beam ==1)
    title('(c) Signal Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
elseif(beam ==2)
    title('(c) Signal Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
elseif(beam ==3)
    title('(c) Signal Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
end
y = ylabel('Range (km)', 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal'); y.FontWeight='Normal';
shading interp;
colormap('jet');
cBarTicks = round(linspace(75, 130,5));
colorbar('Ticks',cBarTicks,"TickDirection","out",'FontSize',12,'FontName','Times New Roman', 'FontWeight', 'Normal');box off;
set(gcf,'renderer','opengl')
subplot(5,2,9);
pcolor(dates,ht,NoiseLevels);
set (gca, 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');
set(gca, 'clim', [-25 12]);
set(gca, 'XTick', xlsp);
set(gca, 'YTick', ylsp);
set(gca,'TickDir','out'); % The only other option is 'in'
xtickformat('HH:mm');
ytickformat('%.1f')
if(beam ==1)
    title('(d) Noise Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
elseif(beam ==2)
    title('(d) Noise Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
elseif(beam ==3)
    title('(d) Noise Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');
end
y = ylabel('Range (km)', 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal'); y.FontWeight='Normal';
shading interp;
colormap('jet');
cBarTicks = round(linspace(-25, 12,5));
colorbar('Ticks',cBarTicks,"TickDirection","out",'FontSize',12,'FontName','Times New Roman', 'FontWeight', 'Normal');box off;

    

subplot(5,2,4);

pcolor(dates,ht,myextractedSNR);

set (gca, 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');

set(gca, 'clim', [-60 60]);

set(gca, 'XTick', xlsp);

set(gca, 'YTick', ylsp);

set(gca,'TickDir','out'); % The only other option is 'in'

set(gca, 'XTickLabels', [{'05:50' '09:43' '13:35' '17:28'}]);

ytickformat('%.1f')

if(beam ==1)

title('(e) Fundamental Signal to Noise Ratio (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

elseif(beam ==2)

title('(e) Fundamental Signal to Noise Ratio (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

elseif(beam ==3)

title('(e) Fundamental Signal to Noise Ratio (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

end

y = ylabel('Range (km)', 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal'); y.FontWeight='Normal';

shading interp;

colormap('jet');

cBarTicks = round(linspace(-60, 60,5));

colorbar('Ticks',cBarTicks,"TickDirection","out",'FontSize',12,'FontName','Times New Roman', 'FontWeight', 'Normal');

box on;

subplot(5,2,6);

pcolor(dates,ht,myextractedDWidth);

set (gca, 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');

set(gca, 'clim', [0 4]);

set(gca, 'XTick', xlsp);

set(gca, 'YTick', ylsp);

set(gca,'TickDir','out'); % The only other option is 'in'

set(gca, 'XTickLabels', [{'05:50' '09:43' '13:35' '17:28'}]);

ytickformat('%.1f')

if(beam ==1)

title('(f) Doppler Width (Hz)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

elseif(beam ==2)

title('(f) Doppler Width (Hz)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

elseif(beam ==3)

title('(f) Doppler Width (Hz)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

end

y = ylabel('Range (km)', 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal'); y.FontWeight='Normal';

shading interp;

colormap('jet');

cBarTicks = round(linspace(0, 4,5));

colorbar('Ticks',cBarTicks,"TickDirection","out",'FontSize',12,'FontName','Times New Roman', 'FontWeight', 'Normal');

box on;

subplot(5,2,8);

pcolor(dates,ht,myextractedPWR);

set(gca, 'clim', [60 160]);

set (gca, 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');

set(gca, 'XTick', xlsp);

set(gca, 'YTick', ylsp);

set(gca,'TickDir','out'); % The only other option is 'in'

set(gca, 'XTickLabels', [{'05:50' '09:43' '13:35' '17:28'}]);

ytickformat('%.1f')

if(beam ==1)

title('(g) Signal Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

elseif(beam ==2)

title('(g) Signal Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

elseif(beam ==3)

title('(g) Signal Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

end

y = ylabel('Range (km)', 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal'); y.FontWeight='Normal';

shading interp;

colormap('jet');

cBarTicks = round(linspace(60, 160,5));

colorbar('Ticks',cBarTicks,"TickDirection","out",'FontSize',12,'FontName','Times New Roman', 'FontWeight', 'Normal');

box on;

subplot(5,2,10);

pcolor(dates,ht,myextractedNoise);

set (gca, 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');

set(gca, 'clim', [40 140]);

set(gca, 'XTick', xlsp);

set(gca, 'YTick', ylsp);

set(gca,'TickDir','out'); % The only other option is 'in'

xtickformat('HH:mm');

ytickformat('%.1f')

if(beam ==1)

title('(h) Noise Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

elseif(beam ==2)

title('(h) Noise Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

elseif(beam ==3)

title('(h) Noise Power (dB)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'bold');

end

% xlabel('Date Time');

y = ylabel('Range (km)', 'FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal'); y.FontWeight='Normal';

str='Comparision of detailed Moments derived using eBayes and FFT';

T = text(min(get(gca, 'xlim')), max(get(gca, 'ylim')), str);

set(T,'position',[-0.201975474640909,242.7541130171058,0],'FontSize', 14,"VerticalAlignment","top","HorizontalAlignment","left",'FontName','Times New Roman', 'FontWeight', 'bold','color','#000000');

y = ylabel('Range (km)','FontSize', 14, 'FontName','Times New Roman', 'FontWeight', 'Normal');%y.FontWeight='bold';

shading interp;

colormap('jet');

cBarTicks = round(linspace(40, 140,5));

colorbar('Ticks',cBarTicks,"TickDirection","out",'FontSize',12,'FontName','Times New Roman', 'FontWeight', 'Normal');

box on;