Adaboost算法和MATLAB實現


 

一、AdaBoost簡介

      Boosting, 也稱為增強學習或提升法,是一種重要的集成學習技術, 能夠將預測精度僅比隨機猜度略高的弱學習器增強為預測精度高的強學習器,這在直接構造強學習器非常困難的情況下,為學習算法的設計提供了一種有效的新思路和新方法。其中最為成功應用的是,Yoav Freund和Robert Schapire在1995年提出的AdaBoost算法。
      AdaBoost是英文"Adaptive Boosting"(自適應增強)的縮寫,它的自適應在於:前一個基本分類器被錯誤分類的樣本的權值會增大,而正確分類的樣本的權值會減小,並再次用來訓練下一個基本分類器。同時,在每一輪迭代中,加入一個新的弱分類器,直到達到某個預定的足夠小的錯誤率或達到預先指定的最大迭代次數才確定最終的強分類器。
Adaboost算法可以簡述為三個步驟:
 (1)首先,是初始化訓練數據的權值分布D1。假設有N個訓練樣本數據,則每一個訓練樣本最開始時,都被賦予相同的權值:w1=1/N。
 (2)然后,訓練弱分類器hi。具體訓練過程中是:如果某個訓練樣本點,被弱分類器hi准確地分類,那么在構造下一個訓練集中,它對應的權值要減小;相反,如果某個訓練樣本點被錯誤分類,那么它的權值就應該增大。權值更新過的樣本集被用於訓練下一個分類器,整個訓練過程如此迭代地進行下去。
 (3)最后,將各個訓練得到的弱分類器組合成一個強分類器。各個弱分類器的訓練過程結束后,加大分類誤差率小的弱分類器的權重,使其在最終的分類函數中起着較大的決定作用,而降低分類誤差率大的弱分類器的權重,使其在最終的分類函數中起着較小的決定作用。
  換而言之,誤差率低的弱分類器在最終分類器中占的權重較大,否則較小。

二、AdaBoost算法過程

    給定訓練數據集:,其中用於表示訓練樣本的類別標簽,i=1,...,N。Adaboost的目的就是從訓練數據中學習一系列弱分類器或基本分類器,然后將這些弱分類器組合成一個強分類器。

相關符號定義:

Adaboost的算法流程如下:

 

 

 

相關說明:

 

綜合上面的推導,可得樣本分錯與分對時,其權值更新的公式為:

三、AdaBoost實例講解

 例:給定如圖所示的訓練樣本,弱分類器采用平行於坐標軸的直線,用Adaboost算法的實現強分類過程。

數據分析:

   將這10個樣本作為訓練數據,根據 X 和Y 的對應關系,可把這10個數據分為兩類,圖中用“+”表示類別1,用“O”表示類別-1。本例使用水平或者垂直的直線作為分類器,圖中已經給出了三個弱分類器,即:

初始化:

   首先需要初始化訓練樣本數據的權值分布,每一個訓練樣本最開始時都被賦予相同的權值:wi=1/N,這樣訓練樣本集的初始權值分布D1(i):

   令每個權值w1i = 1/N = 0.1,其中,N = 10,i = 1,2, ..., 10,然后分別對於t= 1,2,3, ...等值進行迭代(t表示迭代次數,表示第t輪),下表已經給出訓練樣本的權值分布情況:

第1次迭代t=1:

  初試的權值分布D1為1/N(10個數據,每個數據的權值皆初始化為0.1),

D1=[0.1,  0.1, 0.1, 0.1, 0.1, 0.1,0.1, 0.1, 0.1, 0.1]

  在權值分布D1的情況下,取已知的三個弱分類器h1h2h3中誤差率最小的分類器作為第1個基本分類器H1(x)(三個弱分類器的誤差率都是0.3,那就取第1個吧)

在分類器H1(x)=h1情況下,樣本點“5 7 8”被錯分,因此基本分類器H1(x)的誤差率為:

可見,被誤分類樣本的權值之和影響誤差率e,誤差率e影響基本分類器在最終分類器中所占的權重α

然后,更新訓練樣本數據的權值分布,用於下一輪迭代,對於正確分類的訓練樣本“1 2 3 4 6 9 10”(共7個)的權值更新為:

這樣,第1輪迭代后,最后得到各個樣本數據新的權值分布:

D2=[1/14,1/14,1/14,1/14,1/6,1/14,1/6,1/6,1/14,1/14]

  由於樣本數據“5 7 8”被H1(x)分錯了,所以它們的權值由之前的0.1增大到1/6;反之,其它數據皆被分正確,所以它們的權值皆由之前的0.1減小到1/14,下表給出了權值分布的變換情況:

