Categories: 数学建模.

Situation(背景)

某码头设有 c 个泊位供货轮停靠装卸。由于泊位有限,到港货轮需排队等待。保险公司计划推出一款保险产品:货轮到港时缴纳保费 P,若实际等待时间超过承诺时限 T,保险公司支付赔偿金 CC ≫ P)。保费须与承诺时间关联,需为保险公司构建定价模型,探讨盈亏平衡与利润最大化条件。


Task(任务)

  • 将码头作业抽象为数学模型,导出等待时间概率分布
  • 建立保险公司期望利润函数,分析盈亏平衡条件
  • 针对保费与承诺时间的不同耦合关系,构造多种定价模型并求解利润最大化
  • 通过数值实验验证理论推导,分析参数敏感性

Action(行动)

排队论模型:M/M/c 系统

将码头作业抽象为 M/M/c 排队系统——货轮到达服从 Poisson 过程(到达率 λ),每个泊位服务时间服从负指数分布(服务率 μ),FCFS 排队。在稳定性条件 ρ = λ/(cμ) < 1 下,导出稳态概率分布与 Erlang C 公式(排队等待概率):

$$P_Q = \frac{\dfrac{a^c}{c!} \cdot \dfrac{1}{1-\rho}}{\displaystyle\sum_{n=0}^{c-1} \dfrac{a^n}{n!} + \dfrac{a^c}{c!} \cdot \dfrac{1}{1-\rho}}$$

等待时间超出承诺时限 T 的概率为:

P(W > T) = PQ ⋅ e−(cμ − λ)T

保险利润模型

单船期望利润:

π(P, T) = P − C ⋅ PQ ⋅ e−(cμ − λ)T

盈亏平衡条件:P ≥ C ⋅ PQ ⋅ e−(cμ − λ)T

三种定价模型

模型 A:精算公平定价 + 安全附加

P(T) = (1 + θ) ⋅ C ⋅ PQ ⋅ e−(cμ − λ)T

保险公司始终盈利,无约束最优为 T → 0。该模型只考虑保险人视角,未纳入需求侧约束。

模型 B:指数衰减定价(市场驱动型)

P(T) = p0eγT

利润最大化的一阶条件给出最优承诺时间:

$$T^* = \frac{1}{(c\mu - \lambda) - \gamma} \ln\left(\frac{(c\mu - \lambda) C P_Q}{\gamma p_0}\right)$$

极大值条件:保费衰减率 γ 必须低于系统过剩服务率 cμ − λ

模型 C:线性定价

P(T) = p0 − kT

最优承诺时间:

$$T^* = \frac{1}{c\mu - \lambda} \ln\left(\frac{(c\mu - \lambda) C P_Q}{k}\right)$$

二阶条件始终满足。

市场约束

引入两项市场合理性约束: - 最大附加系数 Lmax = 3:保费不得超过期望赔付的 3 倍 - 最低赔付概率 pmin = 0.5%:保单实际赔付概率不得低于 0.5%

数值实验

以基准参数 (λ = 8, μ = 3, c = 3, C = 100000, p0 = 8000, γ = 0.5) 开展系统的数值实验,使用 Python(numpy/scipy/matplotlib)实现排队系统仿真、三种定价模型求解、参数扫描与可视化。


Result(结果)

核心发现

1. 约束收窄策略空间

无约束 FOC 给出的内点解(指数模型 TFOC* = 5.99 天)不满足市场约束。引入 Lmaxpmin 后,可行区间收窄至 [4.61, 5.06] 天(宽度仅 0.45 天),最优解移至 pmin 约束边界。

2. 保费衰减率 γ 的可行区间远比理论预期狭窄

理论区间为 γ ∈ [0.33, 1.0],但实际可行区间仅为 γ ∈ [0.25, 0.519],宽度为理论区间的 27%。γ > 0.519 时盈亏平衡下界超越 pmin 约束上界,可行区间闭合。

3. 利润随泊位数 c 单调递增

c 从 3 增到 6 时,单船利润从 131 元跃升至 4780 元(增幅约 36 倍)。机制:c 增大同时压低 PQ(赔付基准缩小)和抬高 cμ − λ(风险衰减加快),而保费仅以 γ 速率衰减——两者衰减速率之差持续扩大。

4. 到达率 λ 对利润的系统性压制

λ 增大只产生单方向不利影响(PQ 上升),缺乏对冲机制。λ ≥ 8.5 时指数模型丧失可行解,该阈值先于系统稳定性极限 λ = 9 出现——市场约束比系统稳定性设定了更严苛的准入门槛。

5. 利润的数学本质

$$\pi(T) = \underbrace{p_0 e^{-\gamma T}}_{\text{保费端}} - \underbrace{C P_Q e^{-(c\mu-\lambda)T}}_{\text{赔付端}}$$

利润来源于两项指数衰减的速率差。保险公司的盈利空间不取决于保费的绝对水平,而取决于保费与赔付随 T 下降的”速度差”。

三种定价模型对比

模型 最优 T* 单船利润 特点
精算公平 0 天(无约束)/ 5.07 天(约束) 23925 元 / 约束下更低 始终盈利,未考虑需求侧
指数衰减 5.06 天(约束边界) 131 元 可行区间极窄,利润微薄
线性 3.65 天 4097 元 T 较小时保费高,利润显著占优

可视化

数值实验生成 7 张分析图:三种定价模型利润曲线对比、γ 敏感性、λ 系统性影响、(T, λ) 利润热力图、泊位数影响、γ 连续扫描、保费构成分解。

三种定价模型利润曲线
保费衰减率 γ 对比
到达率敏感性
利润热力图
泊位数影响
γ 敏感性
保费构成分解

项目文件

完整数学建模报告(含公式推导、数值实验、7 张可视化图表、Python 代码)位于项目目录。

代码由三个 Python 文件协同完成:insurance_model.py(模型库)、plot_analysis.py(可视化)、compute_tables.py(数值表格)。