《電力系統潮流上機》課程設計報告
院系:電氣與電子工程學院
班級: 電氣1108班
學號: 1111180806
學生姓名: 陳子君
指導教師: 劉寶柱
設計週數: 兩周
成績日期:2023年1月10日
一、課程設計的目的與要求
培養學生的電力系統潮流計算機程式設計能力,掌握計算機潮流計算的相關知識
二、設計正文(詳細內容見附錄)
1.手算: 要求應用牛頓-拉夫遜法或p-q分解法手算求解,要求精度為0.001mw。
節點1為平衡節點,電壓,節點2為pq節點,負荷功率,節點3是pv節點,,兩條支路分別為,,對地支路
2.計算機計算:編寫潮流計算程式,要求如下:
2.1據給定的潮流計算任務書整理潮流計算的基礎資料:節點的分類,線路模型,等值變壓器模型,電壓等級的歸算,標么值的計算;
2.2基礎資料的計算機儲存:節點資料,支路資料(包括變壓器);
2.3用牛頓-拉夫遜法計算;
2.4根據所選潮流計算方法畫流程圖,劃分出功能模組,有資料輸入模組,導納陣形成模組,解線性方程組模組,計算不平衡功率模組,形成雅可比矩陣模組,解修正方程模組,計算線路潮流,網損,pv節點無功功率和平衡節點功率,資料輸出模組;
2.5據上述模組編制程式並上機除錯程式,得出潮流計算結果;
2.6源程式及其程式中的符號說明集、程式流圖
簡單系統如下圖所示,支路資料如下:
,,,,, 節點資料如下:
,,,,
// 2014.cpp : 定義控制台應用程式的入口點。
//#include "stdafx.h"
struct line //線路結構體
;struct bus節點結構體
;#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 awe[100];
int _tmain(int argc, _tchar* argv)
fscanf(fp,"%d,%d,%d",&nb,&nl,&nsh);}
for(i=0;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].type=i2;
};for(i=0;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;ifor(j=1;j yg[i][j]=0.0;
yb[i][j]=0.0;
};for(l=0; 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) //普通支路
else變壓器支路
}/* check the y matrix */
if((fp=fopen("ggbb.txt","w"))==null)
fprintf(fp,"---y matrix---\n");
for(i=1;ifor(j=1;jif(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
else if(sb[i-1].type==1)
u[i][0]=sb[i-1].volt*cos(sb[i-1].phase);
u[i][1]=0.0;
}else if(sb[i-1].type==0)
}fp=fopen("電壓初值.txt","w");
if(fp==null)
fprintf(fp,"---電壓初值---\n");
for(i=1;ifprintf(fp,"e(%d)=%10.5f; f(%d)=%10.5f\n",i,u[i][0],i,u[i][1]);
if(fp!=null) fclose(fp);
printf("賦予的電壓初值為(e+if):\n");
for(i=1;iprintf("%lf+i%lf\n",u[i][0],u[i][1]);
printf("\n");
printf("\n");
迴圈迭代*/
/*求偏移量*/
double p_p[5];
double p_q[5];
double p_uu[5];
int v;
for(v=1;;v++)
{ for(i=1;i { if (sb[i-1].type==2)
{p_p[i]=0.0;
p_q[i]=0.0;
p_uu[i]=0.0;
王超的潮流上機課程設計報告
電力系統潮流上機 課程設計報告 院系 電氣與電子工程學院 班級 電網1102 學號 24 學生姓名 王超 指導教師 王莉麗 設計週數 兩周 成績日期 2013年12月30日 一 課程設計的目的與要求 培養學生的電力系統潮流計算機程式設計能力,掌握計算機潮流計算的相關知識 二 設計正文 詳細內容見附錄...
電力系統潮流上機課程設計報告
課程設計報告 2010 2011年度第一學期 名稱 電力系統潮流上機 院系 電氣與電子工程學院 班級 學號 學生姓名 指導教師 劉崇茹 設計週數 兩周 成績日期 2011年1月4日 一 課程設計的目的與要求 培養學生的電力系統潮流計算機程式設計能力,掌握計算機潮流計算的相關知識 二 設計正文 詳細內...
電力系統潮流上機設計教學指導書
1 c和.h程式可以用記事本編輯。2 將我提供的程式拷貝到d lx目錄下,並將options directories中的輸出目錄改為d lx,就可以執行。我提供的四個檔案,lx0.c導納陣形成檔案 導納陣形成檔案需要的標頭檔案,in.txt,潮流計算的原始資料 g0607.c,用高斯消去法解線方程組...