python案例如何结合士气指数做决策？

wen python案例 2026-08-14 166

本文目录导读：

第一步：构建士气指数（数据准备）
第二步：结合决策场景（示例：是否推行“弹性工作制”）
第三步：进阶决策（多Agent模拟 + 蒙特卡洛模拟）
第四步：可视化决策边界
核心决策逻辑总结
注意事项

将“士气指数”数字化并用于决策，是经典的管理科学与数据科学结合案例，核心逻辑在于：将抽象的“人”的因素（情绪、动力）转化为可量化的KPI，并通过算法模型模拟不同决策对士气的影响，从而辅助管理层选择最优路径。

以下我将从构建指数、建立决策模型、代码实战、应用场景四个维度,用Python给出可运行的案例。

第一步：构建士气指数（数据准备）

士气指数通常由多个可观测的指标加权合成，假设我们关注一个研发团队,每月收集以下数据：

加班时长（小时）
任务完成率（%）
内部培训参与度（%）
员工流失预警（0-1,1为高风险）
即时反馈分数（1-5）

Python代码：生成模拟数据并计算士气指数

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
from sklearn.preprocessing import MinMaxScaler
# 1. 生成模拟数据（3个月，每月30天数据）
np.random.seed(42)
dates = pd.date_range(start='2024-01-01', periods=90, freq='D')
data = {
    '日期': dates,
    '加班时长': np.random.normal(3, 1.5, 90).clip(0, 8),
    '任务完成率': np.random.uniform(0.6, 1.0, 90),
    '培训参与度': np.random.uniform(0.3, 1.0, 90),
    '流失预警': np.random.binomial(1, 0.15, 90),  # 15%概率出预警
    '反馈分数': np.random.uniform(2.5, 5.0, 90)
}
df = pd.DataFrame(data)
# 2. 标准化（注意：流失预警属于反向指标，需反向处理）
scaler = MinMaxScaler()
df['标准加班'] = 1 - scaler.fit_transform(df[['加班时长']])  # 加班越多，分越低
df['标准完成率'] = scaler.fit_transform(df[['任务完成率']])
df['标准培训'] = scaler.fit_transform(df[['培训参与度']])
df['标准流失'] = 1 - df['流失预警']  # 反向
df['标准反馈'] = scaler.fit_transform(df[['反馈分数']])
# 3. 定义权重（可根据业务重要性调整）
weights = {
    '标准加班': 0.25,
    '标准完成率': 0.25,
    '标准培训': 0.15,
    '标准流失': 0.20,
    '标准反馈': 0.15
}
# 4. 计算士气指数（0-1之间，越高越好）
df['士气指数'] = sum(df[col] * weight for col, weight in weights.items())
# 5. 查看前5天的士气
print(df[['日期', '士气指数']].head())

第二步：结合决策场景（示例：是否推行“弹性工作制”）

我们现在面临一个决策：未来一个月是否推行“弹性工作制”（允许员工每周任选2天远程办公）？

决策逻辑：

推行后，预计加班时长将下降20%（士气提升），但任务完成率可能因沟通成本上升下降5%。
如何用Python判断该决策是否有效？

代码实战：模拟决策前后的士气变化

# 假设推行弹性制后：
# 加班时长减少20%
# 任务完成率下降5%
# 培训参与度不变（远程也可参与）
# 流失预警下降50%（因为满意度提升）
# 反馈分数上升10%
def simulate_decision(df, effects):
    """
    effects 为字典：{'加班变化': -0.2, '完成率变化': -0.05, ...}
    返回新的士气指数平均值。
    """
    simulated_df = df.copy()
    # 调整指标
    simulated_df['加班时长'] = (simulated_df['加班时长'] * (1 + effects.get('加班变化', 0))).clip(0, 8)
    simulated_df['任务完成率'] = (simulated_df['任务完成率'] * (1 + effects.get('完成率变化', 0))).clip(0, 1)
    # 流失预警：因为是0/1变量，通过概率变化来模拟（简化处理：直接按比例减少1的个数）
    if effects.get('流失下降率', 0) < 0:
        n_risk = int(len(simulated_df) * (1 + effects['流失下降率']))
        simulated_df['流失预警'] = np.where(simulated_df['流失预警'] == 1, 0, simulated_df['流失预警'])[:n_risk]
    # 反馈分数增加10%
    simulated_df['反馈分数'] = (simulated_df['反馈分数'] * (1 + effects.get('反馈变化', 0))).clip(2.5, 5.0)
    # 重新标准化并计算（为了保证公平，使用原数据的scaler）
    # 由于修改了数据，需要重新标准化（简化：直接使用当前最大值最小值）
    scaler2 = MinMaxScaler()
    simulated_df['标准加班'] = 1 - scaler2.fit_transform(simulated_df[['加班时长']])
    simulated_df['标准完成率'] = scaler2.fit_transform(simulated_df[['任务完成率']])
    simulated_df['标准培训'] = scaler2.fit_transform(simulated_df[['培训参与度']])
    simulated_df['标准流失'] = 1 - simulated_df['流失预警']
    simulated_df['标准反馈'] = scaler2.fit_transform(simulated_df[['反馈分数']])
    simulated_df['士气指数'] = sum(simulated_df[col] * weight for col, weight in weights.items())
    return simulated_df['士气指数'].mean()
