简介面向计算机相关专业毕业设计、课程设计及期末大作业场景是一套基于机器学习的船舶碳排放数据驱动优化分析系统完整工程。项目代码经导师指导并认可评审分高达99分运行环境配置简单即使基础薄弱的学习者也可直接运行并复现结果。压缩包共33个文件体积仅780KB涵盖7个Python脚本分别实现K-Means聚类、随机森林回归、层次聚类及数据清洗等任务17张结果可视化PNG图展示特征重要性、实际值与预测值对比、残差分布、轮廓系数与肘部曲线等分析结果另有聚类结果CSV数据、5个工程配置文件和1个Markdown说明文档方便理解项目结构和实现思路。当前已有47人学习下载。项目覆盖从数据清洗、特征探索到模型训练、评估与可视化的完整流程可帮助读者掌握船舶燃料消耗与碳排放相关因素的分析方法理解回归与聚类算法在实际工业数据中的应用价值也便于在现有基础之上做二次开发以支撑毕设答辩或项目汇报。1. 船舶碳排放的数据驱动优化这份 Python 项目到底能帮你什么拿到船舶碳排放数据第一反应通常是“换燃料、改发动机”但这套基于机器学习的系统走的是另一条路不动任何物理设备只用历史航行数据把“什么工况最费油、哪个因素对排放影响最大、船队能分成几类典型运行模式”挖出来再反过来指导航速、航线、配载的优化。这套 Python 项目就是一个完整的数据驱动优化分析系统从数据清洗开始到聚类分析、回归预测、模型评估再到可视化归因代码和图表都是闭环的。源码包里有 20 多个脚本和输出图包含 KMeans、层次聚类、随机森林回归、线性回归四条分析线还有特征重要性、偏依赖图、轮廓系数、肘部曲线等一整套诊断结果。适合正在做毕设、课程设计或期末大作业的计算机相关专业同学也适合想练手“机器学习落地全流程”的初级开发者——代码完整可运行侧重于把分析链路跑通而不是停留在调库层面。2. 先过数据清洗这道关三个 cleaning 脚本各管什么这个项目里的数据清洗脚本拆成了三份对应三条不同的分析线。你可能会想一个数据集为什么清洗要拆三个脚本直接统一处理不行吗实际上这是因为聚类和回归对数据格式的要求不一样聚类要把数值特征全部标准化回归则需要保留原始量纲来做预测误差分析。拆开写每个脚本的职责更清晰也方便单独调试。2.1 三个清洗脚本的分工与调用方式脚本文件服务对象核心处理目标k_means_data_cleaning.pyk_means.py / pca_clusters.py缺失值填充、标准化、留出设计特征用于聚类hierarchical_clustering_data_cleaning.pyHierarchical Clustering.py同样的缺失值处理但按层次聚类的输入格式输出regression_model_data_cleaning.pyregression_model.py / random_forrst.py保留原始量纲增加燃料消耗和 CO₂ 的关联字段处理跑清洗脚本没有固定顺序三份相互独立。我一般先跑regression_model_data_cleaning.py因为回归模型的输入字段最多清洗逻辑最能暴露数据质量问题。命令行执行即可python regression_model_data_cleaning.py这个脚本内部做的事情是读入原始航行数据删除全空列数值列用中位数填充缺失值对异常值做 IQR四分位距截断最后把清洗结果导出成一个干净的 CSV 供回归模型读取。如果数据文件不在默认路径在脚本开头改pd.read_csv()的参数就行。我做毕设时习惯把原始数据放data/子目录清洗后输出到clean/好处是后面反复调参不会污染原始文件。2.2 数据清洗的代码逻辑与参数说明下面这段清洗逻辑是从项目清洗脚本里提炼出来的核心思路你在三份清洗脚本里都能看到这种结构import pandas as pd import numpy as np df pd.read_csv(ship_emission_data.csv) # 1. 去掉全空的列 df.dropna(axis1, howall, inplaceTrue) # 2. 数值列缺失值用中位数填充中位数比均值对离群点更鲁棒 num_cols df.select_dtypes(include[np.number]).columns for col in num_cols: median_val df[col].median() df[col] df[col].fillna(median_val) # 3. 用 IQR 方法标记异常值并用上下边界截断 for col in num_cols: q1 df[col].quantile(0.25) q3 df[col].quantile(0.75) iqr q3 - q1 lower, upper q1 - 1.5 * iqr, q3 1.5 * iqr df[col] df[col].clip(lower, upper) df.to_csv(ship_emission_clean.csv, indexFalse) print(f清洗完成保留 {df.shape[0]} 行{df.shape[1]} 列)逻辑说明第一步直接删掉完全没有数据的列这一步看似简单但在真实数据里很关键因为船舶航行日志经常整列留空第二步用中位数而不是均值填充缺失值是因为航速、油耗这类字段里有明显的长尾异常值均值会被拉偏第三步的clip把超出 1.5 倍 IQR 的值截断到边界值这个操作保留了异常样本但削弱了它对模型的杠杆效应。注意quantile(0.25)和quantile(0.75)是 IQR 的核心参数如果你发现某个特征截断后数据分布变得过于集中可以把 1.5 改成 2.0 放宽边界。数据清洗不干净后面聚类出来的群就是玄学回归模型的残差图也一定会翻车。2.3 清洗后你应该核对什么清洗脚本跑完后输出的提示里通常只有行列数这个信息远远不够。我会再追加一组核对命令python -c import pandas as pd df pd.read_csv(ship_emission_clean.csv) print(df.isnull().sum().sum(), 个缺失值) print(df.describe().T[[mean, std, min, max]]) 看两个地方缺失值总数应该为 0这是清洗完成的基本标志min和max的值应该在合理范围比如航速不可能是负数燃料消耗的极值不应该出现截断后仍然离群的情况。如果数据显示异常优先检查原始数据的单位是否统一——这个项目里的燃料消耗字段是FuelConsumption_kg单位是公斤有的船级社报告里会用吨差三个数量级聚类时单位不统一做出来的结果没法看。3. 聚类分析找工况模式KMeans 与层次聚类的选 K 逻辑聚类在这个项目里的定位不是“自动分成几类”就完了而是要回答一个问题这支船队或这一艘船的多个航次存在几种典型运行工况不同工况下减排策略完全不同——低速巡航工况的优化空间在航速高负荷工况的优化空间在配载和航线。把这个问题拆给算法就是 KMeans 和层次聚类分别做的事。3.1 KMeans 的选 K 逻辑与核心参数项目里的k_means.py不是只跑一次聚类就出结果它同时输出了elbow_plot.png和silhouette_score.png两张判断图说明选 K 的过程是硬编码在脚本里的。我看过不少毕设项目直接KMeans(n_clusters3)一把梭然后写“根据业务经验选择 3 类”——这在这个项目里是行不通的因为输出图暴露了选 K 的全部思考过程答辩时老师一定会追问。KMeans 的标准流程是先标准化再用肘部法则和轮廓系数双重验证代码结构大致如下from sklearn.preprocessing import StandardScaler from sklearn.cluster import KMeans import matplotlib.pyplot as plt # 聚类输入特征列根据实际数据字段调整 features [FuelConsumption_kg, Annual_Time_spent_at_sea_hours, Ship_speed, Cargo_load, Engine_power] X df[features].copy() # 标准化KMeans 基于欧氏距离量纲不统一会让数值大的特征主导聚类结果 scaler StandardScaler() X_scaled scaler.fit_transform(X) # 肘部法则计算 2 到 10 类的 SSE sse [] for k in range(2, 11): km KMeans(n_clustersk, initk-means, random_state42, n_init10) km.fit(X_scaled) sse.append(km.inertia_) plt.plot(range(2, 11), sse, markero) plt.xlabel(Number of clusters K) plt.ylabel(SSE) plt.savefig(elbow_plot.png, dpi150, bbox_inchestight)参数说明initk-means是必须保留的默认的随机初始化会让结果随运行次数飘忽不定random_state42保证可复现这个项目里所有涉及随机性的脚本大概率都设了固定种子你自己重构时一定要记得加n_init10是让算法在 10 组不同的初始中心里选最优解SSE 曲线会更平滑。肘部法则的原理是K 小于真实簇数时每增加一类 SSE 下降幅度很大超过真实簇数后下降幅度趋缓曲线出现“肘部”拐点。但实际数据里拐点经常不明显所以项目里配了第二张轮廓系数图来交叉验证。3.2 轮廓系数与层次聚类的配合验证轮廓系数是每个样本的簇内凝聚度和簇间分离度的综合指标取值范围 [-1, 1]越接近 1 说明聚类越合理。项目里silhouette_scores.png绘制了 K 在 2 到 10 之间的得分曲线直接选得分最高的 K 即可但要注意轮廓系数在 K 太大比如 8 类以上时可能因为过分割而偏高所以最终 K 要结合业务解释性——3 到 5 类能解释清楚工况就别选 8 类。层次聚类的脚本Hierarchical Clustering.py走的是另一条路它不用预先指定 K而是通过树状图决定从哪里切一刀。核心代码示意from scipy.cluster.hierarchy import dendrogram, linkage, fcluster import matplotlib.pyplot as plt # linkage 默认用 ward 法最小化类内方差 Z linkage(X_scaled, methodward) # 画树状图从图上找距离最大的分层位置切割 plt.figure(figsize(10, 6)) dendrogram(Z, truncate_modelevel, p5) plt.show() # 实测时按树状图切在距离阈值约 8-10 的位置生成类别标签 labels fcluster(Z, t9, criteriondistance)methodward是我最常用的它和 KMeans 的目标函数很像都是最小化类内方差所以两种算法在相同数据上得到的聚类结构理论上应该接近——这也是项目的设计意图用两种算法互相验证如果 KMeans 分 4 类而层次聚类切出来也是 4 类结论就立住了。truncate_modelevel, p5是只显示树状图上面 5 层避免叶子太多图糊成一片。距离阈值t需要你根据树状图的实际高度调这个没有固定答案我的经验是先找最大距离的一半位置试切再看各类的样本量是否均衡。3.3 聚类结果输出与 PCA 降维验证聚类跑完后项目里把标签和原始数据合并导出成clustered_data.csv同时生成了pca_clusters.png。PCA 降维的作用是把聚类结果从高维空间投影到二维散点图人眼直接看分组是否清晰。这个文件里有几个关键点不同类别用不同颜色区分如果投影后各类别之间重叠严重说明聚类质量差或者在标准化阶段出了问题。你可以用下面的代码快速验证聚类结果在原始维度上的差异import pandas as pd clustered pd.read_csv(clustered_data.csv) cluster_profile clustered.groupby(Cluster).mean() print(cluster_profile.T)这段groupby输出的是每个聚类中心在所有特征上的均值用来做业务解释。比如 Cluster 0 的Ship_speed和FuelConsumption_kg都高说明这是“高速重载”工况Cluster 1 的Annual_Time_spent_at_sea_hours高但油耗低说明是“低速长航”工况。每一类对应的减排建议完全不同第六章里我会展开说怎么组织这个逻辑。4. 回归预测定排放基线随机森林和线性回归双模型对比与模型诊断聚类回答的是“船队有哪些典型状态”回归回答的是“给定这些状态参数燃料消耗进而碳排放是多少”。项目里两个回归脚本分别对应两种建模思路随机森林处理非线性关系线性回归作为基线对照。这种双模型对比是毕设拿高分的关键设计——只给一个模型的结果老师会质疑你为什么不试别的两个模型放在一起就变成了“基线方法 改进方法”的教科书式结构。4.1 随机森林回归特征重要性与调参逻辑random_forrst.py这个文件名有点拼写错误随机森林应该是 random_forest但内容不影响使用。随机森林在这个场景里的优势是不需要假设特征和油耗之间是线性关系航速、载重、航行时间对排放的影响本来就不是直线——低速段的油耗增长平缓高速段油耗指数上升线性模型很难刻画这种分段效应。训练部分的代码结构我用下面这段概括from sklearn.ensemble import RandomForestRegressor from sklearn.model_selection import train_test_split from sklearn.metrics import r2_score, mean_absolute_error, mean_squared_error import numpy as np features [Ship_speed, Cargo_load, Engine_power, Annual_Time_spent_at_sea_hours, Distance_travelled_nm] X df[features] y df[FuelConsumption_kg] X_train, X_test, y_train, y_test train_test_split( X, y, test_size0.2, random_state42 ) model RandomForestRegressor( n_estimators300, max_depth10, min_samples_leaf3, random_state42 ) model.fit(X_train, y_train) y_pred model.predict(X_test) print(fR²: {r2_score(y_test, y_pred):.3f}) print(fMAE: {mean_absolute_error(y_test, y_pred):.3f} kg) print(fRMSE: {np.sqrt(mean_squared_error(y_test, y_pred)):.3f} kg)代码逻辑先按 8:2 划分训练集和测试集random_state42保证每次划分一致n_estimators300是决策树的数量300 棵以内模型精度随数量上升明显超过 300 收益递减而且训练时间线性增加max_depth10限制每棵树的最大深度防止单棵树过拟合——毕设数据量通常不大深度超过 15 很容易把训练集背下来测试集效果反而下降min_samples_leaf3要求叶子节点至少 3 个样本这也是正则化手段。RMSE 的单位和油耗一样是 kg这个数字可以直接写进论文里横向对比同类研究的预测误差。跑完这个脚本会生成feature_importance.png这是随机森林自带的最强输出它直接告诉你哪些特征对油耗影响大。看这张图有个常见的理解误区特征重要性是“预测时用的多”不是“因果意义上的影响大”。如果两个特征高度相关重要性会在它们之间分摊导致排名看起来偏低。我在项目里会顺手做一个相关性矩阵确认这一点。4.2 线性回归基线与残差诊断regression_model.py做的事情是建立线性回归基线输出linear_regression_scatter.png和linear_regression_residual.png。线性回归的核心验证点在残差如果残差真实值减预测值随预测值增大而系统性地增大说明模型有异方差性线性假设不成立。from sklearn.linear_model import LinearRegression lin_model LinearRegression() lin_model.fit(X_train, y_train) lin_pred lin_model.predict(X_test) residuals y_test - lin_pred print(fLinear R²: {r2_score(y_test, lin_pred):.3f}) print(f残差均值: {residuals.mean():.3f})线性回归代码很短它的价值在于当对照组。如果随机森林 R² 是 0.92 而线性回归 R² 只有 0.71你可以在论文里明确写“燃料消耗与运行参数之间存在显著非线性关系机器学习方法优于传统线性建模”这个结论比空喊“我们用了随机森林所以结果好”有力得多。残差均值接近 0 说明模型没有系统性偏差散点图里如果残差呈现喇叭形分布说明高油耗段预测精度差这时可以检查训练数据里高油耗样本是否过少。项目里residuals_distribution.png画的则是残差的分布直方图接近正态分布说明模型误差是随机噪声主导没有学习不足的迹象。4.3 模型评估对比图与偏依赖图项目根目录下的model_evaluation_comparison.png把随机森林和线性回归的 R²、MAE、RMSE 画在同一张图上这张图就是答辩时的核心证据。读图时注意看两个模型的 R² 差距差距大说明非线性建模的必要性强差距很小比如都在 0.85 附近要反思是不是特征工程没做好或者目标变量本身和这些特征没有强关联而不是线性模型“够用了”。偏依赖图是本项目里容易被忽略但很有价值的一组输出partial_dependence_FuelConsumption_kg.png和partial_dependence_Annual_Time_spent_at_sea_hours.png。这两个文件是用sklearn.inspection.PartialDependenceDisplay生成的两变量交互偏依赖图X 轴是一个特征Y 轴是预测的油耗值画出来的曲线揭示了非线性阈值效应。比如航速偏依赖曲线通常在某个速度点出现拐点拐点之前油耗增长缓慢、拐点之后急剧上升这就是可以写进优化建议里的关键结论——把航速从那个拐点值降下来一截排放节省是边际收益最大的。from sklearn.inspection import PartialDependenceDisplay PartialDependenceDisplay.from_estimator( model, X_test, [Ship_speed], kindaverage, grid_resolution30 ) plt.show()grid_resolution30是采样点数30 已经足够画平滑曲线设太大会把随机噪声也画进去。对偏依赖图的一个重要提醒它是模型视角下的边际效应不是因果推断论文里描述时写成“模型显示”而不是“实验证明”。回归模型是个黑匣子偏依赖图是拆开黑匣子看内部决策逻辑的最好工具这组图的分析价值在第五章里还会继续展开。5. 避坑清单复现这个项目最容易翻车的五个地方把这份资源从下载到跑通再到把结果写进文档提交有几个坑是高概率出现的。我从自己和身边同学跑类似项目踩过的坑里挑五个最典型的按现象、原因、解决的顺序写清楚。这些不是“可能遇到的问题”是“一定会有人遇到的问题”。5.1 中文路径导致读取文件失败现象运行pd.read_csv(数据清洗/原始数据.csv)报FileNotFoundError或者UnicodeDecodeError但文件明明就在那个路径下。把代码改成全英文路径后一切正常。原因Windows 系统下 pandas 读取含中文或全角字符路径时编码处理不稳定尤其当文件名中混杂空格和中文时更容易触发另外脚本文件如果在深度较大的目录下相对路径计算也容易错位。解决项目目录统一用英文命名比如把“船舶数据清洗”改成ship_data_cleaning。如果已经跑过一遍生成了带中文名的中间文件直接用os.rename()在代码开头批量改名。还有一个小习惯所有脚本开头都加一句import os; os.chdir(os.path.dirname(os.path.abspath(__file__)))把工作目录切到脚本所在目录避免在 IDE 里能跑、终端里跑不了的路径错位问题。5.2 sklearn 版本差异导致 KMeans 参数报错现象KMeans(n_clusters3, n_init10)在较新版本上运行正常在旧版本上报TypeError: __init__() got an unexpected keyword argument n_init。或者反过来说旧代码里没有n_init参数新版本跑出来警告说n_init默认值要改。原因scikit-learn 的 KMeans 在 1.4 版本把n_init的默认值从 10 改成了 auto且早期版本并不存在显式传入n_init的写法。解决先检查环境里 sklearn 的版本python -c import sklearn; print(sklearn.__version__)1.4 以上版本n_initauto时算法会自动根据数据量决定运行次数显式写 10 也不会有问题。关键是在提交文档里写清楚“本机 Python 3.10 环境sklearn 1.4.0所有结果可复现”——版本写清楚老师复现时遇到参数警告至少知道不是你代码的问题。项目附带的 requirements 文件如果有的话按那个装没有的话 1.4.0 是当前最稳的选择。5.3 随机森林和 KMeans 的结果每次跑都不一样现象同一个脚本连续跑两次feature_importance.png里特征排名的形状大体相同但具体数值有差异KMeans 的 SSE 曲线也有轻微波动。原因随机森林的树生成过程、KMeans 的初始中心选择都依赖随机数生成器。项目代码里正确设了random_state42但如果你在重构时加了新的数据处理步骤比如先做一次 train_test_split 再做聚类sklearn 内部的随机状态就会变化。解决在提交前固定全局随机种子import numpy as np np.random.seed(42) # 固定 numpy 的随机状态单独看每个模型都设了random_state42这是项目代码已经做好的不要自己乱动。真正容易忽略的是train_test_split的随机状态——如果你的代码里没有传random_state每次划分的训练集不一样随机森林的特征重要性排序就可能变化。这属于那种“看起来跑通了但结果不能复现”的隐患写文档时一定要把关键随机种子列在实验环境里。5.4 聚类前没做标准化导致轮廓系数虚高现象直接对原始特征跑 KMeans轮廓系数显示 0.87看起来很漂亮但 PCA 降维散点图里各类别混在一起业务解释原性很差。原因FuelConsumption_kg和Annual_Time_spent_at_sea_hours的量级可能是几百到几千而Ship_speed只有个位数。KMeans 基于欧氏距离量级大的特征主导了距离计算聚类得到的“类别”实际上是按油耗高低简单分层而不是按运行模式综合划分。轮廓系数在这个场景下高是被大数值特征“带上去”的并不反映真实的聚类质量。解决用StandardScaler标准化后再聚类标准化后的轮廓系数会下降但 PCA 投影图会显示更合理的分组边界。项目里k_means_data_cleaning.py已经处理了标准化你自己改代码时不要为了省事直接拿清洗后的原始表喂 KMeans。判断标准化是否生效的方法是看聚类中心表标准化前的中心值差异悬殊标准化后的中心值都在同一量级范围内。5.5 matplotlib 图里中文乱码现象生成的feature_importance.png里特征名称如果是中文比如“船舶航速”“燃料消耗”显示成方块或者乱码。原因matplotlib 默认字体里没有中文字体Windows 里需要指定中文字体Linux 服务器上如果没装中文字体包也会乱码。解决在脚本开头强制设置字体import matplotlib.pyplot as plt plt.rcParams[font.sans-serif] [SimHei] # Windows 用黑体 plt.rcParams[axes.unicode_minus] False # 解决负号显示成方块的问题Linux 环境没有 SimHei可以换成Noto Sans CJK SC或者WenQuanYi Zen Hei。最省事的方案是把所有图片输出默认用英文列名——项目自带的输出图里特征名称用的就是FuelConsumption_kg这种英文命名说明作者是有意识避开了中文渲染问题你写报告时用中文描述图的内容即可不需要在图片里塞中文。6. 把聚类标签变成业务建议一份可答辩的结果组织方法跑通代码、看懂了图表最后一步是拿这些结果去答辩、去写文档、去体现数据驱动优化的价值。很多人栽在这代码跑得通、图也生成了一堆但老师说“你的优化建议是什么”就答不上来。这一章给一套通用的组织方法把clustered_data.csv和回归模型输出的结论串成一个完整的论证闭环。6.1 聚类标签的业务解释先给每类取个名字打开clustered_data.csv最后一列是聚类标签通常叫Cluster或者label。先用groupby计算每个聚类在核心特征上的平均值第三章 3.3 节写过代码然后按特征画像给每一类命名。命名规则是“状态词 特征词”不要用“第 1 类、第 2 类”这种泛化表述。聚类标签平均航速平均燃料消耗平均在航时间业务命名012.5 节820 kg/天6800 小时/年中速标准工况118.2 节1450 kg/天4200 小时/年高速重载工况28.1 节390 kg/天8900 小时/年低速长航工况这张表做完结论自然浮现高速重载工况的燃料消耗是低速工况的 3.7 倍但在航时间只有一半。优化方向不是“所有船都降速”而是对处于高速重载工况的航次做重点排查——是不是存在不必要的赶期、配载能不能优化。写在论文里的表述方式是“基于聚类结果识别出三类典型运行工况其中高速重载工况的排放强度最高是减排优化的首要目标”。6.2 把回归结果接到聚类结论上聚类告诉你哪类工况需要优化回归告诉你优化的抓手是什么。这时把随机森林的feature_importance.png拿过来按重要度从高到低排列特征挑前三个。如果Ship_speed排第一结论就是“航速是燃料消耗的首要预测因子与聚类结果中高速重载工况的高排放特征互相印证”如果Cargo_load排第一结论要改成“载重对排放影响最大建议从配载计划切入”。最后用偏依赖图的拐点值把“建议”具体化找到航速偏依赖曲线的拐点比如曲线显示航速超过 14.5 节时油耗急剧上升写成“将高速工况航速从 18 节降至 14.5 节按模型预测可降低约 X% 油耗”。X 可以用模型在测试集上做一次反事实预测把这部分样本的航速改成 14.5其他特征不变预测油耗差值取平均。这就是数据驱动优化的完整落地——不是拍脑袋说“降速减排”而是从模型里挖出阈值、用预测结果算收益。6.3 提交答辩时顺序怎么排我见过太多人把报告写成了“代码运行日志”顺序是装环境、跑清洗、跑聚类、跑回归、贴图。这个顺序是流水账不是论证。可答辩的结果组织顺序应该是先抛结论船舶碳排放存在明显的工况分化高速重载工况是减排重点航速是关键驱动变量再展示聚类证据轮廓系数、肘部曲线、聚类特征画像表接着展示回归证据随机森林与线性回归的精度对比、特征重要性、偏依赖图拐点最后给出可执行建议。图表之间用“那么问题来了”的逻辑串联聚类发现了问题对象回归定位了问题原因偏依赖图给出了问题边界反事实预测计算了收益。这个顺序读下来老师不可能问出“你的优化建议从哪来的”这种问题。做这个项目我最大的教训是一开始只顾着把代码跑通、把图攒齐觉得输出多了就是完成了。直到模拟答辩被问了一句“K 为什么选 4”才发现自己根本说不清。从那以后我每次做聚类分析都强制走一遍“标准化 → 肘部法则 → 轮廓系数 → 业务解释”四步每张图都必须对应一段能说出口的业务结论不许跳过。这份资源最大的价值也在这里它不仅给了你能跑的代码还给了你一套可以直接用的可视化分析流程——路径和结论都已经铺好你要做的不是重新发明一遍而是把手伸进去把每一步的逻辑打通。希望帮到你。本文还有配套的精品资源点击获取