| 查看: 2083 | 回复: 6 | |||
[求助]
就是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 } 真心求助。 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求助翻译德语一段实验操作
已经有2人回复
- 高手请进，帮忙看下我翻译的一段话（中译英），有没有问题 谢谢
已经有2人回复
- 实验室条件下，要想得到几个不同的淬火冷速，可以采用哪些手段呢
已经有23人回复
- 华东理工大学调取初试试卷，需要哪些程序
已经有3人回复
- 【课件】MATLAB 7.0 高级语言程序设计基础及应用
已经有184人回复
- FFT算法的一个问题请教
已经有18人回复
- 【代发帖】结构力学求解
已经有10人回复
- Matlab如何实现这样的计算：赋值——计算——判断结果——重复
已经有10人回复
- 遗传算法，，，什么意思？
已经有15人回复
- 想投EI收录的国内期刊～红外目标检测算法方面的
已经有7人回复
- 求解电化学知识
已经有10人回复
- 跟踪算法：如何判断是否有遮挡？
已经有18人回复
- 如何将取向矩阵中归一化的{hkl}转化成互质化的{HKL}，求算法！！！！
已经有10人回复
- 求助四波混频耦合波方程数值求解
已经有5人回复
- 求解什么是Biomass!!!!
已经有5人回复
2楼2012-05-02 08:44:56
3楼2012-05-02 09:29:19
4楼2012-05-02 09:33:20
gzl9901
铁杆木虫 (文坛精英)
very good
- 应助: 3 (幼儿园)
- 金币: 6279.5
- 红花: 28
- 沙发: 2
- 帖子: 11263
- 在线: 1769.2小时
- 虫号: 1294464
- 注册: 2011-05-13
- 专业: 计算机软件
5楼2012-05-02 17:58:44
6楼2012-05-03 08:37:35
joshuawangn
铁虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 196.4
- 散金: 600
- 红花: 11
- 帖子: 649
- 在线: 87.5小时
- 虫号: 2943254
- 注册: 2014-01-22
- 性别: GG
- 专业: 计算机软件
7楼2014-09-15 20:33:24
10