电力电子仿真

2.Runpf程序流程说明：

载入runpf程式，以‘case9’为例子，跟踪程序流程。

第一步：将断点设在第138行，运行程式，跳过初始化阶段。

初始化阶段主要完成：

①初始化，定义变量名。

②初始化，判断输入参数个数，如果沒有输入，则赋予缺省值。

③初始化，特别情況的判断，如DC、虚功Q限制等。

④初始化，读入数据讯息矩阵。

⑤添加零列。

⑥转换为内部分度。

⑦区分出无限大汇流排，发电机和负载。

初始化以后，

①Mpc为已存的case值

②将不同功能部分分开。

第二步：数据读取和存储完成，进入潮流的计算部分。

①找到发电机端，存在矩阵中。

②t0读入当前时钟，用于计算时间。

③判断程序的冗余程度和是否为dc，如果是则做特殊处理，显然case9跳过这些部分。

DC的计算部分也跳过。

④进入判断潮流计算方法的阶段。这里将分别选用(1) 牛顿-拉夫逊法, (2) 高斯-塞德尔法进行跟踪模拟。

(1)Newton Method

①选用Newton Method

②赋初值（flat start）并归一化。

③计算出导纳矩阵和各bus的复功。

④alg=1，进入newton method的子函数程式newtonpf。

⑤算出Jacobian（雅可比）矩阵，运用Newton Method的算子进行迭代，算出v和angle。

第一次迭代i=1；

Newton Method 完成后：

迭代4次达到精度要求，并跳回runpf函数。

(2)Guass methods

①选用Guass Methods

②初始值和归一化计算和Newton Method 相同，不再叙述。

③alg=4，进入Gauss Methods的子函数程式gausspf。

④运用gauss methods的核心递推公式，可以对解进行逼近，直到误差在允许范围之内。核心程式：

初始值：

迭代一次以后：

Gauss Method完成以后：

迭代212次达到精度要求，返回runpf函数。

⑤用不同方法算出了潮流值以后，将原值更新。

⑥根据算出结果与Q原值限制的比较进行调整，以取得实际的最优化解。

第三步、以一个user-friendly的方式将得到的信息输出。

程序结束。