| 查看: 674 | 回复: 1 | ||
stonechao1金虫 (著名写手)
|
[求助]
关于SVM拉哥朗日乘子对应的的支持向量问题?用cvx_box设计SVM
|
|
load FV600 % FV600是一个数据集，X1对应特征向量，y对应标签，是个两类问题 %下面一段是打乱样本集 [N,m1]=size(X1); [qs,qi]=sort(rand(N,1)); X1=X1(qi,  ;y=y(qi, ;%十倍交叉验证 indices = crossvalind('Kfold',N,10); for F = 1:10 test = (indices == F); train = ~test; Test1=X1(test, ;Train1=X1(train, ;ytrain=y(train);ytest=y(test); [n,m1]=size(Train1); D1=Train1*Train1'; e=ones(n,1); for c=1:10 C=2^c; O=zeros(n,1); cvx_begin variable lambda1(n); maximize(lambda1'*e-0.5*(lambda1.*ytrain)'*D1*(lambda1.*ytrain)); lambda1>=O; lambda1<=C*e; lambda1'*ytrain==0; cvx_end %下面一段是求偏置b，取拉哥朗日乘子最大的10个支持向量来求平均。 [max,ind]=sort(lambda1); for j=1:10 bb(j)=(ytrain(ind(n-j+1))-(lambda1.*ytrain)'*Train1*Train1(ind(n-j+1), ');end b{c}=mean(bb); ypredict1{F,c}=Test1*Train1'*(lambda1.*ytrain)+b{c}; acc1(F,c)=accuracy(ytest,sign(ypredict1{F,c}+10e-10)); aveacc1=mean(acc1); end end |
回复此楼» 猜你喜欢
292求调剂
已经有6人回复
-
化学工程085602 305分求调剂
已经有7人回复
-
332求调剂
已经有5人回复
-
材料求调剂
已经有4人回复
-
303求调剂
已经有5人回复
-
334分 一志愿武理-080500 材料求调剂
已经有4人回复
-
求调剂
已经有3人回复
-
299求调剂
已经有3人回复
-
299求调剂
已经有8人回复
-
316求调剂
已经有8人回复
orcimbalance
木虫 (著名写手)
- 应助: 22 (小学生)
- 金币: 2057.1
- 散金: 937
- 红花: 3
- 帖子: 1306
- 在线: 428.9小时
- 虫号: 1401449
- 注册: 2011-09-14
- 性别: GG
- 专业: 计算机体系结构
你把SVM原理都弄清楚了？强！ 
2楼2012-11-13 15:47:28
6