Python doublet示例

示例#1

文件： test_inputs_generator.py 项目： pablobatun/PyFME

def test_doublet_wrong_not_scalar_offset():
    t_init = 0
    T = 5.
    A = 1.5
    time = np.linspace(0, 10, 11)
    var = np.ones_like(time)
    offset = var

    with pytest.raises(TypeError) as excinfo:
        doublet(t_init, T, A, time, offset=offset)

示例#2

文件： test_inputs_generator.py 项目： AlexS12/PyFME

def test_doublet_wrong_not_scalar_offset():
    t_init = 0
    T = 5.
    A = 1.5
    time = np.linspace(0, 10, 11)
    var = np.ones_like(time)
    offset = var

    with pytest.raises(TypeError) as excinfo:
        doublet(t_init, T, A, time, offset=offset)

示例#3

文件： test_inputs_generator.py 项目： AlexS12/PyFME

def test_ramp_wrong_size_var():
    t_init = 0
    T = 5.
    A = 1.5
    time = np.linspace(0, 10, 11)
    var = np.ones([10])

    with pytest.raises(ValueError) as excinfo:
        doublet(t_init, T, A, time, offset=0, var=var)
    assert ("ValueError: var and time must have the same size"
                in excinfo.exconly())

示例#4

文件： test_inputs_generator.py 项目： pablobatun/PyFME

def test_ramp_wrong_size_var():
    t_init = 0
    T = 5.
    A = 1.5
    time = np.linspace(0, 10, 11)
    var = np.ones([10])

    with pytest.raises(ValueError) as excinfo:
        doublet(t_init, T, A, time, offset=0, var=var)
    assert ("ValueError: var and time must have the same size"
            in excinfo.exconly())

示例#5

文件： test_inputs_generator.py 项目： AlexS12/PyFME

def test_doublet_bounds_not_included():
    t_init = 0.1
    T = 4.8
    A = 3.
    time = np.linspace(0, 10, 11)

    expected_input = np.zeros([11])
    expected_input[1:3] = A/2
    expected_input[3:5] = -A/2

    doublet_input = doublet(t_init, T, A, time, offset=0, var=None)

    assert_almost_equal(doublet_input, expected_input)

示例#6

文件： test_inputs_generator.py 项目： pablobatun/PyFME

def test_doublet():
    t_init = 0
    T = 5.
    A = 3.
    time = np.linspace(0, 10, 11)

    expected_input = np.zeros([11])
    expected_input[0:3] = A / 2
    expected_input[3:6] = -A / 2

    doublet_input = doublet(t_init, T, A, time, offset=0, var=None)

    assert_almost_equal(doublet_input, expected_input)

示例#7

文件： test_inputs_generator.py 项目： pablobatun/PyFME

def test_doublet_bounds_not_included():
    t_init = 0.1
    T = 4.8
    A = 3.
    time = np.linspace(0, 10, 11)

    expected_input = np.zeros([11])
    expected_input[1:3] = A / 2
    expected_input[3:5] = -A / 2

    doublet_input = doublet(t_init, T, A, time, offset=0, var=None)

    assert_almost_equal(doublet_input, expected_input)

示例#8

文件： test_inputs_generator.py 项目： AlexS12/PyFME

def test_doublet():
    t_init = 0
    T = 5.
    A = 3.
    time = np.linspace(0, 10, 11)

    expected_input = np.zeros([11])
    expected_input[0:3] = A/2
    expected_input[3:6] = -A/2

    doublet_input = doublet(t_init, T, A, time, offset=0, var=None)

    assert_almost_equal(doublet_input, expected_input)

示例#9

文件： test_inputs_generator.py 项目： AlexS12/PyFME

def test_doublet_offset():
    t_init = 0
    T = 5
    A = 1.5
    time = np.linspace(0, 10, 11)
    offset = 2.6

    expected_input = np.zeros([11]) + offset
    expected_input[0:3] += A/2
    expected_input[3:6] += -A/2

    doublet_input = doublet(t_init, T, A, time, offset=offset, var=None)

    assert_almost_equal(doublet_input, expected_input)

示例#10

文件： test_inputs_generator.py 项目： pablobatun/PyFME

def test_doublet_offset():
    t_init = 0
    T = 5
    A = 1.5
    time = np.linspace(0, 10, 11)
    offset = 2.6

    expected_input = np.zeros([11]) + offset
    expected_input[0:3] += A / 2
    expected_input[3:6] += -A / 2

    doublet_input = doublet(t_init, T, A, time, offset=offset, var=None)

    assert_almost_equal(doublet_input, expected_input)

示例#11

文件： test_inputs_generator.py 项目： pablobatun/PyFME

