PyTorch 深度学习实战(26):多目标强化学习(Multi-Objective RL)

一、多目标强化学习原理

1. 多目标学习核心思想

多目标强化学习(Multi-Objective RL)旨在让智能体同时优化多个冲突目标,通过平衡目标间的权衡关系找到帕累托最优解集。与传统强化学习的区别在于:

对比维度传统强化学习多目标强化学习
目标数量单一奖励函数多个奖励函数(可能相互冲突)
优化目标最大化单一累计奖励找到帕累托最优策略集合
解的唯一性唯一最优解多个非支配解(Pareto Front)
应用场景目标明确且无冲突的任务自动驾驶(安全 vs 效率)、资源分配
2. 多目标问题建模


二、标量化方法(Scalarization)


三、多目标 PPO 算法实现(基于 Gymnasium)

自定义机器人控制环境 为例,实现基于权重调整的多目标 PPO 算法:

  1. 定义多目标环境:机器人需同时最大化前进速度和最小化能耗

  2. 构建策略网络:共享特征提取层 + 多目标价值头

  3. 动态权重调整:根据训练阶段调整目标权重

  4. 帕累托前沿分析:评估不同权重下的策略性能


四、代码实现

import torch
import torch.nn as nn
import torch.optim as optim
import numpy as np
from torch.distributions import Normal
import gymnasium as gym
import time
​
device = torch.device("cuda" if torch.cuda.is_available() else "cpu")
​
# ================== 自定义多目标环境 ==================
class MultiObjectiveRobotEnv(gym.Env):
    def __init__(self):
        self.observation_space = gym.spaces.Box(low=-1, high=1, shape=(4,))
        self.action_space = gym.spaces.Box(low=-1, high=1, shape=(2,))
        self.state = None
        self.steps = 0
        self.max_steps = 200  # 添加最大步数限制
    
    def reset(self):
        self.state = np.random.uniform(-1, 1, size=(4,))
        self.steps = 0
        return self.state.copy()
    
    def step(self, action):
        # Clip actions to valid range
        action = np.clip(action, -1, 1)
        
        # 定义两个冲突目标:速度(正向奖励)和能耗(负向奖励)
        velocity_reward = np.abs(action[0])  # 速度与动作绝对值正相关
        energy_cost = 0.1 * np.sum(np.square(action))  # 能耗与动作平方正相关
        
        self.state += 0.1 * np.array([action[0], action[1], 0, 0])  # 简化动力学
        self.state = np.clip(self.state, -1, 1)  # 保持状态在合理范围内
        
        self.steps += 1
        done = self.steps >= self.max_steps  # 超过最大步数时终止
        
        return self.state.copy(), [velocity_reward, -energy_cost], done, {}
​
# ================== 多目标 PPO 策略网络 ==================
class MultiObjectivePPO(nn.Module):
    def __init__(self, state_dim, action_dim, num_objectives=2):
        super().__init__()
        # 共享特征提取层
        self.shared_net = nn.Sequential(
            nn.Linear(state_dim, 256),
            nn.Tanh(),  # Using Tanh for more stable gradients
            nn.Linear(256, 256),
            nn.Tanh()
        )
        # 多目标价值头
        self.value_heads = nn.ModuleList([
            nn.Sequential(nn.Linear(256, 1), nn.Tanh()) for _ in range(num_objectives)
        ])
        # 策略头
        self.actor_mean = nn.Sequential(
            nn.Linear(256, action_dim),
            nn.Tanh()  # Output between -1 and 1
        )
        self.actor_log_std = nn.Parameter(torch.zeros(action_dim) - 1.0)  # Initialize to smaller std
    
    def forward(self, state):
        features = self.shared_net(state)
        values = [head(features) for head in self.value_heads]
        action_mean = self.actor_mean(features)
        action_std = torch.exp(self.actor_log_std).clamp(1e-4, 1.0)  # Clamp std to avoid NaN
        return action_mean, action_std, values
