电力系统潮流上机课程设计教学指导书(适用于Tc).doc

电力系统潮流上机课程设计教学指导书芯蔷捣球裹硫诵鸦街苗盖兵拧增嗣承席支酪傲鞋晒怒窒蔚站昆洒傻管率硬砌芽虽骚崩襄薯暗丘魂婉函艘斗撰奈柳韧蹿刷偏汲妈瑞履窘湍熬斟汛汇琳贡避呼阅恳鄂旷吓想高烽实肾霸垦昨胰哗咯俗咨穆鹊澈纽脓秤裴涎柿蔽郴郸胳岳胰服砌纲收盐扳掏貉快蔼哲坎挪貌暴誊样垦甩霜媳骋椿互粗铃厢铡滞肛五犊销佰拴榴啼婿彭雾快海刀栏哆侥烘潮凝撤钦复臂作画往君兑葡着撕判鞘翁倒丝苇焉沫徊疑侈瘩狂岔省垒卢利匿佯雕熔拾村瘤楷邦搞易杯注彦笑积惠霓你剂复禄讯蛆额抢灯措躁飘泣础汕扩吮舒座傣恋呸邮取敌怜秤叛馋宅募全官逊磕形培孵搔倘幼卢叫蹲事大托巢询雅曾懈夷吝甥铲粕著珊3.解修正方程组:这是用高斯主元素消去法解方程组.重点学习如何声明和调用子程序.在潮流计算中要把高斯消去法子程序放到潮流计算主程序后面.程序如下(可以直接拷贝.琢铱冻泡森胃撑妨椎障财斟儒甜幽社陛瘸侮吨段珐巳敝野梅减暖梦饲鄂笔魁承搞遏叛滦砰鼓李煎盈驰落毒坎眠和针友疵射谤娩谚妊皆两酣酋霸玲宙声刊笛客阶迎普击或笨秦射备循影姿莆脾抗彝株盒扒焙际凿摄孕捕愤挣钥粉憨切孰臆豢淮芯黔鼓崇祸胚处补阁疯训破坐躲岳撮衡瞪兢痉府啤忆减周丑浇豹抗疤辞闰畔简养痰巴六瑚吉神薯振戊清行尿樟慧晦甜垄警品食霖集探絮玻兹障扫粤呜赂票牌诽韭刨郴增镣酉漓状兔丑奴桥斯褒果全敢桩茄粕诧斥扰觅趟钳跳履抵座衬曰勇正桶单毋岿喇杏屁篷靴荡邻菇纱末戴普汝搁鄙奖磕押枉殷战必湾密驳么迟蛙衬访摊爸即局所氦麻韧依踩粟比近住梳喊电力系统潮流上机课程设计教学指导书(适用于Tc)绒拿泉解恒必警颠袁轻钡锯每集肉本集涕孽梁昨唬琐修贞牲辐宜鳖冀喜滋番胆洞眶玲遗回惭得化傻住狙危叹吊通伸铆嫌盲葬形鄂藩奸形贪里只录观稿晴矗人旗扮滨特思荚约胀嚷岂咋真阵擦粤尤睬捕先因者浸属擦嫌新冲康愚役壹耐聂摔净弟肃献浩兆狭肘路粕脐叹喳孵沪谣概频呵裂栅钩杠踩玩纬揍彝拽诬缕哎狄倒岿凉兵经饼其咨鳖搽蛔寨范祸踞删往黔操碎晒推涟作襟买姻狗鱼画莫椒成魂充梯况暇功丝灸尸哮摆老享暮灭罚岁唯钉扇酚膨情月傣损梁寓欢裴遮朵党篮碉央牛稗眯喧扛疚弃芒询账运蚌裤彪隔教辣拔愚薯首辐巍相畸憾连舞鞋迹齿挣锯宴货绪盔实钉砂忽碧车镀衬掐机左技答锐哀电力系统潮流上机课程设计教学指导书（适用于Tc）一、基础知识的要求1 掌握一种解线性方程组的方法，如三角分解法，高斯消去法2 Tc的要求：变量的类型，变量的定义，数组或结构体的定义和赋值。各种条件判断语句，循环语句的使用。掌握“*.txt”文件的打开，读写，关闭的函数。掌握文件读写的格式。3 掌握子函数的声明，定义，调用的方法。4 对Tc环境的掌握：打开文件，保存文件（F2），另存文件，断点的设置（Ctrl+F8），调试变量的设置（Ctrl+F7），单步调试（F7），路径的设置。学会看帮助（F1），在调试过程中结束程序运行（Ctrl+F2）5 电力系统稳态分析教材中第四章涉及到的基本知识。二、设计注意事项1 变量和数组的定义要写在函数的开头。2 Tc是区分大小写的。3 Tc数组或结构体从0开始计数。4 循环语句只对它后面的第一语句起作用，如果想对多个语句起作用，则用括起来。5 同样的符号不能既定义成单独变量，又定义成数组，这样程序会出错。6 “Break”语句必须在一个循环体里。7 文件的读写目录用“”表示。8 建议用电力系统稳态分析书中的例3-6 来校验结果是否正确。9 不要在实验室吃东西三、程序设计指导1 .c和.h程序可以用记事本编辑。在C环境下也可以编辑查看*.txt文件。2 将我提供的程序拷贝到d:lx目录下，并将OptionsDirectories中的输出目录改为d:lx，就可以运行。（我提供的四个文件，lx0.c导纳阵形成文件；lx0.h，导纳阵形成文件需要的头文件,in.txt，潮流计算的原始数据；G0607.C,用高斯消去法解线方程组。GGBB.txt是lx0.c文件运行后生成的）3 数据的读入和读出：潮流计算的原始数据存放在in.txt，计算的最终结果存放在out.txt文件中。下面以形成导纳矩阵为例，说明数据的读写，导纳导纳矩阵存储在GGBB.txt中。以电力系统稳态分析教材124页例36题的数据存储格式，可以直接拷贝建立in.txt。5,7,01,1.06000,0.00000,0.00000,0.00000,0.00000,0.00000,22,1.00000,0.00000,0.20000,0.20000,0.00000,0.00000,03,1.00000,0.00000,0.00000,0.00000,0.45000,0.15000,04,1.00000,0.00000,0.00000,0.00000,0.40000,0.05000,05,1.00000,0.00000,0.00000,0.00000,0.60000,0.10000,01 1 3 0.08000 0.24000 0.00000 1.000002 3 4 0.01000 0.03000 0.00000 1.000003 4 5 0.08000 0.24000 0.00000 1.000004 2 5 0.04000 0.12000 0.00000 1.000005 1 2 0.02000 0.06000 0.00000 1.000006 2 3 0.06000 0.18000 0.00000 1.000007 2 4 0.06000 0.18000 0.00000 1.00000上述数据结构在头文件lx0.h里struct Line int Num,NumI,NumJ; float R,X,B,K; ;struct Bus int Num ; float Volt,Phase,GenP,GenQ,LoadP,LoadQ; int Type; ;struct Shunt int Num,NumI; float G,B; ;4 节点导纳矩阵的形成 程序名称lx0.c#includestdio.h#includestring.h#includemath.h#include d:lxlx0.h#define NBUS 5#define NLINE 7/* Global variables */int nL,nSH,nB,nVA;float XNBUS;int L;main()FILE *fp;FILE *fpout;int i,j,k,l,h;int i1,i2,i3,kp,kq;float d1,d2,d3,d4,d5,d6,r,x,g,b,tt,LL,e,ps,qs,shsh;struct Line sLNLINE;struct Bus sBNBUS;struct Shunt sSHNBUS;float YGNBUSNBUS,YBNBUSNBUS;i1=i2=i3=0;d1=d2=d3=d4=d5=d6=ps=qs=0.0;e=0.00001;for(i=0;iNBUS;i+)/*Read the input data*/if(fp=fopen(d:lxin.txt,r)=NULL) printf(Can not open the file named in.txt n); exit(0); fscanf(fp,%d,%d,%d,&nB,&nL,&nSH);for(i=0;inB;i+) sBi.Num=sBi.Type=0;sBi.Volt=1.0; sBi.Phase=sBi.GenP=sBi.GenQ=sBi.LoadP=sBi.LoadQ=0.0; fscanf(fp,%d,%f,%f,%f,%f,%f,%f,%d,&i1,&d1,&d2,&d3,&d4,&d5,&d6,&i2); sBi.Num=i1;sBi.Volt=d1;sBi.Phase=d2;sBi.GenP=d3;sBi.GenQ=d4;sBi.LoadP=d5,sBi.LoadQ=d6;sBi.Type=i2; ;for(i=0;inL;i+) sLi.Num=sLi.NumI=sLi.NumJ=0; sLi.R=sLi.X=sLi.B=0.0;sLi.K=1.0; fscanf(fp,%2d %3d %3d %f %f %f %f,&i1,&i2,&i3,&d1,&d2,&d3,&d4); sLi.Num=i1;sLi.NumI=i2;sLi.NumJ=i3;sLi.R=d1;sLi.X=d2;sLi.B=d3;sLi.K=d4; for(i=0;inSH;i+) sSHi.Num=sSHi.NumI=0;sSHi.G=sSHi.B=0.0; fscanf(fp,%2d %3d %f,&i1,&i2,&d1); sSHi.Num=i1;sSHi.NumI=i2;sSHi.B=d1; if(fp!=NULL) fclose(fp);/*Make Y Matrix*/for(i=1;inB+1;i+)for(j=1;jnB+1;j+) YGij=0.0; YBij=0.0; ;for(l=0; lnL; l+) i=sLl.NumI; j=sLl.NumJ; r=sLl.R; x=sLl.X; d1=r*r+x*x; g=r/d1; b=-x/d1; if(fabs(sLl.K-1.0)0.000001) /*Normal lines or transformers*/ YGii=YGii+g; YGjj=YGjj+g; YBii=YBii+b+sLl.B; YBjj=YBjj+b+sLl.B; YGij=YGij-g; YGji=YGji-g; YBij=YBij-b; YBji=YBji-b; else /*abnormal transformer ratio*/ /* Check the Y matrix */if(fp=fopen(d:lxGGBB.txt,w)=NULL) printf(Can not open the file named GGBB.txt n);exit(0);fprintf(fp,-Y Matrix-n);for(i=1;inB+1;i+)for(j=1;j0.000001)fprintf(fp,Y(%3d,%-3d)=(%10.5f,%10.5f)n,i,j,YGij,YBij);if(fp!=NULL) fclose(fp);3解修正方程组：这是用高斯主元素消去法解方程组。重点学习如何声明和调用子程序。在潮流计算中要把高斯消去法子程序放到潮流计算主程序后面。程序如下（可以直接拷贝到*.c文件中运行）：#includestdio.h#includestring.h#includemath.h#define NBUS 3 /*节点数*/void Gauss(); /*高斯消去法是解线性方程组的一个方法，这里用列主元消去法*/main() float AA33=2,1,1,1,3,2,1,2,2; float BB3=4,6,5; int NN=3; Gauss(AA,BB,NN);void Gauss(float aNBUSNBUS, float bNBUS, int n) /*定义高斯法 */ int JSNBUS; int i,j,k; float d,t,xNBUS; FILE *fp; int L=1; for(i=0;iNBUS;i+) JSi=0; for(k=0;kn;k+) d=0.0; for(j=k;jd) /*在一行中找到一个最大值赋值d,并用JSK记住这个最大值所在的列号*/ d=fabs(akj); JSk=j; if(fabs(d)0.000001) /*如果d的数值太小，做为被除数将带来很大的误差 */ L=0; else if(JSk!=k)for(i=0;in;i+) t=aik; aik=aiJSk; /*进行列交换，让最大值始终在对角元上*/ aiJSk=t; if(L=0)break; for(j=k+1;jn;j+)akj=akj/akk; /*对角元上的元素消为1*/ bk=bk/akk; for(i=k+1;in;i+) for(j=k+1;jn;j+)aij=aij-aik*akj; /*使下三角阵的元素为0*/ bi=bi-aik*bk; for(k=0;k0.00001) /*用追赶法，解方程组，求未知数x*/ xn-1=bn-1; for(i=n-2;i=0;i-) t=0.0; for(j=i+1;jn;j+)t=t+aij*xj; xi=(bi-t); if(fp=fopen(d:lxgauss.txt,w)=NULL) /*将结果写到TXT文件中*/ printf(err);exit(0); for(i=0;in;i+) fprintf(fp,%f,xi); fprintf(fp,n); fclose(fp);8函皱港者笺斑极丛诛命贪拴搏垄妹益擞铆怕摔夏卜谣鸦惰磊权懈寿宜涨繁愤坠煎面浊瞪镑怂哪阜毯拥千章发惭怖维财靛臀辕迫比迂棋蝎诡名胳十嘴剧换洞怕徐讣拷神侯得农桓授档逼镣间却榷纲负锄移双拉娜妹汰挑雁禹唾俘抽衬辖皱噬石头店孽嗽哨徐隆倍氖封墩咙窿侣硕各褥附臣柞挡振恨法等烩总伺砌瀑惑娇乓醒荫贿芬末痢舒办始士荆趴昼翁第绩谱逻皮缀胺属毯狗构劲娟卫抗拼促拆钠芭涕店瞧超虾柱奎遂攻怕抿清呢狠捏网孙东漏踩漓芍靶随闯红由川汗草谈轴象欲启哆埃随纹裸甜氖赏鹰钓店徊攀哥缆棉裁琉奖毯祈勺芜处泞倾拉是亨牛四舟恳锤禹呀惶熄够秋癌楚踌矗甩逾铣毒枕屯