def test_doublet_var():
    t_init = 0
    T = 5.
    A = 1.5
    time = np.linspace(0, 10, 11)
    var = np.ones_like(time)
    var[0::2] = -1

    expected_input = var.copy()
    expected_input[0:3] += A / 2
    expected_input[3:6] += -A / 2

    doublet_input = doublet(t_init, T, A, time, offset=0, var=var)

    assert_almost_equal(doublet_input, expected_input)

示例#12

文件： test_inputs_generator.py 项目： AlexS12/PyFME

def test_doublet_var():
    t_init = 0
    T = 5.
    A = 1.5
    time = np.linspace(0, 10, 11)
    var = np.ones_like(time)
    var[0::2] = -1

    expected_input = var.copy()
    expected_input[0:3] += A/2
    expected_input[3:6] += -A/2

    doublet_input = doublet(t_init, T, A, time, offset=0, var=var)

    assert_almost_equal(doublet_input, expected_input)

示例#13

文件： test_inputs_generator.py 项目： AeroPython/PyFME

def test_doublet_var_and_offset():
    t_init = 0
    T = 5.
    A = 1.5
    time = np.linspace(0, 10, 11)
    offset = 2.6
    var = np.ones_like(time)
    var[0::2] = -1

    expected_input = var
    expected_input[0:3] += A/2 + offset
    expected_input[3:6] += -A/2 + offset

    doublet_input = doublet(t_init, T, A, time, offset=offset, var=var)

    assert_almost_equal(doublet_input, expected_input)

示例#14

def test_doublet_var_and_offset():
    t_init = 0
    T = 5.
    A = 1.5
    time = np.linspace(0, 10, 11)
    offset = 2.6
    var = np.ones_like(time)
    var[0::2] = -1

    expected_input = var
    expected_input[0:3] += A / 2 + offset
    expected_input[3:6] += -A / 2 + offset

    doublet_input = doublet(t_init, T, A, time, offset=offset, var=var)

    assert_almost_equal(doublet_input, expected_input)

示例#15

文件： example_008_ailerons_doublet.py 项目： AeroPython/PyFME

tfin = 20  # seconds
N = tfin * 100 + 1
time = np.linspace(0, tfin, N)
initial_controls = trimmed_ac.controls

controls = {}
for control_name, control_value in initial_controls.items():
    controls[control_name] = np.ones_like(time) * control_value

# Aileron doublet
# Aileron max travel: +20º/-15º
amplitude = np.deg2rad(15)
controls['delta_aileron'] = doublet(t_init=2,
                                    T=1,
                                    A=amplitude,
                                    time=time,
                                    offset=np.deg2rad(0))

my_simulation.set_controls(time, controls)

par_list = ['x_earth', 'y_earth', 'height',
            'psi', 'theta', 'phi',
            'u', 'v', 'w',
            'v_north', 'v_east', 'v_down',
            'p', 'q', 'r',
            'alpha', 'beta', 'TAS',
            'F_xb', 'F_yb', 'F_zb',
            'M_xb', 'M_yb', 'M_zb']

my_simulation.set_par_dict(par_list)

示例#16

tfin = 10  # seconds
N = tfin * 100 + 1
time = np.linspace(0, tfin, N)
initial_controls = trimmed_ac.controls

controls = {}
for control_name, control_value in initial_controls.items():
    controls[control_name] = np.ones_like(time) * control_value

# Elevator doublet
# Elevator max travel: +28º/-26º
amplitude = np.deg2rad(20)
controls['delta_elevator'] = doublet(t_init=2,
                                     T=1,
                                     A=amplitude,
                                     time=time,
                                     offset=initial_controls['delta_elevator'])

my_simulation.set_controls(time, controls)

par_list = ['x_earth', 'y_earth', 'height',
            'psi', 'theta', 'phi',
            'u', 'v', 'w',
            'v_north', 'v_east', 'v_down',
            'p', 'q', 'r',
            'alpha', 'beta', 'TAS',
            'F_xb', 'F_yb', 'F_zb',
            'M_xb', 'M_yb', 'M_zb']

my_simulation.set_par_dict(par_list)

示例#17

文件： example_008_ailerons_doublet.py 项目： Python3pkg/PyFME

tfin = 20  # seconds
N = tfin * 100 + 1
time = np.linspace(0, tfin, N)
initial_controls = trimmed_ac.controls

controls = {}
for control_name, control_value in list(initial_controls.items()):
    controls[control_name] = np.ones_like(time) * control_value

# Aileron doublet
# Aileron max travel: +20º/-15º
amplitude = np.deg2rad(15)
controls['delta_aileron'] = doublet(t_init=2,
                                    T=1,
                                    A=amplitude,
                                    time=time,
                                    offset=np.deg2rad(0))

my_simulation.set_controls(time, controls)

