3 %(4,8) dislocations indicated. Can be used to calculate the defects numbers
with 4 %
no plotting
by setting plotit to 0
5 %
radius defines the distance from 0,0,
for which the defects will be
6 %considered. This can be used to avoid defects that arise naturally from
9 %
DefCount is a vector of the total number of defects counted
for each 19 0.6350 0.0780 0.1840];
21 %Rescale values to be centred
on zero, and in meters
26 DT = delaunayTriangulation([
x y]);
28 triplot(
DT,
'color',[0.6 0.6 0.6],
'LineWidth',2);
32 %
set(gca,
'Color',[0 0 0]);
33 set(gca,
'TickLabelInterpreter',
'latex');
34 set(gca,
'DefaultTextInterpreter',
'Latex')
35 set(gca,
'FontName',
'Latin Modern Roman',
'FontSize',24);
38 xlabel(
'$x$ (m)',
'Interpreter',
'latex');
39 ylabel(
'$y$ (m)',
'Interpreter',
'latex')
46 %Loop over vortices and check the number of
edges per
index, and place
58 if (
length(
DT.vertexAttachments{ii})==6)
62 if (
length(
DT.vertexAttachments{ii})==5)
64 plot(
x(
ii),
y(
ii),
'p',
'MarkerEdgeColor',
'k',
'MarkerSize',
MarkerSize-2,
'MarkerFaceColor',
'w',
'LineWidth',1.5);
67 elseif
length(
DT.vertexAttachments{ii})==7
69 plot(
x(
ii),
y(
ii),
'h',
'MarkerEdgeColor',
'k',
'MarkerSize',
MarkerSize,
'MarkerFaceColor',[0.3 0.3 0.3],
'LineWidth',1.5);
74 if (
length(
DT.vertexAttachments{ii})==3)
76 plot(
x(
ii),
y(
ii),
'^',
'MarkerEdgeColor',
'k',
'MarkerSize',
MarkerSize-2,
'MarkerFaceColor',
RGB(4,:),
'LineWidth',1.5);
79 elseif
length(
DT.vertexAttachments{ii})==9
81 plot(
x(
ii),
y(
ii),
'v',
'MarkerEdgeColor',
'k',
'MarkerSize',
MarkerSize-2,
'MarkerFaceColor',
RGB(5,:),
'LineWidth',1.5);
86 if (
length(
DT.vertexAttachments{ii})==4)
88 plot(
x(
ii),
y(
ii),
'd',
'MarkerEdgeColor',
'k',
'MarkerSize',
MarkerSize-2,
'MarkerFaceColor',
RGB(2,:),
'LineWidth',1.5);
92 elseif
length(
DT.vertexAttachments{ii})==8
94 plot(
x(
ii),
y(
ii),
'o',
'MarkerEdgeColor',
'k',
'MarkerSize',
MarkerSize-2,
'MarkerFaceColor',
RGB(6,:),
'LineWidth',1.5);
% Calculate delaunay triangulation
ylabel(' $y$(m)', 'Interpreter', 'latex') % Defect marker size MarkerSize
and for which the defects will be % considered This can be used to avoid defects that arise naturally from % the edge of the triangulation % h is a handle for generated plot % DefCount is a vector of the total number of defects counted for each % type with the index representing the number of connected edges % some useful colours RGB
% Rescale values to be centred on zero
end % Loop over vortices and check the number of edges per index
and extended by(Dr?) Lee James O 'Riordan. % Ek
elseif length(DT.vertexAttachments{ii})
def delaunay(dataName, dataType, value)
end if sqrt(sum([x(ii), y(ii)].^ 2))< radius %% ignore edges if(length(DT.vertexAttachments
end % Calculate all n choose k pairings of points
if plotit triplot(DT, 'color', [0.6 0.6 0.6], 'LineWidth', 2)
% Indexing needs to % be modified if you wish to use the ordered data sets % Calculate the Voronoi diagram of the resulting data
% Calculate delaunay and setup figure DT
% axis([1e4 1e7 5e-18 1e-10])
% Rescale values to be centred on and in meters x
% If at least one of the indices is
end % Loop over vortices and check the number of edges per and place % appropriate marker for jj
Y Y range of points % radius
% Indexing needs to % be modified if you wish to use the ordered data sets % Calculate the Voronoi diagram of the resulting and plot with the % chosen color scheme The local orientational correlations are returned % as and arr holds the values of the cell color scheme quantity edges
xlabel(' $x$(m)', 'Interpreter', 'latex')
end % idx is for idx is for type for ii