测试sbus、runpf通过

rspower/examples/compare_results.m

@@ -24,7 +24,7 @@ function compare_results(case_name)
         'jac'
         'sdzip'
         'sbus'
-%         'runpf'
+        'runpf'
     };
     cal_makers = {
         @cal_makeybus
@@ -78,17 +78,22 @@ end
 
 function sdzip = cal_makesdzip(case_name)
     mpc = loadcase(case_name);
-    sdzip = full(makeSdzip(mpc));
+    sd = makeSdzip(mpc.baseMVA, mpc.bus);
+    sdzip = [sd.z, sd.i, sd.p];
 end
 
 function sbus = cal_makesbus(case_name)
     mpc = loadcase(case_name);
-    sbus = full(makeSbus(mpc));
+    baseMVA = mpc.baseMVA;
+    bus = mpc.bus;
+    gen = mpc.gen;
+    sbus = makeSbus(baseMVA, bus, gen);
 end
 
-function pf = cal_runpf(case_name)
+function pfv = cal_runpf(case_name)
     mpc = loadcase(case_name);
-    pf = runpf(mpc);
+    bus = runpf(mpc, mpoption('OUT_ALL',0','VERBOSE',0)).bus;
+    pfv = bus(:, 8) .* exp(1j * pi/180 * bus(:, 9));
 end
 
 

rspower/examples/test_make_jac.txt

@@ -7,6 +7,8 @@
 #include ../lib/make_jac.txt
 #include ../lib/bustypes.txt
 
+// 可以同时测试dSbus_dV的计算
+
 jac = make_jac(baseMVA, bus, branch, gen, 1);
 println(jac);
 return jac;

rspower/examples/test_make_sbus.txt

@@ -7,19 +7,19 @@
 r = make_sbus(baseMVA, bus, gen);
 // |a - b| <= max(|a|, |b|) * max_relative(default=1e-8)
 assert_relative_eq(r, [
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)],
-    [c(1,0)]
+    [c(2.324,-0.169)],
+    [c(0.183,0.297)],
+    [c(-0.942,0.044)],
+    [c(-0.478,0.039)],
+    [c(-0.076,-0.016)],
+    [c(-0.112,0.047)],
+    [c(0,0)],
+    [c(0,0.174)],
+    [c(-0.295,-0.166)],
+    [c(-0.090,-0.058)],
+    [c(-0.035,-0.018)],
+    [c(-0.061,-0.016)],
+    [c(-0.135,-0.058)],
+    [c(-0.149,-0.050)]
 ]);
 return r;
\ No newline at end of file