不等式是这样的

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;
    ];

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

FormG.m

@@ -8,6 +8,6 @@ Mat_G=[
     Volt'-mVolt'-bigM*Vbi-0.03;
     Volt'-mVolt'+bigM*Vbi+0.03;
     Vbi;
-   % PDbi;
+    PDbi;
     ];
 end