| 查看: 1599 | 回复: 5 | |||
zlp-lw木虫 (正式写手)
|
[交流]
svm 预测 已有2人参与
|
|
各位大侠： 您好，我用30个案例分析中的第14章的svm回归预测抑制剂的半浓度，并把其程序进行改写，但是出现了一些意想不到的情况。特向您请教：在第14章的例子中，最后的测试为什么不用测试集，而是用原始训练数据做测试，这样的测试效果很好，但是有什么效果呢？在我做的例子中，我用42个数据作为训练，10个数据作为测试集，对测试集我照着您这个例子中的训练集进行了归一化、反归一化处理，最后用训练模型进行测试，效果非常的差，我不知道这是什么原因？是我程序哪里出了问题，您能帮我看一下吗，我也30个案例分析的忠实读者。谢谢您 xx=xlsread('1.xls'); y=xlsread('1-1.xls'); xx=xx'; sh=[y,xx]; %% 数据的提取和预处理 % 数据是一个52*6的double型的矩阵,每一行表示每一天的负荷数据和影响因素 % 提取数据 [m,n] = size(x); ts = x(1:m,1); tsx = x(1:m,2:6); % 数据预处理,将原始数据进行归一化 ts = ts'; tsx = tsx'; % mapminmax为matlab自带的映射函数 % 对ts进行归一化 [TS,TSps] = mapminmax(ts,1,2); % 对TS进行转置,以符合libsvm工具箱的数据格式要求 TS = TS'; % mapminmax为matlab自带的映射函数 % 对tsx进行归一化 [TSX,TSXps] = mapminmax(tsx,1,2); % 对TSX进行转置,以符合libsvm工具箱的数据格式要求 TSX = TSX'; %% 选择回归预测分析最佳的SVM参数c&g bestc=1;bestg=1.6245; bestc,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; xlabel('小时','FontSize',12); ylabel('/Wh','FontSize',12); grid on; figure; error = predict - ts'; plot(error,'rd'); title('误差图(predicted data - original data)','FontSize',12); xlabel('小时','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('小时)','FontSize',12); ylabel('相对误差量','FontSize',12); grid on; |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有173人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
zlp-lw
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 2547.1
- 散金: 85
- 红花: 6
- 帖子: 451
- 在线: 324.1小时
- 虫号: 1502426
- 注册: 2011-11-21
- 专业: 碳素材料与超硬材料
2楼2012-11-08 17:19:58
3楼2012-11-13 17:06:13
zlp-lw
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 2547.1
- 散金: 85
- 红花: 6
- 帖子: 451
- 在线: 324.1小时
- 虫号: 1502426
- 注册: 2011-11-21
- 专业: 碳素材料与超硬材料
|
我用SVMcgForRegress进行参数寻优，可是数据的输入不一样时预测的结果不一样，下面是我的数据输入的程序，帮我看看是什么原因吧？ %% 数据的提取和预处理 % 载入测试数据上证指数(1990.12.19-2009.08.19) % 数据是一个4579*6的double型的矩阵,每一行表示每一天的上证指数 % 6列分别表示当天上证指数的开盘指数,指数最高值,指数最低值,收盘指数,当日交易量,当日交易额. xx=xlsread('TXFZMSHF1.xls'); y=xlsread('ShiyanzhiDuishu1.xls'); xx=xx'; sh=[y,xx]; [m,n] = size(sh); ts = sh(1:m,1); tsx = sh(1:m,  ;% 画出原始上证指数的每日开盘数 %figure; %plot(ts,'LineWidth',2); %title('上证指数的每日开盘数(1990.12.20-2009.08.19)','FontSize',12); %grid on; % 数据预处理,将原始数据进行归一化 ts = ts'; tsx = tsx'; % mapminmax为matlab自带的映射函数 [TS,TSps] = mapminmax(ts); % 将映射函数的范围参数分别置为1和2 TSps.ymin = 1; TSps.ymax = 2; % 对ts进行归一化 [TS,TSps] = mapminmax(ts,TSps); % 画出原始上证指数的每日开盘数归一化后的图像 figure; plot(TS,'LineWidth',2); title('原始上证指数的每日开盘数归一化后的图像','FontSize',12); grid on; % 对TS进行转置,以符合libsvm工具箱的数据格式要求 TS = TS'; TS1 = TS(43:m,1); TS = TS(1:m-10,1); % mapminmax为matlab自带的映射函数 [TSX,TSXps] = mapminmax(tsx); % 将映射函数的范围参数分别置为1和2 TSXps.ymin = 1; TSXps.ymax = 2; % 对tsx进行归一化 [TSX,TSXps] = mapminmax(tsx,TSXps); % 对TSX进行转置,以符合libsvm工具箱的数据格式要求 TSX = TSX'; TSX1 = TSX(43:m,1:n); TSX = TSX(1:m-10,1:n); 这个输入，预测的效果是非常好的，其相关系数都达到了0.9996左右了，而用下面的输入程序，则预测结果不太理想。。 %% 数据的提取和预处理 % 载入测试数据上证指数(1990.12.19-2009.08.19) % 数据是一个4579*6的double型的矩阵,每一行表示每一天的上证指数 % 6列分别表示当天上证指数的开盘指数,指数最高值,指数最低值,收盘指数,当日交易量,当日交易额. xx=xlsread('1.xls'); y=xlsread('1-1.xls'); xx=xx'; sh=[y,xx]; [m,n] = size(sh); ts = sh(1:m,1); tsx = sh(1:m, ;% 画出原始上证指数的每日开盘数 %figure; %plot(ts,'LineWidth',2); %title('上证指数的每日开盘数(1990.12.20-2009.08.19)','FontSize',12); %grid on; % 数据预处理,将原始数据进行归一化 ts = ts'; tsx = tsx'; % mapminmax为matlab自带的映射函数 [TS,TSps] = mapminmax(ts); % 将映射函数的范围参数分别置为1和2 TSps.ymin = 1; TSps.ymax = 2; % 对ts进行归一化 [TS,TSps] = mapminmax(ts,TSps); % 画出原始上证指数的每日开盘数归一化后的图像 figure; plot(TS,'LineWidth',2); title('原始上证指数的每日开盘数归一化后的图像','FontSize',12); grid on; % 对TS进行转置,以符合libsvm工具箱的数据格式要求 TS = TS'; TS1 = TS(43:m,1); TS = TS(1:m-10,1); % mapminmax为matlab自带的映射函数 [TSX,TSXps] = mapminmax(tsx); % 将映射函数的范围参数分别置为1和2 TSXps.ymin = 1; TSXps.ymax = 2; % 对tsx进行归一化 [TSX,TSXps] = mapminmax(tsx,TSXps); % 对TSX进行转置,以符合libsvm工具箱的数据格式要求 TSX = TSX'; TSX1 = TSX(43:m,2:n); TSX = TSX(1:m-10,2:n); 希望您和各位能指点指点。。。 |
4楼2012-11-14 09:19:27
zlp-lw
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 2547.1
- 散金: 85
- 红花: 6
- 帖子: 451
- 在线: 324.1小时
- 虫号: 1502426
- 注册: 2011-11-21
- 专业: 碳素材料与超硬材料
at:5楼2012-11-14 11:45:07
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
6楼2018-05-04 09:34:25
