# TODO: add discrete time option

def _rss_generate(states, inputs, outputs, type):

def _rss_generate(states, inputs, outputs, type, strictly_proper=False):

"""Generate a random state space.

This does the actual random state space generation expected from rss and

# Apply masks.

B = B * Bmask

C = C * Cmask

D = D * Dmask

D = D * Dmask if not strictly_proper else zeros(D.shape)

return StateSpace(A, B, C, D)

raise ValueError("Needs 1 or 2 arguments; received %i." % len(args))

def rss(states=1, outputs=1, inputs=1):

def rss(states=1, outputs=1, inputs=1, strictly_proper=False):

"""

Create a stable *continuous* random state space object.

Number of system inputs

outputs : integer

Number of system outputs

strictly_proper : bool, optional

If set to 'True', returns a proper system (no direct term). Default

value is 'False'.

Returns

-------

"""

return _rss_generate(states, inputs, outputs, 'c')

return _rss_generate(states, inputs, outputs, 'c',

strictly_proper=strictly_proper)

def drss(states=1, outputs=1, inputs=1):

def drss(states=1, outputs=1, inputs=1, strictly_proper=False):

"""

Create a stable *discrete* random state space object.

"""

return _rss_generate(states, inputs, outputs, 'd')

return _rss_generate(states, inputs, outputs, 'd',

strictly_proper=strictly_proper)

def ssdata(sys):