MATLAB实现基于BMA-LSTM 贝叶斯模型平均BMA结合长短期记忆网络LSTM进行股票价格预测的详细项目实例请注意此篇内容只是一个项目介绍 更多详细内容可直接联系博主本人或者访问对应标题的完整博客或者文档下载页面含完整的程序GUI设计和代码详解金融市场价格序列具有显著的非线性、非平稳和高噪声特征单一预测模型往往难以在不同市场阶段和不同波动环境下保持稳定表现。传统时间序列模型如AR、ARIMA、GARCH等在刻画线性依赖结构、条件波动性方面具有一定优势但在捕捉长期记忆、非线性关系和高维特征交互方面存在明显局限。随着深度学习在时序建模中的应用逐步成熟长短期记忆网络LSTM在处理长序列依赖、非线性映射以及复杂模式识别方面展现出独特优势尤其适用于股票价格、指数、期货等金融时间序列的预测任务。单纯依赖LSTM这类深度模型虽然能够极大加强特征提取能力但其结果通常依赖于网络结构、训练超参数、训练初始权重等多种因素的配置。不同超参数组合甚至在相同数据集上会产生明显差异的预测性能同时还会出现一定程度的过拟合和模型不确定性问题。金融市场具有显著的结构性变化和外生冲击对模型稳健性提出更高要求单模型输出缺乏不确定性的量化表达在风险管理、资产配置与监管评估等场景下容易带来决策偏差。贝叶斯模型平均Bayesian Model Averaging, BMA提供了一种系统地综合多模型信息、量化模型不确定性的框架。通过在模型空间上引入先验分布并基于数据得到各候选模型后验概率BMA能够将多个候选模型按其后验概率进行加权从而得到更具鲁棒性的预测结果。与简单平均、加权平均不同BMA的权重具有严格的统计学意义既考虑模型对历史数据的拟合能力又自然惩罚过于复杂的模型在一定程度上平衡了拟合能力和泛化能力。在股票价格预测任务中将LSTM作为候选模型族的核心成员设置不同的网络结构如层数、隐藏单元数、dropout比例、不同的输入特征组合如只用收盘价、加入成交量、技术指标等、不同的时间窗口长度就可以形成一个庞大的候选模型集合。若能结合BMA对这些候选LSTM模型进行系统融合既可以充分利用深度网络对非线性的建模能力又可以通过贝叶斯方法显式引入模型不确定性从而在预测精度和稳定性之间获得更优秀的平衡。MATLAB在数值计算、矩阵操作和可视化方面具有成熟优势同时提供了深度学习工具箱和金融时间序列处理工具对于搭建BMA-LSTM混合框架非常适合。使用MATLAB R2025b环境一方面可以利用其对LSTM网络结构的直接支持另一方面可以方便地构建自定义的BMA权重计算模块对各个LSTM子模型的表现进行度量和更新。在该环境中将股票历史数据导入完成数据预处理、特征工程、训练集/验证集/测试集划分后可以构建多个LSTM模型并批量训练然后基于验证集误差构造近似的BMA权重再在测试集上输出加权预测。金融市场应用要求整个流程能够被稳定复现和调试。采用BMA-LSTM组合方案时可以通过统一的数据读取与预处理接口保证各个子模型共享同一数据基础通过统一的训练控制参数如批量大小、最大迭代次数、优化算法等减小训练过程中的不必要差异通过明确的误差度量指标如均方误差MSE、平均绝对误差MAE、对数似然近似等对各模型的表现进行公平比较。在此基础上将各个模型的误差映射为权重构建出稳健的贝叶斯模型平均预测器。在实践层面该项目不仅关注预测点值的准确性还注重整体预测框架的可扩展性和可维护性。未来可以在同一框架下继续扩展不同结构的深度网络如双向LSTM、GRU、卷积与LSTM混合结构等也可以加入更复杂的先验设定对模型权重的动态更新进行改进。同时该项目也兼顾教育和研究层面的需求通过详细的实现步骤和充分注释的MATLAB代码帮助从业人员以及研究人员系统理解BMA与LSTM融合的思路并在真实的股票数据上开展可控、可重复的实证研究。综上基于BMA-LSTM进行股票价格预测的项目将贝叶斯统计思想与深度学习时序模型紧密结合既回应了金融时间序列预测对精度和鲁棒性的双重需求又为研究者提供了一套可扩展、可深化的实现框架既可用于单一股票的短期预测也可拓展到指数、板块乃至跨市场资产的多维预测任务对金融科技实践具有现实意义和推广价值。项目目标与意义提升股票价格预测精度首要目标是在实际股票价格预测场景中显著提高预测精度和稳定性。传统线性模型难以刻画金融时间序列的非线性状态切换、突发行情和异方差特征而单一深度学习模型对超参数和初始化高度敏感预测效果容易波动。通过构建多组结构不同、输入特征不同的LSTM子模型并使用贝叶斯模型平均对各子模型的预测进行加权能够在模型空间层面进行多样化配置让不同模型在不同市场阶段、不同波动环境下贡献各自优势。表现较好的模型获得较高权重表现一般的模型获得较低权重从而在总体上降低预测误差和方差。该项目将通过在历史股票数据上的系统实验对比单一LSTM、简单平均、多层感知机等多种方案展示BMA-LSTM在均方误差、平均绝对误差和方向预测准确率等指标上带来的改进。同时会通过时间滚动窗口测试展示在不同市场阶段如上升、盘整、下跌中的稳健性力求减少预测模型在极端行情下的失效概率。引入模型不确定性量化与风险视角另一个关键目标是在预测过程中引入模型不确定性量化使预测不仅给出点值还在一定程度上反映预测区间和置信度。纯粹的点值预测无法体现模型对自身判断的信心不利于投资决策中的风险控制。贝叶斯模型平均提供了一种自然的方式可以通过各子模型的后验概率和预测分布组合构建综合预测分布。虽然在实践中会使用近似方法如基于误差转换为权重但仍然能够构造出具有统计意义的权重体系并由此估计预测区间、对尾部风险进行初步评估。在项目框架中可以在测试集阶段输出预测均值和一定置信区间的上下界并与真实价格进行比较帮助投资者更好地理解模型在不同时段的信心程度。这样的结果可以为风险管理提供辅助支撑为仓位控制、止损策略、VaR估计等提供更加丰富的信息。构建可扩展的深度学习与贝叶斯融合框架项目还致力于构建一个可扩展的技术框架使得后续可以方便地引入更多模型、更多特征和更复杂的贝叶斯权重更新策略。在MATLAB R2025b环境中可以通过模块化设计将数据预处理模块、LSTM子模型构建模块、训练模块、BMA权重计算模块、结果可视化模块分离便于替换和扩展。在这个框架中不同LSTM结构可以按配置文件批量生成未来可逐步加入GRU、双向LSTM以及卷积-循环混合结构甚至可以利用attentionLayer构建注意力增强的时序模型。同时在贝叶斯权重计算模块中可以由简单的基于误差的静态权重逐步扩展到时间可变的权重、分市场状态的分段权重等。通过这样的设计项目不仅解决当前单一股票预测任务也为后续研究和实务应用提供了可演化的平台。提升量化研究与MATLAB工程实践能力项目意义还体现在工程实践层面通过完整实现BMA-LSTM股票预测流程帮助从事量化分析和算法交易的人员提升在MATLAB环境下进行深度学习和贝叶斯统计融合建模的能力。项目中将涉及数据读取与清洗、特征工程、网络结构调参与训练过程可视化、模型集成方法实现、以及结果评估与可视化等一系列关键步骤。通过阅读和运行完整的MATLAB脚本可以更直观地理解深度学习工具箱的使用方式、时间序列数据在深度网络中的组织形式、训练选项的设置、以及如何在R2025b版本的约束下规避不兼容函数或参数。这种工程实践能力是把理论方法转化为可部署系统的重要基础对于构建实际投研平台、开发研究型交易策略具有直接推动作用同时也为后续在其他领域如宏观经济预测、能源负荷预测等的迁移应用打下良好基础。项目挑战及解决方案金融时间序列非平稳性与噪声问题股票价格序列通常表现出强烈的非平稳性包含趋势、季节性、结构性突变以及高频噪声等多种成分。价格本身往往带有随机游走特征简单利用历史价格直接预测未来价格容易造成误导。同时在宏观经济政策调整、重大事件、突发新闻等影响下价格序列会出现突发行情和跳变使得基于过去统计规律的模型难以适应。噪声成分会掩盖真正有价值的信号使深度学习模型在训练过程中更容易学习到噪声而非稳定模式从而导致对新数据的泛化能力下降。针对这一挑战项目在数据预处理中采取多层次方法。首要策略是对原始价格序列进行变换例如使用对数收益率、标准化收益率、价格差分等形式弱化趋势与非平稳性同时提升序列的平稳程度。其次在构建输入特征时通过引入技术指标移动平均线、RSI等、成交量以及波动率代理变量等帮助LSTM更好识别市场状态与价格变化之间的关联。还可以在窗口切分时采用重叠滑动窗口通过增加序列样本数量来提升训练稳定性。此外利用BMA框架在模型层面进行平滑不同结构的LSTM对于噪声的敏感程度不同通过贝叶斯加权可以在一定程度上降低单个模型对噪声敏感带来的不利影响从整体上提高对非平稳序列的适应能力。模型选择难度与超参数敏感性LSTM网络包含多种结构超参数如隐藏单元数、层数、dropout比例、序列长度、学习率等每个参数的变化都可能显著影响预测效果。传统做法往往需要大量人工调参或网格搜索既耗费时间又容易过拟合验证集。单一配置即使在某段时间、某只股票上表现优秀也未必在其他股票或其他阶段保持优势。面对这一挑战项目采用“多模型贝叶斯加权”的策略将模型选择问题转化为模型平均问题。具体做法是预先构建多个结构多样的LSTM子模型这些模型覆盖从较浅到较深、从较小到较大隐藏单元数的不同组合使其在模型空间中形成一个有代表性的子集。然后在验证集上评估各模型的性能并将误差指标映射为后验权重使具有更好泛化性能的模型自然获得更高权重表现较差的模型权重被压缩。这样即使没有为某一股票精确寻找单一“最优”模型也可以通过贝叶斯平均获得接近最优甚至更稳健的组合模型。该方案在MATLAB中易于实现通过统一的训练与评估接口批量生成各模型预测结果再基于简单的误差到权重映射公式实现权重计算这在工程实现和理论合理性之间取得平衡。贝叶斯模型平均实现与效率问题在理论上贝叶斯模型平均需要在模型空间上进行积分对每个候选模型计算后验概率需要先验分布、似然函数和边际似然等量的精确或近似计算。在深度学习场景下直接对LSTM这样高维参数空间计算精确后验既不现实也会导致极高的计算成本。项目需要在可接受的时间内完成多模型训练和BMA权重估计尤其考虑到MATLAB R2025b环境可能运行在研究用工作站而非大规模集群上。为此项目在实现策略上作出务实选择采用基于验证集误差的频率学近似方法构建贝叶斯权重。具体做法是使用例如exp(-λ*MSE)形式的函数把各模型的验证误差转换为非负权重再进行归一化使权重之和为1。这一过程等价于在某种特定似然形式和先验设定下的后验比例形式近似既保留了贝叶斯框架的核心思想又避免显式操作高维参数后验。同时通过适度控制候选模型数量例如6到10个和网络规模配合MATLAB的GPU训练如有条件可以在时间和精度之间找到满意的折中。在代码组织方面采用脚本化批量训练、多线程并行如果环境支持以及适当保存中间结果的方式减少重复计算提升整个BMA-LSTM系统的训练与评估效率使其既可用于研究实验也可为实务环境中的定期更新预测提供支撑。项目模型架构数据预处理与特征构建架构模型整体架构的起点是数据预处理与特征构建模块其核心目标是将原始股票历史数据转化为适合作为LSTM输入的标准格式并在此过程中增强特征信息表达。所使用的原始数据通常包括日期、开盘价、最高价、最低价、收盘价、成交量等字段可以从本地文件或数据库读取。首先需要对数据进行时间排序剔除缺失值异常值对连续缺失数据进行适当填补或删除处理。接着会对收盘价进行变换例如计算对数收益率log(C_t/C_{t-1})或者使用收益率和价格共同构成多变量时间序列以减轻非平稳性。特征构建阶段根据项目目标引入技术指标例如N日移动平均、指数移动平均、相对强弱指标、价格振幅等这些指标能够从不同角度刻画趋势、动量和波动状态。为了提高模型训练稳定性需要对各特征进行标准化或归一化处理将其缩放到类似的尺度范围有利于LSTM在反向传播中保持梯度稳定。此外数据还需要组织为监督学习形式将过去T个时间步构成一个序列输入预测T1时刻或未来若干时刻的价格或收益从而形成输入张量和标签向量。该模块的设计直接影响后续LSTM模型能否充分学习有效模式也是BMA-LSTM整体架构的基础层。多LSTM子模型架构在完成基础特征构建后架构的第二层是多LSTM子模型模块。该模块包含若干个结构各异的LSTM网络每个网络在输入维度一致的前提下通过调整网络深度、隐藏单元数、dropout比例、回归层结构等形成差异。其基本原理是LSTM单元在传统RNN的基础上引入输入门、遗忘门和输出门通过门控机制控制信息在时间维度上的传递与遗忘从而有效缓解长序列训练中的梯度消失问题。在具体设计时可以包含诸如“浅层小宽度网络”“浅层大宽度网络”“两层堆叠网络”“带较高dropout的防过拟合网络”等多种结构使这些网络在表达能力和泛化能力上形成互补。在MATLAB中使用sequenceInputLayer、lstmLayer、dropoutLayer、fullyConnectedLayer、regressionLayer等构建序列回归网络每个子模型共享统一的数据输入格式与训练选项便于后续批量训练和比较。通过这样的多模型结构设计使得整个系统不再依赖某一个特定LSTM结构的性能而是在模型空间上建立一个丰富的候选集为后续贝叶斯权重分配提供更多维度的信息来源。贝叶斯模型平均权重计算架构第三个重要层次是贝叶斯模型平均权重计算模块。其核心思想来源于贝叶斯统计学中的模型平均理论即在存在多个候选模型时预测结果应当是各模型预测的加权和而权重由模型的后验概率决定。理论上需要预先给出模型先验概率并基于数据计算边际似然从而得到后验概率。但在深度网络场景中精确计算边际似然非常困难因此采用基于验证集误差的频率学近似方案。具体实现思路为在训练阶段将每个LSTM子模型在训练集上拟合后利用独立的验证集评估模型预测误差采用例如均方误差或平均绝对误差指标进行度量随后把误差通过某种单调递减函数映射为非负权重分数例如exp(-α*MSE)使误差小的模型得到更高分数。该映射过程可以被视为某种隐含似然函数与先验的组合近似参数α控制对误差差异的敏感程度。最终对所有模型的分数进行归一化使其和为1就得到各模型的贝叶斯权重近似。该架构使得在不显式计算高维贝叶斯后验的情况下仍保留贝叶斯模型平均的核心精神为整体输出提供理论合理的权重基础。BMA-LSTM预测集成与输出架构当多LSTM子模型训练完成且权重计算完毕后第四层架构负责将各模型预测结果集成并生成最终预测输出。在测试阶段对每个LSTM子模型分别输入相同的测试集序列数据得到多个预测序列。然后按照预先计算的贝叶斯权重对这些预测进行加权求和得到BMA-LSTM组合预测序列。这个过程本质上是对不同模型预测的线性组合但由于权重来自验证集表现的贝叶斯式估计因此相较于简单平均具有更明确的统计解释。同时若对子模型预测残差的方差进行估计还可以利用加权方差公式构建预测区间进一步体现不确定性信息。该模块还可以扩展为动态权重更新形式例如在滚动窗口环境中定期使用最近一段数据更新权重使模型加权随时间变化更加贴近当前市场状态。输出架构不仅包含数值预测结果还可集成评估指标计算和可视化例如绘制真实价格与预测价格曲线对不同模型与BMA组合进行对比为使用者提供直观的评估依据。整体系统与MATLAB工程组织架构上述各模块需要在MATLAB环境下形成一个连贯一致的工程架构以便于开发、调试和扩展。在文件组织上可以将数据预处理、子模型定义、训练配置、BMA权重计算和预测评估分别构建为相对独立的代码区域或脚本段通过主脚本进行统一调度。MATLAB R2025b对于界面类函数、部分统计学习函数的参数存在一定变更因此在工程架构中尽量采用核心稳定的函数接口避免使用已弃用或行为改变的特性。在训练过程中通过trainingOptions配置训练算法如adam、最大轮数、mini-batch大小、学习率初值等并利用内置的训练记录信息监控损失收敛情况。模型训练结束后通过统一命名规则保存和加载各个子模型便于批量预测与权重计算。整体架构注重清晰、可维护、可扩展为后续加入新的子模型、引入更多特征、调整BMA策略提供便利。这样构建的BMA-LSTM系统不仅完成当前股票价格预测任务也为进一步的研究和应用提供了稳定的技术底座。项目模型描述及代码示例数据读取与基础预处理示例 clear; clc; % 清空工作区变量并清除命令行窗口保证脚本在干净环境下运行 dataTable readtable(stock_data.csv); % 从本地CSV文件中读取股票历史数据表格包含日期和价格等字段 dataTable sortrows(dataTable,Date); % 按日期列对数据进行升序排序确保时间序列顺序正确 dataTable rmmissing(dataTable); % 删除包含缺失值的行避免缺失数据影响后续模型训练 closePrice dataTable.Close; % 提取收盘价列作为主要价格序列用于构建预测目标 closePrice closePrice(:); % 将收盘价转换为列向量格式便于后续矩阵运算和索引 ret diff(log(closePrice)); % 计算相邻时刻对数收益率以减弱价格序列的非平稳特征 ret(isinf(ret) | isnan(ret)) 0; % 将对数收益率中的无穷值与NaN值替换为0避免异常值传播 vol dataTable.Volume(2:end); % 提取对应收益率长度的成交量数据从第二个样本开始对齐 vol vol(:); % 将成交量转换为列向量与收益率保持一致的维度形式 featMat [ret, vol]; % 将收益率和成交量拼接成特征矩阵每行对应一个时间点的多维特征 mu mean(featMat,1); % 计算特征矩阵每一列的均值用于后续标准化中心化操作 sigma std(featMat,[],1) 1e-8; % 计算特征矩阵每一列的标准差并加上微小正数防止除零错误 featNorm (featMat - mu) ./ sigma; % 对特征矩阵进行标准化处理使各特征具有零均值和单位量级 seqLen 20; % 设置LSTM输入序列长度为20个时间步表示以过去20天数据预测未来 numSample size(featNorm,1) - seqLen; % 计算可构建的样本数量扣除形成最后一个序列所需时间步 X cell(numSample,1); % 初始化输入序列单元数组每个单元存储一个时间序列样本 Y cell(numSample,1); % 初始化输出目标单元数组每个单元存储对应的预测值 for i 1:numSample % 遍历所有可用样本索引逐个生成输入序列和预测目标 X{i} featNorm(i:iseqLen-1,:).; % 取连续seqLen行特征并转置为[特征数×时间步]格式作为LSTM输入 Y{i} ret(iseqLen); % 使用序列结束后的第一个收益率作为对应的预测目标 end % 结束样本构建循环完成全部输入输出序列的组织 numTrain floor(0.7 * numSample); % 将70%的样本划为训练集控制训练数据占比 numVal floor(0.15 * numSample); % 将15%的样本划为验证集用于模型选择和权重估计 idxTrain 1:numTrain; % 训练集索引范围从第一个样本到numTrain idxVal numTrain1:numTrainnumVal; % 验证集索引范围紧随训练集样本之后 idxTest numTrainnumVal1:numSample; % 测试集索引范围使用剩余样本进行最终评估 XTrain X(idxTrain); % 从总样本中选取训练集输入序列 YTrain Y(idxTrain); % 从总样本中选取训练集目标序列 XVal X(idxVal); % 从总样本中选取验证集输入序列 YVal Y(idxVal); % 从总样本中选取验证集目标序列 XTest X(idxTest); % 从总样本中选取测试集输入序列 YTest Y(idxTest); % 从总样本中选取测试集目标序列 多个LSTM子模型结构定义示例 numFeatures size(featNorm,2); % 计算输入特征数量等于标准化特征矩阵的列数 numResponses 1; % 设置输出响应维度为1这里预测单一收益率值 layers1 [ ... % 定义第一个LSTM子模型的层结构使用方括号数组表示顺序网络 sequenceInputLayer(numFeatures) ... % 序列输入层输入维度为numFeatures对应特征数 lstmLayer(32,OutputMode,last) ... % LSTM层包含32个隐藏单元仅输出最后时间步的隐藏状态 fullyConnectedLayer(16) ... % 全连接层输出维度为16用于进一步特征映射 reluLayer ... % ReLU激活层引入非线性变换以增强表达能力 fullyConnectedLayer(numResponses) ... % 全连接层将特征映射为单一预测值 regressionLayer]; % 回归输出层用于连续数值预测和损失计算 layers2 [ ... % 定义第二个LSTM子模型的层结构采用不同隐藏单元规模 sequenceInputLayer(numFeatures) ... % 序列输入层与第一个模型保持相同输入维度 lstmLayer(64,OutputMode,last) ... % LSTM层隐藏单元数增加到64以提升表达能力 dropoutLayer(0.3) ... % dropout层随机丢弃30%单元以减轻过拟合 fullyConnectedLayer(32) ... % 全连接层输出维度为32作为中间特征层 reluLayer ... % ReLU激活层增强非线性拟合能力 fullyConnectedLayer(numResponses) ... % 全连接层输出预测值 regressionLayer]; % 回归层定义损失为均方误差 layers3 [ ... % 定义第三个LSTM子模型的层结构采用两层堆叠的LSTM sequenceInputLayer(numFeatures) ... % 序列输入层接受多维特征输入 lstmLayer(32,OutputMode,sequence) ... % 第一层LSTM输出完整序列作为下一层输入 lstmLayer(16,OutputMode,last) ... % 第二层LSTM输出最后时间步的隐藏状态压缩时序信息 fullyConnectedLayer(16) ... % 全连接层将高维特征映射到16维空间 reluLayer ... % ReLU激活层引入非线性变换 fullyConnectedLayer(numResponses) ... % 输出层全连接将特征映射到单值预测 regressionLayer]; % 回归层完成损失定义 lstmModels {layers1, layers2, layers3}; % 将三个模型的层结构存入单元数组以便统一管理 numModels numel(lstmModels); % 记录子模型总数量用于后续循环训练和权重计算 训练选项配置与子模型训练示例 miniBatchSize 64; % 设置每个训练批次包含64个序列样本平衡收敛速度与计算负载 maxEpochs 60; % 设置最大训练轮数为60以保证模型有充分学习机会 options trainingOptions(adam, ... % 选择Adam优化算法适用于非凸问题的自适应学习率方法 MaxEpochs,maxEpochs, ... % 指定最大训练轮数为maxEpochs MiniBatchSize,miniBatchSize, ... % 指定mini-batch大小为miniBatchSize InitialLearnRate,1e-3, ... % 设置初始学习率为0.001在收敛速度和稳定性之间折中 GradientThreshold,1, ... % 限制梯度范数不超过1以避免梯度爆炸 Shuffle,every-epoch, ... % 每个epoch打乱训练数据顺序提高泛化能力 ValidationData,{XVal,YVal}, ... % 指定验证集数据用于监控泛化性能 ValidationFrequency,20, ... % 每训练20个mini-batch进行一次验证评估 Plots,none, ... % 不在训练过程中绘制训练曲线以减少界面开销 Verbose,false); % 关闭详细训练日志输出使命令行更加简洁 trainedNets cell(numModels,1); % 初始化单元数组用于存储训练好的LSTM子模型网络 valMSE zeros(numModels,1); % 初始化向量用于记录每个模型在验证集上的均方误差 for k 1:numModels % 遍历每一个子模型索引对各自结构进行训练和评估 layersK lstmModels{k}; % 取出第k个子模型的层结构定义 netK trainNetwork(XTrain,YTrain,layersK,options); % 使用训练集和指定训练选项训练第k个网络 trainedNets{k} netK; % 将训练完成的网络对象存入数组便于后续预测和集成 YValPred predict(netK,XVal,MiniBatchSize,miniBatchSize); % 使用验证集输入对第k个网络进行预测 YValTrue cell2mat(YVal.); % 将验证集真实目标单元数组拼接为列向量方便计算误差 YValPredVec YValPred(:); % 将预测结果转为列向量确保与真实值维度一致 valMSE(k) mean((YValPredVec - YValTrue).^2); % 计算第k个模型在验证集上的均方误差作为性能指标 end % 完成所有子模型的训练和验证误差计算 基于验证误差的贝叶斯权重计算示例 lambda 50; % 设置误差到权重映射的灵敏度参数lambda控制对MSE差异的放大程度 score exp(-lambda * valMSE); % 使用指数函数将各模型的验证MSE映射为非负得分误差越小得分越大 score(isnan(score) | isinf(score)) 0; % 将得分中的NaN或无穷值置零防止异常值影响权重归一化 if all(score 0) % 检查所有得分是否都为零避免后续归一化出现除零情况 score ones(size(score)); % 若全部为零则将得分统一设为1相当于均等权重起点 end % 结束异常得分检查与修正 weights score / sum(score); % 将得分向量归一化使各模型权重之和为1 weights reshape(weights,[],1); % 将权重向量整理为列向量形式便于矩阵运算 disp(weights); % 在命令行显示各个子模型的贝叶斯平均权重便于观察权重分布情况 测试集预测与BMA集成示例 numTest numel(XTest); % 计算测试集样本数量用于后续预测矩阵维度分配 predAll zeros(numTest,numModels); % 初始化预测结果矩阵每列对应一个子模型在测试集上的预测 for k 1:numModels % 遍历每个子模型对测试集进行预测 netK trainedNets{k}; % 取出第k个训练好的网络 YPredK predict(netK,XTest,MiniBatchSize,miniBatchSize); % 使用测试集输入预测第k个模型的输出 predAll(:,k) YPredK(:); % 将第k个模型的预测结果展平后存入预测矩阵的第k列 end % 完成所有子模型在测试集上的预测 bmaPred predAll * weights; % 使用贝叶斯权重向量对各模型预测进行加权求和得到BMA组合预测 YTestTrue cell2mat(YTest.); % 将测试集真实目标值单元数组拼接为列向量 mseBMA mean((bmaPred - YTestTrue).^2); % 计算BMA组合预测与真实值之间的均方误差 maeBMA mean(abs(bmaPred - YTestTrue)); % 计算BMA组合预测与真实值之间的平均绝对误差 disp(mseBMA); % 显示BMA组合预测的MSE指标用于评价模型精度 disp(maeBMA); % 显示BMA组合预测的MAE指标辅助评估预测误差水平 figure; % 创建新图窗用于可视化测试集真实值与预测值的对比 plot(YTestTrue,b); % 绘制真实对数收益率序列用蓝色线表示 hold on; % 保持当前图像使后续曲线叠加在同一坐标轴上 plot(bmaPred,r); % 绘制BMA组合预测序列用红色线表示 legend({True,BMA-LSTM},Location,best); % 添加图例标识真实序列和BMA预测序列 title(Test Set True vs BMA-LSTM Prediction); % 设置图形标题说明展示内容为测试集真实值与BMA预测对比 xlabel(Time Index); % 设置横轴标签为时间索引表示样本次序 ylabel(Log Return); % 设置纵轴标签为对数收益率直观显示预测目标量纲 colormap(gcf,turbo); % 为当前图窗设置turbo颜色映射兼容R2025b的colormap使用规范 简单预测区间构造示例 residuals predAll - YTestTrue; % 计算每个子模型在测试集上的残差矩阵行对应样本列对应模型 varModel var(residuals,0,1); % 按列计算各模型残差的方差得到每个模型的误差方差估计 bmaVar (weights.^2) * varModel.; % 使用权重平方对方差进行加权求和得到BMA预测方差近似 bmaStd sqrt(bmaVar); % 对BMA预测方差取平方根得到预测标准差 z 1.96; % 设定正态近似的1.96倍标准差对应95%置信区间系数 upperBand bmaPred z * bmaStd; % 计算预测上置信界为预测均值加上z倍标准差 lowerBand bmaPred - z * bmaStd; % 计算预测下置信界为预测均值减去z倍标准差 figure; % 创建新图窗展示预测区间与真实序列 plot(YTestTrue,k); % 绘制真实收益率用黑色线表示 hold on; % 保持当前图形叠加预测均值与区间边界 plot(bmaPred,b); % 绘制BMA预测均值用蓝色线表示 plot(upperBand,r--); % 绘制上置信界用红色虚线表示 plot(lowerBand,r--); % 绘制下置信界同样用红色虚线表示 legend({True,BMA mean,Upper 95%,Lower 95%},Location,best); % 添加图例区分真实值和预测区间 title(BMA-LSTM Prediction with Approximate Confidence Bands); % 设置图形标题说明展示内容为预测均值及置信区间 xlabel(Time Index); % 设置横轴标签表示时间索引 ylabel(Log Return); % 设置纵轴标签表示对数收益率 colormap(gcf,turbo); % 为当前图窗设置turbo颜色映射符合R2025b色图使用规则更多详细内容请访问http://金融预测MATLAB实现基于BMA-LSTM贝叶斯模型平均BMA结合长短期记忆网络LSTM进行股票价格预测的详细项目实例含完整的程序GUI设计和代码详解_随机森林预测模型matlab资源-CSDN下载 https://download.csdn.net/download/xiaoxingkongyuxi/90233666https://download.csdn.net/download/xiaoxingkongyuxi/90233666https://download.csdn.net/download/xiaoxingkongyuxi/90233666