02. LLM Params and FLOPs | 大模型参数量与算力推导
难度: Medium | 环境: CPU-first | 标签: 数学推导, Transformer | 目标人群: 通用基础 (算法/Infra)
🚀 云端运行环境
本章节的实战代码可以点击以下链接在免费 GPU 算力平台上直接运行:
先把参数量的构成拆开,再把前向推理与完整训练的 FLOPs 粗算清楚,这样后面的训练成本、吞吐评估和模型选型才会有统一的底座。
关键词: parameters, FLOPs, MFU
前置阅读
导语: 这一页要把“参数量 -> FLOPs -> 训练成本”这条线讲完整,所以最好先把数据格式和显存推导的直觉接上,再来看算力公式。
- Group 0B: PyTorch Tensors and Autograd | 0B: PyTorch 张量与自动求导
- 01. Data Types and Precision | 大模型的数据格式与混合精度
相关阅读
导语: 如果想继续把“参数量 -> 显存 -> 训练成本 -> 并行策略”这条线补完整,可以继续看:
- 06. VRAM Calculation and ZeRO | 显存计算与 ZeRO 优化:把训练 / 推理阶段的显存拆分接上。
- 22. MoE Parameter and Compute | MoE 模型参数量计算:补充 MoE 的参数量和计算量变化。
- 26. Parallel Strategy Decision Framework | 并行策略决策框架:最后看并行策略怎么影响整体成本。
Q1:假设隐藏层维度为 ,词表大小为 。请推导一个包含 层的标准 Transformer Decoder 的总参数量。
点击展开查看解析
我们把 Transformer 拆解为三大部分(忽略极小的 bias 和 LayerNorm 的权重,它们对百亿参数的占比不到千分之一):
1. 嵌入层 (Embedding Layer) 与 输出层 (LM Head)
- Token Embedding: 形状
,参数量为 。 - LM Head (输出映射): 形状
,参数量为 。 - (注:很多模型如 Gemma/Qwen 会共享这两个权重,参数量减半。这里我们假设不共享)。
2. 注意力机制 (Multi-Head Attention, MHA) 在每个 Decoder 块中:
- 投影 Q, K, V:三个形状为
的矩阵。参数量 。 - 投影 Output (O):一个形状为
的矩阵。参数量 。 - MHA 总参数量 =
。 (如果采用 GQA,K 和 V 的参数量会大幅减少,这里按最原始的 MHA 计算)。
3. 前馈神经网络 (Feed Forward Network, FFN / MLP) 在标准 GPT 架构中,隐藏层会先升维到
- 升维矩阵
: ,参数量 。 - 降维矩阵
: ,参数量 。 - FFN 总参数量 =
。 (如果在 LLaMA 中使用 SwiGLU,维度会变为 ,但有 3 个矩阵,总参数量依然是 )。
综上所述:
- 一个 Block 的参数量 =
(Attn) + (MLP) = 。 - 总参数量
。
带入 LLaMA-7B 感受一下:
实现一个函数,按照 Embedding + Attention + SwiGLU FFN + LayerNorm + LM Head 的方式估算 Transformer 参数量。
def calculate_transformer_params(vocab_size, hidden_dim, num_layers, intermediate_size=None, tie_embeddings=False):
if intermediate_size is None:
intermediate_size = 4 * hidden_dim
embedding_params = vocab_size * hidden_dim
attention_params = num_layers * (4 * hidden_dim * hidden_dim)
ffn_params = num_layers * (3 * hidden_dim * intermediate_size)
layernorm_params = num_layers * (2 * hidden_dim)
lm_head_params = 0 if tie_embeddings else vocab_size * hidden_dim
total_params = embedding_params + attention_params + ffn_params + layernorm_params + lm_head_params
return total_paramsdef test_calculate_transformer_params():
try:
result = calculate_transformer_params(1000, 64, 2, 256, tie_embeddings=True)
assert result == 195328, f"错误:期望 195328,实际 {result}"
result = calculate_transformer_params(1000, 64, 2, 256, tie_embeddings=False)
assert result == 259328, f"错误:期望 259328,实际 {result}"
result = calculate_transformer_params(32000, 4096, 32, 11008, tie_embeddings=False)
assert result > 6000000000, f"错误:LLaMA 级别模型参数量应大于 6B,实际 {result}"
print("✅ 参数量函数测试通过!")
except AssertionError as e:
print(f"❌ 测试失败: {e}")
except Exception as e:
print(f"❌ 运行错误: {e}")
test_calculate_transformer_params()Q1扩展验证:估算 LLaMA-7B 的参数量
使用上面的函数估算一个 7B 级模型的大致参数规模,并观察各模块的占比。
vocab_size = 32000
hidden_dim = 4096
num_layers = 32
intermediate_size = 11008
total_params = calculate_transformer_params(vocab_size, hidden_dim, num_layers, intermediate_size, tie_embeddings=False)
print(f"估算参数量: {total_params / 1e9:.2f}B")
print("提示:不同实现会因为是否共享词嵌入、是否计入偏置而略有差异。")Q2:前向传播 (Inference / Forward Pass) 的 FLOPs 是怎么计算的?
点击展开查看解析
在了解了参数量之后,我们来看大模型在进行推理(前向传播)时需要多少算力。
核心经验法则:1 个参数处理 1 个 Token,大约需要 2 次浮点运算(FLOPs)。 为什么是 2 次?因为在矩阵乘法
推理 FLOPs 公式:
其中:
是前向传播需要的计算量(FLOPs) 是模型的总参数量(Parameters) 是处理的 Token 数量(Tokens)
(注:这里忽略了少量的 Attention 矩阵乘积算力等,因为在大模型中,线性层的矩阵乘法占了绝对大头,通常占 99% 以上)
大模型训练常用近似公式:训练 FLOPs ≈ 6 × 参数量 × token 数。
def calculate_training_flops(num_params_b, num_tokens, flops_per_param_token=6):
return num_params_b * 1_000_000_000 * num_tokens * flops_per_param_tokendef test_calculate_training_flops():
try:
result = calculate_training_flops(7, 1_000_000_000_000)
assert result == 42_000_000_000_000_000_000_000, f"错误:期望 4.2e22,实际 {result}"
result = calculate_training_flops(1, 1_000_000, flops_per_param_token=6)
assert result == 6_000_000_000_000_000, f"错误:期望 6e15,实际 {result}"
print("✅ FLOPs 函数测试通过!")
except AssertionError as e:
print(f"❌ 测试失败: {e}")
except Exception as e:
print(f"❌ 运行错误: {e}")
test_calculate_training_flops()Q2扩展验证:估算训练时间
给定 GPU 的理论算力、GPU 数量和利用率,估算完成训练需要的时间。
def estimate_training_time(num_params_b, num_tokens, gpu_tflops, num_gpus, efficiency=0.35):
total_flops = calculate_training_flops(num_params_b, num_tokens)
effective_flops = gpu_tflops * 1e12 * num_gpus * efficiency
hours = total_flops / effective_flops / 3600
return hoursQ3:训练 (Training) 时包含前向和反向传播,总 FLOPs 是多少?
点击展开查看解析
训练不仅包含前向传播计算损失,还包含反向传播计算梯度。
在反向传播中,我们需要:
- 计算激活值(Activations)的梯度,以便将误差继续向后传(大约需要
FLOPs)。 - 计算权重(Weights)的梯度,用于模型参数更新(大约也需要
FLOPs)。
因此,反向传播的计算量大约是前向传播的 2 倍。
训练 FLOPs 公式:
实战估算: 假设我们要从头预训练一个 LLaMA-7B(70亿参数)模型,训练数据量是 1T(1万亿)个 Tokens。 需要的总理论算力
如果你手里有 1000 张 A100 (每张卡假设实际算力能跑出 150 TFLOPs,即
把参数量、FLOPs 和训练时间串联起来,感受公式在工程中的作用。
scenarios = [
('1x A100 80GB', 312, 1, 0.35),
('8x A100 80GB', 312, 8, 0.35),
('8x H100 80GB', 500, 8, 0.45),
]
print('LLaMA-7B 训练时间粗估(1T tokens):')
print('-' * 70)
for name, gpu_tflops, num_gpus, eff in scenarios:
hours = estimate_training_time(7, 1_000_000_000_000, gpu_tflops, num_gpus, eff)
print(f"{name:<16} {hours:>12.1f} 小时 ({hours/24:>8.1f} 天)")Q4:训练大模型时,什么是算力利用率 (MFU, Model FLOPs Utilization)?
点击展开查看解析
通过前面的 Q3 我们算出了理论所需算力。但在实际工程中,硬件不可能 100% 把所有时间都花在矩阵乘法上。这就引入了 MFU,它是衡量大模型训练工程质量的最核心指标。
- 理论算力 (Peak FLOPs):显卡说明书上写的算力。比如 A100 BF16 理论峰值是 312 TFLOPs(每秒执行 312 万亿次浮点运算)。
- 实际算力 (Observed FLOPs):即我们用
算出的整个训练所需的理论运算量,除以跑完这些步骤所花的实际时间。 - MFU = 实际算力 / 理论算力。
为什么 MFU 很难达到 100%? 因为在真正的训练集群中,存在 Memory-bound (显存墙) 和 Communication (通信瓶颈)。GPU 很多时间在等待数据从内存搬运过来,或者在等其他机器的 All-Reduce 数据传过来,并没有在做有效的乘加运算。
目前顶级的工业界预训练集群,MFU 通常在 40% 到 60% 之间。如果你微调时的 MFU 只有 10%,说明你的代码里存在严重的通信或 IO 阻塞(比如没开梯度累加,或者数据读取成了瓶颈)。
把计算、显存等待和通信等待拆开,看看为什么实际利用率很难接近峰值。
def mfu_from_time_split(compute_ms, memory_wait_ms, comm_wait_ms):
total_ms = compute_ms + memory_wait_ms + comm_wait_ms
if total_ms <= 0:
return {'mfu': 0.0, 'dominant_stall': 'none'}
mfu = compute_ms / total_ms
stalls = {'memory': memory_wait_ms, 'communication': comm_wait_ms}
dominant = max(stalls, key=stalls.get)
return {
'mfu': round(mfu, 3),
'dominant_stall': dominant,
'stall_ratio': round((memory_wait_ms + comm_wait_ms) / total_ms, 3),
}
cases = [
(100, 30, 20),
(100, 60, 40),
(100, 10, 5),
]
for case in cases:
print(case, '->', mfu_from_time_split(*case))
print('MFU is high only when compute dominates the timeline')