24小时热门版块排行榜    

查看: 2035  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

betty920

新虫 (小有名气)

[求助] 就是SVM的BPA求解算法!!!!谁能给我把这段伪代码,写成程序。。。

This  appendix   shows  the  pseudo-code  for the  training is shown  below.   The  algorithm is a
model-trust algorithm  based on the Levenberg-Marquardt algorithm [17].

Input parameters:
out = array of SVM outputs  SVM的输出
target  = array of booleans: is ith example  a positive example?(0,1)
priorl  = number  of positive examples prior0  = number  of negative examples
Outputs:
A, B = parameters of sigmoid

A = 0
B = log((prior0+l)/(priorl+l))
hiTarget = (priorl+l)/(priorl+2) yi为+1时
loTarget = l/(prior0+2)  yi为-1时
lambda  = le-3 olderr  = le300
pp = temp array to store current  estimate of probability of examples
set all pp array elements to (priorl+l)/(prior0+priorl+2)
count = 0
for it = l to l00 {
a = 0, b = 0, c         = 0, d = 0, e = 0
// First, compute  Hessian  & gradient of error function
// with respect  to A & B
for i = l to len {
if (target)
        t = hiTarget else
t = loTarget
dl = pp-t
d2 = pp*(l-pp)
a += out*out*d2 b += d2
c += out*d2
d += out*dl e += dl
}
        // If gradient is really  tiny, then stop if (abs(d)  < 1e-9 && abs(e)  < 1e-9)
break
oldA = A
oldB = B
err = 0
        // Loop until goodness of fit increases while (1) {
det = (a+lambda)*(b+lambda)-c*c
if (det == 0) {         // if determinant of Hessian  is zero,
        // increase stabilizer lambda  *= 10
continue
}
A = oldA + ((b+lambda)*d-c*e)/det
B = oldB + ((a+lambda)*e-c*d)/det
// Now, compute  the goodness of fit err = 0;
for i = 1 to len {
p = 1/(1+exp(out*A+B))
pp = p
// At this step, make sure log(0)  returns  -200 err -= t*log(p)+(1-t)*log(1-p)
}
if (err < olderr*(1+1e-7)) {
lambda  *= 0.1 break
}
// error did not decrease: increase stabilizer by factor  of 10
// & try again lambda  *= 10
        if (lambda  >= 1e6) // something is broken.  Give up break
}
diff = err-olderr
scale = 0.5*(err+olderr+1)
if (diff > -1e-3*scale && diff < 1e-7*scale)
count++
else
        count = 0 olderr  = err
if (count  == 3)
break
}


真心求助。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

joshuawangn

铁虫 (正式写手)

支持向量机。。。
7楼2014-09-15 20:33:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

hengxue

铜虫 (小有名气)

留下个记号,待闲时再过来看看,能不能编。
QQ:474933652
2楼2012-05-02 08:44:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

betty920

新虫 (小有名气)

引用回帖:
2楼: Originally posted by hengxue at 2012-05-02 08:44:56:
留下个记号,待闲时再过来看看,能不能编。

能帮忙编写一下吗
3楼2012-05-02 09:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hengxue

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by betty920 at 2012-05-02 09:29:19:
能帮忙编写一下吗

这周比较忙。
我就搞过一个最小二乘法拟合的算法
现在,偏最小二乘法还没有弄呢。
其它的也没有弄过。
QQ:474933652
4楼2012-05-02 09:33:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0805 316求调剂 +4 大雪深藏 2026-03-18 4/200 2026-03-24 17:34 by dick_runner
[考研] 资源与环境 调剂申请(333分) +7 holy J 2026-03-21 7/350 2026-03-24 17:24 by xiaohai104
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 085600材料与化工调剂 +7 A-哆啦Z梦 2026-03-23 12/600 2026-03-23 23:16 by 星空星月
[考研] 北科281学硕材料求调剂 +8 tcxiaoxx 2026-03-20 9/450 2026-03-23 12:16 by tcxiaoxx
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
信息提示
请填处理意见