simulate201819/PlotAnalyticalModel_v1.m at actinSpeed · SamMachariaPhD/simulate201819 · 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
clear all

Vmax = 7;           % Maximum gliding speed (um/s)
k = 300.0;          % Spring constant (pN/um)
fstall = 9.2;       % Stall force (pN)
frupt = 9.2;        % Rupture force (pN)

Tau1 = 0.025;               % Binding period (s). 1/k_a See Ishigure & Nitta 2015 IEEE Trans. Nanobiosci.
Tau2min = frupt/k/Vmax;

ActiveMotorRatioTemp = 0.01:0.01:0.99;

Tau = Tau1/Tau2min;
f = frupt/fstall;
Gamma = 1.0./ActiveMotorRatioTemp - 1.0;

GammaCritical = ((1.0 + Tau)^2)/2.0/f/Tau
ActiveMotorRatioCritical = 1.0/(GammaCritical + 1.0)

Fid=fopen('Summary.txt','w');
fprintf(Fid,'Critical gamma = %f; Critical ActiveMotorRatio = %f',GammaCritical,ActiveMotorRatioCritical);
fclose(Fid);

Counter = 0;
for I=1:length(Gamma)
    if ActiveMotorRatioTemp(I) > ActiveMotorRatioCritical
        Counter = Counter + 1;
        ActiveMotorRatio(Counter) = ActiveMotorRatioTemp(I);
        Vplus(Counter) = 0.5*Vmax/Tau*(Tau - 1.0 + sqrt((1.0 - Tau)^2 + 4.0*Tau - 2.0*f*Tau.*Gamma(I)));
        Vminus(Counter) = 0.5*Vmax/Tau*(Tau - 1.0 - sqrt((1.0 - Tau)^2 + 4.0*Tau - 2.0*f*Tau.*Gamma(I)));
    end
end


plot(ActiveMotorRatio,Vplus,'r-',ActiveMotorRatio,Vminus,'b.',ActiveMotorRatioCritical,0.0,'ro')
xlabel('Active motor ratio');ylabel('v (um/s)');
axis([0.0 1.0 0.0 8.0])
saveas(gcf,'r-v.fig');
saveas(gcf,'r-v.png');

OutputData = [ActiveMotorRatio;Vplus];
OutputData = OutputData.';
save('r-v.txt','OutputData','-ascii')