rd-model/simulateSpatialTuning.m at master · m3project/rd-model · 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
function varargout = simulateSpatialTuning(args)

tf = 8;

sep = 1; acAngle = 2.5; % degs

duration = 0.5;

sfs = logspace2(1e-3, 1, 100);

[hpf, lpf] = getMantisTemporalFilters();

%% load overrides

if nargin, unpackStruct(args); end

%% body

makePlot = ~nargout;

conds = createTrial(tf, sfs);

args = struct('conds', conds, 'sep', sep, 'acAngle', acAngle, ...
    'hpf', hpf, 'lpf', lpf, 'duration', duration);

ys = simulateMultiSinusoidal(args);

steadyState = mean(ys, 2);

steadyState = steadyState / max(steadyState);

[~, ind] = max(steadyState);

fpref = sfs(ind);

%% output

if makePlot

    clf

    plot(sfs, steadyState);

    xlabel('Spatial Frequency (cpd)');

    title('RD Output');

    set(gca, 'xscale', 'log');

else

    varargout{1} = sfs;

    varargout{2} = steadyState;

    varargout{3} = fpref;

end

end