资源描述：
《浙大电院短学期潮流综合实验实验报告(含源程序).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系统单线图
2、81.2系统参数81.3发电机数据91.4线路参数92.计算结果132.1节点计算结果132.2线路计算结果143.实验结果分析15附录：171.实验源程序171.1主程序——nr.m171.2读取节点与线路数据——fileopen.m201.3对节点重新编号——num.m211.4求节点导纳阵——y.m221.5计算功率偏差——deltaPQ.m231.6求雅克比矩阵——form_jac.m241.7高斯消元法解线性方程组——gs.m261.8节点编号逆回——renum.m271.9计算线路潮流及损耗
3、——flow.m28作业一301.写一个MATLAB函数pifun.m来计算下列级数302.一个平面上的椭圆可以表示成下列方程式：312.1程序——y.m31作业二321.试用牛顿－拉夫逊法求解如下的非线性方程组，要求精度达到10－10以上321.1计算结果321.2程序342.解方程组352.1计算结果362.2程序37一、实验目的和要求电力系统潮流计算目标：取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统各元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗等。研究
4、目的：分析和评价电网运行的安全经济性，改进、规划。二、实验原理与流程图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＝0Y(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,
5、I)=Y(I,J);3）变压器支路（K>0，非标准变比在j侧）4）变压器支路（k<0，非标准变比在i侧）Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+K*K*Yt;Y(I,J)=Y(I,J)-K*Yt;Y(J,I)=Y(I,J);1.2对于初始值和，求出功率偏差和1.3依据公式算出雅可比矩阵设一个系统，含有n个节点，其中m个PQ节点，1个平衡节点，n-m-1个PV节点，则：其中：1.4应用高斯消元法解线性方程组，得到、1.5依据如下公式计算修正量，并形成迭代公式1.6计算线路功率分
6、布，PV节点无功功率和平衡节点注入功率2.程序流程图一、算例及结果1.算例：IEEE-美国新英格兰10机39节点测试系统1.1系统单线图1.2系统参数1）系统容量基准值为100MVA。2)负荷数据见表D-1表D-1负荷数据节点号有功（MW）无功（Mvar）节点号有功（MW）无功（Mvar）BUS-3BUS-4BUS-7BUS-8BUS-12BUS-15BUS-16BUS-18BUS-20322.0500.0233.8522.08.5320.0329.0158.0680.02.4184.084.0176.
7、088.0153.032.330.0103.0BUS-21BUS-23BUS-24BUS-25BUS-26BUS-27BUS-28BUS-29BUS-39274.0247.5308.0224.0139.0281.0206.0283.51104.0115.084.6-92.247.217.075.527.626.91.3发电机数据表D-2发电机数据发电机节点号有功（MW）无功（Mvar）电压（p.u.）PVPQPQPVPVPQPVPVPVBUS-30BUS-32BUS-33BUS-34BUS-35BUS-
8、36BUS-37BUS-38BUS-39250.00650.00632.00508.00650.00560.00540.00830.001000.00175.90103.3596.881.047501.012301.049301.027801.026501.03000平衡节点BUS-310.0(电压幅角)1.0（幅值）1.4线路参数表D-3线路参数序号节点I节点JR(p.u.)X(p.u.)B/2(p.u.)LN1LN2LN3LN4LN5L
