0001 function selMat = selMatrix(selVec)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020 nVar = length(selVec);
0021 if (sum(selVec == -1) == 0)
0022
0023 nSel = sum(selVec);
0024 isel = [1:nSel];
0025 jsel = find(selVec);
0026 selMat = sparse(isel,jsel,ones(nSel,1),nSel,nVar);
0027
0028 else
0029
0030 selFwInd = find(selVec == 1);
0031 selMat = sparse(length(selFwInd),nVar);
0032 for i = 1:length(selFwInd)
0033 selFwID = selFwInd(i);
0034 if (selVec(selFwID+1) == -1)
0035 selMat(i,selFwID) = 1;
0036 selMat(i,selFwID+1) = -1;
0037 else
0038 selMat(i,selFwID) = 1;
0039 end
0040 end
0041 selMat = sparse(selMat);
0042
0043 end