def test_time_data_as_fixed():
d = data.time_data_file(fixed)
d.read(usecols,skiprows=3)
fig = figure(1)
d.plot_input(fig=fig,clear=False,linelabel='Data Input')
d.plot_output(fig=fig,clear=False,linelabel='Data Output')
freq_data = d.experimental_transferfunction(fixed_sine_freq)
fig = figure(2)
freq_data.plot_bode(fig=fig,linelabel='Fixed Freq Bode')
def test_time_data_as_step():
step_d = data.time_data_file(step)
step_d.read(usecols,skiprows=skiprows)
u = step_d.input
y = step_d.output
u = u*5.0
y = y*5.0
step_d.scale_input(5.0)
assert step_d.input.all() == u.all()
assert step_d.output.all() == y.all()
#A = 97385.3465539
#p = 7636.27078234
#c = 94820.4963045
#z = 54.3639192917
num_dict = {'0':'A'}
den_dict = {'2':'1','1':'z','0':'p'}
var_dict = {'A':A,'p':p,'z':z,'B':B}
to_opt_dict = {'A':A,'p':p,'z':z,'B':B}
############################
usecols = (0,2,4)
datafile = 'step.txt'
data = time_data_file(datafile)
data.read(usecols,skiprows=3)
#data.resample(2)
#data.scale_t(1.0/250)
#data.scale_output(1.0/128/4*2*pi)
##############################
# remove steady state stuff
inds = data.t<.7
new_data = time_data(data.t[inds],data.input[inds],data.output[inds])
###############################
model = Model(num_dict,den_dict,var_dict,to_opt_dict)
