Modeling_MS/Draw_Speed.m at master · jpresern/Modeling_MS · 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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Janez Presern, Ales Skorjanc, Tomaz Rodic, Jan Benda 2011-2015
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Function draws dependency of the currents to the stimulus velocity.

%   Function requires:
%       dt..                sampling rate
%       stimamp..           stimulus description (amplitudes changes)
%       stimTime..          stimulus description (time durations of amplitude changes
%       Diagram_x ..        experimentally obtained delay - response curve (velocity
%                           position)
%       Diagram_y ..        experimentally obtained delay - response curve
%                           (currents)
%       var  ..             fitted variables and constants
%       var_names ..        names of variables
%       cmap1,cmap2 ..      color maps
%       fn ..               filename (Results_XX) to identify the model
%                           iteration

%   Function outputs:
%       f             ..    figure handle

function [f, output] = Draw_Speed(dt,stimAmp, stimTime,...
                                Diagram_x, Diagram_y,...
                                var, var_names,...
                                cmap1, cmap2, fn)

f = figure;

s(2) = axes('OuterPosition', [0 0.4 1 0.4]);

hold on;
b = 10;
for a = 1:10
    [~,~,results]=Modeling_DRG_TCM_Engine(horzcat('DRG_TCM_Model_mh','_Report'),...
        stimTime(b,:),stimAmp(b,:),...
        var,var_names,dt);
    p1 = plot(results.t,results.g,'-','Color',cmap2(b,:),'LineWidth',2);
    output.Fig8A.model.Imax(b) = min(results.g);
    output.Fig8A.model.tIMax(b) = results.t(min(results.g)==results.g);
    line ([output.Fig8A.model.tIMax(b) output.Fig8A.model.tIMax(b)],[output.Fig8A.model.Imax(b) 0],'Color','r','LineWidth',1.5);
    b = b-1;
end;
plot (output.Fig8A.model.tIMax, output.Fig8A.model.Imax, '*r','MarkerSize',5);
hold off;
set(gca,'Ylim',[-1.1 0]);
xlim([0 220]);
xlabel ('time [ms]');
ylabel ('I/I_{max}');
grid on;

s(1) = axes('OuterPosition', [0 0.8 1 0.2]);
set(gca,'XTickLabel',[]);
title(horzcat(fn,': ','Stimulus speed encoding'));
for val2 = 1:10
    x0 = 0;
    y0 = 0;
    for val1 = 1:5
        x1 = x0 + stimTime(val2,val1);
        y1 = y0 + stimAmp(val2,val1);
        line([x0 x1],[y0 y1],'Color',cmap2(val2,:),'LineWidth',1.5);
        x0 = x1;
        y0 = y1;
    end;
    line ([output.Fig8A.model.tIMax(val2) output.Fig8A.model.tIMax(val2)],[8.5 0],'Color','r','LineWidth',1.5);
end;
xlim ([0 220]);
ylabel ('[stim amplitude]');
set (gca,'YTick',[8.5]);
set (gca,'YTickLabel',[8.5]);


s(3) = axes('OuterPosition', [0 0.0 1 0.3]);

velocity = (stimAmp(:,2)./stimTime(:,2))*1000;
hold on;
plot (velocity,output.Fig8A.model.Imax./min(output.Fig8A.model.Imax), '*r','LineWidth',1.5,'MarkerSize',10);
plot (Diagram_x, Diagram_y, 's', 'LineWidth',1.5, 'MarkerSize',5,'Color',cmap1(14,:));
set(s(3), 'Xdir','reverse')
xlabel ('stimulus up-ramp velocity [\mum/s]');
ylabel ('I/I_{max}');
grid on;
title(horzcat(fn,': ','Fig8C: Stimulus speed encoding in DRG neuron'));