可得分類函數:f1(x)= α1H1(x) = 0.4236H1(x)。此時,組合一個基本分類器sign(f1(x))作為強分類器在訓練數據集上有3個誤分類點(即5 7 8),此時強分類器的訓練錯誤為:0.3

顯然,H2(x)把樣本“3 4 6”分錯了,根據D2可知它們的權值為D2(3)=1/14,D2(4)=1/14, D2(6)=1/14,所以H2(x)在訓練數據集上的誤差率:

這樣,第2輪迭代后,最后得到各個樣本數據新的權值分布:

D3=[1/22,1/22,1/6,1/6,7/66,1/6,7/66,7/66,1/22,1/22]

  下表給出了權值分布的變換情況:

可得分類函數:f2(x)=0.4236H1(x) + 0.6496H2(x)。此時,組合兩個基本分類器sign(f2(x))作為強分類器在訓練數據集上有3個誤分類點(即3 4 6),此時強分類器的訓練錯誤為:0.3

因此,取當前最小的分類器h3作為第3個基本分類器H3(x):

這樣,第3輪迭代后,得到各個樣本數據新的權值分布為:

D4=[1/6,1/6,11/114,11/114,7/114,11/114,7/114,7/114,1/6,1/38]

  下表給出了權值分布的變換情況:

可得分類函數:f3(x)=0.4236H1(x) + 0.6496H2(x)+0.9229H3(x)。此時,組合三個基本分類器sign(f3(x))作為強分類器,在訓練數據集上有0個誤分類點。至此,整個訓練過程結束。

  整合所有分類器,可得最終的強分類器為:

這個強分類器Hfinal對訓練樣本的錯誤率為0!

    本例Matlab代碼,如下:

 

    先建立Matlab函數文件,定義h1,h2和h3三個弱分類器

function kind = wcH1( X,TH )  
%h1弱分類器  
X1=X(1);  
X2=X(2);  
if X1<TH  
    kind=1;  
else  
    kind=-1;  
end  
end  
function kind = wcH2( X,TH )  
%h2弱分類器  
X1=X(1);  
X2=X(2);  
if X1<TH  
    kind=1;  
else  
    kind=-1;  
end  
end  
function kind = wcH3( X,TH )  
%h3弱分類器  
X1=X(1);  
X2=X(2);  
if X2<TH  
    kind=-1;  
else  
    kind=1;  
end  
end  

主程序Matlab代碼:

clc,clear all;  
%% 訓練樣本數據  
 xData=[1 5;2 2;3 1;4 6;6 8;6 5;7 9;8 7;9 8;10 2] %樣本數據點,對應編號為1,2,...10  
 Y=[1 1 -1 -1 1 -1 1 1 -1 -1]';                             %對應的樣本類別,用1和-1表示  
 xNum=1:10;                                                      %編號  
 format rat  
 %% 繪制樣本分布圖  
 L1=find(Y==1);  
 x=xData(L1,1);y=xData(L1,2);  
plot(x,y,'b+','LineWidth',3,'MarkerSize',12);  
hold on;  
 L2=find(Y==-1);  
x=xData(L2,1);y=xData(L2,2);  
plot(x,y,'ro','LineWidth',3,'MarkerSize',12);  
xlabel('X1');ylabel('X2');axis([0 10 0 10])  
%% ***********************************初試過程************************************  
H1=zeros(10,1);H2=H1;H3=H1  
for i=1:10  
    X=xData(i,:);  
    H1(i) = wcH1( X,2.5 );%弱分類器h1  
    H2(i) = wcH2( X,8.5 );%弱分類器h2  
    H3(i) = wcH3( X,6.5 );%弱分類器h3  
end  
errDataH1=find(H1~=Y);%找到被h1錯分的樣本點的序號  
errDataH2=find(H2~=Y);%找到被h2錯分的樣本點的序號  
errDataH3=find(H3~=Y);%找到被h3錯分的樣本點的序號  
accDataH1=find(H1==Y);%找到被h1正確分的樣本點的序號  
accDataH2=find(H2==Y);%找到被h2正確分的樣本點的序號  
accDataH3=find(H3==Y);%找到被h3正確分的樣本點的序號  
errDataAll=[errDataH1,errDataH2,errDataH3];  
accDataAll=[accDataH1,accDataH2,accDataH3];  
  
