%Chris Olinger
%Prelab 3, exercise 2

clear
clc

%resonator parameters
d = 0.1;
c = 0.5 * 10^-12;
Zo = 50;
epsilon = 5.0;
frequency = (0.3*10^6):(.1*10^6):(3000*10^6);

wavelength = (3*10^8)./(frequency.*sqrt(epsilon));
beta = 1i*2.*pi./wavelength;
omega = 2.*pi.*frequency;
gamma_load = (1./(1i.*omega.*c))./((1./(1i.*omega.*c)) + 2.*Zo);
Zline = Zo.*(1 + gamma_load.*exp(-2.*beta.*d))./(1 - gamma_load.*exp(-2*beta*d));

V_ratio = ((1j.*omega.*c.*Zline)./(1 + 1j.*omega.*c.*Zline)) ...
          .*((1 + gamma_load)./(1 + gamma_load.*exp(-2.*beta.*d))) ...
          .*((1j.*omega.*c.*Zo)./(1 + 1j.*omega.*c.*Zo)).*exp(-beta.*d);

V_ratio_mag = abs(V_ratio);

V_ratio_mag_dbs = 10.*log10(V_ratio_mag);

figure;
plot(frequency, V_ratio_mag_dbs);
title('Prelab exercise 2')
xlabel('Frequency (Hz)')
ylabel('|V_o/V_i| (dB)')


%exercise 3
gamma_L1 = (Zline + 1./(1i.*omega.*c) - Zo)./(Zline + 1./(1i.*omega.*c) + Zo);

transmission = V_ratio .* (1 + gamma_L1);
transmission_db = 10.*log10(transmission);

figure;
subplot(2,1,1)
plot(frequency, transmission);
title('Prelab Exercise 3')
xlabel('Frequency (Hz)')
ylabel('Transmission')
subplot(2,1,2)
plot(frequency, transmission_db);
xlabel('Frequency (Hz)')
ylabel('Transmission (dB)')


%exercise 4

d = 0.1;
c = 0.05 * 10^-12;
Zo = 50;
epsilon = 5.0;
frequency = (0.3*10^6):(.1*10^6):(3000*10^6);

wavelength = (3*10^8)./(frequency.*sqrt(epsilon));
beta = 1i*2.*pi./wavelength;
omega = 2.*pi.*frequency;
gamma_load = (1./(1i.*omega.*c))./((1./(1i.*omega.*c)) + 2.*Zo);
Zline = Zo.*(1 + gamma_load.*exp(-2.*beta.*d))./(1 - gamma_load.*exp(-2*beta*d));

V_ratio = ((1j.*omega.*c.*Zline)./(1 + 1j.*omega.*c.*Zline)) ...
          .*((1 + gamma_load)./(1 + gamma_load.*exp(-2.*beta.*d))) ...
          .*((1j.*omega.*c.*Zo)./(1 + 1j.*omega.*c.*Zo)).*exp(-beta.*d);

      gamma_L1 = (Zline + 1./(1i.*omega.*c) - Zo)./(Zline + 1./(1i.*omega.*c) + Zo);

transmission_1 = V_ratio .* (1 + gamma_L1);
transmission_db_1 = 10.*log10(transmission_1);

figure;
subplot(2,1,1)
plot(frequency, transmission);
title('Prelab Exercise 4')
xlabel('Frequency (Hz)')
ylabel('Transmission')
hold on
plot(frequency, transmission_1, '--r');
legend('C = 0.5 pF', 'C = 0.05 pF', 'Location', 'NorthEastOutside')
legend('hide')
hold off
subplot(2,1,2)
plot(frequency, transmission_db);
xlabel('Frequency (Hz)')
ylabel('Transmission (dB)')
hold on
plot(frequency, transmission_db_1, '--r');
legend('C = 0.5 pF', 'C = 0.05 pF', 'Location', 'NorthEastOutside')
hold off


%exercise 5

tan_delta1 = 0.01;
tan_delta2 = 0.001;
tan_delta3 = 0.0001;

d = 0.1;
c = 0.05 * 10^-12;
Zo = 50;
frequency = (0.3*10^6):(.1*10^6):(3000*10^6);
omega = 2.*pi.*frequency;

