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

%Thz Filter
close all; clear all
Freq=[20 23 24 25 26 27 30 40 50];
Vpp0=[101.6 104.00 108 113.6 122.4 114.4 104.6 87.2 88.]/145.6;
Vpp1=[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];
Vpp0a=[113.6 122.4 118.4 124.0 129.6 124 115.2 98.4 96.0]./[148.8 147.2 145.6 144.0 144.0 143.2 142.4 141.6 142.4];
Vpp1a=[114.4 124.8 116.8 126.4 132.8 124.8 112.4 96.0 92.8]./[148.0 145.6 144.8 144.8 144.0 144.0 143.2 143.2 143.2];
Vpp=[Vpp1; Vpp0a; Vpp1a];

figure
hold on
% plot(Freq,Vpp0,'-ok')
% plot(Freq,Vpp1,'-or','linew',2)
% plot(Freq,Vpp0a,'-ok','linew',2)
% plot(Freq,Vpp1a,'-ob','linew',2)
errorbar(Freq,mean(Vpp),std(Vpp),'-or','linew',2)
set(gca,'fontsize',16)
grid on
xlabel('Frequency, [GHz]')
ylabel('Transmissivity, [1]')
title('GDO-2510, 24 [GHz]')
% legend('Data set 1','Data set 2','Data set 3')
ylim([0.5 1])

% Vpp2=[106.4 104.8 109.6 105.6 104.8 110.4 105.8 114.4 117.6];
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=[53.8 52.8 55.2 48.0 51.2 55.2 45.6 58.4 60.8]./[68 66.4 65.6 65.6 65.6 65.6 65.6 65.6 65.6];
figure
hold on
% plot(Freq,Vpp2)
grid on
plot(Freq,Vpp3,'-or','linew',2)
plot(Freq,Vpp4,'-ok','linew',2)
set(gca,'fontsize',16)
xlabel('Frequency, [GHz]')
ylabel('Transmissivity, [1]')
title('SOL-35310, 35 [GHz]')
legend('With GHA','Without GHA','location','southeast')
ylim([0.5 1])