 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1583 | 回復(fù): 5 | |||
zlp-lw木蟲(chóng) (正式寫(xiě)手)
|
[交流]
svm 預(yù)測(cè) 已有2人參與
|
|
各位大俠： 您好，我用30個(gè)案例分析中的第14章的svm回歸預(yù)測(cè)抑制劑的半濃度，并把其程序進(jìn)行改寫(xiě)，但是出現(xiàn)了一些意想不到的情況。特向您請(qǐng)教：在第14章的例子中，最后的測(cè)試為什么不用測(cè)試集，而是用原始訓(xùn)練數(shù)據(jù)做測(cè)試，這樣的測(cè)試效果很好，但是有什么效果呢？在我做的例子中，我用42個(gè)數(shù)據(jù)作為訓(xùn)練，10個(gè)數(shù)據(jù)作為測(cè)試集，對(duì)測(cè)試集我照著您這個(gè)例子中的訓(xùn)練集進(jìn)行了歸一化、反歸一化處理，最后用訓(xùn)練模型進(jìn)行測(cè)試，效果非常的差，我不知道這是什么原因？是我程序哪里出了問(wèn)題，您能幫我看一下嗎，我也30個(gè)案例分析的忠實(shí)讀者。謝謝您 xx=xlsread('1.xls'); y=xlsread('1-1.xls'); xx=xx'; sh=[y,xx]; %% 數(shù)據(jù)的提取和預(yù)處理 % 數(shù)據(jù)是一個(gè)52*6的double型的矩陣,每一行表示每一天的負(fù)荷數(shù)據(jù)和影響因素 % 提取數(shù)據(jù) [m,n] = size(x); ts = x(1:m,1); tsx = x(1:m,2:6); % 數(shù)據(jù)預(yù)處理,將原始數(shù)據(jù)進(jìn)行歸一化 ts = ts'; tsx = tsx'; % mapminmax為matlab自帶的映射函數(shù) % 對(duì)ts進(jìn)行歸一化 [TS,TSps] = mapminmax(ts,1,2); % 對(duì)TS進(jìn)行轉(zhuǎn)置,以符合libsvm工具箱的數(shù)據(jù)格式要求 TS = TS'; % mapminmax為matlab自帶的映射函數(shù) % 對(duì)tsx進(jìn)行歸一化 [TSX,TSXps] = mapminmax(tsx,1,2); % 對(duì)TSX進(jìn)行轉(zhuǎn)置,以符合libsvm工具箱的數(shù)據(jù)格式要求 TSX = TSX'; %% 選擇回歸預(yù)測(cè)分析最佳的SVM參數(shù)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網(wǎng)絡(luò)回歸預(yù)測(cè) [predict,mse] = svmpredict(TS,TSX,model); predict = mapminmax('reverse',predict',TSps); predict = predict'; % 打印回歸結(jié)果 str = sprintf( '均方誤差 MSE = %g 相關(guān)系數(shù) R = %g%%',mse(2),mse(3)*100); disp(str); %% 結(jié)果分析 figure; hold on; plot(ts,'-o'); plot(predict,'r-^'); legend('原始數(shù)據(jù)','回歸預(yù)測(cè)數(shù)據(jù)'); hold off; xlabel('小時(shí)','FontSize',12); ylabel('/Wh','FontSize',12); grid on; figure; error = predict - ts'; plot(error,'rd'); title('誤差圖(predicted data - original data)','FontSize',12); xlabel('小時(shí)','FontSize',12); ylabel('誤差量','FontSize',12); grid on; figure; error = (predict - ts')./ts'; plot(error,'rd'); title('相對(duì)誤差圖(predicted data - original data)/original data','FontSize',12); xlabel('小時(shí))','FontSize',12); ylabel('相對(duì)誤差量','FontSize',12); grid on; |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有63人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
zlp-lw
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2547.1
- 散金: 85
- 紅花: 6
- 帖子: 451
- 在線: 324.1小時(shí)
- 蟲(chóng)號(hào): 1502426
- 注冊(cè): 2011-11-21
- 專業(yè): 碳素材料與超硬材料
2樓2012-11-08 17:19:58
himyou
銀蟲(chóng) (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 503
- 帖子: 62
- 在線: 23.2小時(shí)
- 蟲(chóng)號(hào): 588041
- 注冊(cè): 2008-08-22
- 性別: GG
- 專業(yè): 生化分析及生物傳感
3樓2012-11-13 17:06:13
zlp-lw
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2547.1
- 散金: 85
- 紅花: 6
- 帖子: 451
- 在線: 324.1小時(shí)
- 蟲(chóng)號(hào): 1502426
- 注冊(cè): 2011-11-21
- 專業(yè): 碳素材料與超硬材料
|
我用SVMcgForRegress進(jìn)行參數(shù)尋優(yōu)，可是數(shù)據(jù)的輸入不一樣時(shí)預(yù)測(cè)的結(jié)果不一樣，下面是我的數(shù)據(jù)輸入的程序，幫我看看是什么原因吧？ %% 數(shù)據(jù)的提取和預(yù)處理 % 載入測(cè)試數(shù)據(jù)上證指數(shù)(1990.12.19-2009.08.19) % 數(shù)據(jù)是一個(gè)4579*6的double型的矩陣,每一行表示每一天的上證指數(shù) % 6列分別表示當(dāng)天上證指數(shù)的開(kāi)盤(pán)指數(shù),指數(shù)最高值,指數(shù)最低值,收盤(pán)指數(shù),當(dāng)日交易量,當(dāng)日交易額. 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,  ;% 畫(huà)出原始上證指數(shù)的每日開(kāi)盤(pán)數(shù) %figure; %plot(ts,'LineWidth',2); %title('上證指數(shù)的每日開(kāi)盤(pán)數(shù)(1990.12.20-2009.08.19)','FontSize',12); %grid on; % 數(shù)據(jù)預(yù)處理,將原始數(shù)據(jù)進(jìn)行歸一化 ts = ts'; tsx = tsx'; % mapminmax為matlab自帶的映射函數(shù) [TS,TSps] = mapminmax(ts); % 將映射函數(shù)的范圍參數(shù)分別置為1和2 TSps.ymin = 1; TSps.ymax = 2; % 對(duì)ts進(jìn)行歸一化 [TS,TSps] = mapminmax(ts,TSps); % 畫(huà)出原始上證指數(shù)的每日開(kāi)盤(pán)數(shù)歸一化后的圖像 figure; plot(TS,'LineWidth',2); title('原始上證指數(shù)的每日開(kāi)盤(pán)數(shù)歸一化后的圖像','FontSize',12); grid on; % 對(duì)TS進(jìn)行轉(zhuǎn)置,以符合libsvm工具箱的數(shù)據(jù)格式要求 TS = TS'; TS1 = TS(43:m,1); TS = TS(1:m-10,1); % mapminmax為matlab自帶的映射函數(shù) [TSX,TSXps] = mapminmax(tsx); % 將映射函數(shù)的范圍參數(shù)分別置為1和2 TSXps.ymin = 1; TSXps.ymax = 2; % 對(duì)tsx進(jìn)行歸一化 [TSX,TSXps] = mapminmax(tsx,TSXps); % 對(duì)TSX進(jìn)行轉(zhuǎn)置,以符合libsvm工具箱的數(shù)據(jù)格式要求 TSX = TSX'; TSX1 = TSX(43:m,1:n); TSX = TSX(1:m-10,1:n); 這個(gè)輸入，預(yù)測(cè)的效果是非常好的，其相關(guān)系數(shù)都達(dá)到了0.9996左右了，而用下面的輸入程序，則預(yù)測(cè)結(jié)果不太理想。。 %% 數(shù)據(jù)的提取和預(yù)處理 % 載入測(cè)試數(shù)據(jù)上證指數(shù)(1990.12.19-2009.08.19) % 數(shù)據(jù)是一個(gè)4579*6的double型的矩陣,每一行表示每一天的上證指數(shù) % 6列分別表示當(dāng)天上證指數(shù)的開(kāi)盤(pán)指數(shù),指數(shù)最高值,指數(shù)最低值,收盤(pán)指數(shù),當(dāng)日交易量,當(dāng)日交易額. 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, ;% 畫(huà)出原始上證指數(shù)的每日開(kāi)盤(pán)數(shù) %figure; %plot(ts,'LineWidth',2); %title('上證指數(shù)的每日開(kāi)盤(pán)數(shù)(1990.12.20-2009.08.19)','FontSize',12); %grid on; % 數(shù)據(jù)預(yù)處理,將原始數(shù)據(jù)進(jìn)行歸一化 ts = ts'; tsx = tsx'; % mapminmax為matlab自帶的映射函數(shù) [TS,TSps] = mapminmax(ts); % 將映射函數(shù)的范圍參數(shù)分別置為1和2 TSps.ymin = 1; TSps.ymax = 2; % 對(duì)ts進(jìn)行歸一化 [TS,TSps] = mapminmax(ts,TSps); % 畫(huà)出原始上證指數(shù)的每日開(kāi)盤(pán)數(shù)歸一化后的圖像 figure; plot(TS,'LineWidth',2); title('原始上證指數(shù)的每日開(kāi)盤(pán)數(shù)歸一化后的圖像','FontSize',12); grid on; % 對(duì)TS進(jìn)行轉(zhuǎn)置,以符合libsvm工具箱的數(shù)據(jù)格式要求 TS = TS'; TS1 = TS(43:m,1); TS = TS(1:m-10,1); % mapminmax為matlab自帶的映射函數(shù) [TSX,TSXps] = mapminmax(tsx); % 將映射函數(shù)的范圍參數(shù)分別置為1和2 TSXps.ymin = 1; TSXps.ymax = 2; % 對(duì)tsx進(jìn)行歸一化 [TSX,TSXps] = mapminmax(tsx,TSXps); % 對(duì)TSX進(jìn)行轉(zhuǎn)置,以符合libsvm工具箱的數(shù)據(jù)格式要求 TSX = TSX'; TSX1 = TSX(43:m,2:n); TSX = TSX(1:m-10,2:n); 希望您和各位能指點(diǎn)指點(diǎn)。。。 |
4樓2012-11-14 09:19:27
zlp-lw
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2547.1
- 散金: 85
- 紅花: 6
- 帖子: 451
- 在線: 324.1小時(shí)
- 蟲(chóng)號(hào): 1502426
- 注冊(cè): 2011-11-21
- 專業(yè): 碳素材料與超硬材料
at:5樓2012-11-14 11:45:07
★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
6樓2018-05-04 09:34:25
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zlp-lw 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +6 | tcxiaoxx 2026-03-20 | 6/300 |
|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +6 | 生物工程調(diào)劑 2026-03-17 | 10/500 |
|
|
|
[考研]
 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|