| 24小時(shí)熱門版塊排行榜 |
| 1 | 1/1 | 返回列表 |
| 查看: 2395 | 回復(fù): 0 | |||
刀光木蟲(chóng) (著名寫手)
鐵蟲(chóng)
|
[求助]
matlab 神經(jīng)網(wǎng)絡(luò)3分類代碼求助!
|
|
本人現(xiàn)在在做畢業(yè)設(shè)計(jì),做了一下支持向量機(jī)3分類和智能算法優(yōu)化支持向量機(jī)參數(shù),還想做一下SVM分類和神經(jīng)網(wǎng)絡(luò),決策樹(shù),KJ近鄰等算法的分類準(zhǔn)確率的比較,由于時(shí)間比較緊急,神經(jīng)網(wǎng)絡(luò)以前沒(méi)接觸過(guò),網(wǎng)上找到一個(gè)參考程序,語(yǔ)音信號(hào)4分類,我的數(shù)據(jù)是放在一個(gè)文本文檔里面的210行9列(第一列是類別標(biāo)簽)的數(shù)據(jù),要實(shí)現(xiàn)我的三分類,代碼要怎么改?求大神相助,感激不盡。! 網(wǎng)上參考代碼如下: web browser http://www.ilovematlab.cn/thread-60056-1-1.html %% 清空環(huán)境變量 clc clear %% 訓(xùn)練數(shù)據(jù)預(yù)測(cè)數(shù)據(jù)提取及歸一化 %下載四類語(yǔ)音信號(hào) load data1 c1 load data2 c2 load data3 c3 load data4 c4 %四個(gè)特征信號(hào)矩陣合成一個(gè)矩陣 data(1:500, =c1(1:500, ;data(501:1000, =c2(1:500, ;data(1001:1500, =c3(1:500, ;data(1501:2000, =c4(1:500, ;%從1到2000間隨機(jī)排序 k=rand(1,2000); [m,n]=sort(k); %輸入輸出數(shù)據(jù) input=data(:,2:25); output1 =data(:,1); %把輸出從1維變成4維 for i=1:2000 switch output1(i) case 1 output(i, =[1 0 0 0];case 2 output(i, =[0 1 0 0];case 3 output(i, =[0 0 1 0];case 4 output(i, =[0 0 0 1];end end %隨機(jī)提取1500個(gè)樣本為訓(xùn)練樣本,500個(gè)樣本為預(yù)測(cè)樣本 input_train=input(n(1:1500), ';output_train=output(n(1:1500), ';input_test=input(n(1501:2000), ';output_test=output(n(1501:2000), ';%輸入數(shù)據(jù)歸一化 [inputn,inputps]=mapminmax(input_train); %% 網(wǎng)絡(luò)結(jié)構(gòu)初始化 innum=24; midnum=25; outnum=4; %權(quán)值初始化 w1=rands(midnum,innum); b1=rands(midnum,1); w2=rands(midnum,outnum); b2=rands(outnum,1); w2_1=w2;w2_2=w2_1; w1_1=w1;w1_2=w1_1; b1_1=b1;b1_2=b1_1; b2_1=b2;b2_2=b2_1; %學(xué)習(xí)率 xite=0.1 alfa=0.01; %% 網(wǎng)絡(luò)訓(xùn)練 for ii=1:10 E(ii)=0; for i=1:1:1500 %% 網(wǎng)絡(luò)預(yù)測(cè)輸出 x=inputn(:,i); % 隱含層輸出 for j=1:1:midnum I(j)=inputn(:,i)'*w1(j, '+b1(j);Iout(j)=1/(1+exp(-I(j))); end % 輸出層輸出 yn=w2'*Iout'+b2; %% 權(quán)值閥值修正 %計(jì)算誤差 e=output_train(:,i)-yn; E(ii)=E(ii)+sum(abs(e)); %計(jì)算權(quán)值變化率 dw2=e*Iout; db2=e'; for j=1:1:midnum S=1/(1+exp(-I(j))); FI(j)=S*(1-S); end for k=1:1:innum for j=1:1:midnum dw1(k,j)=FI(j)*x(k)*(e(1)*w2(j,1)+e(2)*w2(j,2)+e(3)*w2(j,3)+e(4)*w2(j,4)); db1(j)=FI(j)*(e(1)*w2(j,1)+e(2)*w2(j,2)+e(3)*w2(j,3)+e(4)*w2(j,4)); end end w1=w1_1+xite*dw1'; b1=b1_1+xite*db1'; w2=w2_1+xite*dw2'; b2=b2_1+xite*db2'; w1_2=w1_1;w1_1=w1; w2_2=w2_1;w2_1=w2; b1_2=b1_1;b1_1=b1; b2_2=b2_1;b2_1=b2; end end %% 語(yǔ)音特征信號(hào)分類 inputn_test=mapminmax('apply',input_test,inputps); for ii=1:1 for i=1:500%1500 %隱含層輸出 for j=1:1:midnum I(j)=inputn_test(:,i)'*w1(j, '+b1(j);Iout(j)=1/(1+exp(-I(j))); end fore(:,i)=w2'*Iout'+b2; end end %% 結(jié)果分析 %根據(jù)網(wǎng)絡(luò)輸出找出數(shù)據(jù)屬于哪類 for i=1:500 output_fore(i)=find(fore(:,i)==max(fore(:,i))); end %BP網(wǎng)絡(luò)預(yù)測(cè)誤差 error=output_fore-output1(n(1501:2000))'; %畫出預(yù)測(cè)語(yǔ)音種類和實(shí)際語(yǔ)音種類的分類圖 figure(1) plot(output_fore,'r') hold on plot(output1(n(1501:2000))','b') legend('預(yù)測(cè)語(yǔ)音類別','實(shí)際語(yǔ)音類別') %畫出誤差圖 figure(2) plot(error) title('BP網(wǎng)絡(luò)分類誤差','fontsize',12) xlabel('語(yǔ)音信號(hào)','fontsize',12) ylabel('分類誤差','fontsize',12) %print -dtiff -r600 1-4 k=zeros(1,4); %找出判斷錯(cuò)誤的分類屬于哪一類 for i=1:500 if error(i)~=0 [b,c]=max(output_test(:,i)); switch c case 1 k(1)=k(1)+1; case 2 k(2)=k(2)+1; case 3 k(3)=k(3)+1; case 4 k(4)=k(4)+1; end end end %找出每類的個(gè)體和 kk=zeros(1,4); for i=1:500 [b,c]=max(output_test(:,i)); switch c case 1 kk(1)=kk(1)+1; case 2 kk(2)=kk(2)+1; case 3 kk(3)=kk(3)+1; case 4 kk(4)=kk(4)+1; end end %正確率 rightridio=(kk-k)./kk web browser http://www.ilovematlab.cn/thread-60056-1-1.html |
source |

找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 0817調(diào)劑 +3 | 沒(méi)有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|