内容发布更新时间 : 2024/5/4 16:00:27星期一 下面是文章的全部内容请认真阅读。

C 语言程序设计

潮流计算

学 院 自动化学院 专业班级 学 号 姓 名 联系方式

word文档 可自由复制编辑

本程序潮流计算部分采用牛顿拉夫逊极坐标法进行计算，求解一次多元方程采用高斯列主元分解法进行求解。根据工程实际，在存储文件实时记录产生文件时间。此外本程序特意增加了文件查看功能，方便文件的查看。

程序代码：

#include #include #include #include #include #include #include

void dataprepare(void); void initial(void); void Yc(void); void showy(void); double detpqc(void); void showdetav(void); void Jrc(void);

void showdetpq(void); void Gauss(void); void showj(void);

void showsolution(void); char save2file(void); char option(void); int chose(void);

int list(char filename[],int i,char dstn[]); struct PQV {

char name[5]; double vb; double p; double q; double v; double ag; }*pqv; struct Y{ double G; double B; }**y; struct LZ

word文档 可自由复制编辑

{ char name1[5]; double vb1; char name2[5]; double vb2; double r; double x; double b; }*lz; struct TZ { char name1[5]; double vb1; char name2[5]; double vb2; double x; double k1; double k2; }*tz;

char fsource[20];

int pqnum=0,pvnum=0,lznum=0,tznum=0; int temp1=0,temp2=0,temp3=0,temp4=0; int diedai=0;

int numw=1,numofw=0;

double *detpq,*detav,*Pi,*Qi,error,wucha,**J; /* detpq节点功率的误差量，detav修正量，Pi节点的有功功率 */

char c,ch,dstname[20],filename[20]; int main() {

opt:c=option(); if(c=='1')goto opt1; else if(c=='2') { list(\ system(\ while(c!='0') { system(\ printf(\ ################################计算结果文件查看################################\\n\ printf(\ [8]上一个 [2]下一个 [5]查看 [-]删除文件 【0】退出\\n\\n\ list(\ chose();

word文档 可自由复制编辑