




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、/p> int Type;bus9;struct BranchDataStructint ID;int FisrtBus; int EndBus;float R,X,B,K;Branch9;float YG99, YB99;int main(array<System:String > args)ifstream ob("F:yangyuanyuandata.txt"); for(int i=0;i<9;i+) ob >> busi.ID >> busi.U >> busi.A >> busi.PG &g
2、t;> busi.QG >> busi.PL >> /-输入数据- busi.QL >> busi.Type;for(int i=0;i<9;i+) ob.close(); ob >> Branchi.ID >> Branchi.FisrtBus >> Branchi.EndBus >> Branchi.R >> Branchi.X >> Branchi.B >> Branchi.K; 4/-节点导纳矩阵- for(int i=0;i<9;i+) for(i
3、nt j=0;j<9;j+) YGij = YBij = 0.0; for(int i=0;i<9;i+) float yG = Branchi.R / ( Branchi.R* Branchi.R + Branchi.X * Branchi.X); float yB =- Branchi.X / ( Branchi.R* Branchi.R + Branchi.X * Branchi.X); int ii = Branchi.FisrtBus - 1; int jj = Branchi.EndBus - 1; YGiiii += yG;YGjjii += -yG; YBiiii
4、+= yB + Branchi.B;YBjjjj += yB + Branchi.B; YBiijj += -yB; YBjjii += -yB; ofstream ob2("F:yangyuanyuanYmatric.txt"); for(int i=0;i<9;i+)for(int j=0;j<9;j+)if(YBij>=0) ob2<<YGij<<"+"<<YBij<<"j"<<"t" elseob2<<YGij<
5、<YBij<<"j"<<"t" ob2<<"n"ob2.close(); /-设电压初值- float Ue9,Uf9; for (int i=0;i<9;i+) Uei=busi.U; Ufi=0.0; 5/-循环- ofstream ob3("f:yangyuanyuandelta.txt");ofstream ob4("f:yangyuanyuanJacobian.txt");ofstream ob6("f:yangyuanyuan
6、deltaEF.txt");ofstream ob7("f:yangyuanyuanU.txt");for (int step=1;step<5;step+) ob3 <<endl<< endl<< "step" << step << endl; ob4 << endl<< endl<< "step" << step<< endl; ob6 << endl<< endl<
7、< "step" << step<< endl;ob7 <<endl<< endl<< "step" << step<< endl; float deltaP9,deltaQ9,deltaUU9; float aii9=0,bii9=0; float P9,Q9;for (int i=0;i<9;i+) Pi=busi.PG-busi.PL; Qi=busi.QG-busi.QL; for (int i=0;i<9;i+) for (int j=0;j&
8、lt;9;j+)aiii+=YGij*Uej-YBij*Ufj;biii+=YGij*Ufj+YBij*Uej;deltaPi=Pi-Uei*aiii-Ufi*biii;deltaQi=Qi-Ufi*aiii+Uei*biii;deltaUUi=busi.U*busi.U-Uei*Uei-Ufi*Ufi;for(int i=0;i<9;i+) ob3 << "deltaP" << i << "=" << deltaPi << endl; for(int i=0;i<9;i+) ob3
9、 << "deltaQ" << i << "=" << deltaQi << endl; ob3 << "deltaUU" << i << "=" << deltaUUi << endl; for(int i=0;i<9;i+) 6/-形成雅可比矩阵- float H99,N99,J99,L99,R99,S99;for(int i=0;i<9;i+)for(int j=0;j<9;j+) Hij=0; Nij=0
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年教育产业并购整合策略与教育投资并购投资前景研究报告001
- 2025年天然气水合物开采技术地质风险防控策略研究预研报告
- 江西省“红色十校”2025-2026学年高三上学期第一次联考语文试题(含答案)
- 信息科个人年终工作总结(19篇)
- 河北省唐山市玉田县孤树镇方官屯中心小学2024-2025学年五年级下学期3月月考(1-2单元)数学试题
- 2025年第28届中小学学校“推普周”活动总结:以语言为桥 共筑中华民族共同体-
- 应力应变试验培训课件
- 2025年房地产行业投资策略分析报告:恒大碧桂园
- 岩石力学许明课件
- 输电六防课件
- 【MOOC】电工电子实验基础-东南大学 中国大学慕课MOOC答案
- 《专利池与公共利益》课件
- 第1课-远古时期的人类活动【同步练习】
- 演唱会招商方案
- 河北信息技术学业水平考试试题集
- 压力容器使用单位安全总监题库
- 创业合伙人五份协议书模板
- 建筑工程消防查验检查表
- 应征公民政治审查表
- FZ∕T 71006-2021 山羊绒针织绒线
- 慢性创面的治疗及护理课件
评论
0/150
提交评论