source: ETALON/THZ_filters/2017_THz_FILTER/test_param.m @ 733

   clear all; close all; clc
%% Double Cross filter 
n1=1;
n2=1;
n=1;%ref. indx of media between filters
Thetp=90;

a1=1e-4;
A1=0.53;

c=299792458;
Thet=pi/180*0;
SPlam=c/50e+9;
aar=[0.01 0.1 1 2 5]*1e-3;
color=colormap(jet(length(aar)));

figure
hold on
for cc=1:length(aar)

b=aar(cc);
a=1e-4;
g=(SPlam+4*a+2*b)/2.27;
disp(['g=' num2str(round(g*1e+4)/10) ' a=' num2str(a*1e+3) ' b=' num2str(b*1e+3) ' Freq=' num2str(round(c/SPlam*1e-9)) 'GHz' ])

d=g/2;
loss=0;
c=3e+8;
FreqArr=1e+9:1e+8:2e+11;
for cnt=1:length(FreqArr)

freq=FreqArr(cnt);

  % Stot= S(freq, Thet, g, a, b, n1, n2)*T(freq,loss,d,n)* S(freq, Thet, g, a, b, n1, n2);
    Stot=  S(freq, Thet, g, a, b, n1, n2);

TTE(cnt)=abs(1./Stot(2,2))^2;
TTM(cnt)=abs(1./Stot(4,4))^2;
end
 plot(FreqArr*1e-9,TTE,'color',color(cc,:),'linew',2)
legendInfo{cc} = ['b= ' num2str(aar(cc)*1e+3) ' mm']; % or whatever is appropriate

end
%  Ttot=sin(Thetp).^2.*TTE+cos(Thetp).^2.*TTM;
grid on
set(gca,'fontsize',16)
% set(gca,'xscale','log')
xlabel('Frequency, [GHz]')
ylabel('Transmissivity, [1]')
title('Single cross filter 50GHz')
legend(legendInfo,'location','northeast')