Skip to content

Commit fcd71ec

Browse files
emazzochetteemazzochette
committed
Aggregation Files Update
1 parent 4b6f077 commit fcd71ec

6 files changed

Lines changed: 615 additions & 59 deletions

File tree

‎.DS_Store‎

-5 KB
Binary file not shown.

‎DeltaVAggregate.m‎

Lines changed: 318 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -17,7 +17,6 @@
1717
Profiles.(['Force_',num2str(Forces(force))]) = [];
1818
end
1919

20-
2120
%%
2221
for trial = 1:totalTrials
2322
if ismember(trial,[2943:2945])
@@ -40,7 +39,7 @@
4039

4140
if strcmp(trialScore{1,1}{trial,1},'TRUE')
4241

43-
stim = stimNums{1,1}(trial);
42+
stim = stimNums{1,1}(trial);
4443

4544
[ unit.deltaVNorm ] = calculateDeltaVTrace( Stimulus(stim).CurvatureAnalysis.velocitySmoothed, Stimulus(stim).timeData(:,8), Stimulus(stim).StimulusTiming.stimOnFrame );
4645

@@ -62,52 +61,351 @@
6261

6362
end
6463

65-
%%
64+
%% colors:
65+
color_50 = [228/255 26/255 28/255];
66+
color_100 = [55/255 126/255 184/255];
67+
color_500 = [77/255 175/255 74/255];
68+
color_1000 = [152/255 78/255 163/255];
69+
color_5000 = [255/255 127/255 0/255];
70+
color_10000 = [255/255 255/255 51/255];
71+
72+
%% Plot all data points:
73+
6674
markerSize = 8;
6775

