class Agent(object):
def __init__(self, model, env, action_size, args, state):
self.model = model
self.env = env
self.action_size = action_size
self.state = state
self.hx = None
self.cx = None
self.eps_len = 0
self.args = args
self.values = []
self.log_probs = []
self.rewards = []
self.entropies = []
self.done = True # 初始化，可以设置一次新的状态值
self.info = None
self.reward = 0
self.gpu_id = -1
# 参数
self.memory = Experience(history_size=2000)
def fill_experience(self):
prev_state = self.env.last_state
last_action = self.env.last_action
last_reward = self.env.last_reward
last_action_reward = ExperienceFrame.concat_action_and_reward(
last_action, self.action_size, last_reward)
with torch.no_grad():
state = torch.from_numpy(self.env.last_state).unsqueeze(0)
lar = torch.from_numpy(last_action_reward).unsqueeze(0)
_, pi, (self.hx, self.cx) = self.model(task_type='a3c',
states=state,
hx=self.hx,
cx=self.cx,
last_action_rewards=lar)
action_index = pi.max(1)[1].view(1, 1).item()
new_state, reward, terminal, pixel_change = self.env.step(
action_index) # 存储为数组
frame = ExperienceFrame(prev_state, reward, action_index, terminal,
pixel_change, last_action, last_reward)
self.memory.add_frame(frame)
if terminal:
self.env.reset()
if self.memory.is_full():
self.env.reset()
print("Replay buffer filled")
self.done = terminal
def a3c_process(self):
"""
在 on-policy 下运行程序
:return:
"""
states = []
last_action_rewards = []
actions = [] #
rewards = []
values = [] # V
actions_prob = []
terminal_end = False
# t_max times loop
for _ in range(self.args.num_steps):
# Prepare last action reward
last_action = self.env.last_action
last_reward = self.env.last_reward
last_action_reward = ExperienceFrame.concat_action_and_reward(
last_action, self.action_size, last_reward)
state = torch.from_numpy(self.env.last_state).unsqueeze(0)
lar = torch.from_numpy(last_action_reward)
v, pi, (self.hx,
self.cx) = self.model(task_type='a3c',
states=state,
hx=self.hx,
cx=self.cx,
last_action_rewards=lar.unsqueeze(0))
action_index = pi.max(1)[1].view(1, 1).item()
states.append(torch.from_numpy(self.env.last_state))
actions_prob.append(torch.squeeze(pi, dim=0))
last_action_rewards.append(lar)
actions.append(action_index)
values.append(v)
prev_state = self.env.last_state
new_state, reward, terminal, pixel_change = self.env.step(
action_index)
frame = ExperienceFrame(prev_state, reward, action_index, terminal,
pixel_change, last_action, last_reward)
# Store to experience
self.memory.add_frame(frame)
# self.episode_reward += reward
rewards.append(reward)
self.update_lstm_state()
if terminal:
self.env.reset()
break
R = torch.zeros(1, 1)
if not terminal_end:
state = torch.from_numpy(new_state).unsqueeze(0)
lar = torch.from_numpy(frame.get_action_reward(
self.action_size)).unsqueeze(0)
value, _, _ = self.model(task_type='a3c',
states=state,
hx=self.hx,
cx=self.cx,
last_action_rewards=lar)
R = value.data
# 构造误差项
actions.reverse()
rewards.reverse()
values.reverse()
batch_a = []
batch_adv = []
batch_R = []
for (ai, ri, Vi) in zip(actions, rewards, values):
R = ri + self.args.gamma * R
adv = R - Vi
a = np.zeros([self.action_size], dtype=np.float32)
a[ai] = 1.0
batch_a.append(torch.from_numpy(a))
batch_adv.append(adv)
batch_R.append(R)
batch_a.reverse()
batch_adv.reverse()
batch_R.reverse()
# 转换为张量
return batch_a, batch_adv, batch_R, last_action_rewards, states, actions_prob, values
def a3c_loss(self, batch_a, batch_adv, batch_R, last_action_rewards,
states, actions_prob, values):
batch_a = torch.stack(batch_a) # batch, 6
batch_adv = torch.stack(batch_adv) # batch,1,1
last_action_rewards = torch.stack(last_action_rewards) # batch,7
batch_R = torch.stack(batch_R) # batch,1,1
states = torch.stack(states) # batch,3,84,84
actions_prob = torch.stack(actions_prob) # batch,6
values = torch.stack(values)
# 损失函数
log_pi = torch.log(torch.clamp(actions_prob, min=1e-20, max=1.0))
entropy = -torch.sum(log_pi * actions_prob, dim=1)
# 对应的 a_i 的概率
log_pi_a_i = torch.sum(torch.mul(log_pi, batch_a), dim=1)
policy_loss = torch.sum(log_pi_a_i * batch_adv + entropy * 0.001)
# value_loss
value_loss = 0.5 * F.mse_loss(batch_R, values)
return policy_loss + value_loss
def action_train(self):
value, logit, (self.hx, self.cx) = self.model(
(Variable(self.state.unsqueeze(0)), (self.hx, self.cx)))
prob = F.softmax(logit, dim=1)
log_prob = F.log_softmax(logit, dim=1)
entropy = -(log_prob * prob).sum(1)
self.entropies.append(entropy)
action = prob.multinomial(1).data
log_prob = log_prob.gather(1, Variable(action))
state, self.reward, self.done, self.info = self.env.step(
action.cpu().numpy())
self.state = torch.from_numpy(state).float()
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
self.state = self.state.cuda()
self.reward = max(min(self.reward, 1), -1)
self.values.append(value)
self.log_probs.append(log_prob)
self.rewards.append(self.reward)
return self
def action_test(self):
with torch.no_grad():
self.update_lstm_state()
state = torch.from_numpy(self.env.last_state).unsqueeze(0)
last_action = self.env.last_action
last_reward = np.clip(self.env.last_reward, -1, 1)
last_action_reward = ExperienceFrame.concat_action_and_reward(
last_action, self.action_size, last_reward)
lar = torch.from_numpy(last_action_reward)
v, pi, (self.hx,
self.cx) = self.model(task_type='a3c',
states=state,
hx=self.hx,
cx=self.cx,
last_action_rewards=lar.unsqueeze(0))
prob = F.softmax(pi, dim=1)
action = prob.max(1)[1].data.cpu().numpy()
state, self.reward, self.done, pixel_change = self.env.step(action[0])
self.info = 5
self.state = torch.from_numpy(state).float()
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
self.state = self.state.cuda()
self.eps_len += 1
return self
def update_lstm_state(self):
if self.done:
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
self.cx = Variable(torch.zeros(1, 256).cuda())
self.hx = Variable(torch.zeros(1, 256).cuda())
else:
self.cx = Variable(torch.zeros(1, 256))
self.hx = Variable(torch.zeros(1, 256))
else:
self.cx = Variable(self.cx.data)
self.hx = Variable(self.hx.data)
def clear_actions(self):
self.values = []
self.log_probs = []
self.rewards = []
self.entropies = []
return self
class UNREAL(object):
def __init__(self, env, task, visualise):
self.env = env
self.task = task
self.ob_shape = [HEIGHT, WIDTH, CHANNEL]
self.action_n = Environment.get_action_size()
# define the network stored in ps which is used to sync
worker_device = '/job:worker/task:{}'.format(task)
with tf.device(
tf.train.replica_device_setter(1,
worker_device=worker_device)):
with tf.variable_scope('global'):
self.experience = Experience(
EXPERIENCE_HISTORY_SIZE) # exp replay pool
self.network = UnrealModel(self.action_n, self.env,
self.experience)
self.global_step = tf.get_variable('global_step',
dtype=tf.int32,
initializer=tf.constant(
0, dtype=tf.int32),
trainable=False)
# define the local network which is used to calculate the gradient
with tf.device(worker_device):
with tf.variable_scope('local'):
self.local_network = net = UnrealModel(self.action_n, self.env,
self.experience)
net.global_step = self.global_step
# add summaries for losses and norms
self.batch_size = tf.to_float(tf.shape(net.base_input)[0])
base_loss = self.local_network.base_loss
pc_loss = self.local_network.pc_loss
rp_loss = self.local_network.rp_loss
vr_loss = self.local_network.vr_loss
entropy = tf.reduce_sum(self.local_network.entropy)
self.loss = base_loss + pc_loss + rp_loss + vr_loss
grads = tf.gradients(self.loss, net.var_list)
tf.summary.scalar('model/a3c_loss', base_loss / self.batch_size)
tf.summary.scalar('model/pc_loss', pc_loss / self.batch_size)
tf.summary.scalar('model/rp_loss', rp_loss / self.batch_size)
tf.summary.scalar('model/vr_loss', vr_loss / self.batch_size)
tf.summary.scalar('model/grad_global_norm', tf.global_norm(grads))
tf.summary.scalar('model/var_global_norm',
tf.global_norm(net.var_list))
tf.summary.scalar('model/entropy', entropy / self.batch_size)
tf.summary.image('model/state', net.base_input)
self.summary_op = tf.summary.merge_all()
# clip the gradients to avoid gradient explosion
grads, _ = tf.clip_by_global_norm(grads, GRAD_NORM_CLIP)
self.sync = tf.group(*[
v1.assign(v2)
for v1, v2 in zip(net.var_list, self.network.var_list)
])
grads_and_vars = list(zip(grads, self.network.var_list))
inc_step = self.global_step.assign_add(tf.to_int32(
self.batch_size))
lr = log_uniform(LR_LOW, LR_HIGH)
opt = tf.train.RMSPropOptimizer(learning_rate=lr,
decay=RMSP_ALPHA,
momentum=0.0,
epsilon=RMSP_EPSILON)
self.train_op = tf.group(opt.apply_gradients(grads_and_vars),
inc_step)
self.summary_writer = None
self.local_step = 0
def start(self, sess, summary_writer):
self.summary_writer = summary_writer
def fill_experience(self, sess):
"""
Fill experience buffer until buffer is full.
"""
prev_state = self.env.last_state
last_action = self.env.last_action
last_reward = self.env.last_reward
last_action_reward = ExperienceFrame.concat_action_and_reward(
last_action, self.action_n, last_reward)
pi_, _ = self.local_network.run_base_policy_and_value(
sess, prev_state, last_action_reward)
action = choose_action(pi_)
new_state, reward, terminal, pixel_change = self.env.process(action)
frame = ExperienceFrame(prev_state, reward, action, terminal,
pixel_change, last_action, last_reward)
self.experience.add_frame(frame)
if terminal:
self.env.reset()
if self.experience.is_full():
self.env.reset()
print("Replay buffer filled")
def process(self, sess):
sess.run(self.sync) # sync the local network with the global network
if not self.experience.is_full():
self.fill_experience(sess)
return 0
batch_data = self.local_network.get_batch_data(sess)
feed_dict = {
self.local_network.base_input: batch_data[0],
self.local_network.base_last_action_reward_input: batch_data[1],
self.local_network.base_a: batch_data[2],
self.local_network.base_adv: batch_data[3],
self.local_network.base_r: batch_data[4],
self.local_network.base_initial_lstm_state: batch_data[5],
self.local_network.pc_input: batch_data[6],
self.local_network.pc_last_action_reward_input: batch_data[7],
self.local_network.pc_a: batch_data[8],
self.local_network.pc_r: batch_data[9],
self.local_network.vr_input: batch_data[10],
self.local_network.vr_last_action_reward_input: batch_data[11],
self.local_network.vr_r: batch_data[12],
self.local_network.rp_input: batch_data[13],
self.local_network.rp_r: batch_data[14]
}
should_summarize = self.task == 0 and self.local_step % 11 == 0
if should_summarize:
fetches = [self.summary_op, self.train_op, self.global_step]
else:
fetches = [self.train_op, self.global_step]
fetched = sess.run(fetches, feed_dict=feed_dict)
if should_summarize:
self.summary_writer.add_summary(tf.Summary.FromString(fetched[0]),
fetched[-1])
self.summary_writer.flush()
self.local_step += 1
class Trainer(object):
def __init__(self, thread_index, global_network, initial_learning_rate,
learning_rate_input, grad_applier, env_type, env_name,
use_pixel_change, use_value_replay, use_reward_prediction,
pixel_change_lambda, entropy_beta, local_t_max, gamma,
gamma_pc, experience_history_size, max_global_time_step,
device):
self.thread_index = thread_index
self.learning_rate_input = learning_rate_input
self.env_type = env_type
self.env_name = env_name
self.use_pixel_change = use_pixel_change
self.use_value_replay = use_value_replay
self.use_reward_prediction = use_reward_prediction
self.local_t_max = local_t_max
self.gamma = gamma
self.gamma_pc = gamma_pc
self.experience_history_size = experience_history_size
self.max_global_time_step = max_global_time_step
self.action_size = Environment.get_action_size(env_type, env_name)
self.local_network = UnrealModel(self.action_size, thread_index,
use_pixel_change, use_value_replay,
use_reward_prediction,
pixel_change_lambda, entropy_beta,
device)
self.local_network.prepare_loss()
self.apply_gradients = grad_applier.minimize_local(
self.local_network.total_loss, global_network.get_vars(),
self.local_network.get_vars())
self.sync = self.local_network.sync_from(global_network)
self.experience = Experience(self.experience_history_size)
self.local_t = 0
self.initial_learning_rate = initial_learning_rate
self.episode_reward = 0
# For log output
self.prev_local_t = 0
def prepare(self):
print('')
print('trainer creating env...')
print('')
self.environment = Environment.create_environment(
self.env_type, self.env_name)
def stop(self):
self.environment.stop()
def _anneal_learning_rate(self, global_time_step):
learning_rate = self.initial_learning_rate * (
self.max_global_time_step -
global_time_step) / self.max_global_time_step
if learning_rate < 0.0:
learning_rate = 0.0
return learning_rate
def choose_action(self, pi_values):
return np.random.choice(range(len(pi_values)), p=pi_values)
def _record_score(self, sess, summary_writer, summary_op, score_input,
score, global_t):
summary_str = sess.run(summary_op, feed_dict={score_input: score})
summary_writer.add_summary(summary_str, global_t)
summary_writer.flush()
def set_start_time(self, start_time):
self.start_time = start_time
def _fill_experience(self, sess):
"""
Fill experience buffer until buffer is full.
"""
prev_state = self.environment.last_state
last_action = self.environment.last_action
last_reward = self.environment.last_reward
last_action_reward = ExperienceFrame.concat_action_and_reward(
last_action, self.action_size, last_reward)
pi_, _ = self.local_network.run_base_policy_and_value(
sess, self.environment.last_state, last_action_reward)
action = self.choose_action(pi_)
new_state, reward, terminal, pixel_change = self.environment.process(
action)
#print('action:', action, terminal)
frame = ExperienceFrame(prev_state, reward, action, terminal,
pixel_change, last_action, last_reward)
self.experience.add_frame(frame)
if terminal:
self.environment.reset()
if self.experience.is_full():
self.environment.reset()
print("Replay buffer filled")
def _print_log(self, global_t):
if (self.thread_index == 0) and (self.local_t - self.prev_local_t >=
PERFORMANCE_LOG_INTERVAL):
self.prev_local_t += PERFORMANCE_LOG_INTERVAL
elapsed_time = time.time() - self.start_time
steps_per_sec = global_t / elapsed_time
print(
"### Performance : {} STEPS in {:.0f} sec. {:.0f} STEPS/sec. {:.2f}M STEPS/hour"
.format(global_t, elapsed_time, steps_per_sec,
steps_per_sec * 3600 / 1000000.))
def _process_base(self, sess, global_t, summary_writer, summary_op,
score_input):
# [Base A3C]
states = []
last_action_rewards = []
actions = []
rewards = []
values = []
terminal_end = False
start_lstm_state = self.local_network.base_lstm_state_out
# t_max times loop
for _ in range(self.local_t_max):
# Prepare last action reward
last_action = self.environment.last_action
last_reward = self.environment.last_reward
last_action_reward = ExperienceFrame.concat_action_and_reward(
last_action, self.action_size, last_reward)
pi_, value_ = self.local_network.run_base_policy_and_value(
sess, self.environment.last_state, last_action_reward)
action = self.choose_action(pi_)
states.append(self.environment.last_state)
last_action_rewards.append(last_action_reward)
actions.append(action)
values.append(value_)
if (self.thread_index == 0) and (self.local_t % LOG_INTERVAL == 0):
print("pi={}".format(pi_))
print(" V={}".format(value_))
prev_state = self.environment.last_state
# Process game
new_state, reward, terminal, pixel_change = self.environment.process(
action)
frame = ExperienceFrame(prev_state, reward, action, terminal,
pixel_change, last_action, last_reward)
# Store to experience
self.experience.add_frame(frame)
self.episode_reward += reward
rewards.append(reward)
self.local_t += 1
if terminal:
terminal_end = True
print("score={}".format(self.episode_reward))
self._record_score(sess, summary_writer, summary_op,
score_input, self.episode_reward, global_t)
self.episode_reward = 0
self.environment.reset()
self.local_network.reset_state()
break
R = 0.0
if not terminal_end:
R = self.local_network.run_base_value(
sess, new_state, frame.get_action_reward(self.action_size))
actions.reverse()
states.reverse()
rewards.reverse()
values.reverse()
batch_si = []
batch_a = []
batch_adv = []
batch_R = []
for (ai, ri, si, Vi) in zip(actions, rewards, states, values):
R = ri + self.gamma * R
adv = R - Vi
a = np.zeros([self.action_size])
a[ai] = 1.0
batch_si.append(si)
batch_a.append(a)
batch_adv.append(adv)
batch_R.append(R)
batch_si.reverse()
batch_a.reverse()
batch_adv.reverse()
batch_R.reverse()
return batch_si, last_action_rewards, batch_a, batch_adv, batch_R, start_lstm_state
def _process_pc(self, sess):
# [pixel change]
# Sample 20+1 frame (+1 for last next state)
pc_experience_frames = self.experience.sample_sequence(
self.local_t_max + 1)
# Reverse sequence to calculate from the last
pc_experience_frames.reverse()
batch_pc_si = []
batch_pc_a = []
batch_pc_R = []
batch_pc_last_action_reward = []
pc_R = np.zeros([20, 20], dtype=np.float32)
if not pc_experience_frames[1].terminal:
pc_R = self.local_network.run_pc_q_max(
sess, pc_experience_frames[0].state,
pc_experience_frames[0].get_last_action_reward(
self.action_size))
for frame in pc_experience_frames[1:]:
pc_R = frame.pixel_change + self.gamma_pc * pc_R
a = np.zeros([self.action_size])
a[frame.action] = 1.0
last_action_reward = frame.get_last_action_reward(self.action_size)
batch_pc_si.append(frame.state)
batch_pc_a.append(a)
batch_pc_R.append(pc_R)
batch_pc_last_action_reward.append(last_action_reward)
batch_pc_si.reverse()
batch_pc_a.reverse()
batch_pc_R.reverse()
batch_pc_last_action_reward.reverse()
return batch_pc_si, batch_pc_last_action_reward, batch_pc_a, batch_pc_R
def _process_vr(self, sess):
# [Value replay]
# Sample 20+1 frame (+1 for last next state)
vr_experience_frames = self.experience.sample_sequence(
self.local_t_max + 1)
# Reverse sequence to calculate from the last
vr_experience_frames.reverse()
batch_vr_si = []
batch_vr_R = []
batch_vr_last_action_reward = []
vr_R = 0.0
if not vr_experience_frames[1].terminal:
vr_R = self.local_network.run_vr_value(
sess, vr_experience_frames[0].state,
vr_experience_frames[0].get_last_action_reward(
self.action_size))
# t_max times loop
for frame in vr_experience_frames[1:]:
vr_R = frame.reward + self.gamma * vr_R
batch_vr_si.append(frame.state)
batch_vr_R.append(vr_R)
last_action_reward = frame.get_last_action_reward(self.action_size)
batch_vr_last_action_reward.append(last_action_reward)
batch_vr_si.reverse()
batch_vr_R.reverse()
batch_vr_last_action_reward.reverse()
return batch_vr_si, batch_vr_last_action_reward, batch_vr_R
def _process_rp(self):
# [Reward prediction]
rp_experience_frames = self.experience.sample_rp_sequence()
# 4 frames
batch_rp_si = []
batch_rp_c = []
for i in range(3):
batch_rp_si.append(rp_experience_frames[i].state)
# one hot vector for target reward
r = rp_experience_frames[3].reward
rp_c = [0.0, 0.0, 0.0]
if r == 0:
rp_c[0] = 1.0 # zero
elif r > 0:
rp_c[1] = 1.0 # positive
else:
rp_c[2] = 1.0 # negative
batch_rp_c.append(rp_c)
return batch_rp_si, batch_rp_c
def process(self, sess, global_t, summary_writer, summary_op, score_input):
# Fill experience replay buffer
if not self.experience.is_full():
self._fill_experience(sess)
return 0
start_local_t = self.local_t
cur_learning_rate = self._anneal_learning_rate(global_t)
# Copy weights from shared to local
sess.run(self.sync)
# [Base]
batch_si, batch_last_action_rewards, batch_a, batch_adv, batch_R, start_lstm_state = \
self._process_base(sess,
global_t,
summary_writer,
summary_op,
score_input)
feed_dict = {
self.local_network.base_input: batch_si,
self.local_network.base_last_action_reward_input:
batch_last_action_rewards,
self.local_network.base_a: batch_a,
self.local_network.base_adv: batch_adv,
self.local_network.base_r: batch_R,
self.local_network.base_initial_lstm_state: start_lstm_state,
# [common]
self.learning_rate_input: cur_learning_rate
}
# [Pixel change]
if self.use_pixel_change:
batch_pc_si, batch_pc_last_action_reward, batch_pc_a, batch_pc_R = self._process_pc(
sess)
pc_feed_dict = {
self.local_network.pc_input: batch_pc_si,
self.local_network.pc_last_action_reward_input:
batch_pc_last_action_reward,
self.local_network.pc_a: batch_pc_a,
self.local_network.pc_r: batch_pc_R
}
feed_dict.update(pc_feed_dict)
# [Value replay]
if self.use_value_replay:
batch_vr_si, batch_vr_last_action_reward, batch_vr_R = self._process_vr(
sess)
vr_feed_dict = {
self.local_network.vr_input: batch_vr_si,
self.local_network.vr_last_action_reward_input:
batch_vr_last_action_reward,
self.local_network.vr_r: batch_vr_R
}
feed_dict.update(vr_feed_dict)
# [Reward prediction]
if self.use_reward_prediction:
batch_rp_si, batch_rp_c = self._process_rp()
rp_feed_dict = {
self.local_network.rp_input: batch_rp_si,
self.local_network.rp_c_target: batch_rp_c
}
feed_dict.update(rp_feed_dict)
# Calculate gradients and copy them to global network.
sess.run(self.apply_gradients, feed_dict=feed_dict)
self._print_log(global_t)
# Return advanced local step size
diff_local_t = self.local_t - start_local_t
return diff_local_t
class Trainer(object):
def __init__(self, rank, args, shared_model, optimizer, lr):
# CUDA 相关
self.gpu_id = args.gpu_ids[rank % len(args.gpu_ids)]
torch.manual_seed(args.seed + rank)
if self.gpu_id >= 0:
torch.cuda.manual_seed(args.seed + rank)
self.replay_buffer = Experience(history_size=2000)
self.cx = None # todo: 仍然是 一次 step 就前向传播
self.hx = None
self.episodic_score = 0
self.rank = rank
self.args = args
self.shared_model = shared_model
self.optimizer = optimizer
self.local_t = 0
# 初始化
# 初始化环境
print('Training Agent: {}'.format(self.rank))
# todo: 需要给 gym 环境加上 pc 等
# agent 代理对象
model = UNREAL(in_channels=3, action_size=6, enable_pixel_control=True)
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
model = model.cuda()
model.train()
# 学习率
self.initial_learning_rate = lr
self.max_global_time_step = 10 * 10**7
# 记录时间
# For log output
self.prev_local_t = 0
self.model = model
self.env = None
self.reset() # cx hx
def train(self, global_t, summary_writer=None):
t = self.local_t
if not self.replay_buffer.is_full():
self.fill_experience()
return 0 # time_step = 0
# sync
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
self.model.load_state_dict(self.shared_model.state_dict())
else:
self.model.load_state_dict(self.shared_model.state_dict())
loss_a3c, episode_score = self.process_a3c()
# 获取 hx, cx
h0, c0 = self.hx.detach(), self.cx.detach()
loss_pc = self.process_pc(h0=h0, c0=c0)
h0, c0 = self.hx.detach(), self.cx.detach()
loss_vr = self.process_vr(h0, c0)
loss_rp = self.process_rp()
loss = loss_a3c + loss_pc + loss_vr + loss_rp
self.model.zero_grad()
loss.backward()
clip_grad_norm_(self.model.parameters(), 40.0)
ensure_shared_grads(self.model,
self.shared_model,
gpu=self.gpu_id >= 0)
self.adjust_learning_rate(optimizer=self.optimizer,
global_time_step=global_t)
self.optimizer.step()
if summary_writer is not None:
with torch.no_grad():
losses = list(
map(lambda x: float(x.detach().cpu().numpy()),
[loss_a3c, loss_pc, loss_vr, loss_rp, loss]))
tags = dict(
zip(['a3c', 'pc', 'vr', 'rp', 'total_loss'], losses))
summary_writer.add_scalars('losses',
tags,
global_step=global_t)
# 分数
if episode_score:
summary_writer.add_scalars('score',
{'score': episode_score},
global_step=global_t)
self._print_log(global_t)
return self.local_t - t # offset
def init_env(self):
self.env = make_environment(env_type=self.args.env_type,
env_name=self.args.env_name,
channel_first=True)
self.env.reset() # 之后, 每一局结束都会自动 reset
def adjust_learning_rate(self, optimizer, global_time_step):
learning_rate = self.initial_learning_rate * (
self.max_global_time_step -
global_time_step) / self.max_global_time_step
if learning_rate < 0.0:
learning_rate = 0.0
for param_group in optimizer.param_groups:
param_group['lr'] = learning_rate
@property
def action_size(self):
return 6 # self.env.action_space.n
def set_start_time(self, start_time):
self.start_time = start_time
def _print_log(self, global_t):
if (self.rank == 0) and (self.local_t - self.prev_local_t >=
PERFORMANCE_LOG_INTERVAL):
self.prev_local_t += PERFORMANCE_LOG_INTERVAL
elapsed_time = time.time() - self.start_time
steps_per_sec = global_t / elapsed_time
print(
"### Performance : {} STEPS in {:.0f} sec. {:.0f} STEPS/sec. {:.2f}M STEPS/hour"
.format(global_t, elapsed_time, steps_per_sec,
steps_per_sec * 3600 / 1000000.))
def fill_experience(self):
prev_state = self.env.last_state
last_action = self.env.last_action
last_reward = self.env.last_reward
last_action_reward = ExperienceFrame.concat_action_and_reward(
last_action, self.action_size, last_reward)
with torch.no_grad():
state = torch.from_numpy(self.env.last_state['rgb']).unsqueeze(0)
lar = torch.from_numpy(last_action_reward).unsqueeze(0)
# whether to gpu
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
state = state.cuda()
lar = lar.cuda()
_, logits, (self.hx, self.cx) = self.model(task_type='a3c',
states=state,
hx=self.hx,
cx=self.cx,
last_action_reward=lar)
action_index = self.choose_action(
pi_values=F.softmax(logits, 1).cpu().numpy()[0])
obs, reward, terminal, _ = self.env.step(action_index) # 存储为数组
frame = ExperienceFrame(prev_state['rgb'], reward, action_index,
terminal, obs['pixel_change'], last_action,
last_reward)
self.replay_buffer.add_frame(frame)
if terminal:
self.env.reset()
else:
# 更新 LSTM 状态
self.hx = self.hx.detach()
self.cx = self.cx.detach()
if self.replay_buffer.is_full():
self.env.reset()
print("Replay buffer filled")
def choose_action(self, pi_values):
return np.random.choice(range(len(pi_values)), p=pi_values)
def l2_loss(self, x, y):
return F.mse_loss(x, y, reduction='sum') * 0.5
def process_a3c(self):
rewards = []
log_probs = [] # 指定的行为 概率
entropies = []
values = []
action_one_hot = []
# adv = [] # GAE, 采用 advantage 函数
terminal_end = False # 结束采样的时候是否是终止状态
episode_score = None # 决定是否显示 episodic score
for t in range(20):
state = torch.from_numpy(self.env.last_state['rgb']).unsqueeze(
dim=0) # batch = 1
last_action = self.env.last_action
last_reward = self.env.last_reward
last_action_reward = torch.from_numpy(
ExperienceFrame.concat_action_and_reward(
last_action, self.action_size, last_reward)).unsqueeze(0)
# whether to gpu
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
state = state.cuda()
last_action_reward = last_action_reward.cuda()
value, logits, (self.hx, self.cx) = self.model(
'a3c',
state,
hx=self.hx,
cx=self.cx,
last_action_reward=last_action_reward)
prob = F.softmax(logits, dim=1) # batch, 6.
log_prob = torch.log(
prob.clamp(1e-20, 1.0)
) # F.log_softmax(logits, dim=1).clamp(1e-20, 1.0) # batch, 6. NaN
#
entropy = -(log_prob * prob).sum(1)
# 采取行为
with torch.no_grad():
action_index = self.choose_action(
pi_values=F.softmax(logits, 1).cpu().numpy()[0])
prev_state = self.env.last_state['rgb']
observation, reward, terminal, _ = self.env.step(action_index)
# 显示信息
if self.rank == 0 and self.local_t % 100 == 0:
print("pi={}".format(prob.detach().cpu().numpy()))
print(" V={}".format(value.detach().cpu().numpy()))
self.local_t += 1
# 添加到 replay buffer
frame = ExperienceFrame(prev_state, reward, action_index, terminal,
observation['pixel_change'], last_action,
last_reward)
# Store to experience
self.replay_buffer.add_frame(frame)
entropies.append(entropy)
values.append(value)
rewards.append(reward)
log_probs.append(log_prob)
a = torch.zeros(self.action_size, dtype=torch.float32)
# a = np.zeros([self.action_size], dtype=np.float32)
a[action_index] = 1.0
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
a = a.cuda()
action_one_hot.append(a)
self.episodic_score += reward
if terminal:
print('Score: {0}'.format(self.episodic_score))
episode_score = self.episodic_score
self.episodic_score = 0
terminal_end = True
self.env.reset()
self.reset()
break
else:
self.hx = self.hx.detach()
self.cx = self.cx.detach()
# 计算 R
R = torch.zeros(1, 1)
if not terminal_end:
with torch.no_grad():
# 这里进行 bootstrapping
state = torch.from_numpy(observation['rgb']).unsqueeze(0)
lar = torch.from_numpy(
frame.get_action_reward(self.action_size)).unsqueeze(0)
# whether to gpu
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
state = state.cuda()
lar = lar.cuda()
value, _, (_, _) = self.model(task_type='a3c',
states=state,
hx=self.hx,
cx=self.cx,
last_action_reward=lar)
R = value.detach(
) # 这个值为 V(s_t,\theta_v^'), 在计算 actor 的梯度的时候可能会算入, 所以 detach
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
R = R.cuda()
# values.append(R) # 用于 bootstrapping
# 准备计算 loss
policy_loss = 0
value_loss = 0
# 反向计算 loss
for i in reversed(range(len(rewards))): # i=t - 1,...,t_start
R = self.args.gamma * R + rewards[i] # R <- r_t + \gamma * R
adv = R - values[
i] # GAE = advantage, R - V(s_i;\theta^'_v), adv 在反向传播的时候是 \theta_v的梯度, 在 policy 的梯度需要 detach
value_loss += 0.5 * self.l2_loss(
R, values[i]) # 定义为 0.5 * MSELOSS(R - value), 学习率是 actor 一半
log_prob_a = (log_probs[i] * action_one_hot[i]).sum(
1) # log(a_i|s_i;\theta^')
policy_loss += -log_prob_a * adv.detach(
) + entropies[i] * 0.001 # entropy_beta
return value_loss + policy_loss, episode_score
def process_pc(self, h0, c0):
"""
:param h0: PC 起始的 LSTM 的 hidden 状态
:param c0:
:return:
"""
# [pixel change]
# Sample 20+1 frame (+1 for last next state)
pc_experience_frames = self.replay_buffer.sample_sequence(20 + 1)
# Reverse sequence to calculate from the last
pc_experience_frames.reverse()
batch_pc_si = []
batch_pc_a = []
batch_pc_R = []
batch_pc_q = []
batch_pc_last_action_reward = []
pc_R = torch.zeros(20, 20)
if not pc_experience_frames[1].terminal:
with torch.no_grad():
state = torch.from_numpy(
pc_experience_frames[0].state).unsqueeze(
dim=0) # batch = 1
last_action_reward = torch.from_numpy(
pc_experience_frames[0].get_last_action_reward(
self.action_size)).unsqueeze(0)
# whether to gpu
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
state = state.cuda()
last_action_reward = last_action_reward.cuda()
# 每次forward传播的时候需要使用全零的 hx, cx
(hx, cx) = self.get_zero_hx_cx()
_, pc_q_max, (_, _) = self.model(
'pc',
state,
hx=hx,
cx=cx,
last_action_reward=last_action_reward) # 不更新 lstm 的状态
pc_R = pc_q_max.detach() # 共享 LSTM 的参数
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
pc_R.cuda()
for frame in pc_experience_frames[1:]:
pc = torch.from_numpy(frame.pixel_change)
a = torch.zeros(self.action_size, dtype=torch.float32)
# a = np.zeros([self.action_size], dtype=np.float32)
a[frame.action] = 1.0
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
a = a.cuda()
pc = pc.cuda()
pc_R = pc + 0.9 * pc_R # gamma_pc = 0.9
last_action_reward = frame.get_last_action_reward(self.action_size)
batch_pc_si.append(frame.state)
batch_pc_a.append(a)
batch_pc_R.append(pc_R)
batch_pc_last_action_reward.append(last_action_reward)
batch_pc_si.reverse()
batch_pc_a.reverse()
batch_pc_R.reverse()
batch_pc_last_action_reward.reverse()
# 计算输出的
for (state, last_action_reward) in zip(batch_pc_si,
batch_pc_last_action_reward):
# 获取输出的 pc_q, 需要在网络中计算一遍
state = torch.from_numpy(state).unsqueeze(dim=0) # batch = 1
last_action_reward = torch.from_numpy(
last_action_reward).unsqueeze(0)
# whether to gpu
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
state = state.cuda()
last_action_reward = last_action_reward.cuda()
# 在反向传播的时候, LSTM 的Cell和 Hidden 与 A3C 最后一次的输出有关, 在 tf 的版本中, feed_dict设置的 初始LSTM 状态就是 A3C 运行之前的状态
pc_q, _, (h0, c0) = self.model(
'pc',
state,
hx=h0,
cx=c0, # 使用 A3C 初始的状态
last_action_reward=last_action_reward) # 不更新 lstm 的状态
# 修改 LSTM cell 和 hidden 的状态
h0 = h0.detach()
c0 = c0.detach()
batch_pc_q.append(pc_q)
# 损失, 使用 n-step Qlearning
batch_pc_R = torch.cat(batch_pc_R, dim=0)
batch_pc_q = torch.cat(batch_pc_q, dim=0) # 按照第一个维度连接
pc_a_reshape = torch.stack(batch_pc_a).view(-1, self.action_size, 1, 1)
pc_qa_ = torch.mul(batch_pc_q, pc_a_reshape)
pc_qa = pc_qa_.sum(dim=1, keepdim=False) # -1, 20, 20
pc_loss = 0.05 * self.l2_loss(pc_qa,
batch_pc_R) # pixel_change_lambda = 0.05
return pc_loss
def process_vr(self, h0, c0):
# [Value replay]
# Sample 20+1 frame (+1 for last next state)
vr_experience_frames = self.replay_buffer.sample_sequence(20 + 1)
# Reverse sequence to calculate from the last
vr_experience_frames.reverse()
batch_vr_si = []
batch_vr_R = []
batch_vr_last_action_reward = []
batch_values = []
vr_R = torch.zeros(1, 1)
if not vr_experience_frames[1].terminal:
with torch.no_grad():
state = torch.from_numpy(
vr_experience_frames[0].state).unsqueeze(
dim=0) # batch = 1
last_action_reward = torch.from_numpy(
vr_experience_frames[0].get_last_action_reward(
self.action_size)).unsqueeze(0)
# whether to gpu
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
state = state.cuda()
last_action_reward = last_action_reward.cuda()
#
(hx, cx) = self.get_zero_hx_cx()
vr_R, (_, _) = self.model(
'vr',
state,
hx=hx,
cx=cx,
last_action_reward=last_action_reward) # 不更新 lstm 的状态
vr_R = vr_R.detach() # value 的参数, \theta^'_v
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
vr_R.cuda()
# t_max times loop
for frame in vr_experience_frames[1:]:
vr_R = frame.reward + self.args.gamma * vr_R
batch_vr_si.append(frame.state)
batch_vr_R.append(vr_R)
last_action_reward = frame.get_last_action_reward(self.action_size)
batch_vr_last_action_reward.append(last_action_reward)
batch_vr_si.reverse()
batch_vr_R.reverse()
batch_vr_last_action_reward.reverse()
for (state, last_action_reward) in zip(batch_vr_si,
batch_vr_last_action_reward):
# 获取输出的 pc_q, 需要在网络中计算一遍
state = torch.from_numpy(state).unsqueeze(dim=0) # batch = 1
last_action_reward = torch.from_numpy(
last_action_reward).unsqueeze(0)
# whether to gpu
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
state = state.cuda()
last_action_reward = last_action_reward.cuda()
value, (h0, c0) = self.model(
'vr',
state,
hx=h0,
cx=c0,
last_action_reward=last_action_reward) # 不更新 lstm 的状态
h0 = h0.detach()
c0 = c0.detach()
batch_values.append(value)
batch_vr_R = torch.cat(batch_vr_R, dim=0)
batch_values = torch.cat(batch_values, dim=0)
vr_loss = self.l2_loss(batch_vr_R, batch_values)
return vr_loss
def process_rp(self):
# [Reward prediction]
rp_experience_frames = self.replay_buffer.sample_rp_sequence()
# 4 frames
batch_rp_si = []
# batch_rp_c = []
for i in range(3):
batch_rp_si.append(rp_experience_frames[i].state)
# 求输出的 标签 0 + -
states = torch.from_numpy(np.stack(batch_rp_si)) # batch = 3
# whether to gpu
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
states = states.cuda()
logits = self.model('rp',
states,
hx=None,
cx=None,
last_action_reward=None) # 不更新 lstm 的状态, conv的参数
# one hot vector for target reward
r = rp_experience_frames[3].reward
rp_c_target = torch.zeros(3, dtype=torch.float32)
if r == 0:
rp_c_target[0] = 1.0 # zero
elif r > 0:
rp_c_target[1] = 1.0 # positive
else:
rp_c_target[2] = 1.0 # negative
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
rp_c_target = rp_c_target.cuda()
rp_c = F.softmax(logits, 1).clamp(1e-20, 1.0)
rp_loss = -(rp_c_target * torch.log(rp_c)).sum()
return rp_loss
def reset(self):
self.cx, self.hx = self.get_zero_hx_cx()
def get_zero_hx_cx(self, batch_size=1):
if self.gpu_id >= 0:
with torch.cuda.device(self.gpu_id):
cx = Variable(torch.zeros(1, batch_size, 256).cuda())
hx = Variable(torch.zeros(1, batch_size, 256).cuda())
else:
cx = Variable(torch.zeros(1, batch_size,
256)) # todo: 仍然是 一次 step 就前向传播
hx = Variable(torch.zeros(1, batch_size, 256))
return (cx, hx)
def close(self):
self.env.close()