支持向量機(jī)非線性回歸通用matlab源碼(編輯修改稿)
2025-07-23 19:42 本頁面
【文章內(nèi)容簡(jiǎn)介】 OPT=optimset。=39。off39。=39。off39。%%%%%整理輸出回歸方程的系數(shù)Alpha1=(Gamma(1:l,1))39。Alpha2=(Gamma((l+1):end,1))39。Alpha=Alpha1Alpha2。Flag=2*ones(1,l)。%%%%%支持向量的分類Err=。for i=1:lAA=Alpha1(i)。BB=Alpha2(i)。if (abs(AA0)=Err)amp。amp。(abs(BB0)=Err)Flag(i)=0。%非支持向量endif (AAErr)amp。amp。(AA Flag(i)=2。%標(biāo)準(zhǔn)支持向量endif (abs(AA0)=Err)amp。amp。(BBErr)amp。amp。(BB Flag(i)=2。%標(biāo)準(zhǔn)支持向量endif (abs(AAC)=Err)amp。amp。(abs(BB0)=Err)Flag(i)=1。%邊界支持向量endif (abs(AA0)=Err)amp。amp。(abs(BBC)=Err)Flag(i)=1。%邊界支持向量endend%%%%%計(jì)算回歸方程中的常數(shù)項(xiàng)BB=0。counter=0。for i=1:lAA=Alpha1(i)。BB=Alpha2(i)。if (AAErr)amp。amp。(AA %計(jì)算支持向量加權(quán)值SUM=0。for j=1:lif Flag(j)0switch TKFcase 1SUM=SUM+Alpha(j)*sum(X(:,j).*X(:,i))。case 2SUM=SUM+Alpha(j)*(sum(X(:,j).*X(:,i))+c)^p。case 3SUM=SUM+Alpha(j)*exp((norm(X(:,j)X(:,i)))^2/(2*sigma^2))。case 4SUM=SUM+Alpha(j)*exp(norm(X(:,j)X(:,i))/(2*sigma^2))。case 5SUM=SUM+Alpha(j)*1/(1+exp(v*sum(X(:,j).*X(:,i))+c))。otherwiseSUM=SUM+Alpha(j)*exp((sum((X(:,j)X(:,i)).^2)/(2*sigma^2)))。endendendb=Y(i
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1