%tan_delta = 0.01
epsilon = 5*(1 - tan_delta1);
wavelength = (3*10^8)./(frequency.*sqrt(epsilon));
beta = 1i*2.*pi./wavelength;

gamma_load = (1./(1i.*omega.*c))./((1./(1i.*omega.*c)) + 2.*Zo);
Zline = Zo.*(1 + gamma_load.*exp(-2.*beta.*d))./(1 - gamma_load.*exp(-2*beta*d));

V_ratio = ((1j.*omega.*c.*Zline)./(1 + 1j.*omega.*c.*Zline)) ...
          .*((1 + gamma_load)./(1 + gamma_load.*exp(-2.*beta.*d))) ...
          .*((1j.*omega.*c.*Zo)./(1 + 1j.*omega.*c.*Zo)).*exp(-beta.*d);

      gamma_L1 = (Zline + 1./(1i.*omega.*c) - Zo)./(Zline + 1./(1i.*omega.*c) + Zo);

transmission_2 = V_ratio .* (1 + gamma_L1);
transmission_db_2 = 10.*log10(transmission_2);

%tan_delta = 0.001
epsilon = 5*(1 - tan_delta2);
wavelength = (3*10^8)./(frequency.*sqrt(epsilon));
beta = 1i*2.*pi./wavelength;

gamma_load = (1./(1i.*omega.*c))./((1./(1i.*omega.*c)) + 2.*Zo);
Zline = Zo.*(1 + gamma_load.*exp(-2.*beta.*d))./(1 - gamma_load.*exp(-2*beta*d));

V_ratio = ((1j.*omega.*c.*Zline)./(1 + 1j.*omega.*c.*Zline)) ...
          .*((1 + gamma_load)./(1 + gamma_load.*exp(-2.*beta.*d))) ...
          .*((1j.*omega.*c.*Zo)./(1 + 1j.*omega.*c.*Zo)).*exp(-beta.*d);

      gamma_L1 = (Zline + 1./(1i.*omega.*c) - Zo)./(Zline + 1./(1i.*omega.*c) + Zo);

transmission_1 = V_ratio .* (1 + gamma_L1);
transmission_db_1 = 10.*log10(transmission_1);

%tan_delta = 0.0001
epsilon = 5*(1 - tan_delta3);
wavelength = (3*10^8)./(frequency.*sqrt(epsilon));
beta = 1i*2.*pi./wavelength;

gamma_load = (1./(1i.*omega.*c))./((1./(1i.*omega.*c)) + 2.*Zo);
Zline = Zo.*(1 + gamma_load.*exp(-2.*beta.*d))./(1 - gamma_load.*exp(-2*beta*d));

V_ratio = ((1j.*omega.*c.*Zline)./(1 + 1j.*omega.*c.*Zline)) ...
          .*((1 + gamma_load)./(1 + gamma_load.*exp(-2.*beta.*d))) ...
          .*((1j.*omega.*c.*Zo)./(1 + 1j.*omega.*c.*Zo)).*exp(-beta.*d);

      gamma_L1 = (Zline + 1./(1i.*omega.*c) - Zo)./(Zline + 1./(1i.*omega.*c) + Zo);

transmission_3 = V_ratio .* (1 + gamma_L1);
transmission_db_3 = 10.*log10(transmission_3);

%plotting
figure;
subplot(2,1,1)
plot(frequency, transmission_1);
title('Prelab Exercise 5')
xlabel('Frequency (Hz)')
ylabel('Transmission')
hold on
plot(frequency, transmission_2, 'r');
plot(frequency, transmission_3, 'g');
legend('tan(delta) = 0.01', 'tan(delta) = 0.001', 'tan(delta) = 0.0001', 'Location', 'NorthEastOutside')
legend('hide')
hold off
subplot(2,1,2)
plot(frequency, transmission_db_1);
xlabel('Frequency (Hz)')
ylabel('Transmission (dB)')
hold on
plot(frequency, transmission_db_2, 'r');
plot(frequency, transmission_db_3, 'g');
legend('tan(delta) = 0.01', 'tan(delta) = 0.001', 'tan(delta) = 0.0001', 'Location', 'NorthEastOutside')
hold off