LinDiscStimulusEnsemble.m

classdef LinDiscStimulusEnsemble < DiscreteStimulusEnsemble
    
    properties
        circular = false;
    end
    
    methods
        
        function obj = LinDiscStimulusEnsemble(varargin)
            % superclass constructor
            obj = obj@DiscreteStimulusEnsemble();
            
            switch nargin
            case 0
                % do nothing
                
            case 1
                obj.ensemble = double(varargin{1});
                obj.pS = 1.0 ./ double(obj.n) .* ones(1, obj.n);
                
            case 2
                obj.ensemble = double(varargin{1});
                obj.pS = double(varargin{2});
                
                assert(all(size(obj.ensemble) == size(obj.pS)), 'Ensemble value and probability arrays must be the same size')
                
            case 4
                bottom = double(varargin{1});
                top = double(varargin{2});
                interval = double(varargin{3});
                shift = double(varargin{4});

                assert(top > bottom, 'Upper limit must be greater than lower limit')
                assert(interval > 0 && interval <= (top - bottom), 'Invalid interval value')
                assert(mod(top - bottom, interval) == 0, 'Interval must be a divisor of the range')
                assert(shift >= 0 && shift < interval, 'Invalid shift value')
                
                if shift == 0.0
                    obj.ensemble = bottom : interval : (top - interval);
                else
                    obj.ensemble = (bottom + shift) : interval : top;
                end
                
                obj.pS = 1.0 ./ double(obj.n) .* ones(1, obj.n);

            otherwise
                error('Wrong number of arguments')
            end                

        end
        
	end
end