欢迎来到天天文库
浏览记录
ID:79082145
大小:1.23 MB
页数:27页
时间:2022-02-09
《电力系统》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、《潮流计算》实验轮机工程学院2013.04姚玉斌实验目的和内容目的:初步掌握潮流计算程序设计内容:第5周潮流计算:数据输入第7周潮流计算:极坐标牛顿法迭代主程序1建立目录在D盘上建立一个目录(以姓名的拼音命名)在C盘已建目录:电力系统自动化实验/PowerFlowPowerFlow.cpp源文件PowerFlow.h头文件PowerFlowLib.lib库文件book5.txt输入数据文件PowerFlow.dsw工程文件双击,打开该工程2使用VisaulC++6.0开始/程序/MicrosoftViualStudio6.0/Micros
2、oftViualC++6.0建立单个文件建立C++源文件选择目录输入文件名文件/新建2使用VisaulC++6.0编写、编译、链接、执行程序在本窗口输入文件内容编译、链接、执行3数据结构1)数据结构structCLine//支路数据{inti;//首端母线号(k侧为负号)intj;//末端母线号(k侧为负号)doubler;//电阻(在1侧)doublex;//电抗(在1侧)doubleyk;//B/2(或变比)};3数据结构1)数据结构structCBus//节点数据{inti;//母线号intType;//节点类型doublepg;/
3、/有功出力doubleqg;//无功出力doublepl;//有功负荷(负值)doubleql;//无功负荷(负值)doublev0;//电压给定值(PQ节点给定值为1.0)};3数据结构1)数据结构structCCp//电容器(或电抗器){inti;//所在母线号doubleb;//电纳值(电容器为正值,电抗器为负值)};structCYij//导纳{doubleg;//实部doubleb;//虚部};3数据结构1)数据结构structCVp//电压{doublev;//幅值doubleva;//相角};structCVxy//直角坐标
4、电压{doubleve;//实部doublevf;//虚部};3数据结构2)常量constintNS=1000;//节点数据数组大小constintZS=2000;//支路数据数组大小constintNS2=NS+NS;//节点数据数组大小2constintVA=-2;//节点类型:平衡节点constintPV=-1;//节点类型:PV节点constintPQ=0;//节点类型:PQ节点3数据结构3)数组CLineLine[ZS];CBusBus[NS];CCpCp[NS];CVpVp[NS];CVxyVxy[NS];CYijYij[NS
5、][NS];doubleJacb[NS2][NS2];//雅可比矩阵doublew[NS2];//节点功率数组doubleb[NS2];//功率不平衡量数组、解向量3数据结构4)变量intn;//节点数intn2;//intzls;//支路数intncp;//电容器个数intdsd;//最大不平衡量的节点号intniter;//迭代次数intmaxiter;//最大迭代次数doubledsm;//最大不平衡量doubleeps;//收敛精度FILEfp1,fp2;3数据结构5)库函数voiddataio(void);//从文件中读取数据v
6、oidvini();//平启动电压初值voidPr_nod(void);//输出节点计算结果voidbranch(void);//生成支路计算结果voidy2(void);//形成导纳矩阵voidsolution(void);//牛顿法解方程voidnewtonp(void);//极坐标牛顿法voidjcbp(void);//雅可比矩阵voidbbhl_p(void);//不平衡量及平衡节点功率、PV节点无功4源文件1)包含文件4源文件2)主函数4源文件3)数据输入函数4源文件3)数据输入函数4源文件3)数据输入函数4源文件3)数据输入函
7、数4源文件3)数据输入函数4源文件3)数据输入函数4源文件3)数据输入函数5极坐标牛顿法迭代主程序框图电压初始化计算残差P、Q修正量、V,修正、V计算雅可比矩阵P、Q<形成导纳矩阵开始结束原始数据输入结果输出YN5比较工具Windiff选择比较文件红色第1个文件,黄色第2个文件6使用VisaulC++6.0在工程中添加静态库工程/设置选择Win32Debug添加PowerFlowLib.libTheEnd
此文档下载收益归作者所有