clc; clear all;clf;close all;
EQualMomWithRTI=figure(1); set(EQualMomWithRTI, 'PaperType', 'a4'); set(EQualMomWithRTI, 'PaperOrientation', 'portrait');set(EQualMomWithRTI, 'Resize', 'on');set(EQualMomWithRTI, 'Position',[0 0 2480*(95/100) 3508*(95/100)]);
[A,B,C]=xlsread('MSTSample.xlsx',1);
ht1= A(1,(1:44))';
nrgs = numel(ht1);
beam=2;
if(beam==1)
[A,B,C]=xlsread('MSTSample.xlsx',1);
ht= A(1,(1:44))';
SNR_5 =A((2:10),(1:44));SNR_6 =A((11:63),(1:44));SNR_7 =A((64:115),(1:44));SNR_8 =A((116:167),(1:44));SNR_9 =A((168:220),(1:44));SNR_10 =A((221:272),(1:44));SNR_11 =A((273:324),(1:44));SNR_12 =A((325:376),(1:44));SNR_13 =A((377:429),(1:44));SNR_14 =A((430:481),(1:44));SNR_15 =A((482:533),(1:44));SNR_16 =A((534:586),(1:44));SNR_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',2);
tSNR_5 =A((2:10),(1:44));tSNR_6 =A((11:63),(1:44));tSNR_7 =A((64:115),(1:44));tSNR_8 =A((116:167),(1:44));tSNR_9 =A((168:220),(1:44));tSNR_10 =A((221:272),(1:44));tSNR_11 =A((273:324),(1:44));tSNR_12 =A((325:376),(1:44));tSNR_13 =A((377:429),(1:44));tSNR_14 =A((430:481),(1:44));tSNR_15 =A((482:533),(1:44));tSNR_16 =A((534:586),(1:44));tSNR_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',3);
BandPower_5 =A((2:10),(1:44));BandPower_6 =A((11:63),(1:44));BandPower_7 =A((64:115),(1:44));BandPower_8 =A((116:167),(1:44));BandPower_9 =A((168:220),(1:44));BandPower_10 =A((221:272),(1:44));BandPower_11 =A((273:324),(1:44));BandPower_12 =A((325:376),(1:44));BandPower_13 =A((377:429),(1:44));BandPower_14 =A((430:481),(1:44));BandPower_15 =A((482:533),(1:44));BandPower_16 =A((534:586),(1:44));BandPower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',4);
tBandPower_5 =A((2:10),(1:44));tBandPower_6 =A((11:63),(1:44));tBandPower_7 =A((64:115),(1:44));tBandPower_8 =A((116:167),(1:44));tBandPower_9 =A((168:220),(1:44));tBandPower_10 =A((221:272),(1:44));tBandPower_11 =A((273:324),(1:44));tBandPower_12 =A((325:376),(1:44));tBandPower_13 =A((377:429),(1:44));tBandPower_14 =A((430:481),(1:44));tBandPower_15 =A((482:533),(1:44));tBandPower_16 =A((534:586),(1:44));tBandPower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',5);
NoisePower_5 =A((2:10),(1:44));NoisePower_6 =A((11:63),(1:44));NoisePower_7 =A((64:115),(1:44));NoisePower_8 =A((116:167),(1:44));NoisePower_9 =A((168:220),(1:44));NoisePower_10 =A((221:272),(1:44));NoisePower_11 =A((273:324),(1:44));NoisePower_12 =A((325:376),(1:44));NoisePower_13 =A((377:429),(1:44));NoisePower_14 =A((430:481),(1:44));NoisePower_15 =A((482:533),(1:44));NoisePower_16 =A((534:586),(1:44));NoisePower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',6);
tNoisePower_5 =A((2:10),(1:44));tNoisePower_6 =A((11:63),(1:44));tNoisePower_7 =A((64:115),(1:44));tNoisePower_8 =A((116:167),(1:44));tNoisePower_9 =A((168:220),(1:44));tNoisePower_10 =A((221:272),(1:44));tNoisePower_11 =A((273:324),(1:44));tNoisePower_12 =A((325:376),(1:44));tNoisePower_13 =A((377:429),(1:44));tNoisePower_14 =A((430:481),(1:44));tNoisePower_15 =A((482:533),(1:44));tNoisePower_16 =A((534:586),(1:44));tNoisePower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',7);
SinadRatio_5 =A((2:10),(1:44));SinadRatio_6 =A((11:63),(1:44));SinadRatio_7 =A((64:115),(1:44));SinadRatio_8 =A((116:167),(1:44));SinadRatio_9 =A((168:220),(1:44));SinadRatio_10 =A((221:272),(1:44));SinadRatio_11 =A((273:324),(1:44));SinadRatio_12 =A((325:376),(1:44));SinadRatio_13 =A((377:429),(1:44));SinadRatio_14 =A((430:481),(1:44));SinadRatio_15 =A((482:533),(1:44));SinadRatio_16 =A((534:586),(1:44));SinadRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',8);
tSinadRatio_5 =A((2:10),(1:44));tSinadRatio_6 =A((11:63),(1:44));tSinadRatio_7 =A((64:115),(1:44));tSinadRatio_8 =A((116:167),(1:44));tSinadRatio_9 =A((168:220),(1:44));tSinadRatio_10 =A((221:272),(1:44));tSinadRatio_11 =A((273:324),(1:44));tSinadRatio_12 =A((325:376),(1:44));tSinadRatio_13 =A((377:429),(1:44));tSinadRatio_14 =A((430:481),(1:44));tSinadRatio_15 =A((482:533),(1:44));tSinadRatio_16 =A((534:586),(1:44));tSinadRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',9);
SFDRRatio_5 =A((2:10),(1:44));SFDRRatio_6 =A((11:63),(1:44));SFDRRatio_7 =A((64:115),(1:44));SFDRRatio_8 =A((116:167),(1:44));SFDRRatio_9 =A((168:220),(1:44));SFDRRatio_10 =A((221:272),(1:44));SFDRRatio_11 =A((273:324),(1:44));SFDRRatio_12 =A((325:376),(1:44));SFDRRatio_13 =A((377:429),(1:44));SFDRRatio_14 =A((430:481),(1:44));SFDRRatio_15 =A((482:533),(1:44));SFDRRatio_16 =A((534:586),(1:44));SFDRRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',10);
tSFDRRatio_5 =A((2:10),(1:44));tSFDRRatio_6 =A((11:63),(1:44));tSFDRRatio_7 =A((64:115),(1:44));tSFDRRatio_8 =A((116:167),(1:44));tSFDRRatio_9 =A((168:220),(1:44));tSFDRRatio_10 =A((221:272),(1:44));tSFDRRatio_11 =A((273:324),(1:44));tSFDRRatio_12 =A((325:376),(1:44));tSFDRRatio_13 =A((377:429),(1:44));tSFDRRatio_14 =A((430:481),(1:44));tSFDRRatio_15 =A((482:533),(1:44));tSFDRRatio_16 =A((534:586),(1:44));tSFDRRatio_17 =A((587:611),(1:44));
elseif(beam==2)
[A,B,C]=xlsread('MSTSample.xlsx',11);
ht= A(1,(1:44))';
SNR_5 =A((2:10),(1:44));SNR_6 =A((11:63),(1:44));SNR_7 =A((64:115),(1:44));SNR_8 =A((116:167),(1:44));SNR_9 =A((168:220),(1:44));SNR_10 =A((221:272),(1:44));SNR_11 =A((273:324),(1:44));SNR_12 =A((325:376),(1:44));SNR_13 =A((377:429),(1:44));SNR_14 =A((430:481),(1:44));SNR_15 =A((482:533),(1:44));SNR_16 =A((534:586),(1:44));SNR_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',12);
tSNR_5 =A((2:10),(1:44));tSNR_6 =A((11:63),(1:44));tSNR_7 =A((64:115),(1:44));tSNR_8 =A((116:167),(1:44));tSNR_9 =A((168:220),(1:44));tSNR_10 =A((221:272),(1:44));tSNR_11 =A((273:324),(1:44));tSNR_12 =A((325:376),(1:44));tSNR_13 =A((377:429),(1:44));tSNR_14 =A((430:481),(1:44));tSNR_15 =A((482:533),(1:44));tSNR_16 =A((534:586),(1:44));tSNR_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',13);
BandPower_5 =A((2:10),(1:44));BandPower_6 =A((11:63),(1:44));BandPower_7 =A((64:115),(1:44));BandPower_8 =A((116:167),(1:44));BandPower_9 =A((168:220),(1:44));BandPower_10 =A((221:272),(1:44));BandPower_11 =A((273:324),(1:44));BandPower_12 =A((325:376),(1:44));BandPower_13 =A((377:429),(1:44));BandPower_14 =A((430:481),(1:44));BandPower_15 =A((482:533),(1:44));BandPower_16 =A((534:586),(1:44));BandPower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',14);
tBandPower_5 =A((2:10),(1:44));tBandPower_6 =A((11:63),(1:44));tBandPower_7 =A((64:115),(1:44));tBandPower_8 =A((116:167),(1:44));tBandPower_9 =A((168:220),(1:44));tBandPower_10 =A((221:272),(1:44));tBandPower_11 =A((273:324),(1:44));tBandPower_12 =A((325:376),(1:44));tBandPower_13 =A((377:429),(1:44));tBandPower_14 =A((430:481),(1:44));tBandPower_15 =A((482:533),(1:44));tBandPower_16 =A((534:586),(1:44));tBandPower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',15);
NoisePower_5 =A((2:10),(1:44));NoisePower_6 =A((11:63),(1:44));NoisePower_7 =A((64:115),(1:44));NoisePower_8 =A((116:167),(1:44));NoisePower_9 =A((168:220),(1:44));NoisePower_10 =A((221:272),(1:44));NoisePower_11 =A((273:324),(1:44));NoisePower_12 =A((325:376),(1:44));NoisePower_13 =A((377:429),(1:44));NoisePower_14 =A((430:481),(1:44));NoisePower_15 =A((482:533),(1:44));NoisePower_16 =A((534:586),(1:44));NoisePower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',16);
tNoisePower_5 =A((2:10),(1:44));tNoisePower_6 =A((11:63),(1:44));tNoisePower_7 =A((64:115),(1:44));tNoisePower_8 =A((116:167),(1:44));tNoisePower_9 =A((168:220),(1:44));tNoisePower_10 =A((221:272),(1:44));tNoisePower_11 =A((273:324),(1:44));tNoisePower_12 =A((325:376),(1:44));tNoisePower_13 =A((377:429),(1:44));tNoisePower_14 =A((430:481),(1:44));tNoisePower_15 =A((482:533),(1:44));tNoisePower_16 =A((534:586),(1:44));tNoisePower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',17);
SinadRatio_5 =A((2:10),(1:44));SinadRatio_6 =A((11:63),(1:44));SinadRatio_7 =A((64:115),(1:44));SinadRatio_8 =A((116:167),(1:44));SinadRatio_9 =A((168:220),(1:44));SinadRatio_10 =A((221:272),(1:44));SinadRatio_11 =A((273:324),(1:44));SinadRatio_12 =A((325:376),(1:44));SinadRatio_13 =A((377:429),(1:44));SinadRatio_14 =A((430:481),(1:44));SinadRatio_15 =A((482:533),(1:44));SinadRatio_16 =A((534:586),(1:44));SinadRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',18);
tSinadRatio_5 =A((2:10),(1:44));tSinadRatio_6 =A((11:63),(1:44));tSinadRatio_7 =A((64:115),(1:44));tSinadRatio_8 =A((116:167),(1:44));tSinadRatio_9 =A((168:220),(1:44));tSinadRatio_10 =A((221:272),(1:44));tSinadRatio_11 =A((273:324),(1:44));tSinadRatio_12 =A((325:376),(1:44));tSinadRatio_13 =A((377:429),(1:44));tSinadRatio_14 =A((430:481),(1:44));tSinadRatio_15 =A((482:533),(1:44));tSinadRatio_16 =A((534:586),(1:44));tSinadRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',19);
SFDRRatio_5 =A((2:10),(1:44));SFDRRatio_6 =A((11:63),(1:44));SFDRRatio_7 =A((64:115),(1:44));SFDRRatio_8 =A((116:167),(1:44));SFDRRatio_9 =A((168:220),(1:44));SFDRRatio_10 =A((221:272),(1:44));SFDRRatio_11 =A((273:324),(1:44));SFDRRatio_12 =A((325:376),(1:44));SFDRRatio_13 =A((377:429),(1:44));SFDRRatio_14 =A((430:481),(1:44));SFDRRatio_15 =A((482:533),(1:44));SFDRRatio_16 =A((534:586),(1:44));SFDRRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',20);
tSFDRRatio_5 =A((2:10),(1:44));tSFDRRatio_6 =A((11:63),(1:44));tSFDRRatio_7 =A((64:115),(1:44));tSFDRRatio_8 =A((116:167),(1:44));tSFDRRatio_9 =A((168:220),(1:44));tSFDRRatio_10 =A((221:272),(1:44));tSFDRRatio_11 =A((273:324),(1:44));tSFDRRatio_12 =A((325:376),(1:44));tSFDRRatio_13 =A((377:429),(1:44));tSFDRRatio_14 =A((430:481),(1:44));tSFDRRatio_15 =A((482:533),(1:44));tSFDRRatio_16 =A((534:586),(1:44));tSFDRRatio_17 =A((587:611),(1:44));
elseif(beam==3)
[A,B,C]=xlsread('MSTSample.xlsx',21);
ht= A(1,(1:44))';
SNR_5 =A((2:10),(1:44));SNR_6 =A((11:63),(1:44));SNR_7 =A((64:115),(1:44));SNR_8 =A((116:167),(1:44));SNR_9 =A((168:220),(1:44));SNR_10 =A((221:272),(1:44));SNR_11 =A((273:324),(1:44));SNR_12 =A((325:376),(1:44));SNR_13 =A((377:429),(1:44));SNR_14 =A((430:481),(1:44));SNR_15 =A((482:533),(1:44));SNR_16 =A((534:586),(1:44));SNR_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',22);
tSNR_5 =A((2:10),(1:44));tSNR_6 =A((11:63),(1:44));tSNR_7 =A((64:115),(1:44));tSNR_8 =A((116:167),(1:44));tSNR_9 =A((168:220),(1:44));tSNR_10 =A((221:272),(1:44));tSNR_11 =A((273:324),(1:44));tSNR_12 =A((325:376),(1:44));tSNR_13 =A((377:429),(1:44));tSNR_14 =A((430:481),(1:44));tSNR_15 =A((482:533),(1:44));tSNR_16 =A((534:586),(1:44));tSNR_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',23);
BandPower_5 =A((2:10),(1:44));BandPower_6 =A((11:63),(1:44));BandPower_7 =A((64:115),(1:44));BandPower_8 =A((116:167),(1:44));BandPower_9 =A((168:220),(1:44));BandPower_10 =A((221:272),(1:44));BandPower_11 =A((273:324),(1:44));BandPower_12 =A((325:376),(1:44));BandPower_13 =A((377:429),(1:44));BandPower_14 =A((430:481),(1:44));BandPower_15 =A((482:533),(1:44));BandPower_16 =A((534:586),(1:44));BandPower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',24);
tBandPower_5 =A((2:10),(1:44));tBandPower_6 =A((11:63),(1:44));tBandPower_7 =A((64:115),(1:44));tBandPower_8 =A((116:167),(1:44));tBandPower_9 =A((168:220),(1:44));tBandPower_10 =A((221:272),(1:44));tBandPower_11 =A((273:324),(1:44));tBandPower_12 =A((325:376),(1:44));tBandPower_13 =A((377:429),(1:44));tBandPower_14 =A((430:481),(1:44));tBandPower_15 =A((482:533),(1:44));tBandPower_16 =A((534:586),(1:44));tBandPower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',25);
NoisePower_5 =A((2:10),(1:44));NoisePower_6 =A((11:63),(1:44));NoisePower_7 =A((64:115),(1:44));NoisePower_8 =A((116:167),(1:44));NoisePower_9 =A((168:220),(1:44));NoisePower_10 =A((221:272),(1:44));NoisePower_11 =A((273:324),(1:44));NoisePower_12 =A((325:376),(1:44));NoisePower_13 =A((377:429),(1:44));NoisePower_14 =A((430:481),(1:44));NoisePower_15 =A((482:533),(1:44));NoisePower_16 =A((534:586),(1:44));NoisePower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',26);
tNoisePower_5 =A((2:10),(1:44));tNoisePower_6 =A((11:63),(1:44));tNoisePower_7 =A((64:115),(1:44));tNoisePower_8 =A((116:167),(1:44));tNoisePower_9 =A((168:220),(1:44));tNoisePower_10 =A((221:272),(1:44));tNoisePower_11 =A((273:324),(1:44));tNoisePower_12 =A((325:376),(1:44));tNoisePower_13 =A((377:429),(1:44));tNoisePower_14 =A((430:481),(1:44));tNoisePower_15 =A((482:533),(1:44));tNoisePower_16 =A((534:586),(1:44));tNoisePower_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',27);
SinadRatio_5 =A((2:10),(1:44)); SinadRatio_6 =A((11:63),(1:44));SinadRatio_7 =A((64:115),(1:44));SinadRatio_8 =A((116:167),(1:44));SinadRatio_9 =A((168:220),(1:44));SinadRatio_10 =A((221:272),(1:44));SinadRatio_11 =A((273:324),(1:44));SinadRatio_12 =A((325:376),(1:44));SinadRatio_13 =A((377:429),(1:44));SinadRatio_14 =A((430:481),(1:44));SinadRatio_15 =A((482:533),(1:44));SinadRatio_16 =A((534:586),(1:44));SinadRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',28);
tSinadRatio_5 =A((2:10),(1:44));tSinadRatio_6 =A((11:63),(1:44));tSinadRatio_7 =A((64:115),(1:44));tSinadRatio_8 =A((116:167),(1:44));tSinadRatio_9 =A((168:220),(1:44));tSinadRatio_10 =A((221:272),(1:44));tSinadRatio_11 =A((273:324),(1:44));tSinadRatio_12 =A((325:376),(1:44));tSinadRatio_13 =A((377:429),(1:44));tSinadRatio_14 =A((430:481),(1:44));tSinadRatio_15 =A((482:533),(1:44));tSinadRatio_16 =A((534:586),(1:44));tSinadRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',29);
SFDRRatio_5 =A((2:10),(1:44));SFDRRatio_6 =A((11:63),(1:44));SFDRRatio_7 =A((64:115),(1:44));SFDRRatio_8 =A((116:167),(1:44));SFDRRatio_9 =A((168:220),(1:44));SFDRRatio_10 =A((221:272),(1:44));SFDRRatio_11 =A((273:324),(1:44));SFDRRatio_12 =A((325:376),(1:44));SFDRRatio_13 =A((377:429),(1:44));SFDRRatio_14 =A((430:481),(1:44));SFDRRatio_15 =A((482:533),(1:44));SFDRRatio_16 =A((534:586),(1:44));SFDRRatio_17 =A((587:611),(1:44));
[A,B,C]=xlsread('MSTSample.xlsx',30);
tSFDRRatio_5 =A((2:10),(1:44));tSFDRRatio_6 =A((11:63),(1:44));tSFDRRatio_7 =A((64:115),(1:44));tSFDRRatio_8 =A((116:167),(1:44));tSFDRRatio_9 =A((168:220),(1:44));tSFDRRatio_10 =A((221:272),(1:44));tSFDRRatio_11 =A((273:324),(1:44));tSFDRRatio_12 =A((325:376),(1:44));tSFDRRatio_13 =A((377:429),(1:44));tSFDRRatio_14 =A((430:481),(1:44));tSFDRRatio_15 =A((482:533),(1:44));tSFDRRatio_16 =A((534:586),(1:44));tSFDRRatio_17 =A((587:611),(1:44));
end
MySNR=[SNR_5; SNR_6; SNR_7;SNR_8;SNR_9;SNR_10;SNR_11;SNR_12;SNR_13;SNR_14;SNR_15;SNR_16;SNR_17 ];
MytSNR=[tSNR_5; tSNR_6; tSNR_7;tSNR_8;tSNR_9;tSNR_10;tSNR_11;tSNR_12;tSNR_13;tSNR_14;tSNR_15;tSNR_16;tSNR_17 ];
MyBandPower=[BandPower_5; BandPower_6; BandPower_7;BandPower_8;BandPower_9;BandPower_10;BandPower_11;BandPower_12;BandPower_13;BandPower_14;BandPower_15;BandPower_16;BandPower_17 ];
MytBandPower=[tBandPower_5; tBandPower_6; tBandPower_7;tBandPower_8;tBandPower_9;tBandPower_10;tBandPower_11;tBandPower_12;tBandPower_13;tBandPower_14;tBandPower_15;tBandPower_16;tBandPower_17 ];
MyNoisePower=[NoisePower_5; NoisePower_6; NoisePower_7;NoisePower_8;NoisePower_9;NoisePower_10;NoisePower_11;NoisePower_12;NoisePower_13;NoisePower_14;NoisePower_15;NoisePower_16;NoisePower_17 ];
MytNoisePower=[tNoisePower_5; tNoisePower_6; tNoisePower_7;tNoisePower_8;tNoisePower_9;tNoisePower_10;tNoisePower_11;tNoisePower_12;tNoisePower_13;tNoisePower_14;tNoisePower_15;tNoisePower_16;tNoisePower_17 ];
MySinadRatio=[SinadRatio_5; SinadRatio_6; SinadRatio_7;SinadRatio_8;SinadRatio_9;SinadRatio_10;SinadRatio_11;SinadRatio_12;SinadRatio_13;SinadRatio_14;SinadRatio_15;SinadRatio_16;SinadRatio_17 ];
MytSinadRatio=[tSinadRatio_5; tSinadRatio_6; tSinadRatio_7;tSinadRatio_8;tSinadRatio_9;tSinadRatio_10;tSinadRatio_11;tSinadRatio_12;tSinadRatio_13;tSinadRatio_14;tSinadRatio_15;tSinadRatio_16;tSinadRatio_17 ];
MySFDRRatio=[SFDRRatio_5; SFDRRatio_6; SFDRRatio_7;SFDRRatio_8;SFDRRatio_9;SFDRRatio_10;SFDRRatio_11;SFDRRatio_12;SFDRRatio_13;SFDRRatio_14;SFDRRatio_15;SFDRRatio_16;SFDRRatio_17 ];
MytSFDRRatio=[tSFDRRatio_5; tSFDRRatio_6; tSFDRRatio_7;tSFDRRatio_8;tSFDRRatio_9;tSFDRRatio_10;tSFDRRatio_11;tSFDRRatio_12;tSFDRRatio_13;tSFDRRatio_14;tSFDRRatio_15;tSFDRRatio_16;tSFDRRatio_17 ];
SNR_t_avg = nanmean(MySNR,1)';BandPower_t_avg = nanmean(MyBandPower,1)';NoisePower_t_avg = nanmean(MyNoisePower,1)';SinadRatio_t_avg = nanmean(MySinadRatio,1)';SFDRRatio_t_avg = nanmean(MySFDRRatio,1)';
tSNR_t_avg = nanmean(MytSNR,1)';tBandPower_t_avg = nanmean(MytBandPower,1)';tNoisePower_t_avg = nanmean(MytNoisePower,1)';tSinadRatio_t_avg = nanmean(MytSinadRatio,1)';tSFDRRatio_t_avg = nanmean(MytSFDRRatio,1)';
SNR_t_avg_std = 1*nanstd(MySNR,1)';BandPower_t_avg_std = 1*nanstd(MyBandPower,1)';NoisePower_t_avg_std = 1*nanstd(MyNoisePower,1)';SinadRatio_t_avg_std = 1*nanstd(MySinadRatio,1)';SFDRRatio_t_avg_std = 1*nanstd(MySFDRRatio,1)';
tSNR_t_avg_std = 1*nanstd(MytSNR,1)';tBandPower_t_avg_std = 1*nanstd(MytBandPower,1)';tNoisePower_t_avg_std = 1*nanstd(MytNoisePower,1)';tSinadRatio_t_avg_std = 1*nanstd(MytSinadRatio,1)';tSFDRRatio_t_avg_std = 1*nanstd(MytSFDRRatio,1)';
subplot(4,5,6);
arrr= SNR_t_avg; barr= ht1; x = [arrr barr];
arrr= tSNR_t_avg; barr= ht1; y = [arrr barr]; z=[SNR_t_avg tSNR_t_avg];
out=pdist2(x, y);
out1= diag(out);
datas=[SNR_t_avg ;tSNR_t_avg;out1;SNR_t_avg_std;tSNR_t_avg_std];
arrx = linspace(-16,25,4);
AvgGain = mean(out1);
[MaxVal MaxInd]=max(SNR_t_avg);x = [MaxVal ht1(MaxInd)];
[MaxVal MaxInd]=max(tSNR_t_avg);y = [MaxVal ht1(MaxInd)];
MaxGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
[MinVal MinInd]=min(SNR_t_avg);x = [MinVal ht1(MinInd)];
[MinVal MinInd]=min(tSNR_t_avg);y = [MinVal ht1(MinInd)];
MinGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
hold on;
errorbar(SNR_t_avg, ht1, SNR_t_avg_std, 'horizontal', 'color', '#9696FA', 'LineWidth', .01);
errorbar(tSNR_t_avg, ht1, tSNR_t_avg_std, 'horizontal', 'color', '#96FAFF', 'LineWidth', .01);
plot(tSNR_t_avg,ht1,'color','#01AAAA', 'LineWidth', 0.1,"LineStyle","-");
plot(SNR_t_avg,ht1,'color','#1E28FA', 'LineWidth', 0.1,"LineStyle","-");
[maxv maxi] = max(out1);
plot(out1(maxi),ht1(maxi),'o','MarkerSize',6,'MarkerFaceColor', 'magenta');
xxx=out1';yyy=ht1';
mYVal=AvgGain;
xind = find(xxx<=mYVal,1)
xind = 16
xxx = [xxx(1,1:xind-1 ) mYVal (xxx(1,xind:end))];
yyy = [yyy(1,1:xind-1 ) nan (yyy(1, xind: end))];
adash=[xxx(1,xind-1) yyy(1,xind-1)]'; bdash=[mYVal nan]';cdash =[xxx(1,xind+1) yyy(1,xind+1)]';
p=((cdash(1)-bdash(1))/(bdash(1)-adash(1)));
bdash(2)= (cdash(2) + (adash(2)*p))/(1+p);
yyy(1,xind) = bdash(2);
plot(xxx, yyy,'Color','#D95319','LineWidth', 1.25,"LineStyle","-");
plot(bdash(1),bdash(2),'o','MarkerSize',6,'MarkerFaceColor', 'cyan');
vals = [SNR_t_avg; tSNR_t_avg;out1];mystds=[SNR_t_avg_std;tSNR_t_avg_std;zeros(44,1)];
myacts1 = vals-mystds;
myacts2 = vals+mystds;
datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];
plot((min(datas)),(min(ht1)-1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
plot((max(datas)),(max(ht1)+1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
if(beam==1)
myTitle = strcat('(a) SNR');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
elseif(beam==2)
myTitle = strcat('(a) SNR');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
else
myTitle = strcat('(a) SNR');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
end
set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));
x = xlabel('SNR (dB)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
y = ylabel('Range (km)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
y.FontSize = 12;
ytickformat('%.1f');
set(gca, 'XTick',linspace(min(datas),max(datas),3));
set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
xtickformat('%.0f');
grid minor; grid on;
axis tight;
box on;
SNRAvgGain = AvgGain
SNRAvgGain = 9.8702
SNRMaxGain = MaxGain
SNRMaxGain = 18.5240
subplot(4,5,7);
a=BandPower_t_avg;
b=tBandPower_t_avg;
arrr= a; barr= ht1; x = [arrr barr];
arrr= b; barr= ht1; y = [arrr barr];
out=pdist2(x, y);
out1= diag(out);
datas=[a ;b;out1;BandPower_t_avg_std;tBandPower_t_avg_std];
AvgGain = mean(out1);
[MaxVal MaxInd]=max(a);x = [MaxVal ht1(MaxInd)];
[MaxVal MaxInd]=max(b);y = [MaxVal ht1(MaxInd)];
MaxGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
[MinVal MinInd]=min(a);x = [MinVal ht1(MinInd)];
[MinVal MinInd]=min(b);y = [MinVal ht1(MinInd)];
MinGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
hold on;
errorbar(a, ht1, BandPower_t_avg_std, 'horizontal', 'color', '#9696FA', 'LineWidth', .01);
errorbar(b, ht1, tBandPower_t_avg_std, 'horizontal', 'color', '#96FAFF', 'LineWidth', .01);
plot(b,ht1,'color','#01AAAA', 'LineWidth', 0.1,"LineStyle","-");
plot(a,ht1,'color','#1E28FA', 'LineWidth', 0.1,"LineStyle","-");
[maxv maxi] = max(out1);
plot(out1(maxi),ht1(maxi),'o','MarkerSize',6,'MarkerFaceColor', 'magenta');
xxx=out1';yyy=ht1';
mYVal=AvgGain;
xind = find(xxx<=mYVal,1)
xind = 15
xxx = [xxx(1,1:xind-1 ) mYVal (xxx(1,xind:end))];
yyy = [yyy(1,1:xind-1 ) nan (yyy(1, xind: end))];
adash=[xxx(1,xind-1) yyy(1,xind-1)]'; bdash=[mYVal nan]';cdash =[xxx(1,xind+1) yyy(1,xind+1)]';
p=((cdash(1)-bdash(1))/(bdash(1)-adash(1)));
bdash(2)= (cdash(2) + (adash(2)*p))/(1+p);
yyy(1,xind) = bdash(2);
plot(xxx, yyy,'Color','#D95319','LineWidth', 1.25,"LineStyle","-");
plot(bdash(1),bdash(2),'o','MarkerSize',6,'MarkerFaceColor', 'cyan');
vals = [BandPower_t_avg; tBandPower_t_avg;out1];
mystds=[BandPower_t_avg_std;tBandPower_t_avg_std;zeros(44,1)];
myacts1 = vals-mystds;
myacts2 = vals+mystds;
datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];
plot((min(datas)),(min(ht1)-1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
plot((max(datas)),(max(ht1)+1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
if(beam==1)
myTitle = strcat('(b) Band Power');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
elseif(beam==2)
myTitle = strcat('(b) Band Power');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
else
myTitle = strcat('(b) Band Power');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
end
x = xlabel('Band Power (dB)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
x.FontSize = 12;
set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));
set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
ytickformat('%.1f');
set(gca, 'XTick', linspace(min(datas),max(datas),3));
set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
xtickformat('%.0f');
grid minor; grid on;
box on;
axis tight;
BandPowerAvgGain = AvgGain
BandPowerAvgGain = 31.5664
BandPowerMaxGain = MaxGain
BandPowerMaxGain = 39.3689
subplot(4,5,8);
a=NoisePower_t_avg;b=tNoisePower_t_avg;
arrr= a; barr= ht1; x = [arrr barr];
arrr= b; barr= ht1; y = [arrr barr];
out=pdist2(x, y);
out1= diag(out);
datas=[a ;b;out1;NoisePower_t_avg_std;tNoisePower_t_avg_std];
AvgGain = mean(out1);
[MaxVal MaxInd]=max(a);x = [MaxVal ht1(MaxInd)];
[MaxVal MaxInd]=max(b);y = [MaxVal ht1(MaxInd)];
MaxGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
[MinVal MinInd]=min(a);x = [MinVal ht1(MinInd)];
[MinVal MinInd]=min(b);y = [MinVal ht1(MinInd)];
MinGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
hold on;
errorbar(a, ht1, NoisePower_t_avg_std, 'horizontal', 'color', '#9696FA', 'LineWidth', .01);
errorbar(b, ht1, tNoisePower_t_avg_std, 'horizontal', 'color', '#96FAFF', 'LineWidth', .01);
plot(b,ht1,'color','#01AAAA', 'LineWidth', 0.1,"LineStyle","-");
plot(a,ht1,'color','#1E28FA', 'LineWidth', 0.1,"LineStyle","-");
[maxv maxi] = max(out1);
plot(out1(maxi),ht1(maxi),'o','MarkerSize',6,'MarkerFaceColor', 'magenta');
xxx=out1';yyy=ht1';
mYVal=AvgGain;
xind = find(xxx<=mYVal,1)
xind = 12
xxx = [xxx(1,1:xind-1 ) mYVal (xxx(1,xind:end))];
yyy = [yyy(1,1:xind-1 ) nan (yyy(1, xind: end))];
adash=[xxx(1,xind-1) yyy(1,xind-1)]'; bdash=[mYVal nan]';cdash =[xxx(1,xind+1) yyy(1,xind+1)]';
p=((cdash(1)-bdash(1))/(bdash(1)-adash(1)))
p = 2.2593
bdash(2)= (cdash(2) + (adash(2)*p))/(1+p)
yyy(1,xind) = bdash(2);
plot(xxx, yyy,'Color','#D95319','LineWidth', 1.25,"LineStyle","-");
plot(bdash(1),bdash(2),'o','MarkerSize',6,'MarkerFaceColor', 'cyan');
vals = [NoisePower_t_avg; tNoisePower_t_avg;out1];mystds=[NoisePower_t_avg_std;tNoisePower_t_avg_std;zeros(44,1)];
myacts1 = vals-mystds;
myacts2 = vals+mystds;
datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];
plot((min(datas)),(min(ht1)-1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
plot((max(datas)),(max(ht1)+1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
if(beam==1)
myTitle = strcat('(c) Noise Power');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
elseif(beam==2)
myTitle = strcat('(c) Noise Power');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
else
myTitle = strcat('(c) Noise Power');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
end
x = xlabel('Noise Power (dB)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
x.FontSize = 12;
set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));
ytickformat('%.1f');
set(gca, 'XTick', linspace(min(datas),max(datas),3));
set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
xtickformat('%.0f');
grid minor; grid on; box on;
axis tight;
NoisePowerAvgGain = AvgGain
NoisePowerAvgGain = 9.0623
NoisePowerMaxGain = MaxGain
NoisePowerMaxGain = 24.7362
subplot(4,5,9);
a=SinadRatio_t_avg;b=tSinadRatio_t_avg;
arrr= a; barr= ht1; x = [arrr barr];
arrr= b; barr= ht1; y = [arrr barr];
out=pdist2(x, y);
out1= diag(out);
datas=[a ;b;out1;SinadRatio_t_avg_std;tSinadRatio_t_avg_std];
AvgGain = mean(out1);
[MaxVal MaxInd]=max(a);x = [MaxVal ht1(MaxInd)];
[MaxVal MaxInd]=max(b);y = [MaxVal ht1(MaxInd)];
MaxGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
[MinVal MinInd]=min(a);x = [MinVal ht1(MinInd)];
[MinVal MinInd]=min(b);y = [MinVal ht1(MinInd)];
MinGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
hold on;
errorbar(a, ht1, SinadRatio_t_avg_std, 'horizontal', 'color', '#9696FA', 'LineWidth', .01);
errorbar(b, ht1, tSinadRatio_t_avg_std, 'horizontal', 'color', '#96FAFF', 'LineWidth', .01);
plot(b,ht1,'color','#01AAAA', 'LineWidth', 0.1,"LineStyle","-");
plot(a,ht1,'color','#1E28FA', 'LineWidth', 0.1,"LineStyle","-");
[maxv maxi] = max(out1);
plot(out1(maxi),ht1(maxi),'o','MarkerSize',6,'MarkerFaceColor', 'magenta');
xxx=out1';yyy=ht1';
mYVal=AvgGain;
xind = find(xxx<=mYVal,1)
xind = 16
xxx = [xxx(1,1:xind-1 ) mYVal (xxx(1,xind:end))];
yyy = [yyy(1,1:xind-1 ) nan (yyy(1, xind: end))];
adash=[xxx(1,xind-1) yyy(1,xind-1)]'; bdash=[mYVal nan]';cdash =[xxx(1,xind+1) yyy(1,xind+1)]';
p=((cdash(1)-bdash(1))/(bdash(1)-adash(1)))
p = 0.4377
bdash(2)= (cdash(2) + (adash(2)*p))/(1+p)
yyy(1,xind) = bdash(2);
plot(xxx, yyy,'Color','#D95319','LineWidth', 1.25,"LineStyle","-");
plot(bdash(1),bdash(2),'o','MarkerSize',6,'MarkerFaceColor', 'cyan');
vals = [SinadRatio_t_avg; tSinadRatio_t_avg;out1];mystds=[SinadRatio_t_avg_std;tSinadRatio_t_avg_std;zeros(44,1)];
myacts1 = vals-mystds;
myacts2 = vals+mystds;
datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];
plot((min(datas)),(min(ht1)-1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
plot((max(datas)),(max(ht1)+1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
if(beam==1)
myTitle = strcat('(d) Sinad Ratio');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
elseif(beam==2)
myTitle = strcat('(d) Sinad Ratio');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
else
myTitle = strcat('(d) Sinad Ratio');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
end
x = xlabel('Sinad Ratio (dBc)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
x.FontSize = 12;
set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));
set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
ytickformat('%.1f');
set(gca, 'XTick', linspace(min(datas),max(datas),3));
xtickformat('%.0f');
grid minor; grid on; box on;
axis tight;
SinadRatioAvgGain = AvgGain
SinadRatioAvgGain = 11.0939
SinadRatioMaxGain = MaxGain
SinadRatioMaxGain = 20.4998
subplot(4,5,10);
a=SFDRRatio_t_avg;b=tSFDRRatio_t_avg;
%datas=[a;b];
arrr= a; barr= ht1; x = [arrr barr];
arrr= b; barr= ht1; y = [arrr barr];
out=pdist2(x, y);
out1= diag(out);
datas=[a ;b;out1;SFDRRatio_t_avg_std;tSFDRRatio_t_avg_std];
arrx = linspace(0,30,4);
AvgGain = mean(out1);
[MaxVal MaxInd]=max(a);x = [MaxVal ht1(MaxInd)];
[MaxVal MaxInd]=max(b);y = [MaxVal ht1(MaxInd)];
MaxGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
[MinVal MinInd]=min(a);x = [MinVal ht1(MinInd)];
[MinVal MinInd]=min(b);y = [MinVal ht1(MinInd)];
MinGain = sqrt(((y(1,1)-x(1,1))^2)+((y(1,2)-x(1,2))^2));
hold on;
errorbar(a, ht1, SFDRRatio_t_avg_std, 'horizontal', 'color', '#9696FA', 'LineWidth', .01);
errorbar(b, ht1, tSFDRRatio_t_avg_std, 'horizontal', 'color', '#96FAFF', 'LineWidth', .01);
plot(b,ht1,'color','#01AAAA', 'LineWidth', 0.1,"LineStyle","-");
plot(a,ht1,'color','#1E28FA', 'LineWidth', 0.1,"LineStyle","-");
[maxv maxi] = max(out1);
plot(out1(maxi),ht1(maxi),'o','MarkerSize',6,'MarkerFaceColor', 'magenta');
xxx=out1';yyy=ht1';
mYVal=AvgGain;
xind = find(xxx<=mYVal,1)
xind = 15
xxx = [xxx(1,1:xind-1 ) mYVal (xxx(1,xind:end))];
yyy = [yyy(1,1:xind-1 ) nan (yyy(1, xind: end))];
adash=[xxx(1,xind-1) yyy(1,xind-1)]'; bdash=[mYVal nan]';cdash =[xxx(1,xind+1) yyy(1,xind+1)]';
p=((cdash(1)-bdash(1))/(bdash(1)-adash(1)))
p = 1.4256
bdash(2)= (cdash(2) + (adash(2)*p))/(1+p)
yyy(1,xind) = bdash(2);
plot(xxx, yyy,'Color','#D95319','LineWidth', 1.25,"LineStyle","-");
plot(bdash(1),bdash(2),'o','MarkerSize',6,'MarkerFaceColor', 'cyan');
vals = [SFDRRatio_t_avg; tSFDRRatio_t_avg;out1];mystds=[SFDRRatio_t_avg_std;tSFDRRatio_t_avg_std;zeros(44,1)];
myacts1 = vals-mystds;
myacts2 = vals+mystds;
datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];
plot((min(datas)),(min(ht1)-1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
plot((max(datas)),(max(ht1)+1),'o','MarkerSize',0.1,'MarkerFaceColor', 'white');
if(beam==1)
myTitle = strcat('(e) SFDR Ratio');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
elseif(beam==2)
myTitle = strcat('(e) SFDR Ratio');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
else
myTitle = strcat('(e) SFDR Ratio');
t= title(myTitle,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'bold','HorizontalAlignment',"center");
end
x = xlabel('SFDRRatio (dBc)','FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
x.FontSize = 12;
set(gca, 'YTick',linspace( min(ht1), max(ht1), 7));
set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
ytickformat('%.1f');
datas=[adash ;bdash;out1;xxx';vals((1:44),1);vals((45:88),1);];
set(gca, 'XTick', linspace(min(datas),max(datas),3));
set(gca,'FontSize',11,'FontName','Times New Roman', 'FontWeight', 'Normal');
xtickformat('%.0f');
grid minor; grid on; box on;
axis tight;
legend ({'eBayes','FFT','','','','Difference'},"Location","none","Position",[0.553717928430614 0.72805509250522 0.356545279603305 0.0192592596301325],"Orientation","Horizontal","Box","off","FontSize",11,'FontName','Times New Roman', 'FontWeight', 'Bold');
SFDRRatioAvgGain = AvgGain
SFDRRatioAvgGain = 2.7881
SFDRRatioMaxGain = MaxGain
SFDRRatioMaxGain = 8.2955