《浙大电院短学期潮流综合实验实验报告(含源程序).doc》由会员分享，提供在线免费全文阅读可下载，此文档格式为doc，更多相关《浙大电院短学期潮流综合实验实验报告(含源程序).doc》文档请在天天文库搜索。

1、基于matlab的牛顿拉夫逊法潮流计算　　　　　 姓名 学号 指导教师 专业与班级 所在学院 目录一、 实验目的和要求 4二、 实验原理与流程图 41.实验原理 41.1 输入初始数据，节点编号，形成节点导纳矩阵 41.2 对于初始值和，求出功率偏差和 41.3 依据公式算出雅可比矩阵 41.4 应用高斯消元法解线性方程组，得到 、 51.5 依据如下公式计算修正量，并形成迭代公式 51.6 计算线路功率分布，PV节点无功功率和平衡节点注入功率 62.程序流程图 6三、 算例及结果 81.算例：IEEE-美国新英格兰10机39节点测试系统 81.1 系统单线图 81.2 系统参数 81.3发电机数据 91.4线路参数 92.计算结果 132.1 节点计算结果 132.2线路计算结果 143. 实验结果分析 15附录： 171.实验源程序 171.1主程序——nr.m 171.2 读取节点。

2、与线路数据——fileopen.m 201.3 对节点重新编号——num.m 211.4 求节点导纳阵——y.m 221.5 计算功率偏差——deltaPQ.m 231.6求雅克比矩阵——form_jac.m 241.7 高斯消元法解线性方程组——gs.m 261.8 节点编号逆回——renum.m 271.9 计算线路潮流及损耗——flow.m 28作业一 301.写一个 MATLAB 函数pifun.m 来计算下列级数 302.一个平面上的椭圆可以表示成下列方程式： 312.1 程序——y.m 31作业二 321.试用牛顿－拉夫逊法求解如下的非线性方程组，要求精度达到10－10以上 321.1 计算结果 321.2 程序 342. 解方程组 352.1 计算结果 362.2 程序 37一、 实验目的和要求电力系统潮流计算目标：取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统。

3、各元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗等。研究目的：分析和评价电网运行的安全经济性，改进、规划。二、 实验原理与流程图1.实验原理程序流程：1.1 输入初始数据，节点编号，形成节点导纳矩阵1)普通线路(k＝0，j不等于0)Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(I,J)=Y(I,J)-Yt;Y(J,I)=Y(I,J)； Y(I,I)=Y(I,I)+Ym; 2)对于对地支路：k＝0，j＝0 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+Yt/K/K; Y(I,J)=Y(I,J)-Yt/K; Y(J,I)=Y(I,J); 3)变压器支路(K>0，非标准变比在j侧)4)变压器支路(k> NRans =第1次迭代结果：ans = Columns 1 through 2 0.1650。

4、 0.8350 0.2209 0.7791 -0.0231 1.0231 Column 3 1.8511 -24.5779 0.5995ans =第2次迭代结果：ans = Columns 1 through 2 0.1215 0.0435 0.6301 0.5908 -0.8709 -0.1522 Column 3 0.4868 -5.0477 0.2525ans =第3次迭代结果：ans = Columns 1 through 2 0.1405 0.9810 0.7892 0.8410 -0.9423 -0.9286 Column 3 0.2268 -1.1974 0.7536ans =第4次迭代结果：ans = Columns 1 through 2 0.6742 0.4663 0.2922 0.4970 -0.6242 -0.3181 Column 3 0.1465 -0.3639。

5、 0.7536ans =第5次迭代结果：ans = Columns 1 through 2 0.6405 0.0337 0.0793 0.2129 -0.1101 -0.5141 Column 3 0.7927 -0.2851 0.4784ans =第6次迭代结果：ans = Columns 1 through 2 0.6294 0.0111 0.3292 0.7501 -0.0676 -0.0425 Column 3 0.9317 -0.5500 0.5862ans 第7次迭代结果：ans = Columns 1 through 2 0.0145 0.6149 0.5860 0.7432 -0.8688 -0.1988 Column 3 0.0435 -0.7129 0.0515ans =第8次迭代结果：ans = Columns 1 through 2 0.0000 0.0145 -0。

6、.0000 0.5860 -0.8299 -0.0389 Column 3 0 0 -0.00021.2 程序1)主程序——NR.m2)编写函数表示题目中的非线性方程组——F.m3)编写一个函数表示题目中的雅可比矩阵——Fd.m2. 解方程组2.1 计算结果>> NRans =第1次迭代结果：ans = 1.0e+002 * 0.6393 -0.6393 1.4234 0.1457 0.8543 -0.2519 -0.0873 0.0873 0.9354ans =第2次迭代结果：ans = 0.7504 0.1804 20.0979 4.5312 0.0425 -0.3128 -4.9539 -1.7776 -1.5003ans =第3次迭代结果：ans = 1.3757 -0.6254 0.5641 5.8448 -0.3136 0.8733 -4.3082 -0.6457 0.4217ans =第4次迭代结果ans = 1.2147 0.1610 0.6068 5.3148 0.5300 0.2008 -4.8956 -0.4126 -0.6697ans =第5次迭代结果：ans = 1.8884 0.3263 0.9779 5.0221 0.2927 0.1521 -4.9142 -0.9814 -0.7833ans =第6次迭代结果：ans = 1.0000 0.8884 0 5.0000 0.0221 0 -4.0000 -0.9142 02.2 程序1)主程序——NR.m2)编写函数表示题目中的非线性方程组——F.m3)编写一个函数表示题目中的雅可比矩阵——Fd.m。