实验四(华电版)
华中科技大学电力系统综合实验报告
电力系统综合实验报告第一部分综合实验台1、实验目的通过实验加深对电力系统暂态稳定性内容的理解，在对不同类型短路数据的分析中锻炼独立思考的能力，进一步了解不同短路故障对电力系统的危害。
实验方式为在理想实验台上模拟最简单的电力系统暂态稳定性问题，以期巩固学生在前一阶段的学习中对相关内容的掌握.2、实验原理与接线电力系统中不同类型的短路故障引起的最大短路电流可由下式得到,推导过程可参见《电力系统分析》一书相关章节内容。
单相接地短路：I f=3E aj(X ff1+X ff2+X ff0)两相相间短路:I f=√3E aj(X ff1+X ff2)两相接地短路：I f=√3√1−X ff2//X ff0X ff2+X ff0•E aj(X ff1+X ff2//X ff0)三相对称短路:I f=E aj(X ff1)如下图1实验接线模拟了单机无穷大系统.图1单机无穷大系统3、实验结果与数据分析⑴、不同故障类型对短路电流影响在下表中根据QF1～QF6的开断来选择单机无穷大系统的运行方式。
XL1接入双回线运行、XL2接入双回线运行.表格 1 短路切除时间0.5s单相接地短路实验数据根据以上表格得出以下结论：1）在各种不同类型的短路中,系统以双回线运行时短路电流较系统单回线运行时短路电流更大，与序网分析结果一致。
2)在各种不同类型的短路中，XL1接入时短路电流相对XL2接入时短路电流更小,以单相短路为例进行分析，可知接入XL阻抗越小，短路电流越大.判断实验台设置中XL1大于XL2。
3）对比各组实验数据，发现短路电流大小在不同短路类型中呈现有以下关系.两相相间短路>两相接地短路>三相接地短路>单相接地短路根据相关实验原理分析，由于X ff1≈X ff2≫X ff0,一般有三相接地短路>两相接地短路>单相接地短路>两相相间短路。
实验台中元件参数不可测量，经分析判断本次实验中负序阻抗偏小导致了两相短路电流偏大的现象发生。
华中科技大学电力系统综合实验报告
电力系统综合实验报告第一部分综合实验台1、实验目的通过实验加深对电力系统暂态稳定性内容的理解，在对不同类型短路数据的分析中锻炼独立思考的能力，进一步了解不同短路故障对电力系统的危害。
实验方式为在理想实验台上模拟最简单的电力系统暂态稳定性问题，以期巩固学生在前一阶段的学习中对相关内容的掌握。
2、实验原理与接线电力系统中不同类型的短路故障引起的最大短路电流可由下式得到，推导过程可参见《电力系统分析》一书相关章节内容。
单相接地短路：I f=3E aj(X ff1+X ff2+X ff0)两相相间短路：I f=√3E aj(X ff1+X ff2)两相接地短路：I f=√3√1−X ff2//X ff0X ff2+X ff0•E aj(X ff1+X ff2//X ff0)三相对称短路：I f=E aj(X ff1)如下图1实验接线模拟了单机无穷大系统。
图1单机无穷大系统3、实验结果与数据分析⑴、不同故障类型对短路电流影响在下表中根据QF1~QF6的开断来选择单机无穷大系统的运行方式。
XL1接入双回线运行、XL2接入双回线运行。
表格 1 短路切除时间0.5s单相接地短路实验数据根据以上表格得出以下结论：1）在各种不同类型的短路中，系统以双回线运行时短路电流较系统单回线运行时短路电流更大，与序网分析结果一致。
2）在各种不同类型的短路中，XL1接入时短路电流相对XL2接入时短路电流更小，以单相短路为例进行分析，可知接入XL阻抗越小，短路电流越大。
判断实验台设置中XL1大于XL2。
3）对比各组实验数据，发现短路电流大小在不同短路类型中呈现有以下关系。
两相相间短路>两相接地短路>三相接地短路>单相接地短路根据相关实验原理分析，由于X ff1≈X ff2≫X ff0，一般有三相接地短路>两相接地短路>单相接地短路>两相相间短路。
实验台中元件参数不可测量，经分析判断本次实验中负序阻抗偏小导致了两相短路电流偏大的现象发生。
数学模型实验四 综合实验
实验四 综合实验一、 实验目的：通过实验小结，布置小型研究问题（经过数学处理），使学生在练习过程中进一步熟悉MATLAB的使用，以及深入理解数学模型的建模思想。
为后续课程设计教学环节构筑基础。
二、 预备知识：1.具备数学分析、常微分方程、运筹学和概率论的学科知识基础；2.相关学科知识的简单求解方法以及辅助MATLAB求解相关问题。
三、 实验内容及要求(任选一题完成)：1、黄河小浪底调水调沙问题2004 年6 月至7 月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功。
整个试验期为20 多天，小浪底从6 月19 日开始预泄放水，直到7 月13 日恢复正常供水结束。
小浪底水利工程按设计拦沙量为75.5 亿m3，在这之前，小浪底共积泥沙达14.15 亿t。
这次调水调沙试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙，在小浪底水库开闸泄洪以后，从6 月27 日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7 月3 日达到最大流量2700m3/s，使小浪底水库的排沙量也不断地增加。
表7 是由小浪底观测站从6 月29 日到7 月10 检测到的试验数据。
现在，根据试验数据建立数学模型研究下面的问题：（1）给出估计任意时刻的排沙量及总排沙量的方法；（2）确定排沙量与水流量的关系。
2、炼油厂将A, B, C三种原油加工成甲、乙、丙三种汽油。
一桶原油加工成一桶汽油的费用为4元，每天至多能加工汽油14000桶。
原油的买入价、买入量、辛烷值、硫含量，及汽油的卖出价、需求量、辛烷值、硫含量由下表给出。
问如何安排生产计划，在满足需求的条件下使利润最大？一般说来，作广告可以增加销售，估计一天向一种汽油投入一元广告费，可使这种汽油日销量增加10桶，问如何安排生产和广告计划使利润最大？原油类别 买入价（元/桶） 买入量（桶/天）辛烷值硫含量（%）A 45 ≤5000 12 0.5B 35 ≤5000 6 2.0 C25≤50008 3.0汽油类别 卖出价（元/桶） 需求量（桶/天）辛烷值硫含量（%）甲 70 3000 ≥10 ≤1.0 乙 60 2000 ≥8 ≤2.0 丙501000≥6≤1.03、合金的强度y 与其中的碳含量x 有比较密切的关系，今从生产中收集了一批 数据如下表1。
三表法测电路参数实验报告
竭诚为您提供优质文档/双击可除三表法测电路参数实验报告篇一：用三表法测量电路等效参数实验报告(含数据处理)实验七用三表法测量电路等效参数一、实验目的1.学会用交流电压表、交流电流表和功率表测量元件的交流等效参数的方法。
2.学会功率表的接法和使用。
二、原理说明1.正弦交流信号激励下的元件的阻抗值，可以用交流电压表、交流电流表及功率表分别测量出元件两端的电压u、流过该元件的电流I和它所消耗的功率p，然后通过计算得到元件的参数值，这种方法称为三表法。
计算的基本公式为：up，电路的功率因数cos??IuIp等效电阻R＝2＝│Z│cosφ，等效电抗x＝│Z│sinφI阻抗的模Z?2.阻抗性质的判别方法可用在被测元件两端并联电容的方法来判别,若串接在电路中电流表的读数增大，则被测阻抗为容性，电流减小则为感性。
其原理可通过电压、电流的相量图来表示：图7-1并联电容测量法图7-2相量图(:三表法测电路参数实验报告) 3.本实验所用的功率表为智能交流功率表，其电压接线端应与负载并联，电流接线端应与负载串联。
三、实验设备DgJ-1型电工实验装置：交流电压表、交流电流表、功率表、自耦调压器、白炽灯、镇流器、电容器。
四、实验内容测试线路如图7-3所示，根据以下步骤完成表格7-1。
1.按图7-3接线，将调压器调到表1中的规定值。
2.分别测量15w白炽灯(R)、镇流器(L)和4.7μF电容器(c)的电流和功率以及功率因数。
3.测量L、c串联与并联后的电流和功率以及功率因数。
4.如图7-4，用并联电容法判断以上负载的性质。
Z图7-3图7-4五、实验数据的计算和分析根据表格7-1的测量结果，分别计算每个负载的等效参数。
up=2386.6，cos??=1IuIup镇流器L：Z?=551.7，cos??=0.172IuIup1电容器c：Z?=647.2，cos??=0，??2?f,|Z|?，f=50hz，因此c=4.9?FIuI?cupL和c串联：Z?=180.9，cos??=0.35；并联1?F电容后，电流增大，所以是容IuI白炽灯：Z?性负载L和c并联：Z?性负载由以上数据计算等效电阻R＝│Z│cosφ，等效电抗x ＝│Z│sinφ，填入表7-1中。
华北电力大学自动化系硕士研究生导师详细信息
本人为华电在校学生，见证了2011年华电考生如何通过内部题，轻而易举的搞定专业课。
如果你看过往届的华电题型，你会发现华电自动控制原理考试有一个明显的特点，那就是高度的重复，所以说如果你把近三年的题搞定，一般上135是一点问题没有的（除非华电换题型，可能性很小，从近几年出题方式上可以看出来）。
华电考生一般都在12月份才开始看专业课，凭借华电本校生的优势，搞到了内部的复习题，也就是常说的出题老师为照顾本校学生，根据真题出的复习题（根据往届真题和复习题对比发现，基本就是换一个数，或者已知变求解，求解换已知）。
说的再明白一点就是大家期末考试前老师给同学们划的范围，不过华电的这个范围更小，一共16~18到大题，原因很简单，华电专业课没有选择，填空，问答，全部是大题！！！所以说搞到内部复习题和真题比什么都重要。
华电自控专用书是胡寿松的四版自动控制原理，一共600多页，如果你想都看完了估计数学，政治，英语肯定一点复习时间没有，即使看完了你能保证每一个知识点都把握的牢固？！我不相信。
如果你有这样的实力，你可以清华了，果断放弃华电吧，呵呵。
不说废话，本人提供08~10三年内部复习题和08,10年真题，这些题一般情况下华电内部学生是不会透露给外校生的（除非是你华电内部有人，找关系搞到题，那得花多少钱啊，况且老师也不敢给你，出了事老师可就下岗了），所以说市面上是没有的，你可以随便查，网上肯定没有！至少现阶段我发现只有我独家提供，呵呵。
原因2点：1，货真价实！！！2，保证至少节约三个月的时间，这三个月复习数学，政治，英语会有质的突破。
都是学生，不会骗你们，花100~200买来全都是03~07年真题有什么用，还没有答案，连对错都不知道，本人提供详细答案（我自己做的，有时间我会用数学软件编辑出来，没时间发照片），具体到公式在专业书的第几页，什么思路，怎么来的！除此之外，本人免费提供03~07的真题,绝对货真价实！！！不多说了，联系方式QQ：1@1@8@2@06@69@64（实在抱歉去掉@符号）. 想要联系这个号，备注：华电。
华电潮流上机计算实验报告全解
院系：电气与电子工程学院班级：电气1205学号：1121181223学生姓名：王城指导教师：孙英云设计周数：两周成绩：日期：2015年7月7日一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、设计正文（详细内容见附录）1．手算2．计算机计算3．思考题三、课程设计总结或结论（详细内容见附录）四、参考文献1.《电力系统计算：电子数字计算机的应用》，西安交通大学等合编。
北京：水利电力出版社；2.《现代电力系统分析》，王锡凡主编，科学出版社；3.《电力系统稳态分析》，陈珩，中国电力出版社，1995年，第三版；附录（设计流程图、程序、表格、数据等）4.机算潮流程序及结果// dierti.cpp : 定义控制台应用程序的入口点。
//#include "stdafx.h"struct Line //线路结构体{int Num,NumI,NumJ; //线路号左节点名右节点名float R,X,B,K; //电阻电抗电纳变比（K等于1为普通支路，不等于1为变压器支路的变比） };struct Bus //节点结构体{int Num ;float Volt,Phase,GenP,GenQ,LoadP,LoadQ;int Type;};#include"stdio.h"#include"string.h"#include"math.h"#include"stdlib.h"#define NBUS 4#define NLINE 4/* Global variables */int nL,nB,nVA,nSH;float X[NBUS];int L;double def[2*NBUS];double mn[50];void Gauss(double a[50][50],double b[50], int n) /*定义高斯法 */{int JS[50];int i,j,k;float d,t,x[50];FILE *fp;int L=1;for(i=0;i<50;i++) JS[i]=0;for(k=0;k<n;k++){d=0.0;for(j=k;j<n;j++)if(fabs(a[k][j])>d){ /*在一行中找到一个最大值赋值d,并用JS[K]记住这个最大值所在的列号*/ d=fabs(a[k][j]);JS[k]=j;}if(fabs(d)<0.000001) /*如果d的数值太小，做为被除数将带来很大的误差 */L=0;else {if(JS[k]!=k)for(i=0;i<n;i++){t=a[i][k];a[i][k]=a[i][JS[k]]; /*进行列交换，让最大值始终在对角元上*/a[i][JS[k]]=t;}}if(L==0)break;for(j=k+1;j<n;j++)a[k][j]=a[k][j]/a[k][k]; /*对角元上的元素消为1*/b[k]=b[k]/a[k][k];for(i=k+1;i<n;i++){for(j=k+1;j<n;j++)a[i][j]=a[i][j]-a[i][k]*a[k][j]; /*使下三角阵的元素为0*/b[i]=b[i]-a[i][k]*b[k];}}if(fabs(a[n-1][n-1])>0.00001){ /*用追赶法，解方程组，求未知数x*/ x[n-1]=b[n-1];for(i=n-2;i>=0;i--){t=0.0;for(j=i+1;j<n;j++)t=t+a[i][j]*x[j];x[i]=(b[i]-t);}}if((fp=fopen("gauss.txt","w"))==NULL) /*将结果写到TXT文件中*/{printf("err");exit(0);}for(i=0;i<n;i++){fprintf(fp,"%lf",x[i]);mn[i]=x[i];fprintf(fp,"\n");}fclose(fp);if(fp!=NULL) fclose(fp);}int _tmain(int argc, _TCHAR* argv[]){FILE *fp;FILE *fpout;int i,j,k,l,h,n,v;int i1,i2,i3,kp,kq;float d1,d2,d3,d4,d5,d6,r,x,g,b,tt,LL,e,ps,qs,shsh,m;struct Line sL[NLINE];struct Bus sB[NBUS];float YG[NBUS+1][NBUS+1],YB[NBUS+1][NBUS+1];double u[50][2];i1=i2=i3=0;d1=d2=d3=d4=d5=d6=ps=qs=0.0;for(i=0;i<NBUS;i++)if((fp=fopen("in.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;i<nB;i++){sB[i].Num=sB[i].Type=0;sB[i].Volt=1.0;sB[i].Phase=sB[i].GenP=sB[i].GenQ=sB[i].LoadP=sB[i].LoadQ=0.0;fscanf(fp,"%d,%f,%f,%f,%f,%f,%f,%d",&i1,&d1,&d2,&d3,&d4,&d5,&d6,&i2);sB[i].Num=i1;sB[i].Volt=d1;sB[i].Phase=d2;sB[i].GenP=d3;sB[i].GenQ=d4;sB[i].LoadP=d5,sB[i].LoadQ=d6;sB[i].T ype=i2;};for(i=0;i<nL;i++){sL[i].Num=sL[i].NumI=sL[i].NumJ=0;sL[i].R=sL[i].X=sL[i].B=0.0;sL[i].K=1.0;fscanf(fp,"%2d %3d %3d %f %f %f %f",&i1,&i2,&i3,&d1,&d2,&d3,&d4);sL[i].Num=i1;sL[i].NumI=i2;sL[i].NumJ=i3;sL[i].R=d1;sL[i].X=d2;sL[i].B=d3;sL[i].K=d4;}if(fp!=NULL) fclose(fp);/*Make Y Matrix*/for(i=1;i<nB+1;i++)for(j=1;j<nB+1;j++){YG[i][j]=0.0;YB[i][j]=0.0;};for(l=0; l<nL; l++){i=sL[l].NumI;j=sL[l].NumJ;r=sL[l].R;x=sL[l].X;d1=r*r+x*x;g=r/d1;b=-x/d1;m=sL[l].K;if(fabs(sL[l].K-1.0)<0.000001) //普通支路 {YG[i][i]=YG[i][i]+g;YG[j][j]=YG[j][j]+g;YB[i][i]=YB[i][i]+b+sL[l].B;YB[j][j]=YB[j][j]+b+sL[l].B;YG[i][j]=YG[i][j]-g;YG[j][i]=YG[j][i]-g;YB[i][j]=YB[i][j]-b;YB[j][i]=YB[j][i]-b;}else //变压器支路{YG[i][i]=YG[i][i]+g/m+g*(m-1)/m;YG[j][j]=YG[j][j]+g/m+g*(1-m)/m/m;YB[i][i]=YB[i][i]+b/m+b*(m-1)/m;YB[j][j]=YB[j][j]+b/m+b*(1-m)/m/m;YG[i][j]=YG[i][j]-g/m;YG[j][i]=YG[j][i]-g/m;YB[i][j]=YB[i][j]-b/m;YB[j][i]=YB[j][i]-b/m; }}/* Check the Y matrix */if((fp=fopen("GGBB.txt","w"))==NULL){printf("Can not open the file named 'GGBB.txt' \n");exit(0);}fprintf(fp,"---Y Matrix---\n");for(i=1;i<nB+1;i++)for(j=1;j<nB+1;j++)if(fabs(YB[i][j]-0.0)>0.000001) fprintf(fp,"Y(%3d,%-3d)=(%10.5f,%10.5f)\n",i,j,YG[i][j],YB[i][j]);if(fp!=NULL) fclose(fp);/* 节点电压附初值 */for(i=1;i<nB+1;i++){if(sB[i-1].Type==0){u[i][0]=0.0;u[i][1]=1.0;}else if(sB[i-1].Type==1){u[i][1]=sB[i-1].Volt;u[i][0]=0.0;}else if(sB[i-1].Type==2){u[i][1]=sB[i-1].Volt;u[i][0]= sB[i-1].Phase;}}for(v=1;;v++)/* 迭代次数可以无限大 */{/* 节点电压附初值 */printf("迭代第%d次赋予的电压初值为e+jf:\n",v); for(i=1;i<nB+1;i++)printf("%lf,%lf\n",u[i][1],u[i][0]);printf("\n");printf("\n");/* 求偏移量 */double P_P[10];double P_Q[10];double P_UU[10];for(i=1;i<nB+1;i++){if(sB[i-1].Type==2){P_P[i]=0.0;P_Q[i]=0.0;P_UU[i]=1.05;}if(sB[i-1].Type==0){double tempP=0.0;double tempQ=0.0;for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0];tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];}P_P[i]=(sB[i-1].GenP-sB[i-1].LoadP)-tempP*u[i][1]-tempQ*u[i][0]; P_Q[i]=(sB[i-1].GenQ-sB[i-1].LoadQ)-tempP*u[i][0]+tempQ*u[i][1]; P_UU[i]=0.0;}if(sB[i-1].Type==1){double tempP=0.0;double tempQ=0.0;for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0];tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];P_P[i]=(sB[i-1].GenP-sB[i-1].LoadP)-tempP*u[i][1]-tempQ*u[i][0]; }P_UU[i]=sB[i-1].Volt*sB[i-1].Volt-u[i][1]*u[i][1]-u[i][0]*u[i][0]; P_Q[i]=0.0;}}/* 偏移量阵 */double P_PQ[6];int a=0;for(i=1;i<3;i++){P_PQ[a]=P_P[i];a=a+2;}a=1;for(i=1;i<3;i++){P_PQ[a]=P_Q[i];a=a+2;P_PQ[4]=P_P[3];P_PQ[5]=P_UU[3];printf("迭代第%d次的偏移量为：\n",v);for(i=0;i<6;i++){printf("%f",P_PQ[i]);printf("\n");}printf("\n");printf("\n");/* 雅可比矩阵 */double H[6][6],N[6][6],J[6][6],L[6][6],R[6][6],S[6][6],aa[6],bb[6]; for(i=1;i<5;i++){ if(fabs(sB[i-1].Type-2.0)<0.000001)continue;else{for(j=1;j<5;j++)if(i!=j){H[i][j]=-YB[i][j]*u[i][1]+YG[i][j]*u[i][0];N[i][j]=YG[i][j]*u[i][1]+YB[i][j]*u[i][0];J[i][j]=-N[i][j];L[i][j]=H[i][j];R[i][j]=0;S[i][j]=0;}else{aa[i]=bb[i]=0.0;aa[i]+=YG[i][n]*u[n][1]-YB[i][n]*u[n][0];bb[i]+=YG[i][n]*u[n][0]+YB[i][n]*u[n][1];}H[i][i]=-YB[i][i]*u[i][1]+YG[i][i]*u[i][0]+bb[i]; N[i][i]=YG[i][i]*u[i][1]+YB[i][i]*u[i][0]+aa[i]; J[i][i]=-YG[i][i]*u[i][1]-YB[i][i]*u[i][0]+aa[i]; L[i][i]=YG[i][i]*u[i][0]-YB[i][i]*u[i][1]-bb[i]; R[i][i]=2*u[i][0];S[i][i]=2*u[i][1];}}}double ss[50][50];for(i=0;i<6;i++)for(j=0;j<6;j++)ss[i][j]=0.0;for(i=1;i<3;i++)for(j=1;j<4;j++){ss[2*i-2][2*j-2]=H[i][j];ss[2*i-2][2*j-1]=N[i][j];ss[2*i-1][2*j-2]=J[i][j];ss[2*i-1][2*j-1]=L[i][j];}i=3;for(j=1;j<4;j++){ss[2*i-2][2*j-2]=H[i][j];ss[2*i-2][2*j-1]=N[i][j];ss[2*i-1][2*j-2]=R[i][j];ss[2*i-1][2*j-1]=S[i][j];}printf("迭代第%d次的雅可比矩阵为：\n",v);for(i=0;i<6;i++){for(j=0;j<6;j++)printf("%10f",ss[i][j]);printf("\n");}printf("\n");printf("\n");Gauss(ss,P_PQ,6);for(i=1;i<nB;i++){u[i][0]=u[i][0]+mn[2*(i-1)];u[i][1]=u[i][1]+mn[2*i-1];}double max;max=fabs(P_PQ[0]);for(i=0;i<=5;i++)if (max<fabs(P_PQ[i]))max=fabs(P_PQ[i]);if(fabs(max)<0.0001){printf("满足精度要求，迭代终止，迭代次数为%d\n",v); printf("\n");}/* 叠代循环的括号 */printf("最终求得的节点电压值为e+jf:\n");for(i=1;i<nB+1;i++)printf("%lf,%lf\n",u[i][1],u[i][0]);printf("\n");printf("\n");double uu[5],Phase[5];for(i=1;i<nB+1;i++){uu[i]=sqrt(u[i][1]*u[i][1]+u[i][0]*u[i][0]); Phase[i]=atan(u[i][0]/u[i][1]);}for(i=1;i<nB+1;i++)printf("%lf,%lf\n",uu[i],Phase[i]);*计算线路功率和平衡节点 PV节点功率*/double P[5],Q[5];double tempP=0.0;double tempQ=0.0;for(i=1;i<nB+1;i++){for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0]; tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];}P[i]=tempP*u[i][1]+tempQ*u[i][0];Q[i]=tempP*u[i][0]-tempQ*u[i][1];}for(i=1;i<nB+1;i++)printf("节点%d注入功率为%lf,%lf\n",i,P[i],Q[i]);/* 支路功率 */double V[4][2];for(i=1;i<5;i++)for(j=0;j<3;j++)V[i][j]=u[i][j];double sP[5][5],sQ[5][5];double dsq,dsp,dp,sumgen;for(i=1;i<NBUS+1;i++){for(j=1;j<NBUS+1;j++){sP[i][j]=0.0;sQ[i][j]=0.0;}}for(l=0; l<nL; l++){i=sL[l].NumI;j=sL[l].NumJ;r=sL[l].R;x=sL[l].X;d1=r*r+x*x;if(fabs(sL[l].K-1.0)<0.000001){/*Normal lines or transformers*/sP[i][j]=V[i][1]*V[i][1]*g-V[i][1]*V[j][1]*(g*cos(V[i][0]-V[j][0])+b*sin(V[i][0]-V[j][0]));sQ[i][j]=-(V[i][1]*V[i][1]*sL[l].B+V[i][1]*V[i][1]*b+V[i][1]*V[j][1]*(g*sin(V[i][0]-V[j][0])-b*cos(V[i ][0]-V[j][0])));sP[j][i]=V[j][1]*V[j][1]*g-V[i][1]*V[j][1]*(g*cos(V[j][0]-V[i][0])+b*sin(V[j][0]-V[i][0]));sQ[j][i]=-(V[j][1]*V[j][1]*sL[l].B+V[j][1]*V[j][1]*b+V[i][1]*V[j][1]*(g*sin(V[j][0]-V[i][0])-b*cos(V[j ][0]-V[i][0])));}else{/*abnormal transformer ratio*/sP[i][j]=V[i][1]*V[i][1]*g/sL[l].B/sL[l].B-V[i][1]*V[j][1]*(g*cos(V[i][0]-V[j][0])/sL[l].B+b*sin(V[i][ 0]-V[j][0])/sL[l].B);sQ[i][j]=-(V[i][1]*V[i][1]*b/sL[l].B/sL[l].B+V[i][1]*V[j][1]*(g*sin(V[i][0]-V[j][0])/sL[l].B-b*cos(V[i ][0]-V[j][0])/sL[l].B));sP[j][i]=V[j][1]*V[j][1]*g-V[i][1]*V[j][1]*(g*cos(V[j][0]-V[i][0])/sL[l].B+b*sin(V[j][0]-V[i][0])/sL[l ].B);sQ[j][i]=-(V[i][1]*V[i][1]*b+V[i][1]*V[j][1]*(g*sin(V[j][0]-V[i][0])/sL[l].B-b*cos(V[j][0]-V[i][0])/sL [l].B));}}/* 输电效率 */dsp=P[4];sumgen=P[4];for(i=0;i<NBUS;i++){dsp+=sB[i].GenP-sB[i].LoadP;dsq+=sB[i].GenQ-sB[i].LoadQ;sumgen+=sB[i].GenP;}dp=dsp/sumgen*100;/* 输出功率情况 */if((fp=fopen("功率情况.txt","w"))==NULL){printf("Can not open the file named '功率情况.txt' \n");exit(0);}fprintf(fp,"---功率情况---\n");fprintf(fp,"平衡节点功率S=%10.5f+ j%10.5f\n",P[4],Q[4]);for(i=1;i<NBUS+1;i++)for(j=1;j<NBUS+1;j++)if(fabs(sP[i][j]-0.0)>0.000001)fprintf(fp,"S(%3d,%-3d)=(%10.5f,j%10.5f)\n",i,j,sP[i][j],sQ[i][j]); fprintf(fp,"网损为%10.5f+j%10.3f,输电效率为%10.3f\n",dsp,dsq,100-dp);if(fp!=NULL) fclose(fp);return 0;}结果：1．导纳阵Y( 1,1 )=( 1.01534, -8.19201) Y( 1,2 )=( -0.56148, 2.30208) Y( 1,3 )=( 0.00000, 3.66667) Y( 1,4 )=( -0.45386, 1.89107) Y( 2,1 )=( -0.56148, 2.30208) Y( 2,2 )=( 1.04225, -4.67651) Y( 2,4 )=( -0.48077, 2.40385) Y( 3,1 )=( 0.00000, 3.66667) Y( 3,3 )=( 0.00000, -3.33333) Y( 4,1 )=( -0.45386, 1.89107) Y( 4,2 )=( -0.48077, 2.40385) Y( 4,4 )=( 0.93463, -4.26159)2．设定电压初值01.1;01;01)0(3)0(3)0(2)0(2)0(1)0(1j jf e j jf e j jf e +=++=++=+ 3．计算功率和电压偏移；27731.0])()([41)0(11)0(1)0(141)0(11)0(1)0(11)0(11)0(1-=++--=-=∆∑∑==j j j j j j jj s s e B f G ff B eG e P P P P0.05097])()([41)0(11)0(1)0(141)0(11)0(1)0(11)0(11)0(1-=----=-=∆∑∑==jj j j jj jj s s e B f G ef B e G f Q Q Q Q同理可算出52596.0)0(22)0(2-=-=∆P P P s ,0196.0)0(22)0(2=-=∆Q Q Q s 5.0)0(33)0(3=-=∆P P P s ,0.02)0(3232)0(3=-=∆U U U s 4．根据求的第一次迭代时雅可比矩阵各元素的公式计算雅可比矩阵各个元素的具体值：⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤-----------⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡20000.200066667.300003334.40052691.406629.130208.256148.00001821.182612.456148.030208.266667.3030208.256148.006298.803803.1066667.356148.030208.299265.032104.85．求高斯计算后的修正量：⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆-0.0000000.1276520.023566-0.108546-0.006511-0.007919-2)0(3)0(3)0(2)0(2)0(1)0(11)0()0(3)0(3)0(2)0(2)0(1)0(1U P Q P Q P e f e f e f J 6．计算各节点电压的一次近似值：12765.010855.000792.010000.197643.099349.0)0(3)0(3)1(3)0(2)0(2)1(2)0(1)0(1)1(1)0(3)0(3)1(3)0(2)0(2)1(2)0(1)0(1)1(1=∆+=-=∆+=-=∆+==∆+==∆+==∆+=f f f f f f fffe e e e e e e e e返回第三步重新迭代，并校验收敛与否，令410-=ε。
发电厂实验报告完整版
广西大学电气工程学院发电机运行实验报告同步发电机运行与控制专业班级：姓名：学号：实验地点：综合楼电机实验室一、实验目的同步发电机是电力系统最重要又最复杂的电气设备，在电力系统运行中起着十分重要的作用。
通过实验，使学生掌握和巩固同步发电机及其运行的基本概念和基本原理，培养学生的实践能力、分析能力和创新能力，加强工程实线训练，提高学生的综合素质。
二、实验装置及接线实验在电力系统监控实验室进行，每套实验装置以7.5KW直流电动机与同轴的5KW 同步发电机为被控对象，配置常规仪表测量控制屏（常规控制）和计算机监视控制屏（计算机监控）。
可实现对发电机组的测量、控制、信号、保护、调节、并列等功能，本次同步发电机运行实验，仅采用常规控制方式。
直流电动机-同步发电机组的参数如下：直流电动机：型号Z2-52，凸极机额定功率7.5kW额定电压DC220V额定电流41A额定转速1500r/min额定励磁电压DC220V额定励磁电流0.98A(5、6、7号机组为0.5A)同步发电机型号T2-54-55额定功率5kW额定电压AC400V(星接)额定电流9.08A额定功率因数0.8空载励磁电流 2.9A额定励磁电流5A直流电动机-同步发电机组接线如图一所示。
发电机通过空气开关2QS和接触器2KM 可与系统并列，发电机机端装有电压互感器1TV和电流互感器1TA，供测量、同期用，系统侧装有单相电压互感器2TV作同期用，两侧电压通过转换开关6SA接入同期表S （MZ-10）。
发电机励磁电源可以取自380V电网（他励方式），也可以取自机端（自励方式），通同步发电机运行实验过4QS进行切换，交流电源经励磁变压器CB降压隔离后，经分立元件整流装置或模块式晶闸管SCR-L变为直流，再通过灭磁开关3KM供电给发电机励磁绕组FLQ，励磁电流通过调压按钮或电位器2WR进行调节。
Rm为灭磁电阻，通过3KM的常闭触点与励磁绕组FLQ并接。
发电机组上面有一台用皮带带动的原作为励磁机用的直流发电机，在其励磁绕组加上恒定的直流电压（从开关稳压电源引来），则电枢上的电压正比于发电机组的转速，故用一只直流电压表即可测量发电机转子转速。
数学建模与数学实验(第4版)课件第18章
模拟的方法
1、物理模拟： 对实际系统及其过程用功能相似的实物系统去模仿。
例如，军事演习、船艇实验、沙盘作业等。
物理模拟通常花费较大、周期较长，且在物理模 型上改变系统结构和系数都较困难。而且，许多系统 无法进行物理模拟，如社会经济系统、生态系统等。
产生一个［０，１］均匀分布的随机数：rand
例 1的计算机模拟
3.产生 m n 阶均值为 ，方差为 的正态分布的随机数矩阵： normrnd ( , ,m, n)
产生一个均值为 ，方差为 的正态分布的随机数：normrnd( , )
•当研究对象视为大量相互独立的随机变量之和，且其中每 一种变量对总和的影响都很小时，可以认为该对象服从正态 分布。
4，5
6
k1=k1+1 k2=k2+1 k3=k3+1
N k1=k1+1
i＜20? Y N E=(k2+k3)/20 E1=0*k1/20+1*k2/20+2*k3/20
停止
4. 模拟结果
试验 序号
１ ２ ３ ４ ５ ６ ７ ８ ９ １０
投硬币 结 果
正 正 反 正 正 反 正 正 反 反
指示 正确
实验内容
1、模拟的概念。 2、产生随机数的计算机命令。 3、计算机模拟实例。 4、实验作业。
计算机模拟实例 离散系统模拟实例: 排队问题 连续系统模拟实例: 追逐问题 用蒙特卡洛法解非线性规划问题
返回
模拟的概念
模拟就是利用物理的、数学的模型来类比、模仿现实 系统及其演变过程，以寻求过程规律的一种方法。
电力系统分析实验报告四(理工类) (1)
西华大学实验报告（理工类）开课学院及实验室： 实验时间 ： 年 月 日学 生 姓 名学 号成 绩 学生所在学院 年级/专业/班 课 程 名 称课 程 代 码 实验项目名称 电力系统功率特性和功率极限实验项 目 代 码 指 导 教 师项 目 学 分一、实验目的1)初步掌握电力系统物理模拟实验的基本方法。
2)加深理解功率极限的概念，在实验中体会各种提高功率极限措施的作用。
3)通过对实验中各种现象的观察，结合所学的理论知识，培养理论结合实际及分析问题的能力。
二、实验原理所谓简单电力系统，一般是指发电机通过变压器、输电线路与无限大容量母线联接而且不计各元件的电阻和导纳的输电系统。
对于简单系统，如发电机至系统d 轴和g 轴总电抗分别为d X ∑和q X ∑，则发电机的功率特性为2()sin sin 2q d q Eq d d q E UX X U P X X X δδ∑∑∑∑∑-=+⨯⨯当发电机装有励磁调节器时，发电机电势q E 随运行情况而变化，根据一般励磁调节器的性能，可认为保持发电机'q E （或'E ）恒定。
这时发电机的功率特性可表示成''2'''''()sin sin 2Eq q d q d d q E UX X U P X X X δδ∑∑∑∑∑-=+⨯⨯ 或''''sin E q d E UP Xδ∑=这时功率极限为'''Em q d E UP X∑=随着电力系统的发展和扩大，电力系统的稳定性问题更加突出，而提高电力系统稳定性和输送能力的最重要手段之一，就是尽可能提高电力系统的功率极限。
从简单电力系统功率极限的表达式看，要提高功率极限，可以通过发电机装设性能良好的励磁调节器，以提高发电机电势、增加并联运行线路回路数；或通过串联电容补偿等手段，以减少系统电抗，使受端系统维持较高的运行电压水平；或输电线采用中继同步调相机、中继电力系统等手段以稳定系统中继点电压。
华电邹县发电厂四期工程7号机组性能试验工作总结PPT
组织相关技术人员进行培训，确保 他们熟悉并掌握性能试验的标准、 方案及操作规程。
人员准备
01
02
03
组建试验团队
人员分工
安全培训
挑选经验丰富、技能熟练的专业 技术人员组成试验团队，明确各 自的工作职责和任务。
根据试验内容和方案，对试验团 队进行合理分工，确保各项试验 工作的顺利进行。
组织全体试验人员进行安全培训 ，提高他们的安全意识和自我保 护能力。
华电邹县发电厂四期 工程#7号机组性能试 验工作总结
目录
• 试验背景及目的 • 试验前的准备工作 • 试验过程及方法 • 试验结果及讨论 • 结论及建议 • 参考文献及致谢
CHAPTER 01
试验背景及目的
背景介绍
华电邹县发电厂位于山东省邹城市，是华电集团的重要发电 企业之一。该电厂四期工程#7号机组是该电厂近年来新建的 大型燃煤机组，具备先进的性能和工艺。
02
试验过程中发现，机组的热效率、蒸汽参数和排放指标均达到
或超过设计要求，表明机组性能良好。
机组在试验期间未出现重大故障或异常情况，证明机组具有较
03
高的可靠性和稳定性。
建议与展望
建议继续加强机组维护和检修工作，确保设备运 行安全可靠。
展望未来，随着电力市场的变化和技术的不断进 步，应考虑引入更高效、环保的发电技术，以适 应市场需求和绿色发展要求。
试验过程及方法
试验过程描述
试验目的
对华电邹县发电厂四期工程#7号机 组进行性能试验，以评估其性能表 现及确定是否满足设计要求。
试验准备
选取合适的仪器设备，进行试验场 地布置，准备试验所需材料和工具 。
试验操作流程
按照预先设定的操作流程，依次进 行各个性能试验，记录试验数据。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验四 离散系统分析 一、实验目的 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。 二、 实验原理 MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。 例题：离散系统的单位脉冲响应 N=40; a=[1,0.4,-0.12]; b=[1,2]; y=impz(b,a,N); stem(y) xlabel('k');title('h[k]')
N=100; b=[1,2]; a=[1,0.4,-0.12]; x=ones(1,N); y=filter(b,a,x) y =
Columns 1 through 12 1.0000 2.6000 2.0800 2.4800 2.2576 2.3946 2.3131 2.3621 2.3327 2.3504 2.3398 2.3461
Columns 13 through 24 2.3423 2.3446 2.3432 2.3441 2.3436 2.3439 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438
Columns 25 through 36 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438
Columns 37 through 48 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438
Columns 49 through 60 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438
Columns 61 through 72 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3438 2.3438
Columns 73 through 84 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438
Columns 85 through 96 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438
Columns 97 through 100 2.3438 2.3438 2.3438 2.3438 >>
N=100; b=[1,2]; a=[1,0.4,-0.12]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi) y =
Columns 1 through 12 0.8400 2.7840 1.9872 2.5392 2.2228 2.4156 2.3005 2.3697 2.3282 2.3531 2.3381 2.3471
Columns 13 through 24 2.3417 2.3450 2.3430 2.3442 2.3435 2.3439 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438
Columns 25 through 36 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438
Columns 37 through 48 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438
Columns 49 through 60 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438
Columns 61 through 72 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3437 2.3438 2.3438 2.3438 Columns 73 through 84 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438
Columns 85 through 96 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438 2.3438
Columns 97 through 100 2.3438 2.3438 2.3438 2.3438 > 实验内容：
>计算该系统函数的零极点，并画出系统函数零极点分布图。 b=[1,2,0]; a=[1,0.4,-0.12]; z=roots(b) p=roots(a) zplane(b,a)
21112.04.0121)(zzz
zH 分析系统的幅频特性。 b=[1,1]; a=[1,-1,0.5]; [H,w]=freqz(b,a); plot(w,abs(H)) xlabel('Frequency(rad)'); ylabel('Magnitude'); title('Magnitude response');
2115.011)(zzz
zH 三、 实验内容 1.已知某LTI系统的差分方程为：
（1）初始状态 ， 输入 计算系统的完全响应。 （2）当以下三个信号分别通过系统时，分别计算离散系统的零状态响应： （3）该系统具有什么特性？ N=100; b=[0.0675,0.1349,0.0675]; a=[1,-1.143,0.412]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi)
y = Columns 1 through 12 0.3865 0.2322 0.3760 0.6041 0.8054 0.9416 1.0143 1.0413 1.0422 1.0322 1.0202 1.0108
]2[0675.0]1[1349.0][0675.0]2[412.0]1[143.1][kxkxkxkykyky2]2[,1]1[yy][][kukx Columns 13 through 24 1.0049 1.0020 1.0012 1.0015 1.0021 1.0027 1.0031 1.0033 1.0034 1.0035 1.0034 1.0034
Columns 25 through 36 1.0034 1.0034 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033
Columns 37 through 48 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033
Columns 49 through 60 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033
Columns 61 through 72 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033
Columns 73 through 84 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033
Columns 85 through 96 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033 1.0033
Columns 97 through 100 1.0033 1.0033 1.0033 1.0033
N=40; b=[0.0675,0.1349,0.0675];