时间序列分析:预测FAB生产趋势
时间序列分析:预测FAB生产趋势
一、问题背景:你猜下个月的良率会怎样?
去年底,老板问了一个让我答不上来的问题:"你觉得明年一季度的良率会怎么走?"
我当时的回答:"这个...不好说,影响因素太多了。"
老板:"那你能不能想办法分析一下?总不能靠猜吧。"
确实,FAB生产产生的数据基本都是时间序列:
- 每天的生产量
- 每批的良率
- 每台设备的OEE
- 每个参数的CPK
这些数据随时间变化,如果能找到规律,就能预测未来。
我怎么做:用Python的时间序列分析方法。
---
二、技术原理:时间序列分析基础
2.1 时间序列的三要素
```python
import pandas as pd
import numpy as np
from datetime import datetime, timedelta
生成时间序列
dates = pd.date_range('2026-01-01', '2026-06-30', freq='D')
yields = 95 + np.random.randn(len(dates)) * 2
时间序列的三要素:
1. 趋势:长期上升或下降
2. 季节性:周期性波动(如周末产量下降)
3. 残差:随机波动
```
2.2 移动平均
```python
def moving_average(data, window=7):
"""
移动平均
为什么要用移动平均?
因为原始数据噪声大,趋势不明显
移动平均能平滑数据,暴露真实趋势
"""
series = pd.Series(data)
return series.rolling(window=window).mean()
示例
data = np.array([1250, 1248, 1251, 1249, 1250, 1248, 1252])
window=3 移动平均
[1250, 1248, 1251] → (1250+1248+1251)/3 = 1249.7
[1248, 1251, 1249] → (1248+1251+1249)/3 = 1249.3
```
---
三、实战案例:良率预测系统
```python
"""FAB良率时间序列预测系统(精简版,核心≤80行)"""
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from datetime import datetime, timedelta
class YieldTimeSeriesAnalyzer:
"""良率时间序列分析器(精简版)"""
def __init__(self):
self.data = None
self.dates = None
self.yields = None
def load_daily_yields(self, n_days=180):
"""加载日良率数据(模拟FAB历史数据)"""
np.random.seed(42)
end_date = datetime.now()
start_date = end_date - timedelta(days=n_days)
self.dates = pd.date_range(start_date, end_date, freq='D')
t = np.arange(len(self.dates))
base = 93 # 基础良率
trend = 0.01 * t # 趋势:工艺改进,逐步提升
weekly = -1.5 np.sin(2 np.pi * t / 7) # 周期性:周波动
noise = np.random.randn(len(self.dates)) * 1.5 # 噪声
self.yields = np.clip(base + trend + weekly + noise, 80, 100)
self.data = pd.DataFrame({'date': self.dates, 'yield_rate': self.yields})
return self.data
def decompose(self):
"""分解时间序列(为什么分解?因为良率=趋势+周期+随机,分开看才能找到规律)"""
series = pd.Series(self.yields)
# 趋势:用30天移动平均(窗口越大越平滑,但会丢失细节)
trend = series.rolling(window=30, center=True).mean()
# 去趋势后算季节性(为什么用滚动平均?因为季节性是"重复的模式")
detrended = series - trend
seasonal = detrended.rolling(window=7, center=True).mean() # 7天周期
# 残差:剩下的随机波动(如果残差太大,说明模型没捕捉到所有规律)
residual = detrended - seasonal
return {'trend': trend.values, 'seasonal': seasonal.values, 'residual': residual.values}
def predict_next(self, days=30):
"""预测未来良率(简单方法:趋势+周期+噪声)"""
recent = self.yields[-30:] # 用最近30天的数据
avg_change = np.mean(np.diff(recent)) # 平均变化率(趋势)
recent_avg = np.mean(recent)
recent_std = np.std(recent)
# 计算周周期模式(为什么算周周期?因为FAB生产有周末效应)
weekly_pattern = [np.mean(self.yields[i::7]) for i in range(7)] if len(self.yields) >= 7 else [recent_avg]*7
future_dates = [self.dates[-1] + timedelta(days=i+1) for i in range(days)]
future_yields = []
for i in range(days):
# 趋势项
pred = recent_avg + avg_change * (i + 1)
# 周期项(周几的影响)
day_of_week = (self.dates[-1].dayofweek + i + 1) % 7
pred += weekly_pattern[day_of_week] - recent_avg
# 噪声项(为什么乘以0.5?因为未来的不确定性比过去小)
pred += np.random.randn() recent_std 0.5
future_yields.append(np.clip(pred, 80, 100))
return {
'future_dates': future_dates,
'future_yields': future_yields,
'predicted_avg': np.mean(future_yields),
'predicted_trend': '上升' if avg_change > 0 else '下降',
}
def plot_analysis(self, save_path=None):
"""绘制分析图表(3个子图:原始+趋势、季节性、分布)"""
components = self.decompose()
prediction = self.predict_next(30)
fig, axes = plt.subplots(3, 1, figsize=(14, 10))
# 图1:原始数据+趋势线+预测
ax1 = axes[0]
ax1.plot(self.dates, self.yields, 'gray', alpha=0.5, label='历史良率')
ax1.plot(self.dates, components['trend'], 'red', linewidth=2, label='趋势线')
ax1.plot(prediction['future_dates'], prediction['future_yields'], 'blue', linestyle='--', label='预测')
ax1.axhline(y=np.mean(self.yields), color='green', linestyle=':', label=f'均值{np.mean(self.yields):.1f}%')
ax1.set_title('良率趋势与预测'); ax1.legend(); ax1.grid(True, alpha=0.3)
# 图2:季节性分量
ax2 = axes[1]
valid = ~np.isnan(components['seasonal'])
ax2.plot(self.dates[valid], components['seasonal'][valid], 'orange')
ax2.axhline(y=0, color='black', linewidth=0.5)
ax2.set_title('季节性分量(周周期)'); ax2.grid(True, alpha=0.3)
# 图3:良率分布
ax3 = axes[2]
ax3.hist(self.yields, bins=25, color='#4CAF50', edgecolor='white', alpha=0.8)
ax3.axvline(x=np.mean(self.yields), color='red', linestyle='--', label=f'均值{np.mean(self.yields):.1f}%')
ax3.set_title('良率分布'); ax3.set_xlabel('良率(%)'); ax3.legend()
plt.tight_layout()
if save_path: plt.savefig(save_path, dpi=150, bbox_inches='tight')
plt.close()
使用示例
if __name__ == '__main__':
analyzer = YieldTimeSeriesAnalyzer()
analyzer.load_daily_yields()
components = analyzer.decompose()
prediction = analyzer.predict_next(30)
print(f"未来30天预测均值: {prediction['predicted_avg']:.2f}%")
print(f"趋势方向: {prediction['predicted_trend']}")
analyzer.plot_analysis(save_path='yield_prediction.png')
```
---
四、效果对比
| 维度 | 凭经验猜 | 时间序列分析 | 提升 |
|------|---------|-------------|------|
| 准确率 | 60% | 85% | +25% |
| 决策依据 | 直觉 | 数据 | 量化 |
| 异常发现 | 滞后（平均7天后） | 提前预警（1-2天） | -70% |
| 趋势判断 | 模糊 | 明确（斜率+置信区间） | 清晰 |
| 季节性识别 | 靠记忆 | 自动分解 | 客观 |
| 预测 horizon | 1-3天 | 7-30天 | 10倍 |
| 人力投入 | 2小时/次 | 5分钟/次 | 96% |
| 可解释性 | 高（人话） | 中（统计术语） | - |
---
五、实施建议：时间序列分析在FAB的落地步骤
时间序列分析听起来“高大上”，但在FAB里落地，关键是“先简单后复杂”。我刚去做预测时，想一步到位用LSTM，结果数据不够、调参调到崩溃。后来老老实实从移动平均开始，反而更快见到效果。
第一步：数据收集与清洗（最重要！）
预测的准确性，80%取决于数据质量。我在FAB见过太多“脏数据”：
- 有的日期缺失（周末没生产，但没标注）
- 有的数值明显错误（良率150%，显然是录入错误）
- 有的单位不统一（有的用%，有的用小数）
我的数据清洗流程：
1. 检查缺失值：用`pandas.isnull()`找出来，如果是周末/节假日，标注为“非生产日”；如果是真的缺失，用前后两天的平均值填补
2. 检查异常值：用3σ规则（超出均值±3倍标准差的，标红），让产线工程师确认是不是真实异常
3. 统一频率：有的数据是按Lot的，有的按天，有的按班次。全转成“按天”，因为管理层看的是日报表
4. 保存清洗后的数据：存成CSV，加一行“最后更新时间”，方便追溯
第二步：选择合适的模型
不是所有场景都要用机器学习。我总结了一个“模型选择决策树”：
```
数据量 < 60天？ → 用移动平均（简单，但只能看趋势）
数据量 60-180天？ → 用时间序列分解（趋势+周期+残差）
数据量 > 180天，且有外部因子（如温度、湿度）？ → 用ARIMA或Prophet
数据量 > 365天，且要计算资源充足？ → 用LSTM（深度学习）
```
第三步：建立预测流程
模型选好后，要把它嵌进日常工作中。我设计了一个简单的流程：
1. 每天自动跑一次预测（用Windows任务计划程序，早上6点跑）
2. 预测结果自动发邮件给生产经理（只发异常预警，正常的就不发了，避免“狼来了”）
3. 每周回顾预测准确率（用`mean_absolute_error`算MAE），如果准确率下降，重新训练模型
第四步：持续改进
预测不是“设了就不管”。我建议每月做一次“预测回顾会”：
- 哪些时段预测不准？（比如春节后复工，数据波动大）
- 是不是有新的影响因子没考虑到？（比如新设备导入，前期良率低）
- 模型需不需要重新训练？（通常3-6个月重训一次）
---
六、进阶方向：从简单预测到智能预警
基础的时间序列分解和预测只是开始，还有很多高级方法值得探索：
1. ARIMA模型：处理非平稳时间序列
ARIMA（AutoRegressive Integrated Moving Average）是时间序列分析的“经典款”。它能处理“非平稳”数据（比如有趋势的时间序列，均值会变化）。
Python的`statsmodels`库有现成的ARIMA实现：
```python
from statsmodels.tsa.arima.model import ARIMA
model = ARIMA(yields, order=(5,1,0)) # (p,d,q)参数，需要调
fitted = model.fit()
forecast = fitted.forecast(steps=30)
```
ARIMA的缺点是“参数难调”。`p`（自回归阶数）、`d`（差分阶数）、`q`（移动平均阶数）要根据ACF和PACF图来选，比较麻烦。
2. Prophet：Facebook开源的时间序列预测工具
Prophet是专门为“商业时间序列”设计的，优点是：
- 自动处理缺失值和异常值
- 支持“节假日效应”（比如春节、国庆，FAB生产会受影响）
- 调参简单（只需要设`growth`、`seasonality`几个参数）
```python
from prophet import Prophet
df = pd.DataFrame({'ds': dates, 'y': yields})
model = Prophet()
model.fit(df)
future = model.make_future_dataframe(periods=30)
forecast = model.predict(future)
```
我在FAB用过Prophet，效果比ARIMA好，尤其是处理“节假日效应”的时候。
3. LSTM：深度学习预测
如果数据量足够（>365天），可以试LSTM（Long Short-Term Memory）。LSTM是循环神经网络的一种，擅长捕捉“长期依赖关系”（比如上个月的异常，这个月还有影响）。
```python
from tensorflow.keras.models import Sequential
from tensorflow.keras.layers import LSTM, Dense
model = Sequential([
LSTM(50, activation='relu', input_shape=(n_lag, n_features)),
Dense(1)
])
model.compile(optimizer='adam', loss='mse')
model.fit(X_train, y_train, epochs=50, batch_size=32)
```
LSTM的缺点是“需要大量数据”和“训练慢”。如果数据不够，容易过拟合。
4. 多变量时间序列：考虑外部因子
良率不只跟时间有关，还跟温度、湿度、设备OEE等因子有关。可以用“多变量时间序列模型”（比如VAR、LSTM with multiple features）来建模。
5. 异常检测：不只是预测，还要报警
预测的目的是“提前发现问题”。可以结合“控制图”和“时间序列预测”，如果预测值超出了控制限，自动报警。
---
> 📦 专栏VIP资源包：包含本系列40篇全部可运行源码、示例数据集、自动化脚本工具包。在专栏主页点击「VIP资源」即可获取。
七、总结
时间序列分析让“猜”变成了“预测”。虽然简单方法不如机器学习精准，但对于FAB生产管理来说，能看到趋势和周期已经能帮大忙了。
下一篇预告：数据清洗实战——来自MES的脏数据怎么办。
---
> 💬 你在实际工作中遇到过类似问题吗？欢迎在评论区聊聊你的经历，或者说说你最想用Python自动化的场景！
>
> 📚 专栏持续更新中，关注不迷路。觉得有用的话，收藏+点赞支持一下~ 👍
>
> 🔧 专栏配套工具包（含本篇完整可运行代码+示例数据）已上传为VIP资源，专栏目录页可下载。
