Working Droplet Code

%HOW TO SET UP PREREQUISITE DATA:

%For the case of RA-013-2a 0.25g:
%Import .csv file "RA-013-2a 0.25g 30Hz (1)" as a Numeric Matrix and rename
%it to "RA0132a025g1".
%Simarly, do the same for the other 2 .csv files.

%Import the .csv file with the Row Start and Row End columns and 4 more
%adjacent empty columns.

widthDifference = zeros(15, 1);
middleCell = zeros(15, 1);
newRowStart = zeros(15, 1);
newRowEnd = zeros(15, 1);

%Finding the middle cell of all the droplets
for i = 1:15 %height of droplet matrix
    differenceValueTemp = droplets15(i,2) - droplets15(i,1);
    droplets15(i,3) = differenceValueTemp;
    middleValueTemp = droplets15(i,1) + (differenceValueTemp/2);
    droplets15(i,4) = middleValueTemp;
end

largestDifference = max(droplets15(:,3));

%Assigns new start and end points for the droplets using the
%largestDifference. This is so that each droplet dataset is the same width
%so we can add datasets together to find the average.
%This also makes sure that all droplets are aligned through their centre.
for i = 1:15
    middleValueTemp = droplets15(i,4);
    newRowStartTemp = ceil(middleValueTemp - (largestDifference/2));
    newRowEndTemp = newRowStartTemp + largestDifference;
    droplets15(i,5) = newRowStartTemp;
    droplets15(i,6) = newRowEndTemp;
end

B = RA0132a0g1;
A = RA0132a0g2;
C = RA0132a0g3;

%Locate droplets
%column 1: x position
%column 2: y height
drop1 = A((droplets15(1,5):droplets15(1,6)),(1:2));
drop2 = A((droplets15(2,5):droplets15(2,6)),(1:2));
drop3 = A((droplets15(3,5):droplets15(3,6)),(1:2));
drop4 = A((droplets15(4,5):droplets15(4,6)),(1:2));
drop5 = A((droplets15(5,5):droplets15(5,6)),(1:2));
drop6 = B((droplets15(6,5):droplets15(6,6)),(1:2));
drop7 = B((droplets15(7,5):droplets15(7,6)),(1:2));
drop8 = B((droplets15(8,5):droplets15(8,6)),(1:2));
drop9 = B((droplets15(9,5):droplets15(9,6)),(1:2));
drop10 = B((droplets15(10,5):droplets15(10,6)),(1:2));
drop11 = C((droplets15(11,5):droplets15(11,6)),(1:2));
drop12 = C((droplets15(12,5):droplets15(12,6)),(1:2));
drop13 = C((droplets15(13,5):droplets15(13,6)),(1:2));
drop14 = C((droplets15(14,5):droplets15(14,6)),(1:2));
drop15 = C((droplets15(15,5):droplets15(15,6)),(1:2));

%Add all the droplets to a single matrix to make operations easier
largeDropletMatrix = [drop1, drop2, drop3, drop4, drop5, drop6, drop7, drop8, drop9, drop10, drop11, drop12, drop13, drop14, drop15];

%On either side of each droplet, make height zero so that adjacent droplet
%data doesn't affect the average.
for k = 1:2:30
    dropletNumber = (k + 1)/2;
    droplet = largeDropletMatrix(:,(k:(k+1)));
    %left side of drop
    for i = 1:(droplets15(dropletNumber,1) - droplets15(dropletNumber,5) - 1)
        if droplets15(dropletNumber,1) ~= droplets15(dropletNumber,5)
            largeDropletMatrix(i,(k+1)) = 0;
        end
    end
    %right side of drop
    for i = (size(droplet,1)-(droplets15(dropletNumber,6)-droplets15(dropletNumber,2))):size(droplet,1)
        if droplets15(dropletNumber,6) ~= droplets15(dropletNumber,2)
            largeDropletMatrix(i,(k+1)) = 0;
        end
    end
end

%Finding average droplet height over x axis.
%Add all the height values of each droplet together and divide by the
%number of droplets.
meanDroplet = zeros((size(largeDropletMatrix,1)),2);
for k = 1:2:30
    meanDroplet = meanDroplet + largeDropletMatrix(:,(k:(k+1)));
end
meanDroplet = meanDroplet/(size(largeDropletMatrix,2)/2);
meanDroplet2 = meanDroplet;

i = largestDifference + 1;
j = 0;

%Resets all x axis data of droplets to range from 0 to the width of the
%largest droplet so that they can be overlapped.
for r = 1:i
    for q = 1:2:30
        meanDroplet2(r,1) = j;
        %Overlap all droplets
        largeDropletMatrix(r,q) = j;
    end
    j = j + 0.00005;
end

figure;

x = meanDroplet2(:,1);
y = meanDroplet2(:,2);
p = plot(x,y, 'r');
title('RA-013-2a Glycerol: 0.25g');
xlabel('Length [mm]')
ylabel('Height [nm]')
p.LineWidth = 3;

hold on
y2 = largeDropletMatrix(:,2);
plot(x,y2)
y3 = largeDropletMatrix(:,4);
plot(x,y3)
y4 = largeDropletMatrix(:,6);
plot(x,y4)
y5 = largeDropletMatrix(:,8);
plot(x,y5)
y6 = largeDropletMatrix(:,10);
plot(x,y6)
y7 = largeDropletMatrix(:,12);
plot(x,y7)
y8 = largeDropletMatrix(:,14);
plot(x,y8)
y9 = largeDropletMatrix(:,16);
plot(x,y9)
y10 = largeDropletMatrix(:,18);
plot(x,y10)
y11 = largeDropletMatrix(:,20);
plot(x,y11)
y12 = largeDropletMatrix(:,22);
plot(x,y12)
y13 = largeDropletMatrix(:,24);
plot(x,y13)
y14 = largeDropletMatrix(:,26);
plot(x,y14)
y15 = largeDropletMatrix(:,28);
plot(x,y15)
y16 = largeDropletMatrix(:,30);
plot(x,y16)
hold off