source: ETALON/CLIO/Analysis/Data2grating/formfactorComparison.m @ 725

clear all; close all; clc


mt1=1/135*0.75/2;
mt2=1/9*1.41/2*1.46;
dd=3:9;
figure
hold on
grid on
load FFX3mm.mat
 plot (Freq(dd), FF(5,dd)*mt1,'-or','linew',2)
load FFX6mm.mat
plot (Freq(dd), FF(8,dd)*mt2,'-ok','linew',2)
legend('3mm 5.44','6mm 8.24')
ylim([0 1])
set(gca,'fontsize',16)
grid on
xlabel('Frequency, [GHz]')
ylabel('FF, [a.u.]')

figure
hold on
grid on
load FFX3mm.mat
 plot (Freq(dd), FF(4,dd)*mt1,'-or','linew',2)
load FFX6mm.mat
plot (Freq(dd), FF(9,dd)*mt2,'-ok','linew',2)
set(gca,'fontsize',16)
grid on
legend('3mm 5.54','6mm 8.34')
ylim([0 1])
xlabel('Frequency, [GHz]')
ylabel('FF, [a.u.]')
figure
hold on
grid on
load FFX3mm.mat
 plot (Freq(dd), FF(3,dd)*mt1,'-or','linew',2)
load FFX6mm.mat
plot (Freq(dd), FF(10,dd)*mt2,'-ok','linew',2)
set(gca,'fontsize',16)
grid on
legend('3mm 5.64','6mm 8.44')
ylim([0 1])
xlabel('Frequency, [GHz]')
ylabel('FF, [a.u.]')
figure
hold on
grid on
load FFX3mm.mat
 plot (Freq(dd), FF(2,dd)*mt1,'-or','linew',2)
load FFX6mm.mat
plot (Freq(dd), FF(11,dd)*mt2,'-ok','linew',2)
legend('5.74','8.54')
ylim([0 1])
xlabel('Frequency, [GHz]')
ylabel('FF, [a.u.]')
figure
hold on
grid on
load FFX3mm.mat
 plot (Freq(dd), FF(1,dd)*mt1,'-or','linew',2)
load FFX6mm.mat
plot (Freq(dd), FF(12,dd)*mt2,'-ok','linew',2)
legend('5.84','8.64')
ylim([0 1])
xlabel('Frequency, [GHz]')
ylabel('FF, [a.u.]')
figure
hold on
grid on
load FFX3mm.mat
 plot (Freq(dd), FF(6,dd)*mt1,'-or','linew',2)
load FFX6mm.mat
plot (Freq(dd), FF(13,dd)*mt2,'-ok','linew',2)
set(gca,'fontsize',16)
grid on
legend('3mm 5.94','6mm 8.74')
ylim([0 1])
xlabel('Frequency, [GHz]')
ylabel('FF, [a.u.]')
figure
hold on
grid on
load FFX3mm.mat
 plot (Freq(dd), FF(7,dd)*mt1,'-or','linew',2)
load FFX6mm.mat
plot (Freq(dd), FF(14,dd)*mt2,'-ok','linew',2)
legend('6.04','8.94')
ylim([0 1])
xlabel('Frequency, [GHz]')
ylabel('FF, [a.u.]')

% mt1=1/135*0.75;
% mt2=1/9;
% 
% figure
% hold on
% grid on
% load FFX3mm.mat
% FF1=FF;
% load FFX6mm.mat
% FF2=FF;
% dd=4:8;
% plot (FF1(5,dd)*mt1, FF2(8,dd)*mt2,'ok','linew',2)
% plot (FF1(4,dd)*mt1, FF2(9,dd)*mt2,'or','linew',2)
% plot (FF1(3,dd)*mt1, FF2(10,dd)*mt2,'og','linew',2)
% plot (FF1(2,dd)*mt1, FF2(11,dd)*mt2,'ob','linew',2)
% plot (FF1(1,dd)*mt1, FF2(12,dd)*mt2,'*b','linew',2)
% plot (FF1(6,dd)*mt1, FF2(13,dd)*mt2,'*k','linew',2)
% plot (FF1(7,dd)*mt1, FF2(14,dd)*mt2,'*r','linew',2)
% 
% xlim([0 1])
% ylim([0 1])

% i=1;
% j=1;
% for mt1=0.001:0.001:0.01;
% 
% for mt2=0.01:0.01:1;
% Xc=[FF1(5,:) FF1(4,:) FF1(3,:) FF1(2,:) FF1(1,:) FF1(6,:) FF1(7,:)]*mt1;
% Yc=[FF2(8,:) FF2(9,:) FF2(10,:) FF2(11,:) FF2(12,:) FF2(13,:) FF2(14,:)]*mt2;
% rez(i,j)=corr(Xc',Yc');
% i=i+1;
% j=j+1;
% end
% end
% figure
% contourf(rez)