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)