NV_Matlab_Code/DrawSequence.m at main · cliulinnaeus/NV_Matlab_Code · 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
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
function DrawSequence(SEQ, hObject, eventdata, ax)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%% Written by Jeronimo Maze, July 2007 %%%%%%%%%%%%%%%%%%
%%%%%%%%%% Harvard University, Cambridge, USA  %%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

global tmin tmax ScaleT bShowTimes bShowTypes ScaleStr Color;

if isempty(bShowTypes)
    bShowTypes = false;
end
if isempty(bShowTimes)
    bShowTimes = false;
end


Color = {'g','r','b','k','m','c','y','g','r','b','k','m','c','y'};

if numel(SEQ.CHN)== 0
    return
end
cla(ax,'reset')
for ichn = 1 :numel(SEQ.CHN)
    atmin(ichn) = min(SEQ.CHN(ichn).T);
    atmax(ichn) = max(SEQ.CHN(ichn).T + SEQ.CHN(ichn).DT);
end
tmin = min(atmin);
tmax = max(atmax);

[ScaleT ScaleStr] = GetScale(tmax);

ymin = -0.5;
ymax = (size(SEQ.CHN,2) -1)*1.5 + 1 + 0.5;

%axes(ax); cla;
for ichn = size(SEQ.CHN,2):-1:1
    PlotCHN(ax,SEQ,ichn,Color, tmin,tmax,ScaleT,ScaleStr);
end

xlim(ax,ScaleT*[tmin tmax]);
ylim(ax,[ymin ymax]);

function PlotCHN(ax,SEQ,ichn,Color,tmin,tmax,ScaleT, ScaleStr)
global bShowTimes bShowTypes;

yLow = (ichn -1)*1.5;
yHigh = (ichn -1)*1.5 + 1;
s = 0:1/99:1;
one = ones(size(s));
t0 = tmin;
hold(ax,'on');
for irise = 1:SEQ.CHN(ichn).NRise %PLOT each rise
    t1 = SEQ.CHN(ichn).T(irise);
    dt = SEQ.CHN(ichn).DT(irise);
    t2 = t1 + dt;
    xL = t0 + (t1-t0)*s;
    yL = yLow*one;
    xLH = t1*one;
    yLH = yLow + (yHigh - yLow)*s;
    xH = t1 + (t2-t1)*s;
    yH = yHigh*one;
    xHL = t2*one;
    yHL = yHigh + (yLow - yHigh)*s;
    plot(ax,ScaleT*[xL xLH xH xHL],[yL yLH yH yHL],Color{ichn});
    t0 = t2;
    if bShowTimes
        text(ax,ScaleT*xLH,(yL+yH)/2,...
            [' (' NiceNotation(xHL(1)-xLH(1)) ', ' NiceNotation(xLH(1)) ')'],'FontSize',8);
    end
    if bShowTypes
        text(ax,ScaleT*xLH,yH+0.2, SEQ.CHN(ichn).Type(irise),'FontSize',8);
    end
end
xL = t2 + (tmax-t2)*s;
yL = yLow*one;
plot(ax, ScaleT*xL,yL,Color{ichn});
hold(ax,'off');
text(ax, ScaleT*(tmin+0.01*(tmax-tmin)),yLow+0.5,sprintf('PB%.0f',SEQ.CHN(ichn).PBN),'Color',Color{ichn});
xlabel(ax, ScaleStr)
set(ax,'FontSize',8);

function [ScaleT, ScaleStr] = GetScale(tmax)
if tmax > 0 && tmax <= 100e-3
    ScaleT = 1e3;
    ScaleStr = 'ps';
elseif tmax > 100e-3 && tmax <= 100
    ScaleT = 1;
    ScaleStr = 'ns';
elseif tmax > 100 && tmax <= 100e3
    ScaleT = 1e-3;
    ScaleStr = '{\mu}s';
elseif tmax > 100e3 && tmax <= 100e6
    ScaleT = 1e-6;
    ScaleStr = 'ms';
elseif tmax > 100e6 && tmax <= 1000e11
    ScaleT = 1e-9;
    ScaleStr = 's';
end