支持向量机回归预测代码

%% 第29章 支持向量机的回归拟合——混凝土抗压强度预测
% 
% 			该案例作者申明：		1：本人长期驻扎在此板块里，对该案例提问，做到有问必答。	2：此案例有配套的教学视频，视频下载请点击http://www.matlabsky.com/forum-91-1.html。 								3：此案例为原创案例，转载请注明出处（《MATLAB智能算法30个案例分析》）。							4：若此案例碰巧与您的研究有关联，我们欢迎您提意见，要求等，我们考虑后可以加在案例里。						5：以下内容为初稿，与实际发行的书籍内容略有出入，请以书籍中的内容为准。		
% 

%% 清空环境变量
clear all
clc

%% 导入数据
load JT.mat
% 随机产生训练集和测试集
n = randperm（size（T2））;
% 训练集——20个样本
p_train = T（:n（1:8））‘;
t_train = J（:n（1:8））‘;
% 测试集——5个样本
p_test = T（:n（9:end））‘;
t_test = J（:n（9:end））‘;

%% 数据归一化

% 训练集
[pn_traininputps] = mapminmax（p_train‘）;
pn_train = pn_train‘;
pn_test = mapminmax（‘apply‘p_test‘inputps）;
pn_test = pn_test‘;
% 测试集
[tn_trainoutputps] = mapminmax（t_train‘）;
tn_train = tn_train‘;
tn_test = mapminmax（‘apply‘t_test‘outputps）;
tn_test = tn_test‘;

%% SVM模型创建/训练

% 寻找最佳c参数/g参数
[cg] = meshgrid（-10:0.5:10-10:0.5:10）;
[mn] = size（c）;
cg = zeros（mn）;
eps = 10^（-4）;
v = 5;
bestc = 0;
bestg = 0;
error = Inf;
for i = 1:m
    for j = 1:n
        cmd = [‘-v ‘num2str（v）‘ -t 2‘‘ -c ‘num2str（2^c（ij））‘ -g ‘num2str（2^g（ij） ）‘ -s 3 -p 0.1‘];
        cg（ij） = svmtrain（tn_trainpn_traincmd）;
        if