我正在使用fmincon求解器和sqp算法解决约束优化问题。
我所拥有的线性等式和不等式矩阵相当稀疏。我很想知道如何以稀疏形式定义以下Aeq和Beq矩阵以加快计算速度？
最好的问候
cons1=[eye(10),-eye(10),zeros(10,5*19),-eye(10),zeros(10,280)];
cons2=[zeros(19,20),eye(19),-eye(19),zeros(19,19*2),-eye(19),zeros(19,290)];
cons3=[zeros(19,96),diag(tan_fid),zeros(19,10),-eye(19),zeros(19,261)];
cons4=[zeros(19,58),-eye(19),-eye(19),eye(19),zeros(19,290)];
% #### linear_eq=Aeq*X'-Beq; ####
Aeq=[cons1;cons2;cons3;cons4];
Beq=[-pga ; -pda ;zeros(19,1); zeros(19,1)];