%%?第29章?支持向量機的回歸擬合——混凝土抗壓強度預測
%?
%?			該案例作者申明：		1：本人長期駐扎在此板塊里，對該案例提問，做到有問必答。	2：此案例有配套的教學視頻，視頻下載請點擊http://www.matlabsky.com/forum-91-1.html。?								3：此案例為原創案例，轉載請注明出處（《MATLAB智能算法30個案例分析》）。							4：若此案例碰巧與您的研究有關聯，我們歡迎您提意見，要求等，我們考慮后可以加在案例里。						5：以下內容為初稿，與實際發行的書籍內容略有出入，請以書籍中的內容為準。		
%?

%%?清空環境變量
clear?all
clc

%%?導入數據
load?concrete_data.mat
%?隨機產生訓練集和測試集
n?=?randperm（size（attributes2））;
%?訓練集——80個樣本
p_train?=?attributes（:n（1:80））‘;
t_train?=?strength（:n（1:80））‘;
%?測試集——23個樣本
p_test?=?attributes（:n（81:end））‘;
t_test?=?strength（:n（81: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?cg（ij）?????????????error?=?cg（ij）;
????????????bestc?=?2^c（ij）;
????????????bestg?=?2^g（ij）;
????????end
????????if?abs（cg（ij）?-?error）?<=?eps?&&?bestc?>?2^c（ij）
????????????error?=?cg（ij）;
????????????bestc?=?2^c（ij）;
????????????bestg?=?2^g（ij）;
????????end
????end
end
%?創建/訓練SVM??
cmd?=?[‘?-t?2‘‘?-c?‘num2str（bestc）‘?-g?‘num2str（bestg）‘?-s?3?-p?0.01‘];
model?=?svmtrain（tn_trainpn_traincmd）;

%%?SVM仿真預測
[Predict_1error_1]?=?svmpredict（tn_trainpn_trainmodel）;
[Predict_2error_2]?=?svmpredict（tn_testpn_testmodel）;
%?反歸一化
predict_1?=?mapminmax（‘reverse‘Predict_1outputps）;
predict_2?=?mapminmax（‘reverse‘Predict_2outputps）;
%?結果對比
result_1?=?[t_train?predict_1];
result_2?=?[t_test?predict_2];

%%?繪圖
figure（1）
plot（1:length（t_train）t_train‘r-*‘1:length（t_train）predict_1‘b:o‘）
grid?on
legend（‘真實值‘‘預測值‘）
xlabel（‘樣本編號‘）
ylabel（‘耐壓強度‘）
string_1?=?{‘訓練集預測結果對比‘;
???????????[‘mse?=?‘?num2str（error_1（2））?‘?R^2?=?‘?num2str（error_1（3））]};
title（string_1）
figure（2）
plot（1:length（t_test）t_test‘r-*‘1:length（t_test）predict_

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2286??2010-11-04?20:07??混凝土抗壓強度預測\concrete_data.mat
?????文件????????4841??2015-06-17?16:06??混凝土抗壓強度預測\main.m
?????目錄???????????0??2018-03-16?18:12??混凝土抗壓強度預測\