matlab30个案例分析案例14-SVM神经网络的回归预测分析
- 格式:docx
- 大小:19.42 KB
- 文档页数:5
在Matlab中实现神经网络的方法与实例神经网络是一种模拟人类大脑神经系统的计算模型,它能够通过学习数据的模式和关联性来解决各种问题。
在计算机科学和人工智能领域,神经网络被广泛应用于图像识别、自然语言处理、预测等任务。
而Matlab作为一种功能强大的科学计算软件,提供了一套完善的工具箱,可以方便地实现神经网络的建模和训练。
本文将介绍在Matlab中实现神经网络的方法与实例。
首先,我们会简要介绍神经网络的基本原理和结构,然后详细讲解在Matlab中如何创建并训练神经网络模型,最后通过几个实例展示神经网络在不同领域的应用。
一、神经网络的原理和结构神经网络模型由神经元和它们之间的连接构成。
每个神经元接收输入信号,并通过权重和偏置进行加权计算,然后使用激活函数对结果进行非线性变换。
这样,神经网络就能够模拟复杂的非线性关系。
常见的神经网络结构包括前馈神经网络(Feedforward Neural Network)和循环神经网络(Recurrent Neural Network)。
前馈神经网络是最基本的结构,信号只能向前传递,输出不对网络进行反馈;而循环神经网络具有反馈连接,可以对自身的输出进行再处理,适用于序列数据的建模。
神经网络的训练是通过最小化损失函数来优化模型的参数。
常用的训练算法包括梯度下降法和反向传播算法。
其中,梯度下降法通过计算损失函数对参数的梯度来更新参数;反向传播算法是梯度下降法在神经网络中的具体应用,通过反向计算梯度来更新网络的权重和偏置。
二、在Matlab中创建神经网络模型在Matlab中,可以通过Neural Network Toolbox来创建和训练神经网络模型。
首先,我们需要定义神经网络的结构,包括输入层、隐藏层和输出层的神经元数量,以及每个神经元之间的连接权重。
例如,我们可以创建一个三层的前馈神经网络模型:```matlabnet = feedforwardnet([10 8]);```其中,`[10 8]`表示隐藏层的神经元数量分别为10和8。
MATLAB 智能算法30 个案例分析(终极版)1 基于遗传算法的TSP算法(王辉)基于遗传算法和非线性规划的函数寻优算法(史峰)基于遗传算法的 BP神经网络优化算法(王辉)设菲尔德大学的 MATLAB遗传算法工具箱(王辉)基于遗传算法的 LQR控制优化算法(胡斐)遗传算法工具箱详解及应用(胡斐)多种群遗传算法的函数优化算法(王辉)基于量子遗传算法的函数寻优算法(王辉)多目标 Pareto 最优解搜索算法(胡斐)基于多目标 Pareto 的二维背包搜索算法(史峰)基于免疫算法的柔性车间调度算法(史峰)基于免疫算法的运输中心规划算法(史峰)基于粒子群算法的函数寻优算法(史峰)基于粒子群算法的 PID 控制优化算法(史峰)基于混合粒子群算法的 TSP寻优算法(史峰)基于动态粒子群算法的动态环境寻优算法(史峰)粒子群算法工具箱(史峰)基于鱼群算法的函数寻优算法(王辉)基于模拟退火算法的 TSP算法(王辉)基于遗传模拟退火算法的聚类算法(王辉)基于模拟退火算法的 HEV能量管理策略参数优化(胡斐)22 蚁群算法的优化计算——旅行商问题(TSP)优化(郁磊)基于蚁群算法的二维路径规划算法(史峰)基于蚁群算法的三维路径规划算法(史峰)有导师学习神经网络的回归拟合——基于近红外光谱的汽油辛烷值预测(郁磊)有导师学习神经网络的分类——鸢尾花种类识别(郁磊)无导师学习神经网络的分类——矿井突水水源判别(郁磊)支持向量机的分类——基于乳腺组织电阻抗特性的乳腺癌诊断(郁磊)支持向量机的回归拟合——混凝土抗压强度预测(郁磊)极限学习机的回归拟合及分类——对比实验研究(郁磊)智能算法是我们在学习中经常遇到的算法,主要包括遗传算法,免疫算法,粒子群算法,神经网络等,智能算法对于很多人来说,既爱又恨,爱是因为熟练的掌握几种智能算法,能够很方便的解决我们的论坛问题,恨是因为智能算法感觉比较“玄乎” ,很难理解,更难用它来解决问题。
Mat1ab技术神经网络优化应用案例随着人工智能和机器学习的快速发展,神经网络成为了解决复杂问题的重要工具。
而在神经网络中,优化算法则是提高网络性能和准确性的关键。
本文将以Mauab技术为基础,介绍神经网络优化的应用案例。
一、神经网络与优化算法神经网络是一种模仿人脑神经系统的计算模型,通过模拟神经元之间的连接和信息传递来实现人工智能。
而优化算法则是指通过调整网络中的权重和偏置等参数,使得网络的输出结果尽可能地接近期望值。
在神经网络的训练过程中,优化算法起到了至关重要的作用。
常见的优化算法包括梯度下降法、遗传算法、粒子群算法等。
而Mauab作为一种强大的数学工具和编程语言,为神经网络的优化提供了丰富的支持和工具。
二、基于MatIab的神经网络优化案例1.图像识别在图像识别领域,神经网络具有广泛的应用。
例如,我们可以利用神经网络识别手写数字。
首先,我们需要提供大量的手写数字样本作为训练集,然后使用MaUab中的神经网络工具箱搭建一个适当的神经网络结构。
接下来,我们可以使用优化算法对网络中的权重和偏置进行调优,以提高识别准确率。
通过迭代优化过程,不断调整网络参数,使得网络能够对手写数字进行准确的判断和识别。
2.金融市场预测神经网络在金融市场的预测也有广泛的应用。
以股票市场为例,我们可以利用神经网络来预测股票价格的变化趋势。
首先,我们需要收集一段时间内的历史股票价格数据,作为训练集。
然后,利用Mat1ab中的神经网络工具箱搭建一个能够接受历史数据并输出价格预测的神经网络结构。
接着,我们可以使用优化算法对网络参数进行调整,使得网络能够更好地拟合历史数据,并对未来的价格变化作出准确的预测。
3.机器人控制在机器人控制领域,神经网络可以用来优化机器人的运动过程。
例如,我们可以利用神经网络来控制机器人的路径规划和运动轨迹。
首先,我们需要将机器人的运动环境建模,并将其转化为输入-输出映射的训练数据。
接着,利用Mauab中的神经网络工具箱搭建一个能够接受环境输入并输出机器人运动指令的神经网络结构。
MATLAB神经网络43个案例分析简介神经网络(Neural Network)是一种模拟人类神经元行为的计算模型,它通过对大量输入数据进行训练,学习到输入和输出之间的复杂关系。
MATLAB是一个强大的数学计算工具,具有丰富的神经网络函数和工具箱,可以用于神经网络的设计、训练和应用。
本文将介绍43个使用MATLAB进行神经网络分析的案例,主要包括神经网络的基本概念、神经网络模型的建立、参数的调整和优化等方面。
二级标题1: 基本概念在开始具体的案例分析之前,首先理解神经网络的基本概念是非常重要的。
三级标题1: 神经元神经网络的基本单元是神经元(Neuron),它模拟了生物神经元的工作原理。
神经元接收多个输入信号,并通过一个激活函数产生输出信号。
常用的激活函数包括Sigmoid函数、ReLU函数等。
三级标题2: 神经网络的结构神经网络由多层神经元组成,通常包括输入层、隐藏层和输出层。
输入层接收外部输入,隐藏层用于处理中间信息,输出层产生最终的输出。
三级标题3: 前向传播和反向传播神经网络的训练过程主要包括前向传播和反向传播两个步骤。
前向传播是通过输入数据按照网络结构进行计算,得到输出结果。
反向传播是根据输出结果与真实结果之间的误差,通过调整网络参数来提高模型性能。
二级标题2: 案例分析三级标题4: 案例1: 图像分类本案例通过使用神经网络和MATLAB工具箱,对手写数字图像进行分类。
首先,将图像数据转化为向量形式,并通过网络进行训练。
然后,将训练好的网络用于分类未知图像,并评估分类性能。
三级标题5: 案例2: 时序预测本案例使用神经网络来预测时间序列数据。
通过对历史数据进行训练,建立时序模型,并利用该模型来预测未来的数据。
通过调整网络结构和参数,提高预测准确性。
三级标题6: 案例3: 异常检测本案例利用神经网络进行异常检测。
通过对正常数据进行训练,建立正常模型,并使用该模型检测异常数据。
通过调整网络的灵敏度和阈值,提高异常检测的性能。
matlab智能算法30个案例分析Matlab智能算法30个案例分析。
Matlab作为一种强大的数学软件,拥有丰富的算法库和强大的编程能力,能够实现各种复杂的智能算法。
本文将针对Matlab智能算法进行30个案例分析,帮助读者深入了解Matlab在智能算法领域的应用和实践。
1. 遗传算法。
遗传算法是一种模拟自然选择和遗传机制的优化算法,能够有效解决复杂的优化问题。
在Matlab中,可以利用遗传算法工具箱快速实现各种优化问题的求解,例如函数最小化、参数优化等。
2. 神经网络。
神经网络是一种模拟人脑神经元网络的计算模型,能够实现复杂的非线性映射和模式识别。
Matlab提供了丰富的神经网络工具箱,可以用于神经网络的建模、训练和应用,例如分类、回归、聚类等任务。
3. 模糊逻辑。
模糊逻辑是一种处理不确定性和模糊信息的逻辑推理方法,能够有效处理模糊规则和模糊数据。
Matlab中的模糊逻辑工具箱提供了丰富的模糊推理方法和工具,可以用于模糊控制、模糊识别等领域。
4. 粒子群算法。
粒子群算法是一种模拟鸟群觅食行为的优化算法,能够有效处理多维优化问题。
在Matlab中,可以利用粒子群算法工具箱快速实现各种优化问题的求解,例如函数最小化、参数优化等。
5. 蚁群算法。
蚁群算法是一种模拟蚂蚁觅食行为的优化算法,能够有效处理离散优化问题和组合优化问题。
Matlab中的蚁群算法工具箱提供了丰富的蚁群优化方法和工具,可以用于解决各种组合优化问题。
6. 遗传规划算法。
遗传规划算法是一种结合遗传算法和规划算法的优化方法,能够有效处理复杂的规划问题。
在Matlab中,可以利用遗传规划算法工具箱快速实现各种规划问题的求解,例如路径规划、资源分配等。
7. 人工免疫算法。
人工免疫算法是一种模拟免疫系统的优化算法,能够有效处理多峰优化问题和动态优化问题。
在Matlab中,可以利用人工免疫算法工具箱快速实现各种复杂的优化问题的求解。
8. 蜂群算法。
MATLAB_智能算法30个案例分析1.线性回归:使用MATLAB的回归工具箱,对给定的数据集进行线性回归分析,获取拟合的直线方程。
2.逻辑回归:使用MATLAB的分类工具箱,对给定的数据集进行逻辑回归分析,建立分类模型。
3.K均值聚类:使用MATLAB的聚类工具箱,对给定的数据集进行K 均值聚类算法,将数据集分为多个簇。
4.支持向量机:使用MATLAB的SVM工具箱,对给定的数据集进行支持向量机算法,建立分类或回归模型。
5.决策树:使用MATLAB的分类工具箱,对给定的数据集进行决策树分析,建立决策模型。
6.随机森林:使用MATLAB的分类和回归工具箱,对给定的数据集进行随机森林算法,集成多个决策树模型。
7. AdaBoost:使用MATLAB的分类工具箱,对给定的数据集进行AdaBoost算法,提升分类性能。
8.遗传算法:使用MATLAB的全局优化工具箱,利用遗传算法进行优化问题的求解。
9.粒子群优化:使用MATLAB的全局优化工具箱,利用粒子群优化算法进行优化问题的求解。
10.模拟退火算法:使用MATLAB的全局优化工具箱,利用模拟退火算法进行优化问题的求解。
11.神经网络:使用MATLAB的神经网络工具箱,构建和训练多层感知机模型。
12.卷积神经网络:使用MATLAB的深度学习工具箱,构建和训练卷积神经网络模型。
13.循环神经网络:使用MATLAB的深度学习工具箱,构建和训练循环神经网络模型。
14.长短期记忆网络:使用MATLAB的深度学习工具箱,构建和训练长短期记忆网络模型。
15.GAN(生成对抗网络):使用MATLAB的深度学习工具箱,构建和训练生成对抗网络模型。
16.自编码器:使用MATLAB的深度学习工具箱,构建和训练自编码器模型。
17.强化学习:使用MATLAB的强化学习工具箱,构建和训练强化学习模型。
18.关联规则挖掘:使用MATLAB的数据挖掘工具箱,发现数据中的关联规则。
MATLAB_智能算法30个案例分析MATLAB是一种强大的数值计算和编程工具,教育和科研领域中广泛应用于数据分析、机器学习和智能算法的研究。
在本文中,我们将介绍30个MATLAB智能算法的案例分析,并探讨其用途和优势。
分析的案例包括分类、回归、聚类、神经网络和遗传算法等不同类型的智能算法。
1. K均值聚类:利用MATLAB中的kmeans函数对一组数据进行聚类分析,得到不同的簇。
2. 随机森林:利用MATLAB中的TreeBagger函数构建一个随机森林分类器,并通过测试数据进行分类预测。
3. 人工神经网络:使用MATLAB中的feedforwardnet函数构建一个人工神经网络,并通过训练集进行预测。
4. 遗传算法:利用MATLAB中的ga函数对一个优化问题进行求解,找到最优解。
5. 支持向量机:使用MATLAB中的svmtrain和svmclassify函数构建一个支持向量机分类器,并进行分类预测。
6. 极限学习机:使用MATLAB中的elmtrain和elmpredict函数构建一个极限学习机分类器,并进行分类预测。
7. 逻辑回归:使用MATLAB中的mnrfit和mnrval函数构建一个逻辑回归模型,并进行预测。
8. 隐马尔可夫模型:使用MATLAB中的hmmtrain和hmmdecode函数构建一个隐马尔可夫模型,对一系列观测数据进行预测。
9. 神经进化算法:利用MATLAB中的ne_train函数构建一个基于神经进化算法的神经网络分类器,并进行分类预测。
10. 朴素贝叶斯分类器:使用MATLAB中的NaiveBayes对象构建一个朴素贝叶斯分类器,并进行分类预测。
11. 高斯过程回归:使用MATLAB中的fitrgp函数构建一个高斯过程回归模型,并进行回归预测。
12. 最小二乘支持向量机:使用MATLAB中的fitcsvm函数构建一个最小二乘支持向量机分类器,并进行分类预测。
13. 遗传网络:利用MATLAB中的ngenetic函数构建一个基于遗传算法和人工神经网络的分类器,并进行分类预测。
MATLAB智能算法30个案例分析1.线性回归:通过拟合数据,预测未知的连续变量。
2.逻辑回归:基于已知输入和输出数据,通过对数斯蒂格回归模型,进行二元分类。
3.决策树:通过对已知数据进行分类预测,构建一棵决策树模型。
4.随机森林:通过构建多个决策树模型,进行分类和回归分析。
5.支持向量机:通过找到一个最优超平面,对数据进行二元分类。
6.高斯混合模型:基于多个高斯分布,对数据进行聚类分析。
7.K均值聚类:通过对数据进行分组,找到数据的簇结构。
8.主成分分析:找到最具有代表性的主成分,实现数据的降维和可视化。
9.独立成分分析:在多变量数据中,找到相互独立的成分。
10.关联规则挖掘:通过分析大规模数据集,找到数据项之间的关联规则。
11.朴素贝叶斯分类器:基于贝叶斯理论,进行分类和预测。
12.遗传算法:通过模拟进化过程,找到最优解。
13.粒子群算法:通过模拟粒子在空间中的移动,优化问题的解。
14.蚁群算法:通过模拟蚂蚁在空间中的行为,解决优化问题。
15.神经网络:通过多层神经元之间的连接,进行模式识别和预测。
16.卷积神经网络:通过卷积层和池化层,进行图像分类和目标检测。
17.循环神经网络:通过循环连接,进行时间序列预测和自然语言处理。
18.支持张量分解的非负矩阵分解:通过分解张量,进行数据降维和特征提取。
19.马尔科夫链:通过状态转移概率,对随机过程进行建模和分析。
20.K最近邻算法:通过找到与未知样本最接近的训练样本,进行分类和回归分析。
21.高斯过程回归:利用高斯过程进行回归分析和置信区间估计。
22.隐马尔科夫模型:通过观测序列推断隐藏状态序列。
23.时序聚类:通过对时间序列数据进行聚类分析,找到相似的模式。
24.大规模机器学习:通过将数据划分为小批量,进行机器学习模型的训练。
25.非线性最小二乘:通过最小化非线性函数的残差平方和,拟合数据。
26.分类集成学习:通过结合多个分类器的预测结果,提高分类准确率。
%% SVM神经网络中的参数优化---如何更好的提升分类器的性能%% 清空环境变量function chapter13_GAclose all;clear;clc;format compact;%% 数据提取% 载入测试数据wine,其中包含的数据为classnumber = 3,wine:178*13的矩阵,wine_labes:178*1的列向量load chapter13_wine.mat;% 画出测试数据的box可视化图figure;boxplot(wine,'orientation','horizontal','labels',categories);title('wine数据的box可视化图','FontSize',12);xlabel('属性值','FontSize',12);grid on;% 画出测试数据的分维可视化图figuresubplot(3,5,1);hold onfor run = 1:178plot(run,wine_labels(run),'*');endxlabel('样本','FontSize',10);ylabel('类别标签','FontSize',10);title('class','FontSize',10);for run = 2:14subplot(3,5,run);hold on;str = ['attrib ',num2str(run-1)];for i = 1:178plot(i,wine(i,run-1),'*');endxlabel('样本','FontSize',10);ylabel('属性值','FontSize',10);title(str,'FontSize',10);end% 选定训练集和测试集% 将第一类的1-30,第二类的60-95,第三类的131-153做为训练集train_wine = [wine(1:30,:);wine(60:95,:);wine(131:153,:)];% 相应的训练集的标签也要分离出来train_wine_labels = [wine_labels(1:30);wine_labels(60:95);wine_labels(131:153)];% 将第一类的31-59,第二类的96-130,第三类的154-178做为测试集test_wine = [wine(31:59,:);wine(96:130,:);wine(154:178,:)];% 相应的测试集的标签也要分离出来test_wine_labels = [wine_labels(31:59);wine_labels(96:130);wine_labels(154:178)];%% 数据预处理% 数据预处理,将训练集和测试集归一化到[0,1]区间[mtrain,ntrain] = size(train_wine);[mtest,ntest] = size(test_wine);dataset = [train_wine;test_wine];% mapminmax为MATLAB自带的归一化函数[dataset_scale,ps] = mapminmax(dataset',0,1);dataset_scale = dataset_scale';train_wine = dataset_scale(1:mtrain,:);test_wine = dataset_scale( (mtrain+1):(mtrain+mtest),: );%% 选择GA最佳的SVM参数c&g% GA的参数选项初始化ga_option.maxgen = 200;ga_option.sizepop = 20;ga_option.cbound = [0,100];ga_option.gbound = [0,100];ga_option.v = 5;ga_option.ggap = 0.9;[bestacc,bestc,bestg] = gaSVMcgForClass(train_wine_labels,train_wine,ga_option);% 打印选择结果disp('打印选择结果');str = sprintf( 'Best Cross Validation Accuracy = %g%% Best c = %g Best g = %g',bestacc,bestc,bestg);disp(str);%% 利用最佳的参数进行SVM网络训练cmd = ['-c ',num2str(bestc),' -g ',num2str(bestg)];model = svmtrain(train_wine_labels,train_wine,cmd);%% SVM网络预测[predict_label,accuracy] = svmpredict(test_wine_labels,test_wine,model);% 打印测试集分类准确率total = length(test_wine_labels);right = sum(predict_label == test_wine_labels);disp('打印测试集分类准确率');str = sprintf( 'Accuracy = %g%% (%d/%d)',accuracy(1),right,total);disp(str);%% 结果分析% 测试集的实际分类和预测分类图figure;hold on;plot(test_wine_labels,'o');plot(predict_label,'r*');xlabel('测试集样本','FontSize',12);ylabel('类别标签','FontSize',12);legend('实际测试集分类','预测测试集分类');title('测试集的实际分类和预测分类图','FontSize',12);grid on;snapnow;%% 子函数gaSVMcgForClass.mfunction [BestCVaccuracy,Bestc,Bestg,ga_option] = gaSVMcgForClass(train_label,train_data,ga_option)% gaSVMcgForClass%% by faruto%Email:patrick.lee@ QQ:516667408 /faruto BNU%last modified 2010.01.17%Super Moderator @ % 若转载请注明:% faruto and liyang , LIBSVM-farutoUltimateVersion% a toolbox with implements for support vector machines based on libsvm, 2009.% Software available at %% Chih-Chung Chang and Chih-Jen Lin, LIBSVM : a library for% support vector machines, 2001. Software available at% .tw/~cjlin/libsvm% 参数初始化if nargin == 2ga_option = struct('maxgen',200,'sizepop',20,'ggap',0.9,...'cbound',[0,100],'gbound',[0,1000],'v',5);end% maxgen:最大的进化代数,默认为200,一般取值范围为[100,500]% sizepop:种群最大数量,默认为20,一般取值范围为[20,100]% cbound = [cmin,cmax],参数c的变化范围,默认为(0,100]% gbound = [gmin,gmax],参数g的变化范围,默认为[0,1000]% v:SVM Cross Validation参数,默认为5%MAXGEN = ga_option.maxgen;NIND = ga_option.sizepop;NVAR = 2;PRECI = 20;GGAP = ga_option.ggap;trace = zeros(MAXGEN,2);FieldID = ...[rep([PRECI],[1,NVAR]);[ga_option.cbound(1),ga_option.gbound(1);ga_option.cbound(2),ga_opti on.gbound(2)]; ...[1,1;0,0;0,1;1,1]];Chrom = crtbp(NIND,NVAR*PRECI);gen = 1;v = ga_option.v;BestCVaccuracy = 0;Bestc = 0;Bestg = 0;%cg = bs2rv(Chrom,FieldID);for nind = 1:NINDcmd = ['-v ',num2str(v),' -c ',num2str(cg(nind,1)),' -g ',num2str(cg(nind,2))];ObjV(nind,1) = svmtrain(train_label,train_data,cmd);end[BestCVaccuracy,I] = max(ObjV);Bestc = cg(I,1);Bestg = cg(I,2);for gen = 1:MAXGENFitnV = ranking(-ObjV);SelCh = select('sus',Chrom,FitnV,GGAP);SelCh = recombin('xovsp',SelCh,0.7);SelCh = mut(SelCh);cg = bs2rv(SelCh,FieldID);for nind = 1:size(SelCh,1)cmd = ['-v ',num2str(v),' -c ',num2str(cg(nind,1)),' -g ',num2str(cg(nind,2))];ObjVSel(nind,1) = svmtrain(train_label,train_data,cmd);end[Chrom,ObjV] = reins(Chrom,SelCh,1,1,ObjV,ObjVSel);if max(ObjV) <= 50continue;end[NewBestCVaccuracy,I] = max(ObjV);cg_temp = bs2rv(Chrom,FieldID);temp_NewBestCVaccuracy = NewBestCVaccuracy;if NewBestCVaccuracy > BestCVaccuracyBestCVaccuracy = NewBestCVaccuracy;Bestc = cg_temp(I,1);Bestg = cg_temp(I,2);endif abs( NewBestCVaccuracy-BestCVaccuracy ) <= 10^(-2) && ...cg_temp(I,1) < BestcBestCVaccuracy = NewBestCVaccuracy;Bestc = cg_temp(I,1);Bestg = cg_temp(I,2);endtrace(gen,1) = max(ObjV);trace(gen,2) = sum(ObjV)/length(ObjV);end%figure;hold on;trace = round(trace*10000)/10000;plot(trace(1:gen,1),'r*-','LineWidth',1.5);plot(trace(1:gen,2),'o-','LineWidth',1.5);legend('最佳适应度','平均适应度',3);xlabel('进化代数','FontSize',12);ylabel('适应度','FontSize',12);axis([0 gen 0 100]);grid on;axis auto;line1 = '适应度曲线Accuracy[GAmethod]';line2 = ['(终止代数=', ...num2str(gen),',种群数量pop=', ...num2str(NIND),')'];line3 = ['Best c=',num2str(Bestc),' g=',num2str(Bestg), ...' CVAccuracy=',num2str(BestCVaccuracy),'%'];title({line1;line2;line3},'FontSize',12);。
%% SVM神经网络的回归预测分析---上证指数开盘指数预测%%% 清空环境变量function chapter14tic;close all;clear;clc;format compact;%% 数据的提取和预处理% 载入测试数据上证指数(1990.12.19-2009.08.19)% 数据是一个4579*6的double型的矩阵,每一行表示每一天的上证指数% 6列分别表示当天上证指数的开盘指数,指数最高值,指数最低值,收盘指数,当日交易量,当日交易额.load chapter14_sh.mat;% 提取数据[m,n] = size(sh);ts = sh(2:m,1);tsx = sh(1:m-1,:);% 画出原始上证指数的每日开盘数figure;plot(ts,'LineWidth',2);title('上证指数的每日开盘数(1990.12.20-2009.08.19)','FontSize',12);xlabel('交易日天数(1990.12.19-2009.08.19)','FontSize',12);ylabel('开盘数','FontSize',12);grid on;% 数据预处理,将原始数据进行归一化ts = ts';tsx = tsx';% mapminmax为matlab自带的映射函数% 对ts进行归一化[TS,TSps] = mapminmax(ts,1,2);% 画出原始上证指数的每日开盘数归一化后的图像figure;plot(TS,'LineWidth',2);title('原始上证指数的每日开盘数归一化后的图像','FontSize',12);xlabel('交易日天数(1990.12.19-2009.08.19)','FontSize',12);ylabel('归一化后的开盘数','FontSize',12);grid on;% 对TS进行转置,以符合libsvm工具箱的数据格式要求TS = TS';% mapminmax为matlab自带的映射函数% 对tsx进行归一化[TSX,TSXps] = mapminmax(tsx,1,2);% 对TSX进行转置,以符合libsvm工具箱的数据格式要求TSX = TSX';%% 选择回归预测分析最佳的SVM参数c&g% 首先进行粗略选择:[bestmse,bestc,bestg] = SVMcgForRegress(TS,TSX,-8,8,-8,8);% 打印粗略选择结果disp('打印粗略选择结果');str = sprintf( 'Best Cross Validation MSE = %g Best c = %g Best g = %g',bestmse,bestc,bestg); disp(str);% 根据粗略选择的结果图再进行精细选择:[bestmse,bestc,bestg] = SVMcgForRegress(TS,TSX,-4,4,-4,4,3,0.5,0.5,0.05);% 打印精细选择结果disp('打印精细选择结果');str = sprintf( 'Best Cross Validation MSE = %g Best c = %g Best g = %g',bestmse,bestc,bestg); disp(str);%% 利用回归预测分析最佳的参数进行SVM网络训练cmd = ['-c ', num2str(bestc), ' -g ', num2str(bestg) , ' -s 3 -p 0.01'];model = svmtrain(TS,TSX,cmd);%% SVM网络回归预测[predict,mse] = svmpredict(TS,TSX,model);predict = mapminmax('reverse',predict',TSps);predict = predict';% 打印回归结果str = sprintf( '均方误差MSE = %g 相关系数R = %g%%',mse(2),mse(3)*100);disp(str);%% 结果分析figure;hold on;plot(ts,'-o');plot(predict,'r-^');legend('原始数据','回归预测数据');hold off;title('原始数据和回归预测数据对比','FontSize',12);xlabel('交易日天数(1990.12.19-2009.08.19)','FontSize',12);ylabel('开盘数','FontSize',12);grid on;figure;error = predict - ts';plot(error,'rd');title('误差图(predicted data - original data)','FontSize',12);xlabel('交易日天数(1990.12.19-2009.08.19)','FontSize',12);ylabel('误差量','FontSize',12);grid on;figure;error = (predict - ts')./ts';plot(error,'rd');title('相对误差图(predicted data - original data)/original data','FontSize',12);xlabel('交易日天数(1990.12.19-2009.08.19)','FontSize',12);ylabel('相对误差量','FontSize',12);grid on;snapnow;toc;%% 子函数SVMcgForRegress.mfunction [mse,bestc,bestg] = SVMcgForRegress(train_label,train,cmin,cmax,gmin,gmax,v,cstep,gstep,msestep)%SVMcg cross validation by faruto%% by faruto%Email:*************************:516667408/farutoBNU%last modified 2010.01.17%******************************.cn% 若转载请注明:% faruto and liyang , LIBSVM-farutoUltimateVersion% a toolbox with implements for support vector machines based on libsvm, 2009.% Software available at %% Chih-Chung Chang and Chih-Jen Lin, LIBSVM : a library for% support vector machines, 2001. Software available at% .tw/~cjlin/libsvm% about the parameters of SVMcgif nargin < 10msestep = 0.06;endif nargin < 8cstep = 0.8;gstep = 0.8;endif nargin < 7v = 5;endif nargin < 5gmax = 8;gmin = -8;endif nargin < 3cmax = 8;cmin = -8;end% X:c Y:g cg:acc[X,Y] = meshgrid(cmin:cstep:cmax,gmin:gstep:gmax);[m,n] = size(X);cg = zeros(m,n);eps = 10^(-4);bestc = 0;bestg = 0;mse = Inf;basenum = 2;for i = 1:mfor j = 1:ncmd = ['-v ',num2str(v),' -c ',num2str( basenum^X(i,j) ),' -g ',num2str( basenum^Y(i,j) ),' -s 3 -p 0.1'];cg(i,j) = svmtrain(train_label, train, cmd);if cg(i,j) < msemse = cg(i,j);bestc = basenum^X(i,j);bestg = basenum^Y(i,j);endif abs( cg(i,j)-mse )<=eps && bestc > basenum^X(i,j)mse = cg(i,j);bestc = basenum^X(i,j);bestg = basenum^Y(i,j);endendend% to draw the acc with different c & g[cg,ps] = mapminmax(cg,0,1);figure;[C,h] = contour(X,Y,cg,0:msestep:0.5);clabel(C,h,'FontSize',10,'Color','r');xlabel('log2c','FontSize',12);ylabel('log2g','FontSize',12);firstline = 'SVR参数选择结果图(等高线图)[GridSearchMethod]'; secondline = ['Best c=',num2str(bestc),' g=',num2str(bestg), ...' CVmse=',num2str(mse)];title({firstline;secondline},'Fontsize',12);grid on;figure;meshc(X,Y,cg);% mesh(X,Y,cg);% surf(X,Y,cg);axis([cmin,cmax,gmin,gmax,0,1]);xlabel('log2c','FontSize',12);ylabel('log2g','FontSize',12);zlabel('MSE','FontSize',12);firstline = 'SVR参数选择结果图(3D视图)[GridSearchMethod]'; secondline = ['Best c=',num2str(bestc),' g=',num2str(bestg), ...' CVmse=',num2str(mse)];title({firstline;secondline},'Fontsize',12);。