source: ETALON/THZ_filters/Experiment results/MakePlots.m @ 733

%%plots
close all; clear all; clc
figure(2)
hold on

load smallFc.mat
plot(nu_filter*1e-9,F24i,'-.ok','linew',2,'MarkerFaceColor','k')
load smallF.mat
plot(nu_filter*1e-9,F24i,'-.ob','linew',2,'MarkerFaceColor','b')
F0=[350 230 180 104 60 24]/438;
% plot(nu_filter*1e-9,F0,'-or','linew',2)

F37=[79.2 119.3 85.7 93.2 79.4 88 89.7 81.8];
F53=[40.7 38.5 51.2 48.3 47.1 34.6 47.5 53.5];
F75=[42.3 23.6 45.6 22.2 30.9 15.2 29.8 25.4];
F124=[18.5 8.18 17.4 16.7 6.15 9.84 7.0 7.5];
F186=[7.8 11.4 9.5 6.4 4.8 4.1 4.9 4.5 ];
F368=[11.03 3.57 5.36 5.7 7.1 2.53 4.5 2.3];
Frame=[95.8 92.5 92.7 92.05 92 92.1];
Fx=[mean(F37) mean(F53) mean(F75) mean(F124) mean(F186) mean(F368)]./mean(Frame);
eFx=[std(F37) std(F53) std(F75) std(F124) std(F186) std(F368)]./mean(Frame);
errorbar(nu_filter*1e-9,Fx,eFx,'^r','linew',2)
grid on
set(gca,'fontsize',16)
xlabel('Frequency, [GHz]')
ylabel('Transmissivity, [1]')
 legend('Calc. Z_0=2.3', 'Calc. Z_0=1.9','Exp. data')
grid on
title('24GHz source')

figure(3)
hold on

load smallFc.mat
plot(nu_filter*1e-9,F35i,'-.ok','MarkerFaceColor','k','linew',2)
load smallF.mat
plot(nu_filter*1e-9,F35i,'-.ob','MarkerFaceColor','b','linew',2)

% F0=[224 122 88 40 19.8 5.12]/215;
% plot(nu_filter*1e-9,F0,'-or','linew',2)

F37=[74.9 74.4 81.33 83.06 75.6 83.3 90.74 80.42 79.78 75.7];
F53=[63.8 58.6 61.2 57.2 69.7 50.2 94.4 52.01 88.2 54.73];
F75=[51.5 35.7 47.13 43.57 53.47 30.16 64.12 41.9 52.5 35.53];
F124=[33.1 10.79 11.06 24.49 11.35 17.22 12.3 25 14.7 21.6];
F186=[11.38 4.3 9.15 21.64 7.89 6.7 10.3 24.4 13 25.3];
F368=[2.6 0.5 0.9 2.12 2.24 2.3 2.8 1.26 7.17 1.68];
Frame=[74.5 72.6 73.1 72.8 73.56 74.2 74.1];
Fx=[mean(F37) mean(F53) mean(F75) mean(F124) mean(F186) mean(F368)]./mean(Frame);
eFx=[std(F37) std(F53) std(F75) std(F124) std(F186) std(F368)]./mean(Frame);

errorbar(nu_filter*1e-9,Fx,eFx,'^r','linew',2)
grid on
set(gca,'fontsize',16)
xlabel('Frequency, [GHz]')
ylabel('Transmissivity, [1]')
 legend('Calc. Z_0=2.3', 'Calc. Z_0=1.9','Exp. data')
title('35GHz source')
load bigF.mat



figure(4)
plot(nu_filter*1e-9,F24i,'-ok','linew',2)
nu_f0=[37 53 75 124 186 368];
F=[416 368 364 390 420 422 382 352 328]/438;
F0=[350 230 180 104 60 24]/438;

% Freq=[20 23 24 25 26 27 30 40 50];
% Vpp0=[116.8 124.0 120.0 126.4 133.6 124.8 115.2 100 97.6]./[153.6 150.4 148.8 147.2 145.2 145.6 144.8 144.8 143.2];
hold on
% plot(Freq,Vpp0,'-ok')
% plot(Freq,F,'-og')
% % plot(nu_f0,F0,'-or')
% ylim([0 1])
F20=[89.1 84.8 83.2 81.2];
F23=[72 71.4 70.3 72];
F24=[67.7 71 70.5 72.7 ];
F25=[74.5 74.6 74.2 74.2 ];
F26=[81.7 82 81 82.4];
F27=[88 84.2 87.8 88.5];
F30=[77.6 77.4 76.6 79.7];
F40=[61.3 60.21 60.9 62];
F50=[66.2 66 67.3 66];
Frame=[94.3 93.8 93 92.3 93.7 93.2 93.1];
Fx=[mean(F20) mean(F23) mean(F24) mean(F25) mean(F26) mean(F27) mean(F30) mean(F40) mean(F50)]./mean(Frame);
eFx=[std(F20) std(F23) std(F24) std(F25) std(F26) std(F27) std(F30) std(F40) std(F50)]./mean(Frame);
errorbar(nu_filter*1e-9,Fx,eFx,'-.^r','linew',2)
grid on
set(gca,'fontsize',16)
xlabel('Frequency, [GHz]')
ylabel('Transmissivity, [1]')
 legend('Prediction', 'Experiment','location','southeast')
title('24GHz source')
load bigF.mat
ylim([0 1])


figure(5)
plot(nu_filter*1e-9,F35i,'-ok','linew',2)
% Vpp3=[112.8 108.8 111.2 107.2 108.0 112.8 107.2 117.6 118.4]./[122.4 121.6 120.8 120.0 120.0 120.0 120.0 119.2 119.2];
% Vpp4=[260 240 240 230 238 212 192 246 252]/226
 hold on
% plot(Freq,Vpp4,'-or')
% plot(Freq,Vpp3,'-ok')
F20=[94.7 95.4 96.2 94.2 ];
F23=[ 90.3 93.3 87 86.6 ];
F24=[ 87.6 91.6 88.4 89.5 ];
F25=[87.7 88 88.3 87.5];
F26=[85.2 83.9 83.1 84.6];
F27=[84.9 83 80.5 77.3];
F30=[82.6 82 83.3 85.5];
F40=[87 89.5 89.2 88];
F50=[91.4 94.2 94.8 90.7];
Frame=[96.3 93.3 92.7 92 91.8 91.6 91 91.2];
Fx=[mean(F20) mean(F23) mean(F24) mean(F25) mean(F26) mean(F27) mean(F30) mean(F40) mean(F50)]./mean(Frame);
eFx=[std(F20) std(F23) std(F24) std(F25) std(F26) std(F27) std(F30) std(F40) std(F50)]./mean(Frame);
errorbar(nu_filter*1e-9,Fx,eFx,'-.^r','linew',2)
grid on
set(gca,'fontsize',16)
xlabel('Frequency, [GHz]')
ylabel('Transmissivity, [1]')
 legend('Prediction', 'Experiment','location','southeast')
title('35GHz source')
ylim([0 1.1])