| 查看: 714 | 回复: 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 |
回复此楼» 猜你喜欢
这样的filecode谁见过
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
orcimbalance
木虫 (著名写手)
- 应助: 22 (小学生)
- 金币: 2057.1
- 散金: 937
- 红花: 3
- 帖子: 1306
- 在线: 428.9小时
- 虫号: 1401449
- 注册: 2011-09-14
- 性别: GG
- 专业: 计算机体系结构
你把SVM原理都弄清楚了？强！ 
2楼2012-11-13 15:47:28
100