简介面向校园信息化建设与数据分析应用场景基于Python的学生校园消费行为分析资源涵盖完整源码、脱敏数据与最终结果集。资源依托校园一卡通消费记录围绕食堂就餐峰值时段、各食堂及三餐占比、不同专业性别消费对比等维度展开分析并对高校食堂运行状况进行可视化呈现对智慧校园建设、贫困生精准援助背景下的消费行为研究具有实用参考价值。包内共22个文件包括7个Jupyter Notebook分析脚本、3个Python辅助脚本、1份docx详细分析报告、8张jpg可视化图表以及md说明文档和配套数据文件压缩包整体19.08MB目录结构按任务模块组织便于按需检索与复现。目前已有143人学习下载适合数据分析初学者、高校课程设计或校园一卡通项目开发者参考实践既能直接运行现有代码也方便替换数据以拓展自定义分析。1. 基于Python的学生校园消费行为分析小额高频数据里藏着的规律很多做学生工作的老师或者校园信息化负责人一开始想分析消费数据脑子里盘算的是“哪个学生家庭困难需要补助”“哪个学生近期状态不对劲”。但真正拿到一卡通流水之后才发现这份数据远比想象中乱同一天多次刷卡、早餐和夜宵混在一起、商户名称五花八门、金额从几毛到几十块都有。用Python做学生校园消费行为分析核心不是跑通一段代码而是把“消费流水”变成“可解释的行为结论”谁在规律吃饭、谁的生活节奏紊乱、谁的消费能力异常波动。这个过程涉及数据清洗、特征工程、聚类分群和异常检测最后落到一份可汇报的标签结果集。适合辅导员、后勤管理人员也想给数据方向的学生练手——一套完整可复现的分析链路比单点算法更值钱。2. 消费数据从哪来校园一卡通流水的字段设计与样本构造2.1 一张能用的消费流水必须包含哪些字段做消费行为分析最怕拿到一张“看上去什么都有、实际什么都用不上”的表。常见翻车现场是商户名称写成了“档口1”“档口2”没有分类消费时间只有日期没有时刻卡号做了脱敏但把学号信息全抹掉了导致无法关联年级和专业。我一般会要求数据至少包含五个字段用户标识脱敏后的学号或卡号、交易时间精确到秒的时间戳、商户名称或编号、消费金额单位统一为元、支付方式刷脸/刷卡/扫码用于判断消费是否本人操作。这五个字段缺一不可。用户标识负责把流水切分成“一个人”的视角交易时间负责还原行为节律商户名称负责还原偏好金额负责度量消费水平支付方式则是排查重复扣款和盗刷的线索。如果只有金额和时间能做的分析非常有限——分群只能用ARPU值硬切食堂的早餐偏好、夜宵习惯这些深度画像全部做不了。有些厂家导出的数据还有补贴字段、剩余余额字段这些属于锦上添花。补贴字段能判断贫困生资格是否生效剩余余额能算出“月初富月末穷”的月光指数。但这类字段经常缺失或不准我建议先按核心五字段做后续有数据再补齐。2.2 没有现成数据用Python生成一份拟真校园消费流水很多时候你手上没有真实数据尤其是想跑通整个分析流程练手的情况。常见做法是用Python自己造一份拟真数据关键点是“拟真”——消费行为要有规律可循否则后面聚类、异常检测都跑不出有意义的结论。import pandas as pd import numpy as np from datetime import datetime, timedelta np.random.seed(42) # 基础参数模拟 2000 名学生30 天的消费记录 student_ids [fS{1000i} for i in range(2000)] days pd.date_range(2024-03-01, periods30, freqD) # 定义商户池食堂位置 商品类型决定消费金额范围 venues { 第一食堂·早餐: (2.0, 6.5), 第一食堂·正餐: (8.0, 18.0), 第二食堂·正餐: (9.0, 22.0), 清真食堂·正餐: (8.0, 20.0), 面包房: (5.0, 15.0), 超市: (3.0, 30.0), 咖啡吧: (12.0, 28.0), 水果店: (4.0, 35.0), } rows [] for student in student_ids: # 每人每天 2~5 笔消费 for day in days: n_trans np.random.randint(2, 6) for _ in range(n_trans): venue np.random.choice(list(venues.keys())) low, high venues[venue] amount round(np.random.uniform(low, high), 2) # 时间集中在 7:00-22:00早中晚各有高峰 hour np.random.choice([7, 8, 8, 11, 12, 12, 13, 17, 18, 18, 19], p[0.1, 0.1, 0.1, 0.1, 0.15, 0.1, 0.05, 0.1, 0.15, 0.03, 0.02]) minute np.random.randint(0, 60) ts datetime(day.year, day.month, day.day, hour, minute, np.random.randint(0, 60)) rows.append([student, ts, venue, amount, np.random.choice([face, card, card])]) df pd.DataFrame(rows, columns[student_id, ts, venue, amount, pay_method]) df df.sort_values([student_id, ts]).reset_index(dropTrue) df.to_csv(campus_consumption.csv, indexFalse) print(df.shape) print(df.head())这段代码生成的流水有student_id、ts、venue、amount、pay_method五列共约 20 万行。生成逻辑上做了两个关键设计金额区间按商户类型收敛早餐不可能出现 30 元的消费时间点按真实就餐高峰分布11 点到 12 点的概率权重明显高于下午三点。这两点是拟真数据的灵魂——如果金额和时间都是均匀分布后面画出来的图就是一片噪声特征工程也学不到任何东西。参数上需要注意np.random.seed(42)固定随机种子保证每次生成结果一致n_trans控制在 2 到 5 笔符合学生日常消费频率p概率向量要和np.random.choice的取值列表严格一一对应否则会报错或者分布错乱。2.3 拿到流水后的第一步清洗与脱敏的三项检查数据清洗是整个项目里最不性感但最决定成败的环节。我的血泪经验是不要一上来就写特征工程先花半小时确认三件事——时间字段能不能解析、有没有重复记录、金额有没有异常值。import pandas as pd df pd.read_csv(campus_consumption.csv, parse_dates[ts]) # 检查缺失值 print(缺失值统计) print(df.isnull().sum()) # 检查重复记录同一人同一秒同一商户刷了两次大概率是重复流水 dup df.duplicated(subset[student_id, ts, venue, amount], keepFalse) print(f重复记录数{dup.sum()}) # 金额边界检查0 元、负数、超过 100 元 amount_error df[(df[amount] 0) | (df[amount] 100)] print(f金额异常记录数{len(amount_error)}) # 时间范围检查是否落在学期内 print(f时间范围{df[ts].min()} 到 {df[ts].max()})这里的逻辑很简单parse_dates把时间字符串解析成 datetime 类型后续所有按小时、按星期几的聚合都依赖这一步duplicated用全字段匹配查重复真实环境里经常因为网络重传或者终端重复上送产生一模一样的记录金额检查重点关注 0 元和负数0 元可能是补贴发放记录负数是退费这两种在特征计算时要单独处理不能直接混进消费金额求均值。清洗完成后要对student_id做脱敏处理。真实项目里一般不允许直接使用学号常见做法是保留学号的年级前缀比如前两位代表入学年份后几位做哈希映射。这样既保留了分群所需的年级信息又不会暴露个人身份。3. 从流水到特征把消费行为变成四个维度的数值指标3.1 消费频次与金额特征RFM 模型在校园场景的变形RFM 是电商常用的用户价值分析模型三个维度分别是最近一次消费时间Recency、消费频率Frequency、消费金额Monetary。校园消费场景不能直接套用最大的差异是时间窗口——学生天天在校内吃饭R 的区分度极低几乎人人都是“最近消费”状态。我一般把 R 改成“晚餐规律性”或“断卡间隔天数”把 F 改成“日均消费笔数”M 改成“日均消费金额”这样更贴合校园实际。import pandas as pd df[date] df[ts].dt.date df[hour] df[ts].dt.hour # 用户维度聚合 user_features df.groupby(student_id).agg( total_transactions(ts, count), total_amount(amount, sum), avg_amount(amount, mean), active_days(date, nunique), unique_venues(venue, nunique), ).reset_index() # 日均指标 user_features[trans_per_day] user_features[total_transactions] / user_features[active_days] user_features[amount_per_day] user_features[total_amount] / user_features[active_days] # 晚餐消费习惯18:00-19:59 之间的消费次数占比 df[is_dinner] df[hour].between(18, 19) dinner_stats df.groupby(student_id)[is_dinner].mean().rename(dinner_ratio) user_features user_features.merge(dinner_stats, onstudent_id) print(user_features.head())active_days比总天数更能反映真实活跃度——如果一个学生 30 天只在校 15 天用总天数做分母会严重低估他的消费强度。trans_per_day和amount_per_day能识别两种典型人群一种是每天三顿规律食堂、日均 3 笔左右另一种是经常校外解决或点外卖校内消费每天只有一两笔。dinner_ratio是校园场景特有的强特征晚餐消费占比高的学生生活规律性明显更好。RFM 的三个指标不是都要用。我做过对比R最近消费时间在校园数据里贡献度极低去掉之后聚类效果反而更干净。保留 F 和 M 的变形再加上后面讲的时间规律和商户偏好特征矩阵的区分度就出来了。3.2 时间规律特征消费节律比消费金额更有解释力消费金额只能反映“花了多少”消费时间分布能反映“怎么花的”。我习惯把一天切成六个时段早餐 7:00-9:00、上午 9:00-11:00、午餐 11:00-13:00、下午 13:00-17:00、晚餐 17:00-19:00、夜宵 19:00-23:00。然后计算每个学生在每个时段的消费次数占比这组特征能直接刻画生活节律。import pandas as pd def time_period(hour): if 7 hour 9: return breakfast elif 9 hour 11: return morning elif 11 hour 13: return lunch elif 13 hour 17: return afternoon elif 17 hour 19: return dinner else: return night df[period] df[hour].apply(time_period) period_pivot df.pivot_table( indexstudent_id, columnsperiod, valuests, aggfunccount, fill_value0 ) # 转为占比特征每个时段的消费次数 / 总消费次数 period_ratio period_pivot.div(period_pivot.sum(axis1), axis0) period_ratio.columns [fperiod_{c}_ratio for c in period_ratio.columns] user_features user_features.merge(period_ratio, onstudent_id) # 周末效应计算周一到周五与周六周日的日均消费差 df[weekend] pd.to_datetime(df[date]).dt.weekday 5 weekday_amount df[~df[weekend]].groupby(student_id)[amount].mean().rename(weekday_avg) weekend_amount df[df[weekend]].groupby(student_id)[amount].mean().rename(weekend_avg) weekend_effect pd.concat([weekday_amount, weekend_amount], axis1).fillna(0) weekend_effect[weekend_gap] weekend_effect[weekend_avg] - weekend_effect[weekday_avg] user_features user_features.merge(weekend_effect[[weekend_gap]], onstudent_id)时段占比特征的好处是天然消除了消费总量的影响。一个每天消费 5 笔的学生和一个每天消费 2 笔的学生在占比特征上可以直接对比。weekend_gap能识别特殊人群——有些学生周末大量在校外消费校内金额骤降这个值会很明显反之考研党周末泡图书馆食堂消费反而上涨。时间特征容易出现一个坑直接用hour做均值算子。早上 7 点和晚上 19 点的小时数均值是 13语义完全错误。必须先映射成时段标签再做统计或者用环形编码处理。这一点在特征工程里特别容易翻车。3.3 商户偏好与消费结构特征食堂依赖度怎么算校园消费分析和普通电商分析一个显著区别是商户类型天然具有“必需消费”和“改善消费”的层级。食堂是刚需超市是补充咖啡吧是改善。计算每个人在“食堂类商户”的消费占比能得到一个食堂依赖度指标——这个指标对贫困生认定和校园商业规划都非常有用。import pandas as pd # 商户分类根据名称关键词打标签 def categorize_venue(venue): if 食堂 in venue: return canteen elif 超市 in venue: return market elif 面包 in venue: return bakery elif 咖啡 in venue: return cafe elif 水果 in venue: return fruit else: return other df[venue_cat] df[venue].apply(categorize_venue) # 每个学生的商户消费结构 cat_pivot df.pivot_table( indexstudent_id, columnsvenue_cat, valuesamount, aggfuncsum, fill_value0 ) for c in cat_pivot.columns: cat_pivot[f{c}_ratio] cat_pivot[c] / cat_pivot.sum(axis1) user_features[f{c}_ratio] cat_pivot[f{c}_ratio] # 早餐食堂依赖如果早餐时段消费主要在食堂说明生活习惯稳定 breakfast df[df[period] breakfast] breakfast_canteen breakfast[breakfast[venue_cat] canteen].groupby(student_id)[amount].count() breakfast_total breakfast.groupby(student_id)[amount].count() breakfast_canteen_ratio (breakfast_canteen / breakfast_total).rename(breakfast_canteen_ratio) user_features user_features.merge(breakfast_canteen_ratio, onstudent_id, howleft).fillna(0)商户分类必须用关键字规则而非人工打标。真实数据里的商户名通常是“一食堂一楼 12 号窗口”“超市东区店”这样带编码的如果不做关键词匹配几十个档口根本没法聚合。ratio类特征用的是“金额占比”而非“次数占比”因为不同商户的客单价差异极大咖啡吧一笔顶食堂三笔用金额占比更能反映真实消费结构。一个容易被忽视的细节早餐食堂依赖度要用次数占比因为早餐金额本身就低用金额算会被极高的咖啡消费稀释掉。不同特征用不同的聚合口径这是特征工程里最考验经验的地方。3.4 特征宽表落地合并、去重、标准化所有特征都算完之后要合并成一张宽表——每行一个学生每列一个特征。然后做标准化和缺失值处理准备喂给聚类模型。from sklearn.preprocessing import StandardScaler # 剔除纯标识列和可能泄漏的列选择用于建模的特征 feature_cols [ trans_per_day, amount_per_day, avg_amount, dinner_ratio, period_breakfast_ratio, period_lunch_ratio, period_dinner_ratio, period_night_ratio, weekend_gap, canteen_ratio, market_ratio, cafe_ratio, fruit_ratio, breakfast_canteen_ratio ] model_data user_features[[student_id] feature_cols].copy() model_data model_data.dropna() # 标准化消除量纲影响均值 0方差 1 scaler StandardScaler() scaled scaler.fit_transform(model_data[feature_cols]) model_data_scaled pd.DataFrame(scaled, columnsfeature_cols, indexmodel_data.index) # 保存宽表供后续分析使用 pd.concat([model_data[[student_id]], model_data_scaled], axis1).to_csv(user_features.csv, indexFalse) print(f特征宽表形状{model_data_scaled.shape})feature_cols的选取有讲究total_amount和amount_per_day高度相关不能同时进模型period_morning_ratio和period_afternoon_ratio区分度低也建议去掉。标准化是把所有特征压到同一个尺度否则金额类特征动辄几十占比类特征只有零点几聚类时金额会主导距离计算。dropna删掉那些没有早餐记录的极少数学生他们对早餐特征的贡献是缺失值不能简单填 0——填 0 相当于告诉模型“这个人从不吃早饭”与实际情况不符。特征宽表的索引用student_id但标准化后的矩阵是 numpy 数组要在合并时恢复索引否则后面聚类结果回贴到宽表时会对不齐。这种低级错误我犯过一次排错排了半小时。4. 聚类分群与异常检测让模型告诉你校园里有几类人4.1 用 KMeans 做学生消费分群轮廓系数决定 K 值特征准备好之后最先做的应该是无监督聚类——把学生分成几个行为模式相近的群体。校园消费行为的聚类 K 值一般是 3 到 5太少区分不出层次太多就碎成了散户。不要拍脑袋定 K用轮廓系数扫一遍。from sklearn.cluster import KMeans from sklearn.metrics import silhouette_score import numpy as np data pd.read_csv(user_features.csv, index_col0) # 用标准化后的特征矩阵注意读进来之后要再标准化一次或直接用模型输入 X data.values silhouette_scores [] K_range range(2, 7) for k in K_range: km KMeans(n_clustersk, random_state42, n_init10) labels km.fit_predict(X) score silhouette_score(X, labels) silhouette_scores.append(score) print(fK{k}, silhouette_score{score:.4f}) best_k K_range[np.argmax(silhouette_scores)] print(f最佳 K 值{best_k})n_init10是 KMeans 的一个关键参数KMeans 的初始质心是随机选的不同初始值可能收敛到不同局部最优解。n_init10表示从 10 次随机初始化中选最优结果避免 “跑两次结果不一样” 的尴尬。轮廓系数的取值范围是 -1 到 1大于 0.25 就说明聚类结构可以接受校园消费数据一般能到 0.3 到 0.4 之间超过 0.5 反而要警惕是不是特征构造出了问题——真实的消费行为不该有这么清晰的边界。聚完类之后每个学生打上一个群体标签。这个标签就是后面结果集的核心内容之一。最佳 K 值对应的聚类模型要重新训练一次并保存方便后续对新增数据预测。4.2 异常消费检测Z-score 与孤立森林的判断边界聚类的输出是“群体画像”异常检测的输出是“偏离者”。两类方法各有边界Z-score 适合检测单一维度的极端值比如日均消费 80 元的人孤立森林适合检测多维度的“组合异常”比如日均消费只有 5 元但咖啡吧消费占比 60%——单看哪个维度都不离谱组合起来就是异常。from sklearn.ensemble import IsolationForest import numpy as np # 方法一用 Z-score 检测日均消费金额异常 data[amount_per_day_z] (data[amount_per_day] - data[amount_per_day].mean()) / data[amount_per_day].std() high_spenders data[data[amount_per_day_z] 2.5] low_spenders data[data[amount_per_day_z] -1.5] print(f高消费异常Z2.5{len(high_spenders)} 人) print(f低消费异常Z-1.5{len(low_spenders)} 人) # 方法二孤立森林contamination 设为 5% iso_forest IsolationForest( n_estimators200, contamination0.05, random_state42 ) data[anomaly_score] iso_forest.fit_predict(data[feature_cols]) # predict 返回 1 表示正常-1 表示异常 data[is_anomaly] data[anomaly_score] -1 anomaly_ids data[data[is_anomaly]].index.tolist() print(f孤立森林识别异常学生数{len(anomaly_ids)})Z-score 的阈值选取高消费用 2.5 而不是 2因为校园消费的正态性比电商数据好2.5 能筛出真正意义上的极端值低消费用 -1.5 是因为低于均值 1.5 个标准差的人已经明显偏离正常生活水平再低下去样本太少。孤立森林的contamination是关键参数它告诉模型“异常比例约 5%”设高了会把正常波动当异常设低了会漏掉真正的风险对象。注意 Z-score 和孤立森林找出来的异常人群大概率不重合。Z-score 偏重“消费规模异常”孤立森林偏重“消费结构异常”。实际项目中两类异常人列表取交集或者分别处理不要混为一谈。4.3 结果集怎么组织把标签、分数和原始画像拼回一张总表分析结果不能存在模型的临时变量里要落成一份可持续使用的文件。我习惯输出两张表一张是“学生标签总表”包含分群标签、异常标记、关键特征数值另一张是“群体画像汇总表”描述每类人的特征均值和业务解释。import pandas as pd # 用最佳 K 值重新拟合得到最终聚类标签 km_final KMeans(n_clustersbest_k, random_state42, n_init10) data[cluster_label] km_final.fit_predict(X) # 结果集学生维度 result pd.DataFrame({ student_id: data.index, cluster: data[cluster_label], is_anomaly: data[is_anomaly], trans_per_day: data[trans_per_day], amount_per_day: data[amount_per_day], canteen_ratio: data[canteen_ratio], night_ratio: data[period_night_ratio].fillna(0), }) result.to_csv(student_analysis_result.csv, indexFalse, encodingutf-8-sig) # 群体画像汇总每个聚类中心点的特征均值 profile data.groupby(cluster_label)[feature_cols].mean().round(3) profile.to_csv(cluster_profile.csv, indexTrue) print(result.head()) print(profile)导出 CSV 时用utf-8-sig编码是个细节——如果用默认的 utf-8Excel 打开带中文的 CSV 会乱码加了 BOM 头就不会。这个坑几乎每个接手项目的人都会踩一次。profile表是聚类模型的解释层业务人员不看聚类中心坐标看的是“第一类人日均消费 35 元、食堂占比 80%、几乎没有夜宵记录”——这就是他们能直接拿去用的结论。结果集的设计要兼顾机器可读和人类可读。student_analysis_result.csv供后续系统读取或二次分析cluster_profile.csv供人写报告做汇报。两个文件一套流程这才叫完整的“源码数据结果集”。5. 避坑校园消费分析的五个高频坑与排查方法5.1 特征泄漏把全周期的汇总值当成单日特征现象聚类结果异常清晰轮廓系数超过 0.6K2 时两个群体几乎完全分离。检查发现“日均消费金额”是用整个月的总金额除以天数算的而“总消费次数”也是全月汇总两个特征本质上是同一个信息。原因特征之间高度共线模型看到的是同一个因素的三个分身不是三个独立维度。解决计算特征相关性矩阵data.corr()之后删掉相关系数超过 0.8 的冗余特征。我在这个项目里删掉了total_amount和total_transactions只保留人均版本模型解释力反而更强。5.2 时间字段解析失败日期格式隐藏的时区陷阱现象pd.to_datetime报错提示Unknown string format仔细看原始数据发现时间是2024/03/01 12:30而不是2024-03-01 12:30:00。原因校园一卡通厂商的系统导出的时间格式五花八门有的用斜杠有的带中文“上午/下午”有的秒字段缺失。解决统一用format参数强制指定格式pd.to_datetime(df[ts], format%Y/%m/%d %H:%M)。不要依赖parse_dates的自动推断自动推断在数据量大的时候速度慢且容易猜错。5.3 金额为负数退费记录被当成消费现象计算日均消费金额时出现负值个别学生日均消费 -3 元聚类结果出现一个“负数群体”。原因一卡通系统里退费记录和消费记录在同一张表退费金额以负数形式存储。解决分离正负数记录负数标记为refund单独统计退费金额不参与消费均值计算。如果退费次数异常多比如一个学生 30 天退费 10 次本身就是一个风险信号值得单独拉出来排查。5.4 KMeans 对初始值敏感两次运行聚类结果不同现象同样的数据和参数跑两次聚类每次得到的群体人数分布都不一致但买过的人并不多。原因random_state没有设置KMeans 初始质心随机收敛到的局部最优解不同。解决所有涉及随机性的模型都要固定random_state42n_init设置大于 5。这不是玄学是 KMeans 这类基于迭代的算法在数据集量不大时初始值对结果影响的真实存在。固定随机种子是保证分析结果可复现的最低要求。5.5 结果无法解释聚类标签没有业务含义现象聚类跑通了输出标签 0、1、2、3但业务方问“这三类人分别有什么特征”回答不上来。原因聚类只做了分群没有做群体画像的解读。解决聚类之后必须立刻计算每个类别在核心特征上的均值和标准差写成人话。0 类是“日均 3 笔、食堂占比 80%、无夜宵”1 类是“咖啡高频、周末校外消费明显”2 类是“低消费、低活跃、疑似数据不全”。标签本身没有意义标签对应的行为描述才有意义。6. 验证分析质量用轮廓系数、画像雷达图和样本回查确认结果可信聚类做完不是终点结果要经得起三个方向的验证否则交付给业务方就是自找麻烦。第一是结构验证再看一眼轮廓系数是否稳定在 0.25 以上如果低于 0.2说明特征工程没有抓住用户差异要回炉重构特征而不是调聚类参数。第二是业务验证把每类人的特征均值做成雷达图比如“消费稳定型”“餐饮依赖型”“改善消费型”“高风险型”四类画像看各类图形的差异是否直观可见。第三是样本回查随机挑 20 个学生回到原始流水里人工核对——聚类说他是“规律三餐型”翻记录应该能看到工作日每天都有早餐和晚餐记录而不是一天三顿全集中在食堂但时间完全随机。雷达图绘制也有一点值得说的技巧画之前要把特征压缩到相同量纲比如都映射到 0 到 1 的范围否则金额均值 25 和食堂占比 0.8 画在同一个坐标轴上毫无可比性。另外每个特征的“好”方向要统一——消费频次是越高越好还是越低越好取决于业务定义雷达图里默认全是“越高越靠外”如果某个特征的业务含义是“越低越好”比如周末消费波动需要先取倒数或反转方向再画。我自己的习惯是每一轮分析都保留一个checkpoint文件夹里面存原始数据摘要、特征宽表、聚类参数和结果集。这个习惯救过我一次项目上线三个月后业务方回查一批学生的异常标记发现判定依据不一致我直接用当时的 checkpoint 复现了整个分析过程定位到是厂商补发了新旧两版数据导致student_id前缀规则变化重新对齐后结论依然成立。做数据分析过程可追溯比结论漂亮更重要。愿你在这个项目上的分析结果也能经得起这样的回查。希望帮到你。本文还有配套的精品资源点击获取