Rev | Line | |
---|
[37] | 1 | function [z,F] = binaryindex(n,m,Fin) |
---|
| 2 | |
---|
| 3 | % Number of bits needed |
---|
| 4 | k=ceil(log2(m)); |
---|
| 5 | |
---|
| 6 | w=2^k; |
---|
| 7 | b=dec2bin(0:w-1); |
---|
| 8 | for i = 1:m |
---|
| 9 | pattern{i} = str2num(strrep(strrep(b(i,:),'1',' 1 '),'0',' 0 ')); |
---|
| 10 | end |
---|
| 11 | |
---|
| 12 | z=sdpvar(n,m,'full'); |
---|
| 13 | delta = binvar(n,k,'full'); |
---|
| 14 | F=set([]); |
---|
| 15 | M=m; |
---|
| 16 | for iz=1:n |
---|
| 17 | for jz=1:m |
---|
| 18 | t=0; |
---|
| 19 | for i = 1:length(pattern{jz}) |
---|
| 20 | if pattern{jz}(i) |
---|
| 21 | F=F+set( -M*delta(iz,i) <= z(iz,jz) <= M*delta(iz,i)); |
---|
| 22 | t=t+1-delta(iz,i); |
---|
| 23 | else |
---|
| 24 | F=F+set( -M*(1-delta(iz,i)) <=z(iz,jz)<= M*(1-delta(iz,i))); |
---|
| 25 | t=t+delta(iz,i); |
---|
| 26 | end |
---|
| 27 | end |
---|
| 28 | F=F+set(-M*t <= z(iz,jz)-1 <= M*t); |
---|
| 29 | end |
---|
| 30 | end |
---|
| 31 | F=F+set(m-1 >= delta*2.^((k-1:-1:0)')); |
---|
| 32 | |
---|
| 33 | if nargin == 3 |
---|
| 34 | F= F+Fin; |
---|
| 35 | end |
---|
Note: See
TracBrowser
for help on using the repository browser.