% ERIC LUNDQUIST
% Boeing Wire Measurements
% Connectors May 2009

clear all
close all
clc

load -ASCII coax.dat
load -ASCII co12.dat
load -ASCII co1.dat
load -ASCII co2.dat

% Eliminate Initial Spikes
coax(1:10) = 0;
co12(1:10) = 0;
co1(1:10) = 0;
co2(1:10) = 0;

% Length Axes
m2ft=3.048;
length = 20 ; % length of reflectometry measurement (in meters)
pts = size(coax); %number of points
pts = pts(1);
dz = length/pts*m2ft;
z = dz.*[1:pts];

figure
subplot(4,1,1)
plot(z,coax)
axis tight
grid on
title('Coax','fontweight','b')
subplot(4,1,2)
plot(z,co12)
title('Coax + Cable 1 + Cable 2','fontweight','b')
axis tight
grid on
subplot(4,1,3)
plot(z,co1)
title('Coax + Cable 1','fontweight','b')
axis tight
grid on
subplot(4,1,4)
plot(z,co2)
axis tight
grid on
title('Coax + Cable 2','fontweight','b')
xlabel('Distance (ft)')

% DESCRIPTION OF PARTS
% Cable 1 ~ 22 ft
% Cable 2 ~ 14 ft
% Connectors ~ 2"
% 
% Parts Legend: 
% 
% Cable 1 = D28213J MW D28213P
% Cable 2 = Smith's Aerospace Cheltenham England Part No. 0802KLM03

Coax_Length_ft = 1 
Cable1_Length_ft = 22/m2ft
Cable2_Length_ft = 14/m2ft

Coax_Length_m = 1*m2ft 
Cable1_Length_m = 22
Cable2_Length_m = 14

figure
plot(z,co12)
axis tight
grid on
title('Coax + Cable 1 + Cable 2','fontweight','b')
xlabel('Distance (ft)')
ylabel('Reflection Coefficient')

figure
plot(z,co12-co1)
axis tight
grid on
title('Difference between measurements with and without Cable 2 attached','fontweight','b')
xlabel('Distance (ft)')
ylabel('Reflection Coefficient')