 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1586 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
zlp-lw木蟲 (正式寫手)
|
[交流]
svm 預(yù)測(cè) 已有2人參與
|
||
|
各位大俠： 您好，我用30個(gè)案例分析中的第14章的svm回歸預(yù)測(cè)抑制劑的半濃度，并把其程序進(jìn)行改寫，但是出現(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è)試，效果非常的差，我不知道這是什么原因？是我程序哪里出了問題，您能幫我看一下嗎，我也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)有145人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(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ù)
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
6樓2018-05-04 09:34:25
zlp-lw
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2547.1
- 散金: 85
- 紅花: 6
- 帖子: 451
- 在線: 324.1小時(shí)
- 蟲號(hào): 1502426
- 注冊(cè): 2011-11-21
- 專業(yè): 碳素材料與超硬材料
2樓2012-11-08 17:19:58
himyou
銀蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 503
- 帖子: 62
- 在線: 23.2小時(shí)
- 蟲號(hào): 588041
- 注冊(cè): 2008-08-22
- 性別: GG
- 專業(yè): 生化分析及生物傳感
3樓2012-11-13 17:06:13
zlp-lw
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2547.1
- 散金: 85
- 紅花: 6
- 帖子: 451
- 在線: 324.1小時(shí)
- 蟲號(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ù)的開盤指數(shù),指數(shù)最高值,指數(shù)最低值,收盤指數(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,  ;% 畫出原始上證指數(shù)的每日開盤數(shù) %figure; %plot(ts,'LineWidth',2); %title('上證指數(shù)的每日開盤數(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); % 畫出原始上證指數(shù)的每日開盤數(shù)歸一化后的圖像 figure; plot(TS,'LineWidth',2); title('原始上證指數(shù)的每日開盤數(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ù)的開盤指數(shù),指數(shù)最高值,指數(shù)最低值,收盤指數(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, ;% 畫出原始上證指數(shù)的每日開盤數(shù) %figure; %plot(ts,'LineWidth',2); %title('上證指數(shù)的每日開盤數(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); % 畫出原始上證指數(shù)的每日開盤數(shù)歸一化后的圖像 figure; plot(TS,'LineWidth',2); title('原始上證指數(shù)的每日開盤數(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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 284求調(diào)劑 +5 | Zhao anqi 2026-03-22 | 5/250 |
|
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|