c++编写潮流上机程序（C++ writing flow program）.docVIP

c++编写潮流上机程序（C++ writing flow program） # includes stdio.h # includes math.h # includes stdlib.h # define m 矩阵阶数 100 / * * / # define n 100 迭代次数 / * * / int i, j, k, a, b, c; / * 循环控制变量 * / int t, l, g 中间变量; / * * / int n; / * 节点数 * / m, 支路数 / * * / fp, pq节点数 / * * / minutes, pv节点数 / * * / duidi 对地支路数; / * * / double 精度 eps; / * * / double aa [m], bb [m], cc [m], dd [m], max, time, rr, 中间变量); / * * / double mo, c1, d1, c2, d2, 复数运算函数的返回值 jiao; / * * / double g [m] [m], b [m] [m], y [m] [m]; / * 节点导纳矩阵中的实部、虚部及其模方值 * / double ykb [m] [m], d [m], du [m]; / * 雅克比矩阵、不平衡量矩阵 * / struct 节点结构体 jd / * * / (int num, d; / * num为节点号, s为节点类型 * / double p, q, s, e, f, u, zkj, dp, dq, du, de, df; / * 节点有功、无功功率, 功率模值, 电压纵、横分量, 电压模值, 阻抗角 牛顿 - 拉夫逊中功率不平衡量、电压不平衡量 * / } jd [m]. struct 支路结构体 zhl / * * / {int numb numb为支路号; / * * / int p1, p2 支路的两个节点; / * * / double r, x; / * 支路的电阻与电抗 * / } zhl [m]. fp1 file *, * fp2; void date (读取数据) / * * / { int h, number; fp1 = fopen (input.txt, r); if (fp1 = = null) { printf (can not open file! \ n ); exit (0); } fscanf (fp1, % d% d% d% d% d% lf \ n, n m, fp minutes, duidi, eps); / * 输入节点数, 支路数, pq节点数, pv节点数和精度 * / j = 1; k = pq + 1; for (i = 1, the = n; i + +) / * 输入节点编号、类型、输入功率和电压初值 * / { fscanf (fp1, % d% d number h). if (a = = 1 = 1是pq节点 类型h) / * * / { fscanf (fp1, %%%, if, if, if, if i% n, jd jd.p [j], [j].q, jd. jd [j], [j].f); jd = [j].num number; jd [j].s = h; j + +; } if (a = = 2) = 2是pv节点 类型h / * * / { fscanf (fp1, %% lf, lf i n, jd jd.p [k], [k].u); jd = [k].num number; jd.s [k] = h; jd.q [k] = 0; c + +; } if (a = = 3) = 3是平衡节点 类型h / * * / { fscanf (fp1, %% lf, lf i n, jd. jd [n], [n].f); jd = [n].num number; jd.s [n] = h; } } for (i = 1; i = m; i + +) 输入支路阻抗 / * * / fscanf (fp1 % d% d% d%,%, lf, lf i n [i], zhl.numb, zhl [i] [i].p1, zhl.p2, zhl [i] [i].r, zhl.x); fclose (fp1); if ((fp2 = fopen (output.txt, w)) = = null) { printf (can not open file! i n ); exit (0); } fprintf (fp2 潮流上机实习 on n, n ); fprintf (fp2, 班