不等式约束

Mat_G=[
    sparse(PD(Loadi));
    sparse(PD(Loadi))-0.001;
    sparse(QD(Loadi));
    sparse(QD(Loadi))-0.001;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>

func_deltF.m

@@ -4,6 +4,7 @@ deltF=[sparse(length(Loadi),1);
     sparse(2*Busnum,1);
     sparse(ones(Busnum,1));
     sparse(ones(length(Loadi),1));
+    sparse(ones(length(Loadi),1))
 ];
 
 end