siffer/run_simple_escape.m at master · mmchenry/siffer · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
function run_simple_escape
% Wrapper to run a single simulation of larval escape response, without
% suction feeding from a predator.


%% Define paths

% Root directory for saving and loading simulation data
% (set up to deal with Matt's 2 computers)
if isdir('/Users/mmchenry/Dropbox/Projects/Holzmann/sims/danio')
    sim_path = '/Users/mmchenry/Dropbox/Projects/Holzmann/sims/danio';
else
    sim_path = '/Volumes/Flow HD/Dropbox/Projects/Holzmann/sims/danio';
end


%% Set parameter values

% Default parameter values
[sim,pred,prey] = params_danio;


%% Generate flow field

% Create flow field structure
fl = flow_noflow(sim,pred,x2,y2,Field1);


%% Run simulation

% Run simulation
tic
r = siffer(sim,prey,fl);

% Animate simulation results
vis_sim(fl,sim,prey,r)

% Report time
t_lapse = toc;
disp(' ')
disp(['Sim completed in ' num2str(t_lapse) 's'])

% Visualize results
figure;

subplot(2,1,1)
plot(r.t,1000.*r.pos(1,:),'k-',r.t,1000.*r.pred_pos(1,:),'r--')
xlabel('time (s)');ylabel('position (mm)');
legend('prey','pred','Location','West')

if r.cap
    title('Captured!');
else
    title('Escaped!');
end

subplot(2,1,2)
plot(r.t,r.D(1,:),'r-',r.t,r.PF(1,:),'b-',r.t,r.AR(1,:),'g-')
legend('D','PF','AR','Location','SouthWest')
xlabel('time (s)');ylabel('Force (N)');


end