2017-02-14 23 views
0

所以我想确定我的SVM分类器的alpha,betagamma值,我训练了我的SVM分类器并确定了alpha和支持向量。我的支持向量具有每行3个变量,下面是从我的支持向量短短十行:添加符号变量后很长的数字?

0.0904235536887480 -0.269325475875919 -0.678528701392414 
-0.321039098061280 -0.507618180664821 -1.42365662798284 
-0.0737761304021185 -0.269641311369441 -0.647521877863172 
0.00105779420640393 -0.311226557946309 -0.667506146498475 
0.0913098589312967 -0.289462325547514 -0.391261050348894 
0.00622693949845773 -0.166248587146820 -0.149546793127464 
-0.292302915842567 -0.564676268888150 -1.60153093563523 
0.112997393643248 -0.310512134534035 -0.725281274142312 
-0.135361511770186 -0.456321702624641 -1.26973221898260 
-0.173160731078767 -0.434439033384469 -1.22687774941370 

,同样下面是我的阿尔法短短十行:

-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 

所以我基本上做的是确定如下系数:

c = (-1)*[0.0904*x1 + (-0.2693)*x2 + (-0.6785)*x3] + (-1)*[(-0.321)*x1 + (-0.5076)*x2 + (-1.4236)*x3]..... 

等等,直到我的尺寸为alpha

所以我编码在Matlab

syms x1, syms x2, syms x3; 

alpha = SVMStruct.Alpha; 
svm_vec = SVMStruct.SupportVectors; 

for i = 1:size(alpha,1) 
    c(i) = alpha(i)*(svm_vec(i,1)*x1 + svm_vec(i,2)*x2 + svm_vec(i,3)*x3); 
end 

sum_it = sum(c); 

下列但是这创造了一个非常奇怪的输出:

(107845064549358722206080751595348329973204613074833920445585562521882937008164658045489239834546021458299139*x1)/50534761550197893278639420198779799540396107395587434771118149413836407509953624874438129483687080755200 + (95720990302914087945142311872326568914380675701489099929103269189530321664249312169660240242394455632803627*x2)/134178504805697854567421908803656709124500009291732154392279224305703564767807900528680550698065697177600 + (90626366614084720573448362168042659133754200934323766866906741825007634289583081638045482944881264585156183*x3)/125521827076297992982426946945356276277758073208394596044390242092432367040852552107475353878835652198400 

这是预期的输出?为什么我将它作为分数?

从输出看,我是否正确实现了我的方程?

回答

1

您的计算输入是符号变量,因此使用这些变量的任何计算都会产生一个确切值,即也是的一个符号变量。如果你想要的结果的近似浮点表示,你要投的总和为double

double(sum_it) 

为了虽然做到这一点,你需要为x1指定实际值,x2,和

double(subs(sum(c), [x1, x2, x3], [1, 2, 3])) 

话虽这么说,我觉得你其实什么是解方程的系统。您可以使用\操作人员的答复做到这一点

variables = svm_vec \ alpha; 
+0

谢谢,我不能使用双,因为我不知道X1,X2和X3的价值,我只对那些变量的系数感兴趣。我不知道'''''实现了什么,我读了它,并没有完全理解它的目的。话虽如此,我的执行是否正确?如果是这样,我可以用''x1''手工分割这些数字来得到这个系数,对吧? – StuckInPhD

+2

@StuckInPhD:使用['vpa'](https://www.mathworks.com/help/symbolic/vpa.html)。不知道为什么你需要这个应用程序的符号或变量精度。 – horchler

+0

@StuckInPhD你不是想解决'x1','x2'和'x3'吗? – Suever