lmse.m

模式识别较全的讲解

function LMSE
clc;
clear;
% x1 = [0 0 0;1 0 0;1 0 1;1 1 0]';
% x2 = [0 0 1;0 1 1;0 1 0;1 1 1]';
% x1 = [0 0;0 1]';
% x2 = [1 0;1 1]';
x1 = [0 1;0 -1]';
x2 = [1 0;-1 0]';
X = [x1';-x2'];
X = [X [ones(size(x1,2),1);-ones(size(x2,2),1)]];
X_sharp = inv(X'*X)*X';
b = ones(size(x1,2)+size(x2,2),1);
while 1
    w = X_sharp*b;
    e = X*w-b
    if e==0 | norm(e)<1e-4
        w
        break;
    elseif e<0
        fprintf('该模式不可分\n');
        break;
    else
%         b = b+0.5*2*e
        for i=1:length(b)
            if e(i)>0
                b(i) = b(i)+2*e(i);
            end
        end
        continue

    end
end