测试runpf

rspower/data/convert_m_file.m

-function convert_m_file3(input_file, output_file)
+function convert_m_file(input_file, output_file)
     % input_file: 输入.m文件路径（例如 'case14.m'）
     % input_file: 输入.m文件路径（例如 'case14.m'）
     % output_file：输出.txt文件路径（例如 'case14_mems.txt'）
     % output_file：输出.txt文件路径（例如 'case14_mems.txt'）
     % 将matpower矩阵格式转换为mems支持的脚本矩阵格式
     % 将matpower矩阵格式转换为mems支持的脚本矩阵格式

rspower/examples/compare_results.m

@@ -3,7 +3,7 @@ function compare_results(case_name)
     % case_name: 算例名称，如 'case14', 'case30', 'case57' 等
     % case_name: 算例名称，如 'case14', 'case30', 'case57' 等
 
 
     if nargin < 1
     if nargin < 1
-        case_name = 'case141';  % 默认使用 case14
+        case_name = 'case145';  % 默认使用 case14
     end
     end
 
 
     fprintf('==============开始为算例 %s 执行测试...==============\n', case_name);
     fprintf('==============开始为算例 %s 执行测试...==============\n', case_name);
@@ -20,17 +20,17 @@ function compare_results(case_name)
 
 
     % 比较的field_name和使用matpower计算函数
     % 比较的field_name和使用matpower计算函数
     field_names = {
     field_names = {
-        'ybus'
-        'jac'
-        'sdzip'
-        'sbus'
+%         'ybus'
+%         'jac'
+%         'sdzip'
+%         'sbus'
         'runpf'
         'runpf'
     };
     };
     cal_makers = {
     cal_makers = {
-        @cal_makeybus
-        @cal_makejac
-        @cal_makesdzip
-        @cal_makesbus
+%         @cal_makeybus
+%         @cal_makejac
+%         @cal_makesdzip
+%         @cal_makesbus
         @cal_runpf
         @cal_runpf
     };
     };
 
 
@@ -111,8 +111,9 @@ end
 
 
 function pfv = cal_runpf(case_name)
 function pfv = cal_runpf(case_name)
     mpc = loadcase(case_name);
     mpc = loadcase(case_name);
-    for i = 1:100
-        r = runpf(mpc, mpoption('OUT_ALL',0','VERBOSE',0));
+    mpopt = mpoption('exp.use_legacy_core', 1, 'out.all',0','verbose',0);
+    for i = 1:2000
+        r = runpf(mpc, mpopt);
     end
     end
     bus = r.bus;
     bus = r.bus;
     pfv = bus(:, 8) .* exp(1j * pi/180 * bus(:, 9));
     pfv = bus(:, 8) .* exp(1j * pi/180 * bus(:, 9));
@@ -159,7 +160,7 @@ function compare_matrices(matrix1, matrix2, matrix_name)
     fprintf('  最大相对误差: %.2e\n', relative_error);
     fprintf('  最大相对误差: %.2e\n', relative_error);
 
 
     % 判断是否相等（使用容差）
     % 判断是否相等（使用容差）
-    tolerance = 1e-10;
+    tolerance = 1e-6;
     if max_error < tolerance
     if max_error < tolerance
         fprintf('  ✅ 矩阵相等 (容差: %.0e)\n', tolerance);
         fprintf('  ✅ 矩阵相等 (容差: %.0e)\n', tolerance);
     else
     else