%%?第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