0001
0002
0003
0004
0005
0006
0007 function [A_mining_d b_mining_d] = deviation_positivity_constraint...
0008 (Blocks120, numOfPeriods)
0009
0010
0011 [numBlocks n] = size([Blocks120.EBV]');
0012
0013 A_mining_x = zeros(numOfPeriods, numBlocks * numOfPeriods);
0014 A_mining_d = eye(numOfPeriods, numOfPeriods)*-1;
0015
0016 A_mining_x = sparse(A_mining_x);
0017 A_mining_d = sparse(A_mining_d);
0018
0019 b_mining_d = sparse(numOfPeriods, 1);
0020
0021 A_mining_d = [A_mining_x, A_mining_x, A_mining_x, A_mining_x...
0022 A_mining_d, A_mining_d, A_mining_d, A_mining_d];
0023 end
0024