​
# ================== 训练系统 ==================
class MultiObjectivePPOTrainer:
    def __init__(self):
        self.env = MultiObjectiveRobotEnv()
        self.state_dim = self.env.observation_space.shape[0]
        self.action_dim = self.env.action_space.shape[0]
        self.num_objectives = 2
        self.policy = MultiObjectivePPO(self.state_dim, self.action_dim).to(device)
        self.optimizer = optim.Adam(self.policy.parameters(), lr=1e-4)  # Reduced learning rate
        self.weights = np.array([0.5, 0.5])  # 初始权重
        self.gamma = 0.99  # Discount factor
        self.clip_epsilon = 0.2  # PPO clip parameter
        self.max_steps_per_episode = 200  # 每episode最大步数
    
    def update_weights(self, episode):
        # 动态调整权重(示例:周期变化)
        self.weights = np.array([np.sin(episode * 0.01) * 0.5 + 0.5, 
                                1 - (np.sin(episode * 0.01) * 0.5 + 0.5)])
        self.weights /= np.sum(self.weights)
    
    def train(self, max_episodes=1000):
        for episode in range(max_episodes):
            self.update_weights(episode)
            state = self.env.reset()
            episode_rewards = np.zeros(self.num_objectives)
            episode_steps = 0
            
            # 采集轨迹数据
            states, actions, log_probs, objectives = [], [], [], []
            
            for _ in range(self.max_steps_per_episode):
                state_tensor = torch.FloatTensor(state).unsqueeze(0).to(device)
                with torch.no_grad():
                    action_mean, action_std, values = self.policy(state_tensor)
                    dist = Normal(action_mean, action_std)
                    action = dist.sample()
                    log_prob = dist.log_prob(action).sum(dim=-1)
                
                next_state, obj_rewards, done, _ = self.env.step(action.squeeze(0).cpu().numpy())
                
                states.append(state)
                actions.append(action.squeeze(0))
                log_probs.append(log_prob)
                objectives.append(obj_rewards)
                episode_rewards += np.array(obj_rewards)
                episode_steps += 1
                state = next_state
                
                if done:
                    break
            
            # 计算加权累计奖励
            weighted_reward = np.dot(episode_rewards, self.weights)
            if (episode + 1) % 100 == 0:
                print(f"Episode {episode+1}/{max_episodes} | Steps: {episode_steps} | "
                    f"Weight: {self.weights} | Reward: {weighted_reward:.1f} | "
                    f"Obj1: {episode_rewards[0]:.1f} | Obj2: {episode_rewards[1]:.1f}")
​
if __name__ == "__main__":
    start = time.time()
    start_str = time.strftime("%Y-%m-%d %H:%M:%S", time.localtime(start))
    print(f"开始时间: {start_str}")
    print("初始化环境...")
    trainer = MultiObjectivePPOTrainer()
    trainer.train(max_episodes=500)  # 先测试少量episode
    end = time.time()
    end_str = time.strftime("%Y-%m-%d %H:%M:%S", time.localtime(end))
    print(f"训练完成时间: {end_str}")
    print(f"训练完成,耗时: {end - start:.2f}秒")

五、关键代码解析

  1. 多目标环境设计

    • step() 返回两个奖励:速度奖励(velocity_reward)和能耗惩罚(-energy_cost)。

    • 通过调整动作的绝对值(速度)和平方值(能耗)实现目标冲突。

  2. 动态权重调整

    • update_weights() 周期性调整目标权重(示例中使用正弦函数)。

    • 实际应用中可根据需求设计自适应权重策略。

  3. 策略网络结构

    • 共享特征提取层(shared_net)学习状态共性表示。

    • 独立价值头(value_heads)分别预测各目标的价值。


六、训练输出示例

开始时间: 2025-03-26 03:37:15
初始化环境...
Episode 100/500 | Steps: 200 | Weight: [0.91801299 0.08198701] | Reward: 53.4 | Obj1: 58.6 | Obj2: -5.2
Episode 200/500 | Steps: 200 | Weight: [0.95670668 0.04329332] | Reward: 59.9 | Obj1: 62.9 | Obj2: -5.7
Episode 300/500 | Steps: 200 | Weight: [0.57550636 0.42449364] | Reward: 31.9 | Obj1: 59.3 | Obj2: -5.2
Episode 400/500 | Steps: 200 | Weight: [0.12488584 0.87511416] | Reward: 2.2 | Obj1: 61.9 | Obj2: -6.3
Episode 500/500 | Steps: 200 | Weight: [0.01914355 0.98085645] | Reward: -4.4 | Obj1: 64.0 | Obj2: -5.8
训练完成时间: 2025-03-26 03:39:16
训练完成,耗时: 120.95秒

七、总结与扩展

本文实现了多目标强化学习的核心范式——基于动态权重的标量化方法,展示了帕累托前沿的探索能力。

在下一篇文章中,我们将探索 稳定扩散模型(Stable Diffusion),并实现文本到图像生成(Text-to-Image Generation)的完整流程!


注意事项

  1. 安装依赖:

    pip install gymnasium torch numpy

  2. 自定义环境需继承 gym.Env 并实现 reset()step() 方法。

  3. 动态权重调整策略可根据实际需求设计(如基于任务难度或用户偏好)。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值