getborders.m

来自「线性M3网络及其示意图」· M 代码 · 共 39 行

function [A,b]=getBorders(C1,C2)
%[A,b]getBorders get border functions of C1 and C2 in form of output [A,b]


    
    
m=size(C1,1);
n=size(C2,1);

A=[];
b=[];

for i=1:m
    for j=1:n
        
       %Z=(C1(i,1)-C2(j,1))*C1(i,1)+(C1(i,2)-C2(j,2))*C1(i,2)+0.5*(norm(C2(j,:))^2-norm(C1(i,:))^2); 
       % Z for the denotation of decision function ! 
       
       %if Z>0
       %    Z=1;
       %else
       %    Z=-1;
       %end
       Z=1;% no mean
       
       temp=[C1(i,1)-C2(j,1),C1(i,2)-C2(j,2)];% this problem is 2D 
       %temp=temp/abs(C1(i,2)-C2(j,2));  % normalize y notice divide by zero
       A=[A;Z*temp];
       
       temp=-0.5*(norm(C2(j,:))^2-norm(C1(i,:))^2);
       %temp=temp/abs(C1(i,2)-C2(j,2));  % normalize y  notice divide by zero
       b=[b;Z*temp];
        
        
    end
end