# 定义两种决策的预期影响
decision_A_no_change = {'加班变化': 0, '完成率变化': 0}  # 保持现状
decision_B_flexible = {'加班变化': -0.2, '完成率变化': -0.05, '流失下降率': -0.5, '反馈变化': 0.1}
# 计算两种决策的预期平均士气
baseline_morale = df['士气指数'].mean()
morale_if_decision_A = baseline_morale  # 不推行，假设不变
morale_if_decision_B = simulate_decision(df, decision_B_flexible)
print(f"现状平均士气: {morale_if_decision_A:.3f}")
print(f"推行弹性制后平均士气: {morale_if_decision_B:.3f}")
# 决策建议
if morale_if_decision_B > morale_if_decision_A:
    print("建议推行弹性工作制，预计士气提升 {:.2%}".format((morale_if_decision_B - morale_if_decision_A)/morale_if_decision_A))
else:
    print("不建议推行，可能弊大于利")

第三步：进阶决策（多Agent模拟 + 蒙特卡洛模拟）

更高级的做法是使用蒙特卡洛模拟来面对不确定性，我们不确定“弹性制”到底让加班减少多少（可能是10%~30%），这时可以通过随机抽样1000次,生成士气提升的概率分布。

import random
def monte_carlo_decision(df, n_sim=1000):
    results = []
    for _ in range(n_sim):
        # 随机抽样不确定的影响参数
        overtime_effect = -random.uniform(0.1, 0.3)  # 10%~30%减少
        completion_effect = -random.uniform(0.01, 0.08)  # 1%~8%下降
        loss_effect = -random.uniform(0.3, 0.7)  # 流失降低30%~70%
        feedback_effect = random.uniform(0.05, 0.15)  # 反馈提升5%~15%
        effects = {
            '加班变化': overtime_effect,
            '完成率变化': completion_effect,
            '流失下降率': loss_effect,
            '反馈变化': feedback_effect
        }
        morale = simulate_decision(df, effects)
        results.append(morale)
    # 统计结果
    mean_morale = np.mean(results)
    p_improve = np.mean([1 if m > baseline_morale else 0 for m in results])
    return mean_morale, p_improve
mean_morale, p_improve = monte_carlo_decision(df, n_sim=500)
print(f"蒙特卡洛模拟：平均士气 {mean_morale:.3f}, 士气提升概率 {p_improve:.1%}")

第四步：可视化决策边界

为了更直观地辅助决策，可以绘制士气指数热力图，横轴为“加班减少幅度”，纵轴为“完成率下降幅度”,颜色代表士气指数。

# 生成网格
x_axis = np.linspace(-0.3, 0, 20)  # 加班变化（-30% 到 0）
y_axis = np.linspace(-0.1, 0, 20)  # 完成率变化（-10% 到 0）
Z = np.zeros((len(x_axis), len(y_axis)))
for i, x in enumerate(x_axis):
    for j, y in enumerate(y_axis):
        effects = {'加班变化': x, '完成率变化': y, '流失下降率': -0.3, '反馈变化': 0.1}
        Z[j, i] = simulate_decision(df, effects)
plt.figure(figsize=(10, 8))
plt.contourf(x_axis, y_axis, Z, levels=20, cmap='RdYlGn')
plt.colorbar(label='预期士气指数')
plt.xlabel('加班时长变化比例')
plt.ylabel('任务完成率变化比例')'士气指数随决策参数变化热力图')
plt.scatter(decision_B_flexible['加班变化'], decision_B_flexible['完成率变化'], color='blue', s=100, label='建议方案')
plt.legend()
plt.show()

核心决策逻辑总结

决策步骤	Python对应手段	业务价值
量化士气	加权指标合成	将主观感受变为可比较数值
情景模拟	`simulate_decision()` 函数	预测不同政策的影响
风险评估	蒙特卡洛随机模拟	给出胜率而非盲目决策
可视化	热力图	快速定位最佳参数区间