1、模型的建立

詳細(xì)代碼：

<span style="font-size:18px;">%% 該代碼為基于BP_Adaboost的強預(yù)測器預(yù)測%% 清空環(huán)境變量 clc clear%% 下載數(shù)據(jù) load data1 input output%% 權(quán)重初始化 k=rand(1,2000); [m,n]=sort(k);%訓(xùn)練樣本 input_train=input(n(1:1900),:)'; output_train=output(n(1:1900),:)';%測試樣本 input_test=input(n(1901:2000),:)'; output_test=output(n(1901:2000),:)';%樣本權(quán)重 [mm,nn]=size(input_train); D(1,:)=ones(1,nn)/nn;%訓(xùn)練樣本歸一化 [inputn,inputps]=mapminmax(input_train); [outputn,outputps]=mapminmax(output_train);K=10; for i=1:K%弱預(yù)測器訓(xùn)練net=newff(inputn,outputn,5);net.trainParam.epochs=20;net.trainParam.lr=0.1;net=train(net,inputn,outputn);%弱預(yù)測器預(yù)測an1=sim(net,inputn);BPoutput=mapminmax('reverse',an1,outputps);%預(yù)測誤差erroryc(i,:)=output_train-BPoutput;%測試數(shù)據(jù)預(yù)測inputn1=mapminmax('apply',input_test,inputps);an2=sim(net,inputn1);test_simu(i,:)=mapminmax('reverse',an2,outputps);%調(diào)整D值Error(i)=0;for j=1:nnif abs(erroryc(i,j))>0.2 %較大誤差Error(i)=Error(i)+D(i,j);D(i+1,j)=D(i,j)*1.1;elseD(i+1,j)=D(i,j);endend%計算弱預(yù)測器權(quán)重at(i)=0.5/exp(abs(Error(i)));%D值歸一化D(i+1,:)=D(i+1,:)/sum(D(i+1,:));end%% 強預(yù)測器預(yù)測 at=at/sum(at);%% 結(jié)果統(tǒng)計 %強分離器效果 output=at*test_simu; error=output_test-output; plot(abs(error),'-*') hold on for i=1:8 error1(i,:)=test_simu(i,:)-output; end plot(mean(abs(error1)),'-or')title('強預(yù)測器預(yù)測誤差絕對值','fontsize',12) xlabel('預(yù)測樣本','fontsize',12) ylabel('誤差絕對值','fontsize',12) legend('強預(yù)測器預(yù)測','弱預(yù)測器預(yù)測')</span>
分類1 -1 ?代碼：

<span style="font-size:18px;">%% 該代碼為基于BP-Adaboost的強分類器分類%% 清空環(huán)境變量 clc clear%% 下載數(shù)據(jù) load data input_train output_train input_test output_test%% 權(quán)重初始化 [mm,nn]=size(input_train); D(1,:)=ones(1,nn)/nn;%% 弱分類器分類 K=10; for i=1:K%訓(xùn)練樣本歸一化[inputn,inputps]=mapminmax(input_train);[outputn,outputps]=mapminmax(output_train);error(i)=0;%BP神經(jīng)網(wǎng)絡(luò)構(gòu)建net=newff(inputn,outputn,6);net.trainParam.epochs=5;net.trainParam.lr=0.1;net.trainParam.goal=0.00004;%BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練net=train(net,inputn,outputn);%訓(xùn)練數(shù)據(jù)預(yù)測an1=sim(net,inputn);test_simu1(i,:)=mapminmax('reverse',an1,outputps);%測試數(shù)據(jù)預(yù)測inputn_test =mapminmax('apply',input_test,inputps);an=sim(net,inputn_test);test_simu(i,:)=mapminmax('reverse',an,outputps);%統(tǒng)計輸出效果kk1=find(test_simu1(i,:)>0);kk2=find(test_simu1(i,:)<0);aa(kk1)=1;aa(kk2)=-1;%統(tǒng)計錯誤樣本數(shù)for j=1:nnif aa(j)~=output_train(j);error(i)=error(i)+D(i,j);endend%弱分類器i權(quán)重at(i)=0.5*log((1-error(i))/error(i));%更新D值for j=1:nnD(i+1,j)=D(i,j)*exp(-at(i)*aa(j)*test_simu1(i,j));end%D值歸一化Dsum=sum(D(i+1,:));D(i+1,:)=D(i+1,:)/Dsum;end%% 強分類器分類結(jié)果 output=sign(at*test_simu);%% 分類結(jié)果統(tǒng)計 %統(tǒng)計強分類器每類分類錯誤個數(shù) kkk1=0; kkk2=0; for j=1:350if output(j)==1if output(j)~=output_test(j)kkk1=kkk1+1;endendif output(j)==-1if output(j)~=output_test(j)kkk2=kkk2+1;endend endkkk1 kkk2 disp('第一類分類錯誤 第二類分類錯誤 總錯誤'); % 窗口顯示 disp([kkk1 kkk2 kkk1+kkk2]);plot(output) hold on plot(output_test,'g')%統(tǒng)計弱分離器效果 for i=1:Kerror1(i)=0;kk1=find(test_simu(i,:)>0);kk2=find(test_simu(i,:)<0);aa(kk1)=1;aa(kk2)=-1;for j=1:350if aa(j)~=output_test(j);error1(i)=error1(i)+1;endend end disp('統(tǒng)計弱分類器分類效果'); error1disp('強分類器分類誤差率') (kkk1+kkk2)/350disp('弱分類器分類誤差率') (sum(error1)/(K*350)) </span>

總結(jié)

以上是生活随笔為你收集整理的基于BP弱分类器用Adaboost的强分类器的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。