PSS/E实验报告
实验名称：潮流计算（一）
班级：XXXXXXX 班
姓名：XXX
学号：XXXXXXX
任课教师：XXXX
实验老师：XXXX
实验时间：2018年11月17日
成绩：
实验二潮流计算（1）
一．实验目的（5分）
通过对9 节点 3 机系统的潮流计算，掌握PSS/E 软件的基本操作方法，提高对电力系统分析的感性认识，学会PSS/E 潮流分布图的绘制方法，并为今后的实验提供数据基础。
二．实验预习要求(5分)
实验前请完成以下预习要求：
1)简述动态仿真和静态仿真的区别；
动态仿真：电力系统遭受一些大的扰动，如短路故障等，系统的结果和参数发生重大改变时，对系统进行潮流计算分析。
静态仿真：静态仿真不涉及元件的动态特性和电力系统动态过程，本质上是系统运行的稳态分析。
2)简述高斯—塞德尔法解潮流和牛顿—拉夫逊法解潮流的步骤；
高斯塞德尔法：
1.设各节点电压初值，并给定迭代误差判据
2.对每一个PQ节点以前一次迭代的电压值来计算新值
3.求平衡节点注入功率
4.求支路功率分部和支路功率损耗
牛顿-拉夫逊法：
1.形成节点导纳矩阵
2.给各节点电压初值
3.求修正方程的常数项向量
4.求雅克比矩阵元素
5.求解修正方程是式，求修正量
6.求取节点电压新值
7.检查是否收敛
8.计算功率分布，PV节点无功功率和平衡节点注入功率
3)简述牛顿—拉夫逊法解潮流问题。
1.稀疏矩阵表示法
2.高斯消去法
3.节点的优化编号
4.牛顿拉夫逊法的收敛性
三．实验内容和步骤(60分)
1）打开已建立的原始数据文件
①打开软件，点击File>>New>>Case Data
②打开数据如图所示
③潮流计算前节点电压和发电机功率记录在表2-1
表2-1 潮流计算前的节点电压和发电机功率
节点电压发电机功率（MW/Mvar）
V1=1.5 V2=1.3 V3=1.025 Pgen1=0 Qgen1=0 Qmin gen1=-9999 V4=1.00 V5=1.00 V6=1.00 Pgen2=163 Qgen2=0 Qmin gen2=-9999 V7=1.00 V8=1.00 V9=1.00 Pgen3=85 Qgen3=0 Qmin gen3=-9999
2）进行潮流计算
Newton-Raphson法
①选择菜单栏>>Power Flow>>Solution>>Solve选项，得到如图的对话框，可以看到PSS/E使用的不同的潮流解法，Newton法和Gauss法。Newton 解法有三种，Gauss-Seidel 解法有两种。
选择潮流计算
牛顿-拉夫逊选项
②选择第一种解法：Fixed slope decoupled Newton-Raphson(固定斜率解耦法)
③点击“Slove”，开始潮流计算，点击“Close”关闭，计算后的结果如图所示
由图可知，使用固定斜率解耦法，潮流收敛，一共迭代了13次，系统总的修正量为0.07MVA。
④检查结果有无越限。选择>>Power Flow>>Repout>>limit checking report 得到如下图窗口，并可选择输出不同母线参数限制结果，如发电机参数限制，支路参数限制等，如图。
检查窗口
检查结果
由检查结果可知，潮流计算结果均无越限，在正常范围内。
④查看潮流计算结果。选择>>Power Flow>>Report>>Bus based reports.选择所
有母线都列表，得到潮流结果报表，如图
⑤记录结果表2-2
表2-2 潮流计算后节点电压和发电机功率
节点电压（KV）发电机功率(MW/MVAR)
V1=150 V2=130 V3=102.5 Pgen1=79.5 Qgen1=211.5 Qmin gen1=-9999.0 V4=141.91 V5=136.34 V6=132.17 Pgen2=163.0 Qgen2=15.8 Qmin gen2=-9999.0 V7=129.48 V8=122.80 V9=115.19 Pgen3=85 Qgen3=-220.3 Qmin gen3=-9999.0 结果分析：通过对比表2-1,2-2,可以看到，1,2,3,节点电压没有改变。
选择第二种解法：Full Newton-Raphson
①选择牛顿拉夫逊解法中的第二项，点击“Slove”，得到以下结果
迭代了2次
②记录表2-3
表2-3 潮流计算后的节点电压和发电机功率
节点电压（KV）发电机功率(MW/MVAR)
V1=150 V2=130 V3=102.5 Pgen1=79.4 Qgen1=211.5 Qmin gen1=-9999.0 V4=141.91 V5=136.34 V6=132.17 Pgen2=163.0 Qgen2=15.8 Qmin gen2=-9999.0 V7=129.48 V8=122.80 V9=115.19 Pgen3=85 Qgen3=-220.3 Qmin gen3=-9999.0
选择第三种解法：Decoupled Newton-Rsphson(P、Q解耦法)
迭代了一次
分析：观察三种解法下得出的数据，可以看到三种方法下计算出的潮流基本相同，并无明显差别。PQ分解法相较于牛顿拉夫逊法迭代次数更少，但结果一致，其原因是系统规模太小。对比三种牛顿拉夫逊解法，可以看到当电力系统网络数据量小，节点数少时，牛顿拉夫逊法下不同解法计算出的潮流基本一致。
Gauss法计算潮流
1.高斯塞德尔法
①在潮流求解窗口选择“Gauss”,求解
②得到结果
迭代18次，系统修正量为0.74MVA,潮流收敛
③记录表2-4
表2-4 潮流计算后节点电压和发电机功率
节点电压（KV）发电机功率(MW/MVAR)
V1=150 V2=130 V3=102.5 Pgen1=79.4 Qgen1=211.4 Qmin gen1=-9999.0 V4=141.91 V5=136.34 V6=132.17 Pgen2=163.0 Qgen2=15.8 Qmin gen2=-9999.0 V7=129.47 V8=122.80 V9=115.19 Pgen3=85 Qgen3=-220.3 Qmin gen3=-9999.0
2.修正高斯塞德尔法
①在潮流求解窗口选择“Modified Gauss-Seildel”,求解
②得到以下结果
迭代23次，系统修正量为0.31MVA,潮流收敛
③记录表格2-5
表2-5 潮流计算后节点电压和发电机功率
节点电压（KV）发电机功率(MW/MVAR)
V1=150 V2=130 V3=102.5 Pgen1=79.4 Qgen1=211.4 Qmin gen1=-9999.0 V4=141.91 V5=136.34 V6=132.17 Pgen2=163.0 Qgen2=15.8 Qmin gen2=-9999.0 V7=129.48 V8=122.80 V9=115.19 Pgen3=85 Qgen3=-220.3 Qmin gen3=-9999.0 分析：通过观察两组数据，可以看出所得结果相差无几，但是修正高斯-赛德尔
法所需更多次迭代，系统修正量更小。综合对比牛顿拉夫逊法和高斯-赛德尔法
可知，在系统较小，节点数很少时，两种方法计算潮流效果相同。
3）计算标准数据下的潮流，进行对比
①在软件中打开标准数据文件
②由于两种方法计算出的潮流结果基本相同，故选取牛顿-拉夫逊法进行计算
牛顿-拉夫逊法
1.选择“Full Newton-Raphson”法进行计算
2.得到以下结果
3.记录表3-1
表3-1 标准数据下计算潮流后节点电压
节点电压（KV）发电机功率(MW/MVAR)
V1=104 V2=102.5 V3=102.5 Pgen1=71.6 Qgen1=27 Qmin gen1=-9999.0 V4=102.58 V5=99.563 V6=101.27 Pgen2=163.0 Qgen2=6.7 Qmin gen2=-9999.0 V7=102.58 V8=101.59 V9=103.24 Pgen3=85.0 Qgen3=-10.9 Qmin gen3=-9999.0
4.对比分析：
在修改后的数据中，将1号母线发电机的输出电压104KV增大为150KV，2号发电机输出电压增大为130KV，将交流线路中4节点至5节点间的电阻
增大为0.02，电抗减小为0.08。对比两次计算结果和表2-3与表3-1，可以看到许多参数发生了大的变化。增大了发电机的端电压，明显的看出对发电机的有功出力影响不大，但是却严重减小发电机的无功出力，其无功输出量级发生了变化，可见端电压对无功输出影响较大。而三个发电机的有功出力并没有太大改变，可以说明发电机有功调节主要依赖于功角，本例中功角没有发生变化，因此有功变化不大。发电机端电压减小导致线路中各节点电压均有所减小。并且系统中某个节点参数发生变化，都会改变系统的潮流分析线路损耗，增大了4,5节点间的电阻，很显然从结构中可以看出，线路4,5之间的有功损耗增大了许多，虽然减小了电抗，但减小程度较小，且由于各节点电压均有增大，故所有线路的有功损耗和无功损耗都增大了。
再次以4,5节点为例，标准数据下，从4节点流向5节点的有功为41（MW），无功为22.9（MVar）；改变数据后，从4流向5节点的有功功率为
38.17,（MW），无功功率为71.6（MVar），因为增加了4、5线路间的电阻，
减小了电抗，显然有功略有减小，无功显著增大，说明高压线路中功率变化主要取决于电抗。并且节点电压幅值的改变引起了无功的显著变化，因此无功功率和电压幅值相关。
4）绘制潮流分布图
①在潮流计算结束后，直接通过 File>>new>>Diagram 新建一个画图文件。
②在界面菜单中选择 Diagram>Auto Draw Style，可设置为 Horizontal buses 或者 vertical buses,根据自己的喜好选择
③绘制潮流分部，如图
④进行保存。直接点击 save 进行保存，得到的文件格式将为*.sld。
四.思考题
1）分别记录完成表2-1 和表2-2 的数据；
2）请比较两表数据，你发现了什么？例如节点1 的电压有变化吗？对于节点9 上的无功出力超出其极限了吗？请给出电压有变化和无功力超出其极限的节点。平衡节点的电压有所变化吗？
3）请回答分别进行Newton-Raphson 和Gauss-Seidel 潮流计算时列出的所有问题。
4）将完成的潮流分布图分别以*.sld 和Word 的格式进行保存到你所建立的文件夹，以备考核。
五.学习体会
《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 2011年10月26日
一、实验题目： 数值积分 二、实验目的： 1．熟悉matlab 编写及运行数值计算程序的方法。 2．进一步理解数值积分的基础理论。 3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。 三、实验内容： 1．分别用复合梯形求积公式及复合辛普森求积公式计算积分xdx x ln 10 ? , 要求计算精度达到410-,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分dx x x ?+3 021,使误差不超过510-. 3．用3=n 的高斯-勒让德公式计算积分?3 1 sin x e x ,给出计算结果. 4．用辛普森公式(取2==M N ) 计算二重积分.5 .00 5 .00 dydx e x y ? ? - 四、实验结果： 1.（1）复合梯形法: 将区间[a,b]划分为n 等份，分点n k n a b h kh a x k ,2,1,0,,=-=+=在每个区间[1,+k k x x ](k=0,1,2,···n-1)上采用梯形公式，则得 )()]()([2)()(1 11 1 f R x f x f h dx x f dx x f I n n k k k b a n k x x k k ++===∑?∑? -=+-=+ 故)]()(2)([21 1 b f x f a f h T n k k n ++=∑-=称为复合梯形公式 计算步长和划分的区间 Eps=1E-4 h1=sqrt(Eps/abs(-(1-0)/12*1/(2+1))) h1 =0.0600 N1=ceil(1/h1) N1 =17 用复合梯形需要计算17个结点。 复合梯形： function T=trap(f,a,b,n) h=(b-a)/n;
课设报告 福建工程学院软件学院 题目：汇编计算器 班级： 1301 姓名 学号： 指导老师： 日期：
目录 1、设计目的 (3) 2、概要设计 (3) 2.1 系统总体分析 (3) 2.2 主模块框图及说明 (3) 3、详细设计 (4) 3.1 主模块及子模块概述 (4) 3.2各模块详运算 (4) 4、程序调试 (7) 4.1 运行界面分析 (7) 算法分析 (7) 4.2 调试过程与分析 (9) 5、心得体会 (11) 5.1 设计体会 (11) 5.2 系统改进 (11) 附录： (11)
1、设计目的 本课程设计是一次程序设计方法及技能的基本训练，通过实际程序的开发及调试，巩固课堂上学到的关于程序设计的基本知识和基本方法，进一步熟悉汇编语言的结构特点和使用，达到能独立阅读、设计编写和调试具有一定规模的汇编程序的水平。 2、概要设计 用8086汇编语言编写一个能实现四则混合运算、带括号功能的整数计算器程序。程序能实现键盘十进制运算表达式的输入和显示（例如输入：“1+2*(3-4)”），按“=”后输出十进制表示的运算结果。 2.1 系统总体分析 在8086的操作环境下，该计算器分成输入，数据存储，运算功能，输出几个大模块，实现了使用者使用该计算器时输入一个算式，能让系统进行计算。此计算器的实现功能是基本的数学的四则运算，结果范围在0~65535。 2.2 主模块框图及说明 此流程图简要的表现出了所要实现的功能以及一些功能的大概算法，同时也是我编写的一个总体的框架。 程序流程图说明：通过流程图，可以看出程序运行时，首先输出提示语气，当用户输入后，程序根据所输入内容进行判断，通过判断的结果来决定调用哪个功能模块，首要先要要判断的是否为0-9，“+”“-”“*”“/”这些字符，若不是就会报错，实则根据运算符号调用其功能模块完成运算。最后将运算的结果显示在主频幕上，返回主程序，使用户可以重新输入。
. / 《计算方法》上机实验报告 班级：XXXXXX 小组成员：XXXXXXX XXXXXXX XXXXXXX XXXXXXX 任课教师：XXX 二〇一八年五月二十五日
前言 通过进行多次的上机实验，我们结合课本上的内容以及老师对我们的指导，能够较为熟练地掌握Newton 迭代法、Jacobi 迭代法、Gauss-Seidel 迭代法、Newton 插值法、Lagrange 插值法和Gauss 求积公式等六种算法的原理和使用方法，并参考课本例题进行了MATLAB 程序的编写。 以下为本次上机实验报告，按照实验内容共分为六部分。 实验一： 一、实验名称及题目： Newton 迭代法 例2.7(P38):应用Newton 迭代法求在附近的数 值解，并使其满足. 二、解题思路： 设'x 是0)(=x f 的根,选取0x 作为'x 初始近似值,过点())(,00x f x 做曲线)(x f y =的切线L ,L 的方程为))((')(000x x x f x f y -+=,求出L 与x 轴交
点的横坐标) (') (0001x f x f x x - =,称1x 为'x 的一次近似值,过点))(,(11x f x 做曲线)(x f y =的切线,求该切线与x 轴的横坐标) (') (1112x f x f x x - =称2x 为'x 的二次近似值,重复以上过程,得'x 的近似值序列{}n x ,把) (') (1n n n n x f x f x x - =+称为'x 的1+n 次近似值，这种求解方法就是牛顿迭代法。 三、Matlab 程序代码： function newton_iteration(x0,tol) syms z %定义自变量 format long %定义精度 f=z*z*z-z-1; f1=diff(f);%求导 y=subs(f,z,x0); y1=subs(f1,z,x0);%向函数中代值 x1=x0-y/y1; k=1; while abs(x1-x0)>=tol x0=x1; y=subs(f,z,x0); y1=subs(f1,z,x0); x1=x0-y/y1;k=k+1; end x=double(x1) K 四、运行结果：
课程设计报告 学生姓名:学号： 学院:电气工程学院 班级: 题目: 电力系统潮流计算 职称: 副教授 指导教师：李翠萍职称: 副教授 2014年 01月10日
1 潮流计算的目的与意义 潮流计算的目的：已知电网的接线方式与参数及运行条件，计算电力系统稳态运行各母线电压、个支路电流与功率及网损。对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。 潮流计算的意义： (1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。 (3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。 (4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 2 潮流计算数学模型 1.变压器的数学模型： 变压器忽略对地支路等值电路：
2.输电线的数学模型： π型等值电路： 3 数值方法与计算流程 利用牛顿拉夫逊法进行求解，用MATLAB 软件编程，可以求解系统潮流分 布根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用matpower 进行潮流分析，将两者进行比较。 牛顿—拉夫逊法 1、牛顿—拉夫逊法概要 首先对一般的牛顿—拉夫逊法作一简单的说明。已知一个变量X 函数为： 0)(=X f 到此方程时，由适当的近似值) 0(X 出发，根据： ,......)2,1() ()() ()() () 1(='-=+n X f X f X X n n n n 反复进行计算，当) (n X 满足适当的收敛条件就是上面方程的根。这样的方 法就是所谓的牛顿—拉夫逊法。 这一方法还可以做下面的解释，设第n 次迭代得到的解语真值之差，即) (n X 的误差为ε时，则： 0)()(=+εn X f 把)() (ε+n X f 在) (n X 附近对ε用泰勒级数展开 0......)(! 2)()()()(2 )() () (=+''+ '+=+n n n n X f X f X f X f εεε 上式省略去2ε以后部分 0)()()()(≈'+n n X f X f ε
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日
y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"
心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。
(此文档为word格式，下载后您可任意编辑修改！) 2012级6班###(学号)计算机数值方法 实验报告成绩册 姓名:宋元台 学号： 成绩：
数值计算方法与算法实验报告 学期： 2014 至 2015 第 1 学期 2014年 12月1日课程名称: 数值计算方法与算法专业:信息与计算科学班级 12级5班 实验编号： 1实验项目Neton插值多项式指导教师：孙峪怀 姓名：宋元台学号：实验成绩： 一、实验目的及要求 实验目的： 掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。 实验要求： 1. 给出Newton插值算法 2. 用C语言实现算法 二、实验内容 三、实验步骤(该部分不够填写.请填写附页)
1.算法分析： 下面用伪码描述Newton插值多项式的算法： Step1 输入插值节点数n，插值点序列{x(i),f(i)},i=1,2,……,n,要计算的插值点x. Step2 形成差商表 for i=0 to n for j=n to i f(j)=((f(j)-f(j-1)(x(j)-x(j-1-i)); Step3 置初始值temp=1,newton=f(0) Step4 for i=1 to n temp=(x-x(i-1))*temp*由temp(k)=(x-x(k-1))*temp(k-1)形成 (x-x(0).....(x-x(i-1)* Newton=newton+temp*f(i); Step5 输出f（x）的近似数值newton(x)=newton. 2.用C语言实现算法的程序代码 #include
课程设计报告 学生:学号： 学院: 班级: 题目: 电力系统潮流计算课程设计
课设题目及要求 一 .题目原始资料 1、系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。 2、发电厂资料： 母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量为（ 100MW ），最大负荷和最小负荷分别为50MW 和20MW ；发电厂二总装机容量为（ 200MW ）。 3、变电所资料： (一) 变电所1、2、3、4低压母线的电压等级分别为：35KV 10KV 35KV 10KV (二) 变电所的负荷分别为： 60MW 40MW 40MW 50MW (三) 每个变电所的功率因数均为cos φ=0.85； 变电所1 变电所母线 电厂一 电厂二
(四) 变电所1和变电所3分别配有两台容量为75MVA 的变压器，短路损 耗414KW ，短路电压（%）=16.7；变电所2和变电所4分别配有两台容 量为63MVA 的变压器，短路损耗为245KW ，短路电压（%）=10.5； 4、输电线路资料： 发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为Ω17.0，单位长度的电抗为Ω0.402，单位长度的电纳为S -610*2.78。 二、 课程设计基本容： 1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。 2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷 情况下的潮流计算，并对计算结果进行分析。 3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。 1) 4个变电所的负荷同时以2%的比例增大； 2) 4个变电所的负荷同时以2%的比例下降 3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的 负荷同时以2%的比例上升； 4. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要 求，进行电压的调整。（变电所低压母线电压10KV 要求调整围在9.5-10.5 之间；电压35KV 要求调整围在35-36之间） 5. 轮流断开支路双回线中的一条，分析潮流的分布。（几条支路断几次） 6. 利用DDRTS 软件，进行绘制系统图进行上述各种情况潮流的分析，并进 行结果的比较。 7. 最终形成课程设计成品说明书。 三、课程设计成品基本要求： 1. 在读懂程序的基础上画出潮流计算基本流程图 2. 通过输入数据，进行潮流计算输出结果 3. 对不同的负荷变化，分析潮流分布，写出分析说明。 4. 对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并 列表表示调节控制的参数变化。 5. 打印利用DDRTS 进行潮流分析绘制的系统图，以及潮流分布图。
简单计算器 姓名: 周吉祥 实验目的：模仿日常生活中所用的计算器，自行设计一个简单的计算器程序，实现简单的计算功能。 实验内容： (1)体系设计： 程序是一个简单的计算器，能正确输入数据，能实现加、减、乘、除等算术运算，运算结果能正确显示，可以清楚数据等。 (2)设计思路： 1)先在Visual C++ 6.0中建立一个MFC工程文件，名为 calculator. 2)在对话框中添加适当的编辑框、按钮、静态文件、复选框和单 选框 3)设计按钮，并修改其相应的ID与Caption. 4)选择和设置各控件的单击鼠标事件。 5)为编辑框添加double类型的关联变量m_edit1. 6)在calculatorDlg.h中添加math.h头文件，然后添加public成 员。 7)打开calculatorDlg.cpp文件，在构造函数中，进行成员初始 化和完善各控件的响应函数代码。 (3)程序清单：
●添加的public成员： double tempvalue; //存储中间变量 double result; //存储显示结果的值 int sort; //判断后面是何种运算：1.加法2.减法3. 乘法 4.除法 int append; //判断后面是否添加数字 ●成员初始化： CCalculatorDlg::CCalculatorDlg(CWnd* pParent /*=NULL*/) : CDialog(CCalculatorDlg::IDD, pParent) { //{{AFX_DATA_INIT(CCalculatorDlg) m_edit1 = 0.0; //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME); tempvalue=0; result=0; sort=0; append=0; }
目录 一、设计任务和性能指标 (1) 1.1设计任务 (2) 1.2性能指标 (2) 二、设计方案 (2) 3 3 4 5 5 6 6 7 7 20 20 20 20 21 参考文献 (21) 附录1、系统硬件电路图 (22) 附录2、硬件实物图 (23) 附录3、器件清单 (24)
一、设计任务和性能指标 1.1设计任务 自制一个单片机最小系统，包括复位电路，采用外部小键盘输入数据，能够实现加法、乘法及一个科学计算，计算结果显示在四位一体的数码管上。 要求用Protel 画出系统的电路原理图（要求以最少组件，实现系统设计所要 显 位 监测模块采用二极管和扬声器（实验室用二极管代替）组成电路。 键盘电路采用4*4矩阵键盘电路。 显示模块采用4枚共阳极数码管和74ls273锁存芯片构成等器件构成。 整个单片机的接口电路： P0用于显示输出； P1用于键扫描输入； P2用于数码管位选控制； P3用于键盘扩展（部分运算符输入）；
三.系统硬件设计 3.1单片机最小系统 单片机最小系统就是支持主芯片正常工作的最小电路部分，包括主控芯片、复位电路和晶振电路。 主控芯片选取STC89C52RC芯片，因其具有良好的性能及稳定性，价格便宜应用方便。 扩展键：“log”，“ln”，“x^2”“小数点”，“开方” 共计25个按键，采用4*4矩阵键盘，键盘的行和列之间都有公共端相连，四行和四列的8个公共端分别接P1.0~P1.7，这样扫描P1口就可以完成对矩阵键盘的扫描，通过对16个按键进行编码，从而得到键盘的口地址，对比P1口德扫描结果和各按键的地址，我们就可以得到是哪个键按下，从而完成键盘的功能。 以下为键盘接口电路的硬件电路图
计算方法实验报告（四） 方程和方程组的迭代解法 一、实验问题 利用简单迭代法，两种加速技术，牛顿法，改进牛顿法，弦割法求解习题5-1，5-2，5-3中的一题，并尽可能准确。 选取5-3：求在x=1.5附近的根。 二、问题的分析（描述算法的步骤等） （1）简单迭代法算法： 给定初始近似值,求的解。 Step 1 令i=0； Step 2 令(计算)； Step 3 如果,则迭代终止，否则重复Step 2。 （2）Aitken加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令-(计算得到一个新的序列，其中k=0,1,2…)；Step 3 如果,则迭代终止，否则重复Step 2。 (3)插值加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令+(计算得到一个新的序列，其中k=1,2,3…)； Step 3 如果,则迭代终止，否则重复Step 2。 (4)牛顿法算法
Step 1给定初始近似值； Step 2令,其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （5）改进牛顿法的算法 Step 1给定初始近似值； Step 2令,其中k迭代计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （6）弦割法算法（双点弦割法） Step 1给定初始近似值，； Step 2令其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 三、程序设计 (1)简单迭代法 利用迭代公式进行迭代运算。 #include
东北电力大学课程设计改革试用任务书: 电力系统潮流计算课程设计任务书 设计名称：电力系统潮流计算课程设计 设计性质：理论计算，计算机仿真与验证 计划学时：两周 一、设计目的 1.培养学生独立分析问题、解决问题的能力； 2.培养学生的工程意识，灵活运用所学知识分析工程问题的能力 3.编制程序或利用电力系统分析计算软件进行电力系统潮流分析。 二、原始资料 1、系统图：IEEE14节点。 2、原始资料：见IEEE14节点标准数据库 三、课程设计基本内容： 1.采用PSAT仿真工具中的潮流计算软件计算系统潮流； 1)熟悉PSAT仿真工具的功能； 2)掌握IEEE标准数据格式内容； 3)将IEEE标准数据转化为PSAT计算数据； 2.分别采用NR法和PQ分解法计算潮流，观察NR法计算潮流中雅可比矩阵的变化情况， 分析两种方法计算潮流的优缺点； 3.分析系统潮流情况，包括电压幅值、相角，线路过载情况以及全网有功损耗情况。
4.选择以下内容之一进行分析： 1)找出系统中有功损耗最大的一条线路，给出减小该线路损耗的措施，比较各种措施 的特点，并仿真验证； 2)找出系统中电压最低的节点，给出调压措施，比较各种措施的特点，并仿真验证； 3)找出系统中流过有功功率最大的一条线路，给出减小该线路有功功率的措施，比较 各种措施的特点，并仿真验证； 5.任选以下内容之一作为深入研究：（不做要求） 1)找出系统中有功功率损耗最大的一条线路，改变发电机有功出力，分析对该线路有 功功率损耗灵敏度最大的发电机有功功率，并进行有效调整，减小该线路的损耗； 2)找出系统中有功功率损耗最大的一条线路，进行无功功率补偿，分析对该线路有功 功率损耗灵敏度最大的负荷无功功率，并进行有效调整，减小该线路的损耗； 3)找出系统中电压最低的节点，分析对该节点电压幅值灵敏度最大的发电机端电压， 并有效调整发电机端电压，提高该节点电压水平； 四、课程设计成品基本要求： 1.绘制系统潮流图，潮流图应包括： 1)系统网络参数 2)节点电压幅值及相角 3)线路和变压器的首末端有功功率和无功功率 2.撰写设计报告，报告内容应包括以下几点： 1)本次设计的目的和设计的任务； 2)电力系统潮流计算的计算机方法原理，分析NR法和PQ分解法计算潮流的特点； 3)对潮流计算结果进行分析，评价该潮流断面的运行方式安全性和经济性； 4)找出系统中运行的薄弱环节，如电压较低点或负载较大线路，给出调整措施； 5)分析各种调整措施的特点并比较它们之间的差异； 6)结论部分以及设计心得； 五、考核形式 1.纪律考核：学生组织出勤情况和工作态度等； 2.书面考核：设计成品的完成质量、撰写水平等； 3.答辩考核：参照设计成品，对计算机方法进行电力系统潮流计算的相关问题等进行答辩； 4.采用五级评分制：优、良、中、及格、不及格五个等级。
计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格
形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 1．理解数值计算稳定性的概念． 2．了解数值计算方法的必要性． 3．体会数值计算的收敛性与收敛速度． 三、实验内容 计算dx x x I n n ? += 1 10 ，1,2,,10n = ． 四、算法描述 由 dx x x I n n ? += 1 10 ，知 dx x x I n n ?+=--101110，则
科学计算器课程设计报告C课程设计 集团标准化工作小组 [Q8QX9QT-X8QQB8Q8-NQ8QJ8-M8QMN]
计算机科学与技术学部 C++课程设计 题目科学计算器 学部计算机科学与技术 班级计科1103 指导教师李军 姓名刘明 学号 2012年6月27日
摘要 计算器的产生和发展是建立在电子计算机基础之上的。硬件方面，自1946年第一台电子计算机诞生以来，计算机技术的发展可谓日新月异，从庞大的只能在实验室里供研究使用的计算机到如今能适应不同环境满足不同需求的各种各样的计算机；运算速度从每秒几千次到每秒几百亿次；处理器从焊有上百万个电子管的大的惊人的电子板到只有指甲大小的集成电路；现在计算机在硬件方面的发展已达到了每三个月更新换代一次的惊人速度。软件方面，也已从机器语言、汇编语言、高级语言发展到现如今的第四代语言——非结构化、面向对象、可视化的语言。 在这个计算器里面，我们实现了简单的四则运算以及更高功能的科学计算，它的外观简洁美观，使人们能快捷简单的操作。能准确的得到计算结果，大大减少了数字计算所需要的时间，为人们的生活带来便利。此系统在Windows 7环境下，使用VC++ 进行编写。 简单计算器包括双目运算和单目运算功能，双目运算符包含基本的四则运算及乘幂功能，单目运算符包含正余弦，对数，开方，阶乘，倒数，进制转换等运算。可对其输入任意操作数，包括小数和整数及正数和负数进行以上的所有运算并能连续运算。并且包含清除，退格功能等。我们所做的计算器其功能较Windows 7下的计算器还是很不够多，没有其菜单的实现功能项，没有其小巧的标准计算器。 关键词：计算器；运算；VC++等
目录 一、设计任务 (1) 1.1 课程设计要求 (1) 1.2 课程设计题目 (1) 1.3 课程设计基本容 (2) 二、问题分析 (3) 2.1 节点设置及分类 (3) 2.2 参数求取 (3) 2.3 计算方法 (4) 三、问题求解 (7) 3.1 等值电路的计算 (7) 3.2画出系统等值电路图： (7) 3.3 潮流计算 (8) 四、误差分析 (29) 五、心得体会及总结 (38) 附录: (39) 参考文献 (39) 程序 (39)
电力系统潮流计算课程设计 一、设计任务 1.1 课程设计要求 1、在读懂程序的基础上画出潮流计算基本流程图 2、通过输入数据，进行潮流计算输出结果 3、对不同的负荷变化，分析潮流分布，写出分析说明。 4、对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并 列表表示调节控制的参数变化。 5、打印利用DDRTS进行潮流分析绘制的系统图，以及潮流分布图。 1.2 课程设计题目 系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。 变电所1 变电所2 母线电厂一电厂二
发电厂资料： 母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量为（ 100MW ），最大负荷和最小负荷分别为50MW 和20MW ；发电厂 二总装机容量为（ 200MW ）。 变电所资料： (一)变电所1、2、3、4低压母线的电压等级分别为：35KV 10KV 35KV 10KV (二)变电所的负荷分别为：60MW 40MW 70MW 50MW (三)每个变电所的功率因数均为cos φ=0.85； (四)变电所1和变电所3分别配有两台容量为75MVA 的变压器，短路损耗414KW ，短路 电压（%）=16.7；变电所2和变电所4分别配有两台容量为63MVA 的变压器，短路损耗为245KW ，短路电压（%）=10.5； 输电线路资料： 发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为 Ω17.0，单位长度的电抗为Ω0.402，单位长度的电纳为S -610*2.78。 1.3 课程设计基本容 1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。 2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮 流计算，并对计算结果进行分析。 3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。 1) 4个变电所的负荷同时以2%的比例增大； 2) 4个变电所的负荷同时以2%的比例下降 3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以 2%的比例上升； 4. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要求，进行电 压的调整。（变电所低压母线电压10KV 要求调整围在9.5-10.5之间；电压35KV 要求调整围在35-36之间） 5. 轮流断开支路双回线中的一条，分析潮流的分布。（几条支路断几次） 6. 利用DDRTS 软件，进行绘制系统图进行上述各种情况潮流的分析，并进行结果的比 较。 7. 最终形成课程设计成品说明书。
高级语言程序（JAVA）课程设计报告 系部名称：商学系专业班级：营销*** 学生姓名：墨璇 墨兰学号： ********** ********** 指导教师：王芬教师职称：讲师 2014年06月26日
目录 一、课程设计目的及意义 .................................... 错误!未定义书签。 二、课程设计任务 .......................................... 错误!未定义书签。 2.1程序设计要求....................................... 错误!未定义书签。 三、课程设计时间 .......................................... 错误!未定义书签。 四、课程设计地点 .......................................... 错误!未定义书签。 五、课程设计内容 .......................................... 错误!未定义书签。 5.1开发工具与平台..................................... 错误!未定义书签。 (1).开发工具 ...................................... 错误!未定义书签。 (2).开发平台 ...................................... 错误!未定义书签。 5.2设计思路........................................... 错误!未定义书签。 5.3 程序测试 .......................................... 错误!未定义书签。 5.4实验总结........................................... 错误!未定义书签。 六、课程设计感想 .......................................... 错误!未定义书签。 七、附录（程序代码） ...................................... 错误!未定义书签。
数学计算方法实验报告 习题二 2.估计用二分法求方程f(x)=x3+4x2-10=0在区间[1,2]内根的近似值，为使方程不超过10时所需的二分次数。f(x k) 程序过程： function two (tolerance) a=1;b=2;counter=0; while (abs(b-a)>tolerance) c=(a+b)/2; fa=a^3+4*a^2-10;
fb=b^3+4*b^2-10; fc=c^3+4*c^2-10; if ((fa==0|fb==0)) disp(counter); elseif (fa*fc<0) b=c;counter=counter+1; elseif (fb*fc<0) a=c;counter=counter+1; elseif (fb==0) disp(counter); end end solution=(a+b)/2; disp(solution); disp(counter); 实验结果： 6.取x0=1.5,用牛顿迭代法求第三中的方程根.f(x)=x3+4x2-10=0的近似值（精确到||x k+1-x k|≦10-5,并将迭代次数与3题比较。 程序过程： function six (g) a=1.5; fa=a^3+4*a^2-10;
ga=3*a^2+8*a; b=a-fa/ga; k=1; while(abs(b-a)>g) a=b; fa=a^3+4*a^2-10; ga=3*a^2+8*a; b=a-fa/ga; k=k+1; end format long; disp(a); disp(k); 实验结果：程序结果计算结果 8.用弦割法求方程f(x)=x3-3x2-x+9=0在区间[-2,-1]内的一个实根近似值x k,|f(x k)|≦10-5. 程序过程： function eight (t) a=-2; b=-1; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=1; while(abs(c-b)>t) a=b; b=c; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=k+1; end
目录 摘要................................................. - 1 - 1.设计意义与要求..................................... - 2 - 1.1设计意义 ...................................... - 2 - 1.2设计要求（具体题目）........................... - 2 - 2.题目解析........................................... - 3 - 2.1设计思路 ...................................... - 3 - 2.2详细设计 ...................................... - 4 - 2.2.1节点类型.................................. - 4 - 2.2.2待求量 ................................... - 4 - 2.2.3导纳矩阵.................................. - 4 - 2.2.4潮流方程.................................. - 5 - 2.2.5牛顿—拉夫逊算法.......................... - 6 - 2.2.5.1牛顿算法数学原理：................... - 6 - 2.2.5.2修正方程............................. - 7 - 2.2.5.3收敛条件............................. - 9 - 3.结果分析.......................................... - 10 - 4.小结.............................................. - 11 - 参考文献............................................ - 12 -
南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用matlab编制以函数为基的多项式最小二乘拟合程序，并用于对下列数据作三次多项式最小二乘拟合(取权函数wi=1) x -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 y -2.30 -1 -0.14 -0.25 0.61 1.03 1.75 2.75 4.42 6.94 给定直线方程为：y=1/4*x3+1/2*x2+x+1 三、实验结论： 最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。 一般地。当测量数据的散布图无明显的规律时，习惯上取n次代数多项式。 程序运行结果为： a = 0.9731 1.1023 0.4862 0.2238 即拟合的三次方程为：y=0.9731+1.1023x+0.4862*x2+0.2238*x3
-2.5 -2-1.5-1-0.5 00.51 1.52 2.5 -4-20246 81012 x 轴 y 轴 拟合图 离散点 y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3 结论： 一般情况下，拟合函数使得所有的残差为零是不可能的。由图形可以看出最小二乘解决了残差的正负相互抵消的问题，使得拟合函数更加密合实验数据。 优点：曲线拟合是使拟合函数和一系列的离散点与观测值的偏差平方和达到最小。 缺点：由于计算方法简单，若要保证数据的精确度，需要大量的数据代入计算。
计算机技术综合课程设计 设计题目锅炉液位控制系统学生姓名史婷艳 专业班级自动化1302班学号20134460203 指导老师洪镇南 2017年1 月3日
目录 前言 (2) 1 锅炉汽包水位控制对象与控制指标 (4) 1.1锅炉汽包水位的特征 (4) 1.2汽包水位动态特性 (4) 1.2.1汽包水位在给水流量W作用下的动态特性 (4) 1.2.2汽包水位在蒸汽流量D扰动下的动态特性 (5) 1.2.3燃料量B扰动下汽包水位的动态特性 (6) 2. 汽包水位控制方案 (7) 2.1单冲量控制方式 (7) 2.2 双冲量控制方式 (8) 2.3 三冲量控制方式 (9) 3. 三冲量串级PID控制 (11) 3.1 串级PID控制 (11) 3.2 智能整定PID控制 (12) 4 汽包水位模糊控制器设计及仿真 (12) 4.1 输入输出变量 (12) 4.2 隶属度函数 (15) 4.3基于MATLAB/Simulink 环境建立的系统仿真分析 (16) 4.3.1 基于MATLAB/Simulink 的系统模型 (16)
4.3.2 仿真结果分析 (18) 总结与体会 (18) 参考文献 (20) 前言 锅炉是典型的复杂热工系统，目前，中国各种类型的锅炉有几十万台，由于设备分散、管理不善或技术原因，使大多数锅炉难以处于良好工况，增加了锅炉的燃料消耗，降低了效率。同时，锅炉工作过程中各项指标的调节难以建立数学模型，具有非线性、不稳定性、时滞等特点，所以如何改善对锅炉的控制，保证其正常工作，提高效率一直是人们关注的焦点。而汽包液位是锅炉安全、稳定运行的重要指标，保证液位在给定范围内，对于高蒸汽品质、减少设备损耗和运行损耗、确保整个网络安全运行具有要意义。 现代锅炉的特点之一就是蒸发量显著提高，汽包容积相对变小，水位变化速度很快，稍不注意就容易造成汽包满水或者烧成干锅，这都对汽包液位控制系统提出了更高的要求。汽包液位过高，会影响汽包内汽液分离效果，使汽包出口的饱和蒸汽带水增多，蒸汽带水会使汽轮机产生水冲击，引起轴封破损、叶片断裂等事故。同时会使饱和蒸汽中含盐量增高，降低过热蒸汽品质，增加在过热器管壁和汽轮机叶片上的结垢。水位过低，则可能破坏自然循环锅炉汽水循环系统中某些薄弱环节，以致局部水冷管壁被烧坏，严重时会造成爆炸事故。 目前，对汽包液位位控制大多采用常规PID控制方式，从控制方式来看，它们要么系统结构简单成本低，不能有效的控制锅炉汽包“虚假水位”现象，要么能够在一定程度上控制“虚假现象”，系统却过于复杂，成本投入过大。常用的蒸汽锅炉液位调节系统有三种基本结构：单冲量调节系统结构、双冲量调节系统结构、串级三冲量调节系统结