最新电力系统潮流分析报告(五篇)
在经济发展迅速的今天,报告不再是罕见的东西,报告中提到的所有信息应该是准确无误的。那么报告应该怎么制定才合适呢?下面我给大家整理了一些优秀的报告范文,希望能够帮助到大家,我们一起来看一看吧。
电力系统潮流分析报告篇一
#include “stdafx.h” #include #include“complex.h” #include“wanjing.h” #include“gauss.h” using namespace std; int _tmain(int argc, _tchar* argv[]){ int i; //i作为整个程序的循环变量 int n=bus::scanfbusno();//输入节点个数 int l=line::scanflineno();//输入支路个数 if((l&&n)==0){return 0;} //如果找不到两个文件中的任意一个,退出 line *line=new line[l];//动态分配支路结构体 line::scanflinedata(line);//输入支路参数 line::printflinedata(line,l);//输出支路参数 bus *bus=new bus[n];//动态分配结点结构体 for(int i=0;i bus[i].=0; bus[i].=0;} bus::scanfbusdata(bus);//输入节点参数 bus::printfbusdata(bus,n);//输出结点参数 complex **x;x=new complex *[n];for(i=0;i bus::jisuannodednz(x,line,bus,l,n);//计算节点导纳矩阵 bus::printfnodednz(x,n);//输出节点导纳矩阵 int nn=(n-1)*2;double **jacaug;jacaug=new double *[nn];for(i=0;i double *x;x=new double[nn];int count=1; loop: bus::jisuannodei(x,bus,n);//计算节点注入电流 bus::jisuannodescal(x,bus,n);//计算节点功率 bus::jisuannodescal(x,bus,n);//计算节点功率 bus::jisuannodesdelta(bus,n);//计算节点功率差值 bus::printfnodescal(x,bus,n);//输出节点功率差值 int icon=wehcon1(bus,n);//whether converbence看迭代是否结束 if(icon==1){ cout<<“icon=”< bus::jisuanjacaug(jacaug,x,bus,n);//计算雅可比增广矩阵 // bus::printfjacaug(jacaug,n); gauss::gauss_slove(jacaug,x,nn);//解方程组求出电压差值 bus::revisenodev(bus,x,n);//修正节点电压 // bus::printfnodev(bus,n); count++; goto loop;} else { for(i=0;i { int statemp,endtemp; complex aa,bb,cc,dd,b; =0; =-line[i].b; statemp=line[i].start; endtemp=line[i].end; aa=complex::productcomplex(complex::getconj(bus[statemp-1].v), b); bb=complex::subcomplex (complex::getconj(bus[statemp-1].v), complex::getconj(bus[endtemp-1].v)); cc=complex::productcomplex(bb , complex::getconj(line[i].y)); dd=complex::caddc(aa,cc); line[i].stoe=complex::productcomplex(bus[statemp-1].v,dd); aa=complex::productcomplex(complex::getconj(bus[endtemp-1].v), b); bb=complex::subcomplex (complex::getconj(bus[endtemp-1].v), complex::getconj(bus[statemp-1].v)); cc=complex::productcomplex(bb , complex::getconj(line[i].y)); dd=complex::caddc(aa,cc); line[i].etos=complex::productcomplex(bus[endtemp-1].v,dd); } cout<<“icon=”< bus::jisuannodescal(x,bus,n);//计算节点功率 for(i=0;i { bus[i]. = bus[i]. + bus[i].;//发电机功率=注入功率+负荷功率 bus[i].= bus[i].+ bus[i].; bus[i].v=complex::rec2polar(bus[i].v); } cout<<“====节点电压===============发电机发出功率======”< for(i=0;i { cout<<“节点”<<(i+1)<<'t'; complex::printfcomplex(bus[i].v); coutt(bus[i].); coutt(bus[i].); cout< } cout<<“======线路传输功率==========”< for(i=0;i { int statemp,endtemp; statemp=line[i].start; endtemp=line[i].end; cout< complex::printfcomplex(complex::comdivrea(line[i].stoe,0.01)); complex::printfcomplex(complex::comdivrea(line[i].etos,0.01)); cout< } } return 0;} #include “stdafx.h” #include #include“complex.h” #include“wanjing.h” #include“gauss.h” using namespace std; int _tmain(int argc, _tchar* argv[]){ int i; //i作为整个程序的循环变量 int n=bus::scanfbusno();//输入节点个数 int l=line::scanflineno();//输入支路个数 if((l&&n)==0){return 0;} //如果找不到两个文件中的任意一个,退出 line *line=new line[l];//动态分配支路结构体 line::scanflinedata(line);//输入支路参数 line::printflinedata(line,l);//输出支路参数 bus *bus=new bus[n];//动态分配结点结构体 for(int i=0;i bus[i].=0; bus[i].=0;} bus::scanfbusdata(bus);//输入节点参数 bus::printfbusdata(bus,n);//输出结点参数 complex **x;x=new complex *[n];for(i=0;i bus::jisuannodednz(x,line,bus,l,n);//计算节点导纳矩阵 bus::printfnodednz(x,n);//输出节点导纳矩阵 int nn=(n-1)*2;double **jacaug;jacaug=new double *[nn];for(i=0;i double *x;x=new double[nn];int count=1; loop: bus::jisuannodei(x,bus,n);//计算节点注入电流 bus::jisuannodescal(x,bus,n);//计算节点功率 bus::jisuannodescal(x,bus,n);//计算节点功率 bus::jisuannodesdelta(bus,n);//计算节点功率差值 bus::printfnodescal(x,bus,n);//输出节点功率差值 int icon=wehcon1(bus,n);//whether converbence看迭代是否结束 if(icon==1){ cout<<“icon=”< bus::jisuanjacaug(jacaug,x,bus,n);//计算雅可比增广矩阵 // bus::printfjacaug(jacaug,n); gauss::gauss_slove(jacaug,x,nn);//解方程组求出电压差值 bus::revisenodev(bus,x,n);//修正节点电压 // bus::printfnodev(bus,n); count++; goto loop;} else { for(i=0;i { int statemp,endtemp; complex aa,bb,cc,dd,b; =0; =-line[i].b; statemp=line[i].start; endtemp=line[i].end; aa=complex::productcomplex(complex::getconj(bus[statemp-1].v), b); bb=complex::subcomplex (complex::getconj(bus[statemp-1].v), complex::getconj(bus[endtemp-1].v)); cc=complex::productcomplex(bb , complex::getconj(line[i].y)); dd=complex::caddc(aa,cc); line[i].stoe=complex::productcomplex(bus[statemp-1].v,dd); aa=complex::productcomplex(complex::getconj(bus[endtemp-1].v), b); bb=complex::subcomplex (complex::getconj(bus[endtemp-1].v), complex::getconj(bus[statemp-1].v)); cc=complex::productcomplex(bb , complex::getconj(line[i].y)); dd=complex::caddc(aa,cc); line[i].etos=complex::productcomplex(bus[endtemp-1].v,dd); } cout<<“icon=”< bus::jisuannodescal(x,bus,n);//计算节点功率 for(i=0;i { bus[i]. = bus[i]. + bus[i].;//发电机功率=注入功率+负荷功率 bus[i].= bus[i].+ bus[i].; bus[i].v=complex::rec2polar(bus[i].v); } cout<<“====节点电压===============发电机发出功率======”< for(i=0;i { cout<<“节点”<<(i+1)<<'t'; complex::printfcomplex(bus[i].v); coutt(bus[i].); coutt(bus[i].); cout< } cout<<“======线路传输功率==========”< for(i=0;i { int statemp,endtemp; statemp=line[i].start; endtemp=line[i].end; cout< complex::printfcomplex(complex::comdivrea(line[i].stoe,0.01)); complex::printfcomplex(complex::comdivrea(line[i].etos,0.01)); cout< } } return 0;} #include
2024六年级毕业评语简短(精选十五篇)
06-11
六年级毕业致同学们的一封信范文2篇
06-11
初中毕业典礼校长发言稿2024(优秀10篇)
06-11
2024年六年级毕业活动策划书(二篇)
06-11
2024初中毕业典礼感谢老师发言稿范文8篇
06-11
2024年大学全民国家安全教育日活动方案(精选9篇)
04-07
党课讲稿2024专题党课(精选5篇)
04-02
2024脱贫攻坚与乡村振兴有效衔接工作规划精选5篇
03-18
关于创三甲医院动员大会讲话稿范文7篇
03-20
2024红领巾奖章四星章事迹材料15篇
06-07