)
1. 從“拍腦袋”到“算數(shù)據(jù)”為什么蔬菜定價補貨值得建模每年國賽C題一出總能引起一陣熱議。2023年的這道“蔬菜類商品的自動定價與補貨決策”乍一看似乎是個超市運營的日常問題離我們熟悉的微分方程、物理仿真有點遠。很多同學(xué)的第一反應(yīng)可能是“這不就是個簡單的預(yù)測和優(yōu)化嗎用個時間序列預(yù)測銷量再套個庫存模型不就行了” 我當(dāng)年剛開始接觸建模時也這么想過結(jié)果被現(xiàn)實狠狠教育了一番。這道題的魅力恰恰在于它用看似平凡的商業(yè)場景包裹了從數(shù)據(jù)處理、不確定性建模、多目標(biāo)決策到方案落地的完整鏈條非?？简灳C合能力。超市里綠葉菜今天賣不完明天就可能蔫掉報廢土豆洋蔥這類耐儲的囤多了又占壓資金。價格定高了顧客不買賬定低了毛利覆蓋不了損耗。補貨更是門學(xué)問補多了損耗和倉儲成本飆升補少了貨架空空錯過銷售機會顧客滿意度下降。傳統(tǒng)上這些決策嚴重依賴店長的個人經(jīng)驗也就是我們常說的“拍腦袋”。今天下雨了感覺菜不好賣就少進點明天有促銷就多備點貨。這種方式在小規(guī)模經(jīng)營時或許有效但在連鎖化、數(shù)字化的當(dāng)下它的粗放和低效就成了致命傷。因此這道題的核心價值是讓我們用數(shù)學(xué)和算法的“望遠鏡”與“顯微鏡”去透視并優(yōu)化這個充滿不確定性的日常商業(yè)過程。它要求我們不再是簡單的“調(diào)包俠”把數(shù)據(jù)扔進模型等結(jié)果而是要深入思考數(shù)據(jù)反映了什么規(guī)律模型背后的假設(shè)是否成立我們的決策建議在真實的、嘈雜的、充滿約束的商業(yè)環(huán)境中是否真的可行接下來我將結(jié)合常見的思路誤區(qū)拆解這道題的幾個核心關(guān)卡并分享一套經(jīng)過實戰(zhàn)檢驗的、可擴展的參考實現(xiàn)框架。2. 破題第一步別急著找算法先理解數(shù)據(jù)和問題邊界看到題目很多人的習(xí)慣是立刻去想“我用什么模型”——LSTM預(yù)測ARIMA還是隨機森林這種思路往往會導(dǎo)致后期陷入僵局。正確的打開方式是像偵探一樣審視題目給出的所有信息定義清楚我們要解決的到底是什么問題。2.1 問題本質(zhì)拆解這不是一個題而是三個題的組合仔細閱讀賽題描述這里基于常見賽題結(jié)構(gòu)進行還原你會發(fā)現(xiàn)它通常包含幾個環(huán)環(huán)相扣的子問題銷量與需求預(yù)測這是所有決策的基礎(chǔ)。給定歷史銷售數(shù)據(jù)、成本、定價等信息預(yù)測未來一段時間如每天、每周各種蔬菜的需求量。注意這里預(yù)測的是“需求”而不是“銷量”。銷量是需求在庫存約束下的實現(xiàn)值。如果貨備少了銷量就等于需求如果貨備多了沒賣完銷量就小于需求。預(yù)測模型必須能區(qū)分這個概念。定價決策在預(yù)測需求的基礎(chǔ)上考慮蔬菜的成本、新鮮度損耗率、市場競爭、顧客價格敏感度等因素制定一個能最大化利潤或兼顧其他目標(biāo)的價格。蔬菜定價的特殊性在于其“易腐性”今天的定價直接影響今天的銷量和明天的剩余庫存質(zhì)量進而影響明天的定價決策這是一個動態(tài)問題。補貨決策在給定定價策略和需求預(yù)測下決定每天或每次的補貨量。這需要平衡多個矛盾的目標(biāo)滿足預(yù)測需求、控制庫存持有成本、最小化因腐爛造成的損耗、利用供應(yīng)商的批量折扣、遵守倉庫容量和貨架空間限制等。這三個問題不是獨立的而是強耦合的。定價影響需求需求預(yù)測指導(dǎo)補貨補貨量又反過來影響庫存成本和損耗進而影響未來的定價空間。因此一個割裂的、先預(yù)測、再定價、最后補貨的流水線式解法往往得不到全局優(yōu)解。2.2 數(shù)據(jù)預(yù)處理清洗與特征工程是模型的“地基”題目通常會提供一段時間內(nèi)如過去一年的詳細數(shù)據(jù)可能包括銷售流水日期、商品編碼、銷售量、銷售單價、銷售金額。商品信息商品編碼、品類葉菜、根莖、茄果等、成本價、供應(yīng)商、保質(zhì)期/損耗特性。庫存記錄每日期初庫存、入庫量、損耗報損量。外部信息可能包含節(jié)假日、天氣溫度、降水量、促銷活動標(biāo)記等。第一步數(shù)據(jù)清洗。檢查并處理缺失值、異常值如某天銷量為0但庫存充足可能是數(shù)據(jù)記錄缺失某天銷量奇高可能是促銷或數(shù)據(jù)錯誤。對于異常值不能簡單刪除要結(jié)合業(yè)務(wù)邏輯判斷。例如國慶節(jié)當(dāng)天銷量暴增是合理的應(yīng)予以保留并作為節(jié)假日特征。第二步構(gòu)建“真實需求”。這是最關(guān)鍵的一步。由于存在缺貨歷史銷量數(shù)據(jù)是“被截斷”的需求。我們可以用一些方法來估算真實需求簡單方法對于缺貨日當(dāng)日銷量等于期初庫存可以將該日需求估算為“銷量 × 一個大于1的系數(shù)”這個系數(shù)可以根據(jù)不缺貨日期的平均銷量與庫存比例來估計。進階方法使用生存分析Survival Analysis中的概念將缺貨視為“刪失”數(shù)據(jù)利用不缺貨日期的數(shù)據(jù)分布來估計整個需求分布。第三步特征工程。將原始數(shù)據(jù)轉(zhuǎn)化為模型能理解的特征。這比選擇模型更重要。對于時間序列預(yù)測有效的特征可能包括時間特征年、月、日、星期幾、是否為周末、是否為節(jié)假日及其前后幾天、季度。滯后特征過去1天、3天、7天、14天、30天的銷量、價格。滾動統(tǒng)計特征過去7天的平均銷量、標(biāo)準(zhǔn)差、最大值、最小值。商品屬性特征品類獨熱編碼、成本價、是否易腐如葉菜為1土豆為0。交互特征價格與品類的交互項不同品類對價格敏感度不同、節(jié)假日與品類的交互項。外部特征天氣數(shù)據(jù)如最高溫、降水量、是否有促銷。注意特征不是越多越好。要防止過擬合可以使用特征重要性評估如基于樹模型或遞歸特征消除來選擇關(guān)鍵特征。3. 模型構(gòu)建預(yù)測、定價與補貨的“鐵三角”理解了問題構(gòu)建了特征我們才進入模型選擇階段。這里提供一個分層、可迭代的建?？蚣?。3.1 需求預(yù)測模型融合時序與回歸的混合思路單純用時間序列模型如ARIMA、Prophet可能忽略了價格、促銷等外部因素。單純用回歸模型如線性回歸、XGBoost可能捕捉不到復(fù)雜的周期性和趨勢。因此混合模型是更穩(wěn)健的選擇。一種有效的架構(gòu)是“殘差學(xué)習(xí)”框架基準(zhǔn)時序預(yù)測首先使用一個穩(wěn)健的時序模型如季節(jié)性分解的STL或Prophet對每個單品的歷史銷量進行擬合得到僅基于時間規(guī)律的預(yù)測值y_temporal。這個模型只使用時間特征。殘差建模計算歷史真實銷量或估算需求與y_temporal的殘差residual y_true - y_temporal。這個殘差包含了時間趨勢無法解釋的部分比如價格變動、促銷、異常天氣的影響。機器學(xué)習(xí)模型擬合殘差使用XGBoost或LightGBM等模型以價格、促銷、天氣、品類等特征為輸入預(yù)測這個殘差值residual_pred。最終預(yù)測最終的預(yù)測需求y_pred y_temporal residual_pred。這種方法的優(yōu)點是時序模型抓住了主體規(guī)律機器學(xué)習(xí)模型負責(zé)精細化修正兩者互補且解釋性較強。在編程實現(xiàn)上我們可以用statsmodels或prophet做時序部分用scikit-learn或xgboost做殘差回歸部分。# 偽代碼示例混合預(yù)測模型框架 import pandas as pd from prophet import Prophet import xgboost as xgb from sklearn.model_selection import train_test_split # 假設(shè) df 是包含日期‘ds’、銷量‘y’、價格‘price’、是否促銷‘promo’等特征的DataFrame df pd.read_csv(sales_data.csv) # 步驟1使用Prophet進行基準(zhǔn)時序預(yù)測 model_prophet Prophet(yearly_seasonalityTrue, weekly_seasonalityTrue, daily_seasonalityFalse) model_prophet.fit(df[[ds, y]]) future model_prophet.make_future_dataframe(periods30) # 預(yù)測未來30天 forecast model_prophet.predict(future) df[y_temporal] forecast[yhat].values[:len(df)] # 將擬合值對齊到歷史數(shù)據(jù) # 步驟2計算殘差 df[residual] df[y] - df[y_temporal] # 步驟3準(zhǔn)備特征預(yù)測殘差 features [price, promo, temperature, is_holiday, day_of_week] # 示例特征 X df[features] y_residual df[residual] X_train, X_test, y_train, y_test train_test_split(X, y_residual, test_size0.2, shuffleFalse) # 注意時間序列避免隨機打亂 model_xgb xgb.XGBRegressor(objectivereg:squarederror, n_estimators100) model_xgb.fit(X_train, y_train) # 步驟4對于未來預(yù)測先有時序預(yù)測再用未來特征預(yù)測殘差 future_features ... # 需要構(gòu)建未來日期的特征如計劃價格、已知節(jié)假日 future_temporal forecast[yhat].values[-30:] # 未來30天的時序預(yù)測 future_residual_pred model_xgb.predict(future_features) final_demand_pred future_temporal future_residual_pred3.2 定價模型在需求彈性與損耗間尋找平衡點定價不是孤立的它建立在需求預(yù)測的基礎(chǔ)上并且需求通常與價格負相關(guān)。我們需要刻畫這種關(guān)系即需求的價格彈性。需求函數(shù)建模對于每種蔬菜或同質(zhì)品類可以假設(shè)一個簡化的需求函數(shù)。最常見的是線性需求函數(shù)和指數(shù)需求函數(shù)。線性需求D(p) a - b * p。其中D是預(yù)測需求量p是定價a是潛在市場規(guī)模b是價格敏感系數(shù)?？梢酝ㄟ^歷史數(shù)據(jù)價格銷量進行線性回歸估計a和b。指數(shù)需求常彈性D(p) A * p^(-e)。取對數(shù)后變?yōu)閘og(D) log(A) - e * log(p)這是一個線性關(guān)系e就是價格彈性系數(shù)通常為正數(shù)。這種形式更符合經(jīng)濟學(xué)常識彈性恒定。單日定價優(yōu)化模型在知道成本c、需求函數(shù)D(p)后單日利潤為π(p) (p - c) * D(p)。通過求導(dǎo)令導(dǎo)數(shù)為零可以解出理論最優(yōu)價格p* (a b*c) / (2b)線性需求或p* (e/(e-1)) * c指數(shù)需求需滿足e1。但蔬菜定價的復(fù)雜性在于損耗蔬菜會腐爛剩余庫存會貶值。假設(shè)每日損耗率為θ如葉菜0.2根莖0.05那么今天的剩余庫存I_tomorrow (I_today Q_today - D_today) * (1-θ)。今天的定價p_today通過影響D_today直接影響I_tomorrow進而影響明天的決策。這變成了一個多期動態(tài)規(guī)劃問題。庫存約束價格不能低到讓需求超過庫存否則會產(chǎn)生無法滿足的需求影響顧客滿意度。多目標(biāo)可能不僅要利潤還要考慮銷量市場份額、損耗率、庫存周轉(zhuǎn)率等。因此更實際的定價模型是一個帶約束的動態(tài)優(yōu)化問題。我們可以將其離散化構(gòu)建一個有限期的決策模型。3.3 補貨模型庫存論與優(yōu)化算法的結(jié)合補貨決策的輸入是預(yù)測需求D和定價策略決定的預(yù)期銷量輸出是補貨量Q。目標(biāo)是平衡多個成本訂貨成本每次下單有固定成本如物流、管理。持有成本庫存占用資金和倉庫空間的成本按天計。缺貨成本需求大于庫存導(dǎo)致的銷售損失和商譽損失。損耗成本因腐爛而丟棄的蔬菜成本。對于不易腐的標(biāo)準(zhǔn)品經(jīng)典的(s, S)或(Q, R)庫存策略很有效。但對于易腐品需要引入損耗率和保質(zhì)期。一種實用的建模方法是構(gòu)建一個有限期的隨機動態(tài)規(guī)劃SDP或近似動態(tài)規(guī)劃ADP模型狀態(tài)變量S_t (I_t, d_t)其中I_t是當(dāng)前庫存水平d_t是當(dāng)天的預(yù)測需求隨機變量有其分布。決策變量補貨量Q_t或定價p_t如果定價也是決策。狀態(tài)轉(zhuǎn)移I_{t1} max(0, I_t Q_t - min(D_t, I_t Q_t)) * (1-θ)。這里min(D_t, I_tQ_t)是實際銷量。即時成本/收益R_t p_t * min(D_t, I_tQ_t) - c * Q_t - h * I_t - v * (I_t * θ)。其中c是進貨成本h是單位持有成本v是單位損耗成本可能等于c。目標(biāo)最大化總期望利潤E[∑ R_t]。直接求解這個SDP在商品種類多、周期長時是“維數(shù)災(zāi)難”。因此我們通常采用近似方法模型預(yù)測控制MPC在每個決策點每天我們只優(yōu)化未來一個較短的時間窗口如未來3-7天的決策并只執(zhí)行第一天的決策。第二天根據(jù)新的實際銷售和庫存狀態(tài)重新滾動優(yōu)化。這大大降低了計算復(fù)雜度且能適應(yīng)不確定性?；诜抡娴膬?yōu)化將補貨策略參數(shù)化例如策略為“當(dāng)庫存低于R時補貨到S水平”然后通過大量隨機需求仿真評估不同(R, S)參數(shù)下的平均利潤選擇最優(yōu)參數(shù)。4. 方案集成與求解從理論模型到可執(zhí)行代碼將預(yù)測、定價、補貨三個模塊集成起來并求解出可執(zhí)行的決策是最后也是最關(guān)鍵的一步。這里容易陷入兩個極端一是追求模型的復(fù)雜和完美導(dǎo)致無法求解二是模型過于簡化脫離實際。4.1 設(shè)計求解框架MPC是實戰(zhàn)中的“瑞士軍刀”我強烈推薦使用模型預(yù)測控制MPC作為整體求解框架。它的思想非常契合這類動態(tài)隨機優(yōu)化問題且易于理解和實現(xiàn)。MPC工作流程如下初始化在第一天開始時擁有當(dāng)前庫存I_0和未來一段時間如7天的外部信息預(yù)測天氣、節(jié)假日。滾動窗口優(yōu)化 a.預(yù)測基于當(dāng)前信息和歷史數(shù)據(jù)使用第3.1節(jié)的混合模型預(yù)測未來N天預(yù)測窗口內(nèi)每天的需求分布可以是一個點估計也可以是多個情景。 b.優(yōu)化建立一個未來N天的確定性或隨機優(yōu)化模型。決策變量是未來N天每天的補貨量Q_t和定價p_t。目標(biāo)函數(shù)是這N天的總期望利潤。約束包括庫存動態(tài)、倉庫容量、補貨上限等。此時我們將未來需求的不確定性用上一步得到的預(yù)測值或幾個代表性情景來近似替代。c.執(zhí)行求解這個優(yōu)化問題得到未來N天的最優(yōu)決策序列{Q_1*, p_1*, Q_2*, p_2*, ..., Q_N*, p_N*}。但只執(zhí)行第一天的決策即按照Q_1*補貨按照p_1*定價。滾動更新進入第二天。觀察到第一天的實際銷量D_1_real和實際損耗更新庫存狀態(tài)I_1。將時間窗口向前滾動一天重復(fù)步驟2。MPC的優(yōu)點在于“反饋校正”。它用最新的實際數(shù)據(jù)來修正模型預(yù)測和優(yōu)化決策對預(yù)測誤差有很強的魯棒性。預(yù)測窗口N的選擇是個權(quán)衡N太短可能目光短淺N太長預(yù)測不準(zhǔn)優(yōu)化問題也難解。對于蔬菜考慮到其易腐性和銷售周期性N3到7天通常是個合理范圍。4.2 優(yōu)化模型求解線性/整數(shù)規(guī)劃與啟發(fā)式算法MPC框架中的核心是第2步的優(yōu)化模型。這個模型通?？梢员硎鰹橐粋€混合整數(shù)線性規(guī)劃MILP問題。模型要素示例集合天數(shù)t in {1,...,N}商品i in {1,...,M}。參數(shù)d_{i,t}商品i在第t天的預(yù)測需求點估計。c_i商品i的進貨成本。p_{i,t}決策變量商品i在第t天的售價或作為參數(shù)如果定價策略已定。h_i單位持有成本。θ_i日損耗率。Cap倉庫總?cè)萘?。Shelf_i商品i的貨架空間。決策變量Q_{i,t}連續(xù)變量商品i在第t天的補貨量。I_{i,t}連續(xù)變量商品i在第t天結(jié)束時的庫存。Sales_{i,t}連續(xù)變量商品i在第t天的實際銷量≤ min(需求 可用庫存)。目標(biāo)函數(shù)最大化總利潤∑_t ∑_i (p_{i,t} * Sales_{i,t} - c_i * Q_{i,t} - h_i * I_{i,t})。約束條件庫存平衡I_{i,t} (I_{i,t-1} Q_{i,t} - Sales_{i,t}) * (1 - θ_i)。銷量限制Sales_{i,t} ≤ d_{i,t}且Sales_{i,t} ≤ I_{i,t-1} Q_{i,t}。容量約束∑_i I_{i,t} ≤ Cap和I_{i,t} ≤ Shelf_i。非負與邏輯約束所有變量≥0。這個MILP模型可以用專業(yè)的優(yōu)化求解器如PuLP調(diào)用CBC、ortools、Gurobi、CPLEX來求解。對于大規(guī)模問題商品數(shù)M很多求解可能較慢這時可以分解按品類或關(guān)聯(lián)性不強的商品分組獨立優(yōu)化。啟發(fā)式規(guī)則先用簡單的經(jīng)驗規(guī)則如報童模型計算出每個單品的大致補貨量再用優(yōu)化模型進行微調(diào)和容量協(xié)調(diào)。遺傳算法/模擬退火如果模型非線性程度高可以嘗試元啟發(fā)式算法尋找滿意解。4.3 參考代碼結(jié)構(gòu)一個模塊化的實現(xiàn)藍圖下面給出一個高度概括的、模塊化的Python代碼結(jié)構(gòu)展示了如何將上述思路串聯(lián)起來。這只是一個框架每個函數(shù)都需要根據(jù)具體數(shù)據(jù)和模型進行填充。import pandas as pd import numpy as np from prophet import Prophet import xgboost as xgb from pulp import LpProblem, LpVariable, lpSum, LpMaximize, LpStatus, PULP_CBC_CMD # 或者使用 ortools.linear_solver import pywraplp class VegetablePricingReplenishmentSystem: def __init__(self, historical_data, cost_data, config): 初始化系統(tǒng)。 historical_data: 歷史銷售、庫存、價格、特征數(shù)據(jù)。 cost_data: 商品成本、持有成本、損耗率等。 config: 配置參數(shù)如預(yù)測窗口、優(yōu)化窗口、倉庫容量等。 self.hist_data historical_data self.cost_data cost_data self.config config self.demand_models {} # 存儲每個商品的需求預(yù)測模型 self.current_inventory self._get_initial_inventory() def _build_demand_forecast_model(self, product_id): 為單個商品構(gòu)建混合需求預(yù)測模型如3.1節(jié)所述。 # 1. 提取該商品歷史數(shù)據(jù) product_data self.hist_data[self.hist_data[product_id] product_id].copy() # 2. 使用時序模型如Prophet擬合基準(zhǔn) # 3. 計算殘差 # 4. 使用XGBoost擬合殘差 # 5. 保存模型到 self.demand_models[product_id] pass def forecast_demand(self, product_id, future_features_df): 給定未來特征預(yù)測商品在未來N天的需求。 if product_id not in self.demand_models: self._build_demand_forecast_model(product_id) model self.demand_models[product_id] # 使用時序部分預(yù)測 temporal_pred model[temporal].predict(future_features_df[[ds]]) # 用機器學(xué)習(xí)部分預(yù)測殘差 residual_pred model[ml].predict(future_features_df[model[feature_cols]]) final_pred temporal_pred residual_pred return final_pred def mpc_optimization(self, current_date, horizon7): 執(zhí)行一次MPC滾動優(yōu)化返回當(dāng)天決策。 decisions {} all_products self.cost_data[product_id].unique() # 步驟1預(yù)測未來horizon天的需求 future_dates pd.date_range(startcurrent_date, periodshorizon, freqD) future_features self._generate_future_features(future_dates) # 構(gòu)建未來特征已知節(jié)假日、計劃天氣等 demand_forecast {} for pid in all_products: demand_forecast[pid] self.forecast_demand(pid, future_features) # 步驟2構(gòu)建并求解未來horizon天的優(yōu)化模型 prob LpProblem(Daily_Replenishment_Pricing, LpMaximize) # 定義決策變量字典 Q LpVariable.dicts(Q, (all_products, range(horizon)), lowBound0) # 補貨量 I LpVariable.dicts(I, (all_products, range(horizon1)), lowBound0) # 庫存多一期用于初始 S LpVariable.dicts(S, (all_products, range(horizon)), lowBound0) # 銷量 # 設(shè)置初始庫存 for pid in all_products: I[pid][0] self.current_inventory.get(pid, 0) # 目標(biāo)函數(shù)總利潤 profit_expr [] for t in range(horizon): for pid in all_products: cost self.cost_data.loc[pid, unit_cost] price ... # 這里價格可以是決策變量也可以是需求函數(shù)反推的結(jié)果簡化起見可以先固定或根據(jù)規(guī)則設(shè)定 holding_cost self.cost_data.loc[pid, holding_cost] profit_expr.append(price * S[pid][t] - cost * Q[pid][t] - holding_cost * I[pid][t]) prob lpSum(profit_expr) # 約束條件 for t in range(horizon): # 庫存平衡約束 for pid in all_products: spoilage_rate self.cost_data.loc[pid, spoilage_rate] prob I[pid][t1] (I[pid][t] Q[pid][t] - S[pid][t]) * (1 - spoilage_rate) # 銷量約束不能超過預(yù)測需求和可用庫存 for pid in all_products: prob S[pid][t] demand_forecast[pid][t] prob S[pid][t] I[pid][t] Q[pid][t] # 倉庫容量約束簡化示例只考慮總庫存體積 total_volume lpSum(I[pid][t1] * self.cost_data.loc[pid, volume_per_unit] for pid in all_products) prob total_volume self.config[warehouse_capacity] # 求解 solver PULP_CBC_CMD(msgFalse) prob.solve(solver) # 步驟3提取并返回第一天的決策 if LpStatus[prob.status] Optimal: for pid in all_products: decisions[pid] { replenish: Q[pid][0].varValue, price: ... # 從解中獲取或根據(jù)規(guī)則計算 } else: # 如果無解啟用備用啟發(fā)式規(guī)則 decisions self._fallback_heuristic(all_products, demand_forecast) return decisions def _fallback_heuristic(self, product_ids, demand_forecast): 備用啟發(fā)式規(guī)則例如采用報童模型計算補貨量價格采用成本加成法。 decisions {} for pid in product_ids: # 簡單報童模型補貨量 預(yù)測需求 安全庫存 forecast_d demand_forecast[pid][0] safety_stock ... # 根據(jù)歷史需求波動計算 decisions[pid] { replenish: max(0, forecast_d safety_stock - self.current_inventory.get(pid, 0)), price: self.cost_data.loc[pid, unit_cost] * 1.3 # 成本加成30% } return decisions def run_daily_operation(self, start_date, end_date): 模擬系統(tǒng)從start_date到end_date的每日運行。 current_date start_date decision_log [] while current_date end_date: # 1. 執(zhí)行MPC優(yōu)化得到今日決策 daily_decisions self.mpc_optimization(current_date) # 2. 模擬執(zhí)行決策在實際系統(tǒng)中這是真實操作 actual_sales self._simulate_sales(current_date, daily_decisions) # 3. 更新庫存狀態(tài)考慮實際銷量和損耗 self._update_inventory(current_date, daily_decisions, actual_sales) # 4. 記錄日志 decision_log.append({date: current_date, decisions: daily_decisions, actual_sales: actual_sales}) # 5. 進入下一天 current_date pd.Timedelta(days1) return pd.DataFrame(decision_log) # 主程序入口 if __name__ __main__: # 加載數(shù)據(jù) hist_df pd.read_csv(historical_sales.csv) cost_df pd.read_csv(product_costs.csv).set_index(product_id) config {warehouse_capacity: 10000, forecast_horizon: 7} # 初始化系統(tǒng) system VegetablePricingReplenishmentSystem(hist_df, cost_df, config) # 運行模擬 start_date pd.Timestamp(2023-06-01) end_date pd.Timestamp(2023-06-30) operation_log system.run_daily_operation(start_date, end_date) print(operation_log.tail())5. 模型檢驗、靈敏度分析與論文呈現(xiàn)要點建好模型和系統(tǒng)后如何驗證其有效性并在論文中清晰呈現(xiàn)是拿高分的關(guān)鍵。5.1 模型檢驗回溯測試與關(guān)鍵績效指標(biāo)KPI絕不能只用訓(xùn)練集上的表現(xiàn)來評價模型。必須進行嚴格的回溯測試。方法將歷史數(shù)據(jù)按時間劃分為訓(xùn)練集和測試集例如用前10個月訓(xùn)練后2個月測試。在測試集上模擬系統(tǒng)逐日運行每天開始時系統(tǒng)只知道過去的信息用它來做預(yù)測和優(yōu)化決策然后將決策與當(dāng)天的“真實”數(shù)據(jù)在測試集中進行比較。需要計算的KPI包括利潤相關(guān)總利潤、平均毛利率、投資回報率。庫存相關(guān)平均庫存水平、庫存周轉(zhuǎn)天數(shù)、缺貨率需求大于庫存的天數(shù)比例。損耗相關(guān)總損耗金額、損耗率損耗量/總進貨量。預(yù)測準(zhǔn)確性使用均方根誤差RMSE、平均絕對百分比誤差MAPE評估需求預(yù)測的精度。特別注意評估在庫存不足可能缺貨的那些日子的預(yù)測誤差這更能反映模型估算“真實需求”的能力。將你的模型與基準(zhǔn)策略對比例如簡單策略每天補貨到固定水平如過去7天平均銷量按成本固定比例加成定價。經(jīng)驗策略模擬有經(jīng)驗的店長的決策規(guī)則可從數(shù)據(jù)中歸納如周末多進20%。 只有顯著、穩(wěn)定地超越了基準(zhǔn)策略你的模型才算有價值。5.2 靈敏度分析找出模型的“脆弱點”和“穩(wěn)健性”模型依賴于許多假設(shè)和參數(shù)如損耗率θ、價格彈性e、持有成本h。靈敏度分析就是觀察當(dāng)這些參數(shù)在合理范圍內(nèi)變動時模型性能如總利潤的變化情況。如何做單因素分析固定其他參數(shù)讓一個參數(shù)如損耗率在[θ-Δθ, θΔθ]范圍內(nèi)變化運行多次模擬觀察利潤的變化曲線。如果利潤對該參數(shù)變化非常敏感說明模型在該參數(shù)上很脆弱需要更精確地估計該參數(shù)或者在策略中加入魯棒性設(shè)計。多因素分析可以使用拉丁超立方抽樣等方法同時擾動多個參數(shù)觀察利潤的分布。這能評估模型在參數(shù)不確定性下的整體穩(wěn)健性。場景分析模擬一些極端但可能發(fā)生的場景例如連續(xù)三天暴雨導(dǎo)致客流減少50%某主要供應(yīng)商突然斷供一周。你的系統(tǒng)決策會產(chǎn)生什么變化利潤會下降多少這能體現(xiàn)系統(tǒng)的抗風(fēng)險能力。在論文中用圖表清晰地展示靈敏度分析的結(jié)果例如畫出利潤隨關(guān)鍵參數(shù)變化的折線圖并給出文字解釋“如圖所示當(dāng)損耗率估計誤差在±20%內(nèi)時系統(tǒng)利潤波動在5%以內(nèi)表明模型對損耗率參數(shù)具有一定的魯棒性?！?.3 論文寫作與可視化講好一個邏輯閉環(huán)的故事國賽論文評審時間緊清晰的邏輯和直觀的可視化至關(guān)重要。行文邏輯建議問題重述與分析不要照抄題目要用自己的話精煉概括并明確指出問題的核心矛盾如利潤最大化與損耗最小化的權(quán)衡和子問題劃分。模型準(zhǔn)備詳細闡述數(shù)據(jù)預(yù)處理過程特別是“真實需求”估算、特征工程方法。這是體現(xiàn)你工作扎實程度的地方。模型建立分節(jié)闡述預(yù)測模型、定價模型、補貨模型。對每個模型都要說明為什么選擇這個模型與替代模型相比的優(yōu)勢以及模型如何與問題結(jié)合公式、變量含義要解釋清楚。模型求解重點介紹MPC框架和優(yōu)化模型的求解方法。如果是用求解器說明是什么求解器如果用了啟發(fā)式算法描述算法步驟。模型檢驗與靈敏度分析展示回溯測試的結(jié)果表格和對比圖用數(shù)據(jù)證明模型的有效性。展示靈敏度分析的圖表和結(jié)論。模型評價與推廣客觀評價模型的優(yōu)點和局限性例如假設(shè)需求函數(shù)形式固定未考慮競爭對手突然降價等。提出可能的改進方向并簡要說明模型思想可以推廣到其他易腐品零售如水果、鮮花、烘焙食品?？梢暬c系統(tǒng)框架圖用清晰的框圖展示“數(shù)據(jù)輸入 - 預(yù)測模塊 - 優(yōu)化模塊 - 決策輸出 - 滾動更新”的整個流程。預(yù)測效果圖對于幾個代表性商品畫出歷史期真實銷量、模型擬合銷量以及測試期預(yù)測銷量與真實銷量的對比折線圖。決策模擬圖選擇一段時間如測試期兩周畫出庫存水平、補貨量、價格、銷量的聯(lián)動變化圖直觀展示模型決策的動態(tài)過程。靈敏度分析圖用折線圖或熱力圖展示關(guān)鍵參數(shù)變化對核心KPI的影響。結(jié)果對比圖用柱狀圖對比你的模型與基準(zhǔn)策略在多個KPI上的表現(xiàn)。記住論文的本質(zhì)是向評委講述你發(fā)現(xiàn)問題、分析問題、解決問題、驗證方案的完整思考過程。代碼是工具模型是手段清晰的邏輯和令人信服的證據(jù)鏈才是獲得高分的關(guān)鍵。這道C題就像一個微縮的商業(yè)智能系統(tǒng)項目完整地走一遍無論比賽結(jié)果如何對解決實際問題的能力都是極大的提升。 基于MPC與混合模型的蔬菜定價補貨決策建模實戰(zhàn)
1. 從“拍腦袋”到“算數(shù)據(jù)”為什么蔬菜定價補貨值得建模每年國賽C題一出總能引起一陣熱議。2023年的這道“蔬菜類商品的自動定價與補貨決策”乍一看似乎是個超市運營的日常問題離我們熟悉的微分方程、物理仿真有點遠。很多同學(xué)的第一反應(yīng)可能是“這不就是個簡單的預(yù)測和優(yōu)化嗎用個時間序列預(yù)測銷量再套個庫存模型不就行了” 我當(dāng)年剛開始接觸建模時也這么想過結(jié)果被現(xiàn)實狠狠教育了一番。這道題的魅力恰恰在于它用看似平凡的商業(yè)場景包裹了從數(shù)據(jù)處理、不確定性建模、多目標(biāo)決策到方案落地的完整鏈條非?？简灳C合能力。超市里綠葉菜今天賣不完明天就可能蔫掉報廢土豆洋蔥這類耐儲的囤多了又占壓資金。價格定高了顧客不買賬定低了毛利覆蓋不了損耗。補貨更是門學(xué)問補多了損耗和倉儲成本飆升補少了貨架空空錯過銷售機會顧客滿意度下降。傳統(tǒng)上這些決策嚴重依賴店長的個人經(jīng)驗也就是我們常說的“拍腦袋”。今天下雨了感覺菜不好賣就少進點明天有促銷就多備點貨。這種方式在小規(guī)模經(jīng)營時或許有效但在連鎖化、數(shù)字化的當(dāng)下它的粗放和低效就成了致命傷。因此這道題的核心價值是讓我們用數(shù)學(xué)和算法的“望遠鏡”與“顯微鏡”去透視并優(yōu)化這個充滿不確定性的日常商業(yè)過程。它要求我們不再是簡單的“調(diào)包俠”把數(shù)據(jù)扔進模型等結(jié)果而是要深入思考數(shù)據(jù)反映了什么規(guī)律模型背后的假設(shè)是否成立我們的決策建議在真實的、嘈雜的、充滿約束的商業(yè)環(huán)境中是否真的可行接下來我將結(jié)合常見的思路誤區(qū)拆解這道題的幾個核心關(guān)卡并分享一套經(jīng)過實戰(zhàn)檢驗的、可擴展的參考實現(xiàn)框架。2. 破題第一步別急著找算法先理解數(shù)據(jù)和問題邊界看到題目很多人的習(xí)慣是立刻去想“我用什么模型”——LSTM預(yù)測ARIMA還是隨機森林這種思路往往會導(dǎo)致后期陷入僵局。正確的打開方式是像偵探一樣審視題目給出的所有信息定義清楚我們要解決的到底是什么問題。2.1 問題本質(zhì)拆解這不是一個題而是三個題的組合仔細閱讀賽題描述這里基于常見賽題結(jié)構(gòu)進行還原你會發(fā)現(xiàn)它通常包含幾個環(huán)環(huán)相扣的子問題銷量與需求預(yù)測這是所有決策的基礎(chǔ)。給定歷史銷售數(shù)據(jù)、成本、定價等信息預(yù)測未來一段時間如每天、每周各種蔬菜的需求量。注意這里預(yù)測的是“需求”而不是“銷量”。銷量是需求在庫存約束下的實現(xiàn)值。如果貨備少了銷量就等于需求如果貨備多了沒賣完銷量就小于需求。預(yù)測模型必須能區(qū)分這個概念。定價決策在預(yù)測需求的基礎(chǔ)上考慮蔬菜的成本、新鮮度損耗率、市場競爭、顧客價格敏感度等因素制定一個能最大化利潤或兼顧其他目標(biāo)的價格。蔬菜定價的特殊性在于其“易腐性”今天的定價直接影響今天的銷量和明天的剩余庫存質(zhì)量進而影響明天的定價決策這是一個動態(tài)問題。補貨決策在給定定價策略和需求預(yù)測下決定每天或每次的補貨量。這需要平衡多個矛盾的目標(biāo)滿足預(yù)測需求、控制庫存持有成本、最小化因腐爛造成的損耗、利用供應(yīng)商的批量折扣、遵守倉庫容量和貨架空間限制等。這三個問題不是獨立的而是強耦合的。定價影響需求需求預(yù)測指導(dǎo)補貨補貨量又反過來影響庫存成本和損耗進而影響未來的定價空間。因此一個割裂的、先預(yù)測、再定價、最后補貨的流水線式解法往往得不到全局優(yōu)解。2.2 數(shù)據(jù)預(yù)處理清洗與特征工程是模型的“地基”題目通常會提供一段時間內(nèi)如過去一年的詳細數(shù)據(jù)可能包括銷售流水日期、商品編碼、銷售量、銷售單價、銷售金額。商品信息商品編碼、品類葉菜、根莖、茄果等、成本價、供應(yīng)商、保質(zhì)期/損耗特性。庫存記錄每日期初庫存、入庫量、損耗報損量。外部信息可能包含節(jié)假日、天氣溫度、降水量、促銷活動標(biāo)記等。第一步數(shù)據(jù)清洗。檢查并處理缺失值、異常值如某天銷量為0但庫存充足可能是數(shù)據(jù)記錄缺失某天銷量奇高可能是促銷或數(shù)據(jù)錯誤。對于異常值不能簡單刪除要結(jié)合業(yè)務(wù)邏輯判斷。例如國慶節(jié)當(dāng)天銷量暴增是合理的應(yīng)予以保留并作為節(jié)假日特征。第二步構(gòu)建“真實需求”。這是最關(guān)鍵的一步。由于存在缺貨歷史銷量數(shù)據(jù)是“被截斷”的需求。我們可以用一些方法來估算真實需求簡單方法對于缺貨日當(dāng)日銷量等于期初庫存可以將該日需求估算為“銷量 × 一個大于1的系數(shù)”這個系數(shù)可以根據(jù)不缺貨日期的平均銷量與庫存比例來估計。進階方法使用生存分析Survival Analysis中的概念將缺貨視為“刪失”數(shù)據(jù)利用不缺貨日期的數(shù)據(jù)分布來估計整個需求分布。第三步特征工程。將原始數(shù)據(jù)轉(zhuǎn)化為模型能理解的特征。這比選擇模型更重要。對于時間序列預(yù)測有效的特征可能包括時間特征年、月、日、星期幾、是否為周末、是否為節(jié)假日及其前后幾天、季度。滯后特征過去1天、3天、7天、14天、30天的銷量、價格。滾動統(tǒng)計特征過去7天的平均銷量、標(biāo)準(zhǔn)差、最大值、最小值。商品屬性特征品類獨熱編碼、成本價、是否易腐如葉菜為1土豆為0。交互特征價格與品類的交互項不同品類對價格敏感度不同、節(jié)假日與品類的交互項。外部特征天氣數(shù)據(jù)如最高溫、降水量、是否有促銷。注意特征不是越多越好。要防止過擬合可以使用特征重要性評估如基于樹模型或遞歸特征消除來選擇關(guān)鍵特征。3. 模型構(gòu)建預(yù)測、定價與補貨的“鐵三角”理解了問題構(gòu)建了特征我們才進入模型選擇階段。這里提供一個分層、可迭代的建?？蚣?。3.1 需求預(yù)測模型融合時序與回歸的混合思路單純用時間序列模型如ARIMA、Prophet可能忽略了價格、促銷等外部因素。單純用回歸模型如線性回歸、XGBoost可能捕捉不到復(fù)雜的周期性和趨勢。因此混合模型是更穩(wěn)健的選擇。一種有效的架構(gòu)是“殘差學(xué)習(xí)”框架基準(zhǔn)時序預(yù)測首先使用一個穩(wěn)健的時序模型如季節(jié)性分解的STL或Prophet對每個單品的歷史銷量進行擬合得到僅基于時間規(guī)律的預(yù)測值y_temporal。這個模型只使用時間特征。殘差建模計算歷史真實銷量或估算需求與y_temporal的殘差residual y_true - y_temporal。這個殘差包含了時間趨勢無法解釋的部分比如價格變動、促銷、異常天氣的影響。機器學(xué)習(xí)模型擬合殘差使用XGBoost或LightGBM等模型以價格、促銷、天氣、品類等特征為輸入預(yù)測這個殘差值residual_pred。最終預(yù)測最終的預(yù)測需求y_pred y_temporal residual_pred。這種方法的優(yōu)點是時序模型抓住了主體規(guī)律機器學(xué)習(xí)模型負責(zé)精細化修正兩者互補且解釋性較強。在編程實現(xiàn)上我們可以用statsmodels或prophet做時序部分用scikit-learn或xgboost做殘差回歸部分。# 偽代碼示例混合預(yù)測模型框架 import pandas as pd from prophet import Prophet import xgboost as xgb from sklearn.model_selection import train_test_split # 假設(shè) df 是包含日期‘ds’、銷量‘y’、價格‘price’、是否促銷‘promo’等特征的DataFrame df pd.read_csv(sales_data.csv) # 步驟1使用Prophet進行基準(zhǔn)時序預(yù)測 model_prophet Prophet(yearly_seasonalityTrue, weekly_seasonalityTrue, daily_seasonalityFalse) model_prophet.fit(df[[ds, y]]) future model_prophet.make_future_dataframe(periods30) # 預(yù)測未來30天 forecast model_prophet.predict(future) df[y_temporal] forecast[yhat].values[:len(df)] # 將擬合值對齊到歷史數(shù)據(jù) # 步驟2計算殘差 df[residual] df[y] - df[y_temporal] # 步驟3準(zhǔn)備特征預(yù)測殘差 features [price, promo, temperature, is_holiday, day_of_week] # 示例特征 X df[features] y_residual df[residual] X_train, X_test, y_train, y_test train_test_split(X, y_residual, test_size0.2, shuffleFalse) # 注意時間序列避免隨機打亂 model_xgb xgb.XGBRegressor(objectivereg:squarederror, n_estimators100) model_xgb.fit(X_train, y_train) # 步驟4對于未來預(yù)測先有時序預(yù)測再用未來特征預(yù)測殘差 future_features ... # 需要構(gòu)建未來日期的特征如計劃價格、已知節(jié)假日 future_temporal forecast[yhat].values[-30:] # 未來30天的時序預(yù)測 future_residual_pred model_xgb.predict(future_features) final_demand_pred future_temporal future_residual_pred3.2 定價模型在需求彈性與損耗間尋找平衡點定價不是孤立的它建立在需求預(yù)測的基礎(chǔ)上并且需求通常與價格負相關(guān)。我們需要刻畫這種關(guān)系即需求的價格彈性。需求函數(shù)建模對于每種蔬菜或同質(zhì)品類可以假設(shè)一個簡化的需求函數(shù)。最常見的是線性需求函數(shù)和指數(shù)需求函數(shù)。線性需求D(p) a - b * p。其中D是預(yù)測需求量p是定價a是潛在市場規(guī)模b是價格敏感系數(shù)?？梢酝ㄟ^歷史數(shù)據(jù)價格銷量進行線性回歸估計a和b。指數(shù)需求常彈性D(p) A * p^(-e)。取對數(shù)后變?yōu)閘og(D) log(A) - e * log(p)這是一個線性關(guān)系e就是價格彈性系數(shù)通常為正數(shù)。這種形式更符合經(jīng)濟學(xué)常識彈性恒定。單日定價優(yōu)化模型在知道成本c、需求函數(shù)D(p)后單日利潤為π(p) (p - c) * D(p)。通過求導(dǎo)令導(dǎo)數(shù)為零可以解出理論最優(yōu)價格p* (a b*c) / (2b)線性需求或p* (e/(e-1)) * c指數(shù)需求需滿足e1。但蔬菜定價的復(fù)雜性在于損耗蔬菜會腐爛剩余庫存會貶值。假設(shè)每日損耗率為θ如葉菜0.2根莖0.05那么今天的剩余庫存I_tomorrow (I_today Q_today - D_today) * (1-θ)。今天的定價p_today通過影響D_today直接影響I_tomorrow進而影響明天的決策。這變成了一個多期動態(tài)規(guī)劃問題。庫存約束價格不能低到讓需求超過庫存否則會產(chǎn)生無法滿足的需求影響顧客滿意度。多目標(biāo)可能不僅要利潤還要考慮銷量市場份額、損耗率、庫存周轉(zhuǎn)率等。因此更實際的定價模型是一個帶約束的動態(tài)優(yōu)化問題。我們可以將其離散化構(gòu)建一個有限期的決策模型。3.3 補貨模型庫存論與優(yōu)化算法的結(jié)合補貨決策的輸入是預(yù)測需求D和定價策略決定的預(yù)期銷量輸出是補貨量Q。目標(biāo)是平衡多個成本訂貨成本每次下單有固定成本如物流、管理。持有成本庫存占用資金和倉庫空間的成本按天計。缺貨成本需求大于庫存導(dǎo)致的銷售損失和商譽損失。損耗成本因腐爛而丟棄的蔬菜成本。對于不易腐的標(biāo)準(zhǔn)品經(jīng)典的(s, S)或(Q, R)庫存策略很有效。但對于易腐品需要引入損耗率和保質(zhì)期。一種實用的建模方法是構(gòu)建一個有限期的隨機動態(tài)規(guī)劃SDP或近似動態(tài)規(guī)劃ADP模型狀態(tài)變量S_t (I_t, d_t)其中I_t是當(dāng)前庫存水平d_t是當(dāng)天的預(yù)測需求隨機變量有其分布。決策變量補貨量Q_t或定價p_t如果定價也是決策。狀態(tài)轉(zhuǎn)移I_{t1} max(0, I_t Q_t - min(D_t, I_t Q_t)) * (1-θ)。這里min(D_t, I_tQ_t)是實際銷量。即時成本/收益R_t p_t * min(D_t, I_tQ_t) - c * Q_t - h * I_t - v * (I_t * θ)。其中c是進貨成本h是單位持有成本v是單位損耗成本可能等于c。目標(biāo)最大化總期望利潤E[∑ R_t]。直接求解這個SDP在商品種類多、周期長時是“維數(shù)災(zāi)難”。因此我們通常采用近似方法模型預(yù)測控制MPC在每個決策點每天我們只優(yōu)化未來一個較短的時間窗口如未來3-7天的決策并只執(zhí)行第一天的決策。第二天根據(jù)新的實際銷售和庫存狀態(tài)重新滾動優(yōu)化。這大大降低了計算復(fù)雜度且能適應(yīng)不確定性?；诜抡娴膬?yōu)化將補貨策略參數(shù)化例如策略為“當(dāng)庫存低于R時補貨到S水平”然后通過大量隨機需求仿真評估不同(R, S)參數(shù)下的平均利潤選擇最優(yōu)參數(shù)。4. 方案集成與求解從理論模型到可執(zhí)行代碼將預(yù)測、定價、補貨三個模塊集成起來并求解出可執(zhí)行的決策是最后也是最關(guān)鍵的一步。這里容易陷入兩個極端一是追求模型的復(fù)雜和完美導(dǎo)致無法求解二是模型過于簡化脫離實際。4.1 設(shè)計求解框架MPC是實戰(zhàn)中的“瑞士軍刀”我強烈推薦使用模型預(yù)測控制MPC作為整體求解框架。它的思想非常契合這類動態(tài)隨機優(yōu)化問題且易于理解和實現(xiàn)。MPC工作流程如下初始化在第一天開始時擁有當(dāng)前庫存I_0和未來一段時間如7天的外部信息預(yù)測天氣、節(jié)假日。滾動窗口優(yōu)化 a.預(yù)測基于當(dāng)前信息和歷史數(shù)據(jù)使用第3.1節(jié)的混合模型預(yù)測未來N天預(yù)測窗口內(nèi)每天的需求分布可以是一個點估計也可以是多個情景。 b.優(yōu)化建立一個未來N天的確定性或隨機優(yōu)化模型。決策變量是未來N天每天的補貨量Q_t和定價p_t。目標(biāo)函數(shù)是這N天的總期望利潤。約束包括庫存動態(tài)、倉庫容量、補貨上限等。此時我們將未來需求的不確定性用上一步得到的預(yù)測值或幾個代表性情景來近似替代。c.執(zhí)行求解這個優(yōu)化問題得到未來N天的最優(yōu)決策序列{Q_1*, p_1*, Q_2*, p_2*, ..., Q_N*, p_N*}。但只執(zhí)行第一天的決策即按照Q_1*補貨按照p_1*定價。滾動更新進入第二天。觀察到第一天的實際銷量D_1_real和實際損耗更新庫存狀態(tài)I_1。將時間窗口向前滾動一天重復(fù)步驟2。MPC的優(yōu)點在于“反饋校正”。它用最新的實際數(shù)據(jù)來修正模型預(yù)測和優(yōu)化決策對預(yù)測誤差有很強的魯棒性。預(yù)測窗口N的選擇是個權(quán)衡N太短可能目光短淺N太長預(yù)測不準(zhǔn)優(yōu)化問題也難解。對于蔬菜考慮到其易腐性和銷售周期性N3到7天通常是個合理范圍。4.2 優(yōu)化模型求解線性/整數(shù)規(guī)劃與啟發(fā)式算法MPC框架中的核心是第2步的優(yōu)化模型。這個模型通?？梢员硎鰹橐粋€混合整數(shù)線性規(guī)劃MILP問題。模型要素示例集合天數(shù)t in {1,...,N}商品i in {1,...,M}。參數(shù)d_{i,t}商品i在第t天的預(yù)測需求點估計。c_i商品i的進貨成本。p_{i,t}決策變量商品i在第t天的售價或作為參數(shù)如果定價策略已定。h_i單位持有成本。θ_i日損耗率。Cap倉庫總?cè)萘?。Shelf_i商品i的貨架空間。決策變量Q_{i,t}連續(xù)變量商品i在第t天的補貨量。I_{i,t}連續(xù)變量商品i在第t天結(jié)束時的庫存。Sales_{i,t}連續(xù)變量商品i在第t天的實際銷量≤ min(需求 可用庫存)。目標(biāo)函數(shù)最大化總利潤∑_t ∑_i (p_{i,t} * Sales_{i,t} - c_i * Q_{i,t} - h_i * I_{i,t})。約束條件庫存平衡I_{i,t} (I_{i,t-1} Q_{i,t} - Sales_{i,t}) * (1 - θ_i)。銷量限制Sales_{i,t} ≤ d_{i,t}且Sales_{i,t} ≤ I_{i,t-1} Q_{i,t}。容量約束∑_i I_{i,t} ≤ Cap和I_{i,t} ≤ Shelf_i。非負與邏輯約束所有變量≥0。這個MILP模型可以用專業(yè)的優(yōu)化求解器如PuLP調(diào)用CBC、ortools、Gurobi、CPLEX來求解。對于大規(guī)模問題商品數(shù)M很多求解可能較慢這時可以分解按品類或關(guān)聯(lián)性不強的商品分組獨立優(yōu)化。啟發(fā)式規(guī)則先用簡單的經(jīng)驗規(guī)則如報童模型計算出每個單品的大致補貨量再用優(yōu)化模型進行微調(diào)和容量協(xié)調(diào)。遺傳算法/模擬退火如果模型非線性程度高可以嘗試元啟發(fā)式算法尋找滿意解。4.3 參考代碼結(jié)構(gòu)一個模塊化的實現(xiàn)藍圖下面給出一個高度概括的、模塊化的Python代碼結(jié)構(gòu)展示了如何將上述思路串聯(lián)起來。這只是一個框架每個函數(shù)都需要根據(jù)具體數(shù)據(jù)和模型進行填充。import pandas as pd import numpy as np from prophet import Prophet import xgboost as xgb from pulp import LpProblem, LpVariable, lpSum, LpMaximize, LpStatus, PULP_CBC_CMD # 或者使用 ortools.linear_solver import pywraplp class VegetablePricingReplenishmentSystem: def __init__(self, historical_data, cost_data, config): 初始化系統(tǒng)。 historical_data: 歷史銷售、庫存、價格、特征數(shù)據(jù)。 cost_data: 商品成本、持有成本、損耗率等。 config: 配置參數(shù)如預(yù)測窗口、優(yōu)化窗口、倉庫容量等。 self.hist_data historical_data self.cost_data cost_data self.config config self.demand_models {} # 存儲每個商品的需求預(yù)測模型 self.current_inventory self._get_initial_inventory() def _build_demand_forecast_model(self, product_id): 為單個商品構(gòu)建混合需求預(yù)測模型如3.1節(jié)所述。 # 1. 提取該商品歷史數(shù)據(jù) product_data self.hist_data[self.hist_data[product_id] product_id].copy() # 2. 使用時序模型如Prophet擬合基準(zhǔn) # 3. 計算殘差 # 4. 使用XGBoost擬合殘差 # 5. 保存模型到 self.demand_models[product_id] pass def forecast_demand(self, product_id, future_features_df): 給定未來特征預(yù)測商品在未來N天的需求。 if product_id not in self.demand_models: self._build_demand_forecast_model(product_id) model self.demand_models[product_id] # 使用時序部分預(yù)測 temporal_pred model[temporal].predict(future_features_df[[ds]]) # 用機器學(xué)習(xí)部分預(yù)測殘差 residual_pred model[ml].predict(future_features_df[model[feature_cols]]) final_pred temporal_pred residual_pred return final_pred def mpc_optimization(self, current_date, horizon7): 執(zhí)行一次MPC滾動優(yōu)化返回當(dāng)天決策。 decisions {} all_products self.cost_data[product_id].unique() # 步驟1預(yù)測未來horizon天的需求 future_dates pd.date_range(startcurrent_date, periodshorizon, freqD) future_features self._generate_future_features(future_dates) # 構(gòu)建未來特征已知節(jié)假日、計劃天氣等 demand_forecast {} for pid in all_products: demand_forecast[pid] self.forecast_demand(pid, future_features) # 步驟2構(gòu)建并求解未來horizon天的優(yōu)化模型 prob LpProblem(Daily_Replenishment_Pricing, LpMaximize) # 定義決策變量字典 Q LpVariable.dicts(Q, (all_products, range(horizon)), lowBound0) # 補貨量 I LpVariable.dicts(I, (all_products, range(horizon1)), lowBound0) # 庫存多一期用于初始 S LpVariable.dicts(S, (all_products, range(horizon)), lowBound0) # 銷量 # 設(shè)置初始庫存 for pid in all_products: I[pid][0] self.current_inventory.get(pid, 0) # 目標(biāo)函數(shù)總利潤 profit_expr [] for t in range(horizon): for pid in all_products: cost self.cost_data.loc[pid, unit_cost] price ... # 這里價格可以是決策變量也可以是需求函數(shù)反推的結(jié)果簡化起見可以先固定或根據(jù)規(guī)則設(shè)定 holding_cost self.cost_data.loc[pid, holding_cost] profit_expr.append(price * S[pid][t] - cost * Q[pid][t] - holding_cost * I[pid][t]) prob lpSum(profit_expr) # 約束條件 for t in range(horizon): # 庫存平衡約束 for pid in all_products: spoilage_rate self.cost_data.loc[pid, spoilage_rate] prob I[pid][t1] (I[pid][t] Q[pid][t] - S[pid][t]) * (1 - spoilage_rate) # 銷量約束不能超過預(yù)測需求和可用庫存 for pid in all_products: prob S[pid][t] demand_forecast[pid][t] prob S[pid][t] I[pid][t] Q[pid][t] # 倉庫容量約束簡化示例只考慮總庫存體積 total_volume lpSum(I[pid][t1] * self.cost_data.loc[pid, volume_per_unit] for pid in all_products) prob total_volume self.config[warehouse_capacity] # 求解 solver PULP_CBC_CMD(msgFalse) prob.solve(solver) # 步驟3提取并返回第一天的決策 if LpStatus[prob.status] Optimal: for pid in all_products: decisions[pid] { replenish: Q[pid][0].varValue, price: ... # 從解中獲取或根據(jù)規(guī)則計算 } else: # 如果無解啟用備用啟發(fā)式規(guī)則 decisions self._fallback_heuristic(all_products, demand_forecast) return decisions def _fallback_heuristic(self, product_ids, demand_forecast): 備用啟發(fā)式規(guī)則例如采用報童模型計算補貨量價格采用成本加成法。 decisions {} for pid in product_ids: # 簡單報童模型補貨量 預(yù)測需求 安全庫存 forecast_d demand_forecast[pid][0] safety_stock ... # 根據(jù)歷史需求波動計算 decisions[pid] { replenish: max(0, forecast_d safety_stock - self.current_inventory.get(pid, 0)), price: self.cost_data.loc[pid, unit_cost] * 1.3 # 成本加成30% } return decisions def run_daily_operation(self, start_date, end_date): 模擬系統(tǒng)從start_date到end_date的每日運行。 current_date start_date decision_log [] while current_date end_date: # 1. 執(zhí)行MPC優(yōu)化得到今日決策 daily_decisions self.mpc_optimization(current_date) # 2. 模擬執(zhí)行決策在實際系統(tǒng)中這是真實操作 actual_sales self._simulate_sales(current_date, daily_decisions) # 3. 更新庫存狀態(tài)考慮實際銷量和損耗 self._update_inventory(current_date, daily_decisions, actual_sales) # 4. 記錄日志 decision_log.append({date: current_date, decisions: daily_decisions, actual_sales: actual_sales}) # 5. 進入下一天 current_date pd.Timedelta(days1) return pd.DataFrame(decision_log) # 主程序入口 if __name__ __main__: # 加載數(shù)據(jù) hist_df pd.read_csv(historical_sales.csv) cost_df pd.read_csv(product_costs.csv).set_index(product_id) config {warehouse_capacity: 10000, forecast_horizon: 7} # 初始化系統(tǒng) system VegetablePricingReplenishmentSystem(hist_df, cost_df, config) # 運行模擬 start_date pd.Timestamp(2023-06-01) end_date pd.Timestamp(2023-06-30) operation_log system.run_daily_operation(start_date, end_date) print(operation_log.tail())5. 模型檢驗、靈敏度分析與論文呈現(xiàn)要點建好模型和系統(tǒng)后如何驗證其有效性并在論文中清晰呈現(xiàn)是拿高分的關(guān)鍵。5.1 模型檢驗回溯測試與關(guān)鍵績效指標(biāo)KPI絕不能只用訓(xùn)練集上的表現(xiàn)來評價模型。必須進行嚴格的回溯測試。方法將歷史數(shù)據(jù)按時間劃分為訓(xùn)練集和測試集例如用前10個月訓(xùn)練后2個月測試。在測試集上模擬系統(tǒng)逐日運行每天開始時系統(tǒng)只知道過去的信息用它來做預(yù)測和優(yōu)化決策然后將決策與當(dāng)天的“真實”數(shù)據(jù)在測試集中進行比較。需要計算的KPI包括利潤相關(guān)總利潤、平均毛利率、投資回報率。庫存相關(guān)平均庫存水平、庫存周轉(zhuǎn)天數(shù)、缺貨率需求大于庫存的天數(shù)比例。損耗相關(guān)總損耗金額、損耗率損耗量/總進貨量。預(yù)測準(zhǔn)確性使用均方根誤差RMSE、平均絕對百分比誤差MAPE評估需求預(yù)測的精度。特別注意評估在庫存不足可能缺貨的那些日子的預(yù)測誤差這更能反映模型估算“真實需求”的能力。將你的模型與基準(zhǔn)策略對比例如簡單策略每天補貨到固定水平如過去7天平均銷量按成本固定比例加成定價。經(jīng)驗策略模擬有經(jīng)驗的店長的決策規(guī)則可從數(shù)據(jù)中歸納如周末多進20%。 只有顯著、穩(wěn)定地超越了基準(zhǔn)策略你的模型才算有價值。5.2 靈敏度分析找出模型的“脆弱點”和“穩(wěn)健性”模型依賴于許多假設(shè)和參數(shù)如損耗率θ、價格彈性e、持有成本h。靈敏度分析就是觀察當(dāng)這些參數(shù)在合理范圍內(nèi)變動時模型性能如總利潤的變化情況。如何做單因素分析固定其他參數(shù)讓一個參數(shù)如損耗率在[θ-Δθ, θΔθ]范圍內(nèi)變化運行多次模擬觀察利潤的變化曲線。如果利潤對該參數(shù)變化非常敏感說明模型在該參數(shù)上很脆弱需要更精確地估計該參數(shù)或者在策略中加入魯棒性設(shè)計。多因素分析可以使用拉丁超立方抽樣等方法同時擾動多個參數(shù)觀察利潤的分布。這能評估模型在參數(shù)不確定性下的整體穩(wěn)健性。場景分析模擬一些極端但可能發(fā)生的場景例如連續(xù)三天暴雨導(dǎo)致客流減少50%某主要供應(yīng)商突然斷供一周。你的系統(tǒng)決策會產(chǎn)生什么變化利潤會下降多少這能體現(xiàn)系統(tǒng)的抗風(fēng)險能力。在論文中用圖表清晰地展示靈敏度分析的結(jié)果例如畫出利潤隨關(guān)鍵參數(shù)變化的折線圖并給出文字解釋“如圖所示當(dāng)損耗率估計誤差在±20%內(nèi)時系統(tǒng)利潤波動在5%以內(nèi)表明模型對損耗率參數(shù)具有一定的魯棒性?！?.3 論文寫作與可視化講好一個邏輯閉環(huán)的故事國賽論文評審時間緊清晰的邏輯和直觀的可視化至關(guān)重要。行文邏輯建議問題重述與分析不要照抄題目要用自己的話精煉概括并明確指出問題的核心矛盾如利潤最大化與損耗最小化的權(quán)衡和子問題劃分。模型準(zhǔn)備詳細闡述數(shù)據(jù)預(yù)處理過程特別是“真實需求”估算、特征工程方法。這是體現(xiàn)你工作扎實程度的地方。模型建立分節(jié)闡述預(yù)測模型、定價模型、補貨模型。對每個模型都要說明為什么選擇這個模型與替代模型相比的優(yōu)勢以及模型如何與問題結(jié)合公式、變量含義要解釋清楚。模型求解重點介紹MPC框架和優(yōu)化模型的求解方法。如果是用求解器說明是什么求解器如果用了啟發(fā)式算法描述算法步驟。模型檢驗與靈敏度分析展示回溯測試的結(jié)果表格和對比圖用數(shù)據(jù)證明模型的有效性。展示靈敏度分析的圖表和結(jié)論。模型評價與推廣客觀評價模型的優(yōu)點和局限性例如假設(shè)需求函數(shù)形式固定未考慮競爭對手突然降價等。提出可能的改進方向并簡要說明模型思想可以推廣到其他易腐品零售如水果、鮮花、烘焙食品?？梢暬c系統(tǒng)框架圖用清晰的框圖展示“數(shù)據(jù)輸入 - 預(yù)測模塊 - 優(yōu)化模塊 - 決策輸出 - 滾動更新”的整個流程。預(yù)測效果圖對于幾個代表性商品畫出歷史期真實銷量、模型擬合銷量以及測試期預(yù)測銷量與真實銷量的對比折線圖。決策模擬圖選擇一段時間如測試期兩周畫出庫存水平、補貨量、價格、銷量的聯(lián)動變化圖直觀展示模型決策的動態(tài)過程。靈敏度分析圖用折線圖或熱力圖展示關(guān)鍵參數(shù)變化對核心KPI的影響。結(jié)果對比圖用柱狀圖對比你的模型與基準(zhǔn)策略在多個KPI上的表現(xiàn)。記住論文的本質(zhì)是向評委講述你發(fā)現(xiàn)問題、分析問題、解決問題、驗證方案的完整思考過程。代碼是工具模型是手段清晰的邏輯和令人信服的證據(jù)鏈才是獲得高分的關(guān)鍵。這道C題就像一個微縮的商業(yè)智能系統(tǒng)項目完整地走一遍無論比賽結(jié)果如何對解決實際問題的能力都是極大的提升。