par_list = ['x_earth', 'y_earth', 'height',
            'psi', 'theta', 'phi',
            'u', 'v', 'w',
            'v_north', 'v_east', 'v_down',
            'p', 'q', 'r',
            'alpha', 'beta', 'TAS',
            'F_xb', 'F_yb', 'F_zb',
            'M_xb', 'M_yb', 'M_zb']

my_simulation.set_par_dict(par_list)

示例#18

文件： example_007_rudder_doublet.py 项目： AeroPython/PyFME

tfin = 20  # seconds
N = tfin * 100 + 1
time = np.linspace(0, tfin, N)
initial_controls = trimmed_ac.controls

controls = {}
for control_name, control_value in initial_controls.items():
    controls[control_name] = np.ones_like(time) * control_value

# Rudder doublet
# Rudder max travel: +16º/-16º
amplitude = np.deg2rad(32)

controls['delta_rudder'] = doublet(t_init=2,
                                   T=2,
                                   A=amplitude,
                                   time=time,
                                   offset=np.deg2rad(0.0))

my_simulation.set_controls(time, controls)

par_list = ['x_earth', 'y_earth', 'height',
            'psi', 'theta', 'phi',
            'u', 'v', 'w',
            'v_north', 'v_east', 'v_down',
            'p', 'q', 'r',
            'alpha', 'beta', 'TAS',
            'F_xb', 'F_yb', 'F_zb',
            'M_xb', 'M_yb', 'M_zb']

my_simulation.set_par_dict(par_list)

示例#19

文件： example_007_rudder_doublet.py 项目： sanguinariojoe/PyFME

tfin = 20  # seconds
N = tfin * 100 + 1
time = np.linspace(0, tfin, N)
initial_controls = trimmed_ac.controls

controls = {}
for control_name, control_value in initial_controls.items():
    controls[control_name] = np.ones_like(time) * control_value

# Rudder doublet
# Rudder max travel: +16º/-16º
amplitude = np.deg2rad(32)

controls['delta_rudder'] = doublet(t_init=2,
                                   T=2,
                                   A=amplitude,
                                   time=time,
                                   offset=np.deg2rad(0.0))

my_simulation.set_controls(time, controls)

par_list = [
    'x_earth', 'y_earth', 'height', 'psi', 'theta', 'phi', 'u', 'v', 'w',
    'v_north', 'v_east', 'v_down', 'p', 'q', 'r', 'alpha', 'beta', 'TAS',
    'F_xb', 'F_yb', 'F_zb', 'M_xb', 'M_yb', 'M_zb'
]

my_simulation.set_par_dict(par_list)
my_simulation.run_simulation()

# print(my_simulation.par_dict)

示例#20

文件： example_009_elev_aile_rud_doublet.py 项目： AeroPython/PyFME

initial_controls = trimmed_ac.controls

controls = {}
for control_name, control_value in initial_controls.items():
    controls[control_name] = np.ones_like(time) * control_value

# Elevator, aileron and rudder doublets
# Elevator max travel: +28º/-26º
# Aileron max travel: +20º/-15º
# Rudder max travel: +16º/-16º

amplitude_elev = np.deg2rad(26)
controls['delta_elevator'] = initial_controls['delta_elevator'] + \
                             doublet(t_init=2,
                                     T=2,
                                     A=amplitude_elev,
                                     time=time,
                                     offset=np.deg2rad(0.0))

amplitude_aile = np.deg2rad(15)
controls['delta_aileron'] = doublet(t_init=5,
                                    T=2,
                                    A=amplitude_aile,
                                    time=time,
                                    offset=np.deg2rad(0.0))

amplitude_rud = np.deg2rad(32)
controls['delta_rudder'] = doublet(t_init=10,
                                   T=2,
                                   A=amplitude_rud,
                                   time=time,

示例#21

文件： example_009_elev_aile_rud_doublet.py 项目： Python3pkg/PyFME

initial_controls = trimmed_ac.controls

controls = {}
for control_name, control_value in list(initial_controls.items()):
    controls[control_name] = np.ones_like(time) * control_value

# Elevator, aileron and rudder doublets
# Elevator max travel: +28º/-26º
# Aileron max travel: +20º/-15º
# Rudder max travel: +16º/-16º

amplitude_elev = np.deg2rad(26)
controls['delta_elevator'] = initial_controls['delta_elevator'] + \
                             doublet(t_init=2,
                                     T=2,
                                     A=amplitude_elev,
                                     time=time,
                                     offset=np.deg2rad(0.0))

amplitude_aile = np.deg2rad(15)
controls['delta_aileron'] = doublet(t_init=5,
                                    T=2,
                                    A=amplitude_aile,
                                    time=time,
                                    offset=np.deg2rad(0.0))

amplitude_rud = np.deg2rad(32)
controls['delta_rudder'] = doublet(t_init=10,
                                   T=2,
                                   A=amplitude_rud,
                                   time=time,