机器学习 Softmax classifier (无隐含层)

程序实现 Softmax classifer, 没有隐含层，

f=wx+b

y=efi∑jefj

%% Softmax classifier

function  Out=Softmax_Classifier(train_x,  train_y, opts)

% setting learning parameters

step_size=opts.step_size;

reg=opts.reg;

batchsize = opts.batchsize;

numepochs = opts.numepochs;

K=opts.class;

%% learn on the training set

% initialize the parameters

D=size(train_x, 2);

W=0.01*randn(D,K);

b=zeros(1,K);

loss(1 : numepochs)=0;

num_examples=size(train_x, 1);

numbatches = num_examples / batchsize;

for epoch=1:numepochs

% %         tic;

% %

% %         sprintf('epoch %d:  \n, ' , epoch)

        kk = randperm(num_examples);

        loss(epoch)=0;

    for bat=1:numbatches

        batch_x = train_x(kk((bat - 1) * batchsize + 1 : bat * batchsize), :);

        batch_y = train_y(kk((bat - 1) * batchsize + 1 : bat * batchsize), :);

        cc=repmat(b,  batchsize  , 1);

        scores=batch_x*W+cc;

        exp_scores=exp(scores);

        dd=repmat(sum(exp_scores, 2), 1, K);

        probs=exp_scores./dd;

        correct_logprobs=-log(sum(probs.*batch_y, 2));

        data_loss=sum(correct_logprobs)/batchsize;

        reg_loss=0.5*reg*sum(sum(W.*W));

        loss(epoch) =loss(epoch)+ data_loss + reg_loss;

        dscores = probs-batch_y;

        dscores=dscores/batchsize;

        dW=batch_x'*dscores;

        db=sum(dscores);

        dW=dW+reg*W;

        W=W-step_size*dW;

        b=b-step_size*db;

    end

    loss(epoch)=loss(epoch)/numbatches;

    if (mod(epoch, 10)==0)

        sprintf('epoch: %d, training loss is  %f:  \n', epoch, loss(epoch))

    end

% %     toc;

end

Out.W=W;

Out.b=b;

Out.loss=loss;

end