Matlab牛拉法代码.docVIP

牛顿-拉夫逊法潮流计算的Matlab实现 (2009-09-26 13:23:26)转载 标签： 电力系统 潮流计算 牛顿-拉夫逊 牛顿拉夫逊 matlab newton-raphson 写了近一个礼拜的程序，藏着可惜，拿出来给需要的人 转载请注明出处 有不少朋友反映程序报错： ??? Output argument bus (and maybe others) not assigned during call to E:\MAT\best\OpDF_.mOpDF_. 是因为输入数据的格式问题 输入数据中是包含bus,line的，并非只有两个矩阵，如： bus = [ 1 1.00 0.00 -0.30 -0.18 1; 2 1.00 0.00 -0.55 -0.13 1; 3 1.10 0.00 0.50 0.00 2; 4 1.05 0.00 0.00 0.00 3]; line = [ 1 2 0.10 0.40 0.0 0.01528 0; 4 2 0.08 0.40 0.0 0.01413 0; 1 4 0.12 0.50 0.0 0.0192 0; 3 1 0.00 0.3 0.0 0.0 -1.1]; 所有蓝色文字都为数据文件中的内容 ---------------------分---------------------割---------------------线--------------------- 【注意】 一、此程序只适用于求解节点电压以极坐标形式表示的潮流方程，没有考虑节点优化编号 二、程序在Matlab 6.5上测试通过，应该适用于目前其所有后续版本 三、程序主要通过文件方式输入输出（同时程序也返回结果向量），对输入文件格式有严格要求，具体如下： 1.输入文件可以直接在Matlab中新建m文件编写，也内容可以以文本方式编写，但最后必须改后缀名成为“.m”文件。文件名第一个字符必须是字母，后面可以跟字母、数字和下划线的任何组合，但不能和已有文件和函数冲突，不能含中文。 2.输入文件内容格式： 内容中包括bus（节点数据）和line（线路数据）两个矩阵 其中bus的格式为： bus = [节点编号 节点电压 节点相角（弧度制） 有功注入 无功注入 节点类型]; line的格式为： line = [节点i 节点j 线路电阻 电抗 电导 电纳 变压器变比（普通线路为零）]; 四、程序分主程序和各子程序，需全部放入Matlab当前目录下调用，计算开始后会在当前目录下生成文件Result.m保存运算结果。 ---------------------分---------------------割---------------------线--------------------- 以下为程序代码（红色部分，“%”后为注释）： 主程序 PowerFlow_NR.m function [bus_res,S_res] = PowerFlow_NR_2 % 牛顿-拉夫逊法解潮流方程的主程序 %%%%%%%%%%%%%%%%%%%%% by longdinhohe %%%%%%%%%%%%%%%%%%%%% /longdinhohe %%%%%%%%%%%%%%%%%%%%% [bus,line] = OpDF_; % 打开数据文件的子程序，返回bus（节点数据）和line（线路数据）回主程序 [nb,mb]=size(bus); [nl,ml]=size(line); % 计算bus和line矩阵的行数和列数 [bus,line,nPQ,nPV,nodenum] = Num_(bus,line); % 对节点重新排序的子程序 Y = y_(bus,line); % 计算节点导纳矩阵的子程序 myf = fopen(Result.m,w); fprintf(myf,-------