N=10;  
D1=zeros(10,1)+1/N       % 初始化權值分布  
%% ***********************************第一次迭代***********************************  
err1=sum(D1(errDataH1,:));%所有被錯分類的樣本點的權值之和即為誤差率  
err2=sum(D1(errDataH2,:));%所有被錯分類的樣本點的權值之和即為誤差率  
err3=sum(D1(errDataH3,:));%所有被錯分類的樣本點的權值之和即為誤差率  
errAll=[err1,err2,err3];  
[minErr,minIndex]=min(errAll);  
%根據誤差率e1計算H1的系數:  
a1=0.5*log((1-minErr)/minErr)  
minErrData=errDataAll(:,minIndex);  
minAccData=accDataAll(:,minIndex);  
D2=D1;  
for i=minAccData'  
    D2(i)=D2(i)/(2*(1-minErr));  
end  
for i=minErrData'  
     D2(i)=D2(i)/(2*minErr);  
end  
D2  
%分類函數  
f1=a1.*H1;  
kindFinal=sign(f1)%此時強分類器的分類結果  
  
%% ***********************************第二次迭代***********************************  
err1=sum(D2(errDataH1,:));%所有被錯分類的樣本點的權值之和即為誤差率  
err2=sum(D2(errDataH2,:));%所有被錯分類的樣本點的權值之和即為誤差率  
err3=sum(D2(errDataH3,:));%所有被錯分類的樣本點的權值之和即為誤差率  
errAll=[err1,err2,err3];  
[minErr,minIndex]=min(errAll);  
% 根據誤差率e2計算H2的系數:  
a2=0.5*log((1-minErr)/minErr)  
minErrData=errDataAll(:,minIndex);  
minAccData=accDataAll(:,minIndex);  
D3=D2;  
for i=minAccData'  
    D3(i)=D3(i)/(2*(1-minErr));  
end  
for i=minErrData'  
     D3(i)=D3(i)/(2*minErr);  
end  
D3  
% 分類函數  
f2=a1.*H1+a2*H2;  
kindFinal=sign(f2)%此時強分類器的分類結果  
  
%% ***********************************第三次迭代***********************************  
err1=sum(D3(errDataH1,:));%所有被錯分類的樣本點的權值之和即為誤差率  
err2=sum(D3(errDataH2,:));%所有被錯分類的樣本點的權值之和即為誤差率  
err3=sum(D3(errDataH3,:));%所有被錯分類的樣本點的權值之和即為誤差率  
errAll=[err1,err2,err3];  
[minErr,minIndex]=min(errAll);  
% 根據誤差率e3計算G3的系數:  
a3=0.5*log((1-minErr)/minErr)  
minErrData=errDataAll(:,minIndex);  
minAccData=accDataAll(:,minIndex);  
D4=D3;  
for i=minAccData'  
    D4(i)=D4(i)/(2*(1-minErr));  
end  
for i=minErrData'  
     D4(i)=D4(i)/(2*minErr);  
end  
D4  
% 分類函數  
f3=a1.*H1+a2*H2+a3*H3;  
kindFinal=sign(f3)%此時強分類器的分類結果  
%%  

四、AdaBoost的優點和缺點

優點

     (1)Adaboost提供一種框架,在框架內可以使用各種方法構建子分類器。可以使用簡單的弱分類器,不用對特征進行篩選,也不存在過擬合的現象。

     (2)Adaboost算法不需要弱分類器的先驗知識,最后得到的強分類器的分類精度依賴於所有弱分類器。無論是應用於人造數據還是真實數據,Adaboost都能顯著的提高學習精度。

     (3)Adaboost算法不需要預先知道弱分類器的錯誤率上限,且最后得到的強分類器的分類精度依賴於所有弱分類器的分類精度,可以深挖分類器的能力。Adaboost可以根據弱分類器的反饋,自適應地調整假定的錯誤率,執行的效率高。

     (4)Adaboost對同一個訓練樣本集訓練不同的弱分類器,按照一定的方法把這些弱分類器集合起來,構造一個分類能力很強的強分類器,即“三個臭皮匠賽過一個諸葛亮”。

缺點:

     在Adaboost訓練過程中,Adaboost會使得難於分類樣本的權值呈指數增長,訓練將會過於偏向這類困難的樣本,導致Adaboost算法易受噪聲干擾。此外,Adaboost依賴於弱分類器,而弱分類器的訓練時間往往很長。

 轉載自https://blog.csdn.net/guyuealian/article/details/70995333

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM