source: ETALON/Filters/Experiment results/Experiment2translationStage.m @ 713

clear all; close all; clc;

F20=[91.2 96 99.2 100.8 99.2 96 92 89.6 88 90.4 88 89.6 90.4 93.6 98.4 104 105.6 104.8 100 96  92.8 92 92.8 92 92.8 93.6];
F23=[100.8 104 105.6 104.8 104 100 97.6 96.8 94.4 92.8 91.2 95.2 100 105.6 108.8 111.2 111.2 108.8 104.8 103.2 100.8 97.6 95.2 96 99.2 103.2];
F24=[96.8 100.8 104.8 105.6 105.6 104 101.6 99.2 96 95.2 96 96.8 98.4 102.4 108 111.2 112.8 112 110.4 105.6 100.8 99.2 98.4 98.4 100 102.4];
F25=[107.2 111.2 114.4 117.6 116.8 115.2 110.4 107.2 102.4 100.8 103.2 105.6 108 112.8 117.6 122.4 124.8 124 121.6 116.8 110.4 107.2 106.4 106.4 109.6 111.2];
F26=[104.8 108 110.4 109.6 108 104.8 103.2 102.4 100.8 100 99.2 100 104.8 111.2 115.2 117.6 118.4 114.4 111.2 110.4 109.6 108.0 104.8 103.2 106.4 110.4];
F27=[107.2 113.6 116.8 116.0 113.6 110.4 106.4 101.6 98.4 96.8 97.6 100 104.8 112 117.6 122.4 122.4 120 114.4 108.8 104.8 103.2 102.4 103.2 104.8 108.8];
F30=[104 110.4 113.6 115.2 113.6 111.2 106.2 101.6 97.6 94.4 96.0 100 105.6 111.2 117.6 120.8 120.8 119.2 115.2 109.6 104.8 102.4 101.6 102.4 104.0 109.6];
F40=[109.6 117.6 120.0 116.8 112.8 106.8 99.2 95.2 91.2 88.8 91.2 97.6 106.4 116 124.2 124.8 121.6 117.6 110.4 104 98.4 95.2 95.2 98.4 104.8 112.0];
F50=[99.2 107.2 109.6 106.4 100.8 92.0 84.0 80.8 78.4 77.6 78.4 86.4 96.0 107.2 112.8 115.2 110.4 104.0 96.0 88.8 84.0 84.0 84.8 88.8 96.0 103.2];
Block=[23 17.8 17 21.4 26.6 34.6 42.0 46.8 46 44.8 41.2 35.6 28.8 20 15.4 16.0 20.6 27.2 35.6 43.8 46.8 46 42.8 37.2 28.4 20];
Frame=[122.4 121.3 126.4 132.0 136 136.8 135.2 132.8 127.2 123.2 121.6 124.0 126.4 125.6 127.2 133.6 140.8 143.2 141.6 140.0 136.0 131.2 125.6 122.4 123.2 124];

F=[F20; F23; F24; F25; F26; F27; F30; F40; F50];
F=F-repmat(Block,[9 1]);
Frame=Frame-Block;
F=F./repmat(Frame,[9 1]);
% F=F./repmat(F(4,:),[9 1]);
nu_f=[20 23 24 25 26 27 30 40 50]
Len=(0:25)*1/2;
figure
[hc hc]=contourf(Len, nu_f,F)
set(hc,'edgecolor','none')
colorbar
set(gca,'fontsize',16)
ylabel('Frequency, [GHz]')
xlabel('Frame position, [mm]')
title('Transmissivity, 24GHz')


figure
hold on
cm=colormap(jet(8));
for ii=1:8
plot(nu_f,F(:,ii+1),'-^','color',cm(9-ii,:),'linew',2)
end
grid on
set(gca,'fontsize',16)
xlabel('Frequency, [GHz]')
ylabel('Transmissivity, [1]')
 legend('0.5 mm','1 mm','1.5 mm','2 mm','2.5 mm','3 mm','3.5 mm','4 mm','Location','northeast')
title('24GHz source')
xlim([19 66])


Fi=pchip(1:26,F,1:0.1:26);
for k=1:9
for c=-50:50
   R(c+51)=corr(Fi(k,1:(126-c+1))',Fi(k,(126+c-1):251)');
end
    [~,Rk(k)]=max(R);
end
% plot(nu_f,(Rk-50+126)/10)
% hold on
% plot(nu_f,ones(size(nu_f))*299/24)
% ylim([12 13.5])
299.792/mean((Rk-50+126)/10)% GHz
299.792/mean((Rk-50+126)/10)^2*std((Rk-50+126)/10)% Error GHz

% figure
% errorbar(nu_f,mean(F,2),std(F'))
% xlim([19 51])