1.一种MATLAB与BPA联合仿真计算含IPFC电力系统潮流的方法，其特征在于，包括如下步骤：S1：BPA潮流计算程序对DAT文件进行初始化潮流计算，生成潮流计算PFO文件，并令迭代次数n＝1；S2：在MATLAB中设定IPFC主控线路有功功率控制目标P ijref 与无功功率控制目标Q ijref 、辅控线路有功功率控制目标P ikref ，并由数据接口将其传给MATLAB中的IPFC计算程序；S3：数据接口读取PFO文件，获得第n次迭代时的输入信息 并将其传给MATLAB中的IPFC计算程序；其中 为第n次迭代时i节点电压幅值、 为第n次迭代时i节点电压相角、 为第n次迭代时j节点电压幅值、 为第n次迭代时j节点电压相角、 为第n次迭代时k节点电压幅值、 为第n次迭代时k节点电压相角、 为第n次迭代时线路ij有功功率、 为第n次迭代时线路ij无功功率、 为第n次迭代时线路ik有功功率；S4：IPFC计算程序将控制目标与输入信息代入式(1)～(4)进行计算，求得 其中， 为第n次迭代时IPFC注入节点j的有功功率、 为第n次迭代时IPFC注入节点j的无功功率、 为第n次迭代时IPFC注入节点k的有功功率、 为第n次迭代时IPFC主控线路等效电压源在p轴的分量、 为第n次迭代时IPFC主控线路等效电压源在q轴的分量、 为第n次迭代时IPFC辅控线路等效电压源在p轴的分量、 为第n次迭代时IPFC辅控线路等效电压源在q轴的分量， 为第n次迭代时线路ij不含IPFC时的自然有功潮流、 为第n次迭代时线路ij不含IPFC时的自然无功潮流、 为第n次迭代时线路ik不含IPFC时的自然有功潮流；X seij 为IPFC主控线路串联耦合变压器的电抗、X seik 为IPFC辅控线路串联耦合变压器的电抗； S5：根据上一步计算出的 进一步计算，得到第n次迭代时IPFC注入节点i的有功功率 和无功功率 注入节点k的无功功率 S6：将输出信息 传给数据接口，并通过数据接口将输出信息写入DAT文件中；S7：判断是否满足下式的收敛条件：若不满足则返回步骤S3，并令n＝n+1；若满足则进入下一步；式中，ε为收敛精度；S8：将求得的V seijp 、V seijq 、V seikp 、V seikq 代入下式，计算出IPFC控制参数V seij 、θ seij 、V seik 、θ seik ，计算结束；其中，V seijp 为IPFC主控线路等效电压源在p轴的分量、V seijq 为IPFC主控线路等效电压源在q轴的分量、V seikp 为IPFC辅控线路等效电压源在p轴的分量、V seikq 为IPFC辅控线路等效电压源在q轴的分量V seij 为IPFC主控线路等效电压源幅值、θ seij 为IPFC主控线路等效电压源相角、V seik 为IPFC辅控线路等效电压源幅值、θ seik 为辅控线路等效电压源相角。
2.根据权利要求1所述的一种MATLAB与BPA联合仿真计算含IPFC电力系统潮流的方法，其特征在于：步骤S5中，所述 由下面三个方程得到：
