-
Notifications
You must be signed in to change notification settings - Fork 0
/
main_dI1.m
65 lines (55 loc) · 1.39 KB
/
main_dI1.m
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
64
65
main_init;
addpath('C:\Users\bai_f\Documents\research\Axon_growth\DelayedFeedback\mine_dde23');
%% Initial values of bifurcation parameters
% tau = N/nu/(1-rho_bulk)
% N==1, nu==1, rho_bulk==1/2, then we have the following initial values.
rhobulkK = 0.5;
rhobulkD = 0.5;
pE = 6;
wE = 5;
pI = 6;
wI = 5;
vK = 1;
vD = 1;
dE1 = 1;
dE2 = 1;
dI1 = 1.7;
dI2 = 1;
contpar=ind_dI1;
namepar='dI1';
N = 10;
% file_id can be used to generate a new file name to save data.
file_id=4;
% no_y specifies which signal will be set as y-axis in the bifurcation diagram.
% E1, E2, I1, I2 -> 1,2,3,4
no_y = 3;
% label_y specify the ylabel in the bifurcation diagram.
label_y = 'I1';
% min max and step for extending 1par branch
min_bound = 0.1;
max_bound = 2;
max_step = 0.002;
% number of attempts when extending 1par branch
bf_attempts = 4000;
% time instants for extracting frequencies
t_start = 500;
t_end = 1000;
tspan = [0, t_end];
dt = 0.1;
% Need to refine near dI1=0
t_start_near0 = 8000;
t_end_near0 = 10000;
tspan_near0 = [0, t_end_near0];
dt_near0 = 0.1;
%% Setup the first fixed point
main_setup_first_fixed_pt;
%% Initialize the branch
main_initialize_branch;
%% Continuation of the branch
main_cont_branch;
%% Plot bifurcation diagram
main_plot_bf;
%% Extract frequency and amplitude of oscillation
main_extract_freq_amp;
%% output
main_save_data_1par;