68-
for trialPlot = 1:length(Profiles.Force_10000)
69-
plot(Profiles.Force_10000(trialPlot).positionApplied,Profiles.Force_10000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [255/255 255/255 51/255],'MarkerEdgeColor','none')
76+
maxTrials = max([length(Profiles.Force_50),length(Profiles.Force_100)]);
77+
maxTrials2 = max([length(Profiles.Force_500),length(Profiles.Force_1000),length(Profiles.Force_5000),length(Profiles.Force_10000)]);
78+
for trialPlot = maxTrials2+1:maxTrials
79+
% for trialPlot = 1:length(Profiles.Force_100)
80+
try plot(Profiles.Force_100(trialPlot).positionApplied,Profiles.Force_100(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_100,'MarkerEdgeColor','none'); end
7081
hold on
7182

83+
% end
84+
%
85+
% for trialPlot = 1:length(Profiles.Force_50)
86+
try plot(Profiles.Force_50(trialPlot).positionApplied,Profiles.Force_50(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_50,'MarkerEdgeColor','none'); end
87+
hold on
7288
end
7389

74-
for trialPlot = 1:length(Profiles.Force_5000)
75-
plot(Profiles.Force_5000(trialPlot).positionApplied,Profiles.Force_5000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [255/255 127/255 0/255],'MarkerEdgeColor','none')
90+
for trialPlot = 1:maxTrials2
91+
92+
% for trialPlot = 1:length(Profiles.Force_100)
93+
try plot(Profiles.Force_100(trialPlot).positionApplied,Profiles.Force_100(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_100,'MarkerEdgeColor','none'); end
7694
hold on
7795

78-
end
96+
% end
97+
%
98+
% for trialPlot = 1:length(Profiles.Force_50)
99+
try plot(Profiles.Force_50(trialPlot).positionApplied,Profiles.Force_50(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_50,'MarkerEdgeColor','none'); end
100+
hold on
101+
102+
try plot(Profiles.Force_10000(trialPlot).positionApplied,Profiles.Force_10000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_10000,'MarkerEdgeColor','none'); end
103+
hold on
79104

80-
for trialPlot = 1:length(Profiles.Force_1000)
81-
plot(Profiles.Force_1000(trialPlot).positionApplied,Profiles.Force_1000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [152/255 78/255 163/255],'MarkerEdgeColor','none')
105+
% end
106+
%
107+
% for trialPlot = 1:length(Profiles.Force_5000)
108+
try plot(Profiles.Force_5000(trialPlot).positionApplied,Profiles.Force_5000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_5000,'MarkerEdgeColor','none'); end
82109
hold on
83110

84-
end
111+
% end
112+
%
113+
% for trialPlot = 1:length(Profiles.Force_1000)
114+
try plot(Profiles.Force_1000(trialPlot).positionApplied,Profiles.Force_1000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_1000,'MarkerEdgeColor','none'); end
115+
hold on
85116

86-
for trialPlot = 1:length(Profiles.Force_500)
87-
plot(Profiles.Force_500(trialPlot).positionApplied,Profiles.Force_500(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [77/255 175/255 74/255],'MarkerEdgeColor','none')
117+
% end
118+
%
119+
% for trialPlot = 1:length(Profiles.Force_500)
120+
try plot(Profiles.Force_500(trialPlot).positionApplied,Profiles.Force_500(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_500,'MarkerEdgeColor','none'); end
88121
hold on
89122

123+
% end
124+
%
125+
90126
end
91127

92-
for trialPlot = 1:length(Profiles.Force_100)
93-
plot(Profiles.Force_100(trialPlot).positionApplied,Profiles.Force_100(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [55/255 126/255 184/255],'MarkerEdgeColor','none')
94-
hold on
128+
%%
129+
130+
axis([0 100 -6 12])
131+
set(gca,'FontSize',20)
132+
[xNew yNew] = MiriamAxes(gca,'xy');
95133

134+
%% rebin by target
135+
markerSize = 10;
136+
lineWidth = 4;
137+
targetBins = [10:10:90];
138+
for bin = 1:length(targetBins)
139+
Bins_50.(['TargetBin_',num2str(bin)]) = [];
140+
Bins_100.(['TargetBin_',num2str(bin)]) = [];
141+
Bins_500.(['TargetBin_',num2str(bin)]) = [];
142+
Bins_1000.(['TargetBin_',num2str(bin)]) = [];
143+
Bins_5000.(['TargetBin_',num2str(bin)]) = [];
144+
Bins_10000.(['TargetBin_',num2str(bin)]) = [];
145+
Bins_50.average = [];
146+
Bins_100.average = [];
147+
Bins_500.average = [];
148+
Bins_1000.average = [];
149+
Bins_5000.average = [];
150+
Bins_10000.average = [];
96151
end
97152

98153
for trialPlot = 1:length(Profiles.Force_50)
99-
plot(Profiles.Force_50(trialPlot).positionApplied,Profiles.Force_50(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [228/255 26/255 28/255],'MarkerEdgeColor','none')
100-
hold on
154+
bin = find(targetBins>Profiles.Force_50(trialPlot).positionApplied,1,'first')-1;
155+
Bins_50.(['TargetBin_',num2str(bin)]) = [Bins_50.(['TargetBin_',num2str(bin)]); Profiles.Force_50(trialPlot).positionApplied Profiles.Force_50(trialPlot).maxDeltaVNorm];
101156
end
102157

158+
for trialPlot = 1:length(Profiles.Force_100)
159+
bin = find(targetBins>Profiles.Force_100(trialPlot).positionApplied,1,'first')-1;
160+
Bins_100.(['TargetBin_',num2str(bin)]) = [Bins_100.(['TargetBin_',num2str(bin)]); Profiles.Force_100(trialPlot).positionApplied Profiles.Force_100(trialPlot).maxDeltaVNorm];
161+
end
162+
163+
for trialPlot = 1:length(Profiles.Force_500)
164+
bin = find(targetBins>Profiles.Force_500(trialPlot).positionApplied,1,'first')-1;
165+
Bins_500.(['TargetBin_',num2str(bin)]) = [Bins_500.(['TargetBin_',num2str(bin)]); Profiles.Force_500(trialPlot).positionApplied Profiles.Force_500(trialPlot).maxDeltaVNorm];
166+
end
103167

168+
for trialPlot = 1:length(Profiles.Force_1000)
169+
bin = find(targetBins>Profiles.Force_1000(trialPlot).positionApplied,1,'first')-1;
170+
Bins_1000.(['TargetBin_',num2str(bin)]) = [Bins_1000.(['TargetBin_',num2str(bin)]); Profiles.Force_1000(trialPlot).positionApplied Profiles.Force_1000(trialPlot).maxDeltaVNorm];
171+
end
104172

173+
for trialPlot = 1:length(Profiles.Force_5000)
174+
bin = find(targetBins>Profiles.Force_5000(trialPlot).positionApplied,1,'first')-1;
175+
Bins_5000.(['TargetBin_',num2str(bin)]) = [Bins_5000.(['TargetBin_',num2str(bin)]); Profiles.Force_5000(trialPlot).positionApplied Profiles.Force_5000(trialPlot).maxDeltaVNorm];
176+
end
105177

178+
for trialPlot = 1:length(Profiles.Force_10000)
179+
bin = find(targetBins>Profiles.Force_10000(trialPlot).positionApplied,1,'first')-1;
180+
Bins_10000.(['TargetBin_',num2str(bin)]) = [Bins_10000.(['TargetBin_',num2str(bin)]); Profiles.Force_10000(trialPlot).positionApplied Profiles.Force_10000(trialPlot).maxDeltaVNorm];
181+
end
182+
minTrials = 5;
183+
for bin = 1:length(targetBins)
106184

107185

186+
if size(Bins_50.(['TargetBin_',num2str(bin)]),1)> minTrials
187+
Bins_50.average = [Bins_50.average; mean( Bins_50.(['TargetBin_',num2str(bin)]),1) std(Bins_50.(['TargetBin_',num2str(bin)]),1), length(Bins_50.(['TargetBin_',num2str(bin)]))];
188+
plot(Bins_50.average(:,1),Bins_50.average(:,2), 'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [228/255 26/255 28/255],'MarkerEdgeColor','none','Color', [228/255 26/255 28/255],'LineWidth',lineWidth);
189+
hold on
190+
end
191+
192+
193+
if size(Bins_100.(['TargetBin_',num2str(bin)]),1)> minTrials
194+
Bins_100.average = [Bins_100.average; mean( Bins_100.(['TargetBin_',num2str(bin)]),1) std(Bins_100.(['TargetBin_',num2str(bin)]),1), length(Bins_100.(['TargetBin_',num2str(bin)]))];
195+
plot(Bins_100.average(:,1),Bins_100.average(:,2), 'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [55/255 126/255 184/255],'MarkerEdgeColor','none','Color', [55/255 126/255 184/255],'LineWidth',lineWidth);
196+
hold on
197+
end
198+
108199

200+
if size(Bins_500.(['TargetBin_',num2str(bin)]),1)> minTrials
201+
Bins_500.average = [Bins_500.average; mean( Bins_500.(['TargetBin_',num2str(bin)]),1) std(Bins_500.(['TargetBin_',num2str(bin)]),1), length(Bins_500.(['TargetBin_',num2str(bin)]))];
202+
plot(Bins_500.average(:,1),Bins_500.average(:,2), 'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [77/255 175/255 74/255],'MarkerEdgeColor','none','Color', [77/255 175/255 74/255],'LineWidth',lineWidth);
203+
hold on
204+
end
205+
206+
207+
if size(Bins_1000.(['TargetBin_',num2str(bin)]),1)> minTrials
208+
Bins_1000.average = [Bins_1000.average; mean( Bins_1000.(['TargetBin_',num2str(bin)]),1) std(Bins_1000.(['TargetBin_',num2str(bin)]),1), length(Bins_1000.(['TargetBin_',num2str(bin)]))];
209+
plot(Bins_1000.average(:,1),Bins_1000.average(:,2), 'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [152/255 78/255 163/255],'MarkerEdgeColor','none','Color', [152/255 78/255 163/255],'LineWidth',lineWidth);
210+
hold on
211+
end
212+
109213

110-
%%
111-
axis([0 100 -6 12])
214+
if size(Bins_5000.(['TargetBin_',num2str(bin)]),1)> minTrials
215+
Bins_5000.average = [Bins_5000.average; mean( Bins_5000.(['TargetBin_',num2str(bin)]),1) std(Bins_5000.(['TargetBin_',num2str(bin)]),1), length(Bins_5000.(['TargetBin_',num2str(bin)]))];
216+
plot(Bins_5000.average(:,1),Bins_5000.average(:,2), 'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [255/255 127/255 0/255],'MarkerEdgeColor','none','Color', [255/255 127/255 0/255],'LineWidth',lineWidth);
217+
hold on
218+
end
219+
220+
221+
if size(Bins_10000.(['TargetBin_',num2str(bin)]),1)> minTrials
222+
Bins_10000.average = [Bins_10000.average; mean( Bins_10000.(['TargetBin_',num2str(bin)]),1) std(Bins_10000.(['TargetBin_',num2str(bin)]),1), length(Bins_10000.(['TargetBin_',num2str(bin)]))];
223+
plot(Bins_10000.average(:,1),Bins_10000.average(:,2), 'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', [255/255 255/255 51/255],'MarkerEdgeColor','none','Color', [255/255 255/255 51/255],'LineWidth',lineWidth);
224+
hold on
225+
end
226+
227+
end
228+
229+
axis([0 100 -3 3])
112230
set(gca,'FontSize',20)
113231
[xNew yNew] = MiriamAxes(gca,'xy');
232+
233+
234+
%% Organize by Force, Create histograms:
235+
236+
for trial = 1:length(Profiles.Force_50)
237+
data_50(trial) = Profiles.Force_50(trial).maxDeltaVNorm;
238+
end
239+
for trial = 1:length(Profiles.Force_100)
240+
data_100(trial) = Profiles.Force_100(trial).maxDeltaVNorm;
241+
end
242+
243+
for trial = 1:length(Profiles.Force_500)
244+
data_500(trial) = Profiles.Force_500(trial).maxDeltaVNorm;
245+
end
246+
247+
for trial = 1:length(Profiles.Force_1000)
248+
data_1000(trial) = Profiles.Force_1000(trial).maxDeltaVNorm;
249+
end
250+
251+
for trial = 1:length(Profiles.Force_5000)
252+
data_5000(trial) = Profiles.Force_5000(trial).maxDeltaVNorm;
253+
end
254+
255+
for trial = 1:length(Profiles.Force_10000)
256+
data_10000(trial) = Profiles.Force_10000(trial).maxDeltaVNorm;
257+
end
258+
259+
bins = [-6:0.5:8];
260+
reversalThreshold = 0.75;
261+
reversalBound = [reversalThreshold reversalThreshold];
262+
speedUpThreshold = -1;
263+
speedUpBound = [speedUpThreshold speedUpThreshold];
264+
265+
close all
266+
xLim = [-5 7];
267+
yLim = [0 0.6];
268+
fontSize = 27;
269+
lineWidth = 4;
270+
lineStyle = '--';
271+
displayStyle = {'bar'};
272+
reversalLineColor = [0 0 0];%{'k'};
273+
speedUpLineColor = [0 0 1];%{'b'};
274+
275+
276+
figure(1)
277+
histogram(data_50, bins, 'Normalization', 'probability' ,'DisplayStyle',displayStyle,'FaceColor',color_50, 'EdgeColor','none');
278+
hold on
279+
plot(reversalBound,yLim,'LineWidth',lineWidth,'Color',reversalLineColor,'LineStyle',lineStyle);
280+
plot(speedUpBound,yLim,'LineWidth',lineWidth,'Color',speedUpLineColor,'LineStyle',lineStyle);
281+
set(gca,'FontSize',fontSize)
282+
set(gca,'ylim', yLim);
283+
set(gca,'xlim', xLim);
284+
[xNew yNew] = MiriamAxes(gca,'xy');
285+
286+
figure(2)
287+
histogram(data_100, bins, 'Normalization', 'probability','DisplayStyle',displayStyle,'FaceColor',color_100, 'EdgeColor','none')
288+
hold on
289+
plot(reversalBound,yLim,'LineWidth',lineWidth,'Color',reversalLineColor,'LineStyle',lineStyle);
290+
plot(speedUpBound,yLim,'LineWidth',lineWidth,'Color',speedUpLineColor,'LineStyle',lineStyle);
291+
set(gca,'FontSize',fontSize)
292+
set(gca,'ylim', yLim);
293+
set(gca,'xlim', xLim);
294+
[xNew yNew] = MiriamAxes(gca,'xy');
295+
296+
figure(3)
297+
histogram(data_500, bins, 'Normalization', 'probability','DisplayStyle',displayStyle,'FaceColor',color_500, 'EdgeColor','none')
298+
hold on
299+
plot(reversalBound,yLim,'LineWidth',lineWidth,'Color',reversalLineColor,'LineStyle',lineStyle);
300+
plot(speedUpBound,yLim,'LineWidth',lineWidth,'Color',speedUpLineColor,'LineStyle',lineStyle);
301+
set(gca,'FontSize',fontSize)
302+
set(gca,'ylim', yLim);
303+
set(gca,'xlim', xLim);
304+
[xNew yNew] = MiriamAxes(gca,'xy');
305+
306+
figure(4)
307+
histogram(data_1000, bins, 'Normalization', 'probability','DisplayStyle',displayStyle,'FaceColor',color_1000, 'EdgeColor','none')
308+
hold on
309+
plot(reversalBound,yLim,'LineWidth',lineWidth,'Color',reversalLineColor,'LineStyle',lineStyle);
310+
plot(speedUpBound,yLim,'LineWidth',lineWidth,'Color',speedUpLineColor,'LineStyle',lineStyle);
311+
set(gca,'FontSize',fontSize)
312+
set(gca,'ylim', yLim);
313+
set(gca,'xlim', xLim);
314+
[xNew yNew] = MiriamAxes(gca,'xy');
315+
316+
figure(5)
317+
histogram(data_5000, bins, 'Normalization', 'probability','DisplayStyle',displayStyle,'FaceColor',color_5000, 'EdgeColor','none')
318+
hold on
319+
plot(reversalBound,yLim,'LineWidth',lineWidth,'Color',reversalLineColor,'LineStyle',lineStyle);
320+
plot(speedUpBound,yLim,'LineWidth',lineWidth,'Color',speedUpLineColor,'LineStyle',lineStyle);
321+
set(gca,'FontSize',fontSize)
322+
set(gca,'ylim', yLim);
323+
set(gca,'xlim', xLim);
324+
[xNew yNew] = MiriamAxes(gca,'xy');
325+
326+
figure(6)
327+
histogram(data_10000, bins, 'Normalization', 'probability','DisplayStyle',displayStyle,'FaceColor',color_10000, 'EdgeColor','none')
328+
hold on
329+
plot(reversalBound,yLim,'LineWidth',lineWidth,'Color',reversalLineColor,'LineStyle',lineStyle);
330+
plot(speedUpBound,yLim,'LineWidth',lineWidth,'Color',speedUpLineColor,'LineStyle',lineStyle);
331+
set(gca,'FontSize',fontSize)
332+
set(gca,'ylim', yLim);
333+
set(gca,'xlim', xLim);
334+
[xNew yNew] = MiriamAxes(gca,'xy');
335+
336+
%% Plot filtered out points based on thresholds, plot:
337+
338+
figure(7)
339+
markerSize = 8;
340+
341+
maxTrials = max([length(Profiles.Force_50),length(Profiles.Force_100)]);
342+
maxTrials2 = max([length(Profiles.Force_500),length(Profiles.Force_1000),length(Profiles.Force_5000),length(Profiles.Force_10000)]);
343+
for trialPlot = maxTrials2+1:maxTrials
344+
345+
try
346+
if (Profiles.Force_100(trialPlot).maxDeltaVNorm > reversalBound | Profiles.Force_100(trialPlot).maxDeltaVNorm < speedUpBound)
347+
plot(Profiles.Force_100(trialPlot).positionApplied,Profiles.Force_100(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_100,'MarkerEdgeColor','none');
348+
end
349+
end
350+
hold on
351+
352+
try
353+
if (Profiles.Force_50(trialPlot).maxDeltaVNorm > reversalBound | Profiles.Force_50(trialPlot).maxDeltaVNorm < speedUpBound)
354+
plot(Profiles.Force_50(trialPlot).positionApplied,Profiles.Force_50(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_50,'MarkerEdgeColor','none');
355+
end
356+
end
357+
hold on
358+
end
359+
360+
for trialPlot = 1:maxTrials2
361+
362+
try
363+
if (Profiles.Force_100(trialPlot).maxDeltaVNorm > reversalBound | Profiles.Force_100(trialPlot).maxDeltaVNorm < speedUpBound)
364+
plot(Profiles.Force_100(trialPlot).positionApplied,Profiles.Force_100(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_100,'MarkerEdgeColor','none');
365+
end
366+
end
367+
hold on
368+
369+
370+
try
371+
if (Profiles.Force_50(trialPlot).maxDeltaVNorm > reversalBound | Profiles.Force_50(trialPlot).maxDeltaVNorm < speedUpBound)
372+
plot(Profiles.Force_50(trialPlot).positionApplied,Profiles.Force_50(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_50,'MarkerEdgeColor','none');
373+
end
374+
end
375+
hold on
376+
377+
try
378+
if (Profiles.Force_10000(trialPlot).maxDeltaVNorm > reversalBound | Profiles.Force_10000(trialPlot).maxDeltaVNorm < speedUpBound)
379+
plot(Profiles.Force_10000(trialPlot).positionApplied,Profiles.Force_10000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_10000,'MarkerEdgeColor','none');
380+
end
381+
end
382+
hold on
383+
384+
385+
try
386+
if (Profiles.Force_5000(trialPlot).maxDeltaVNorm > reversalBound | Profiles.Force_5000(trialPlot).maxDeltaVNorm < speedUpBound)
387+
plot(Profiles.Force_5000(trialPlot).positionApplied,Profiles.Force_5000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_5000,'MarkerEdgeColor','none');
388+
end
389+
end
390+
hold on
391+
392+
try
393+
if (Profiles.Force_1000(trialPlot).maxDeltaVNorm > reversalBound | Profiles.Force_1000(trialPlot).maxDeltaVNorm < speedUpBound)
394+
plot(Profiles.Force_1000(trialPlot).positionApplied,Profiles.Force_1000(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_1000,'MarkerEdgeColor','none');
395+
end
396+
end
397+
hold on
398+
399+
try
400+
if (Profiles.Force_500(trialPlot).maxDeltaVNorm > reversalBound | Profiles.Force_500(trialPlot).maxDeltaVNorm < speedUpBound)
401+
plot(Profiles.Force_500(trialPlot).positionApplied,Profiles.Force_500(trialPlot).maxDeltaVNorm,'Marker','o','MarkerSize',markerSize,'MarkerFaceColor', color_500,'MarkerEdgeColor','none');
402+
end
403+
end
404+
hold on
405+
406+
end
407+
408+
set(gca,'FontSize',fontSize)
409+
set(gca,'ylim', [-6 12]);
410+
set(gca,'xlim', [0 100]);
411+
[xNew yNew] = MiriamAxes(gca,'xy');

0 commit comments

Comments
 (0)