MATLAB数值计算绘图模拟仿真以及使用总结

Work1 1-1

1-2

1-3

1-4

2-1

2-1-（1）

2-2

2-3

2-3-(1)

2-4

2-4-(1)

2-5

2-6和2-7

3-1

3-1-（1）

4-1

4-2

4-2-（2）

4-2-3

5-1-（1）

5-1-2

5-1-（3）

5-2简述MATLAB命令窗的主要作用？

（1）命令窗口(Command Window)位于MATLAB 操作桌面的右方，用于输入命令并显示除图形以外的所有执行结果，是MATLAB 的主要交互窗口。

（2）Matlab既可以运行命令也可以执行程序，在命令窗口中可以运行单独的命令也可以调用程序，相当方便，而编辑调试窗口和图像窗口都是程序运行结果展示窗口，可以很直观的对程序运行过程中出现的矩阵或者是变量等等进行监视。

（3）在MATLAB 命令窗口中可以看到有一个“>>”，该符号为命令提示符，表示MATLAB 正在处于准备状态。在命令提示符后输入命令并按回车键后，MATLAB 就会解释执行所输入的命令，并在命令后面给出计算结果。

5-3简述MATLAB绘制二维图形的一般步骤

MATLAB绘制图形一般采取以下7个步骤：

（1）准备数据

（2）设置当前绘图区

（3）绘制图形

（4）设置图形中曲线和标记点格式

（5）设置坐标轴和网格线

（6）标注图形

（7）保存和导出图形

5-4启动Simulink的方式有几种？

1. 启动Simulink

启动Simulink通常有三种方式：

1)直接从Matlab指令窗口选取菜单File| New| Modal命令，Matlab将会打开Simulink库浏览器和名为untitled的模型窗口。 2)在Matlab命令窗口中键人Simulink命令，Matlab将会打开Simulink库浏览器。

3)点击Matlab命令窗口工具条的图标，启动Simulink库浏览器。

由启动Simulink的三种方式，要新建一个模型文件，至少可以采用两种方式：

1)直接从Matlab指令窗口选取菜单File|New|Modal命令。 2)先启动Simulink库浏览器，然

后点击Simulink库浏览器的工具条中的“新建模型”图标，建立新的模型文件。

如果模型文件已经存在，至少有三种方法打开模型文件： 1)从Matlab指令窗口选取菜单File|Open命令。

2)先启动Simulink库浏览器，然后点击Simulink库浏览器的工具条中的“打开模型”图标，打开已经存在的模型文件，对它进行编辑、修改和仿真。

3)在Matlab命令窗口中键人模型文件名称，不需要．mdl后缀。 Simulink用不同的窗口显示模块库、模型号和仿真输出图形结果．这些窗口不是Matlab图形窗口，不能用句柄图形命令来操作。

5-5

5-5-1

5-6

5-6-1

5-7

matlab电力系统潮流计算

华中科技大学信息工程学院课程设计报告书题目: 电力系统潮流计算专业：电气工程及其自动化班级：学号：学生姓名：指导教师： 2015年 11 月 10 日

2015年11月12日

信息工程学院课程设计成绩评定表

Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation

潮流计算(matlab)实例计算

潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。一、设计原理 1．牛顿-拉夫逊原理牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新

的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。牛顿—拉夫逊迭代法的一般步骤：（1）形成各节点导纳矩阵Y。（2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。（3）计算各个节点的功率不平衡量。（4）根据收敛条件判断是否满足，若不满足则向下进行。（5）计算雅可比矩阵中的各元素。（6）修正方程式个节点电压（7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。（8）计算平衡节点输出功率和各线路功率 2．网络节点的优化１）静态地按最少出线支路数编号这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。二、设计内容 1.设计流程图

MATLAB绘图函数

一、MA TLAB通用图形函数命令有关命令行环境的一些操作: (1) clc 擦去一页命令窗口,光标回屏幕左上角 (2) clear 从工作空间清除所有变量 (3) clf 清除图形窗口内容命令1 figure 功能创建一个新的图形对象。图形对象为在屏幕上单独的窗口，在窗口中可以输出图形。用法figure 用缺省的属性值创建一个新的图形对象。命令2 subplot 功能生成与控制多个坐标轴。把当前图形窗口分隔成几个矩形部分，不同的部分是按行方向以数字进行标号的。每一部分有一坐标轴，后面的图形输出于当前的部分中。用法subplot(m,n,p) 将一图形窗口分成m*n个小窗口，在第p个小窗口中创建一坐标轴。则新的坐标轴成为当前坐标轴。若p为一向量，则创建一坐标轴，包含所有罗列在p 中的小窗口。命令3 hold 功能保持当前图形窗口中的图形。该命令是决定是否在当前坐标轴中只能增加新的图形对象还是覆盖原有图形对象。用法hold on 保留当前图形与当前坐标轴的属性值，后面的图形命令只能在当前存在的坐标轴中增加图形。但是，当新图形的数据范围超出了当前坐标轴的范围，则命令会自动地改变坐标轴的范围，以适应新图形。 hold off 在画新图形之前，重新设置坐标轴的属性为缺省值。命令4 axis 功能坐标轴的刻度与外在显示用法axis([xmin xmax ymin ymax]) 设置当前坐标轴的x-轴与y-轴的范围。命令5 close 功能关闭指定的图形窗口。用法close 关闭当前的图形窗口。二、MA TLAB绘图参数控制命令1 plot 功能这是最基本、最常用的绘图函数，用于绘制线性二维图。有多条曲线时，循环使用由坐标轴颜色顺序属性定义的颜色，以区别不同的曲线；之后再循环使用由坐标轴线型顺序属性定义的线型，以区别不同的曲线。

matlab潮流计算

附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 % 牛拉法计算潮流程序 %----------------------------------------------------------------------- % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点； %------------------------------------------------------------------------ clear all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance='); pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %------------------------------------------------------------------ for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4); %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元1侧 end %求导纳矩阵 disp('导纳矩阵 Y='); disp(Y) %------------------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部

实验三 MATLAB绘图(含实验报告)

实验三 MATLAB 绘图一、实验目的 1．掌握二维图形的绘制。 2．掌握图形的标注 3．了解三维曲线和曲面图形的绘制。二、实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。设计提示 1．Matlab 允许在一个图形中画多条曲线：plot(x1，y1，x2，y2，……) 指令绘制y 1 = f 1(x 1), y 2 = f 2 (x 2 )等多条曲线。Matlab 自动给这些曲线以不同颜色。标注可用text 函数。 2．绘图时可以考虑极坐标和直角坐标的转换。 3．三维曲线绘图函数为plot3，注意参考帮助中的示例。三、实验内容 1．生成1×10 维的随机数向量a ，分别用红、黄、蓝、绿色绘出其连线图、杆图、阶梯图和条形图，并分别标出标题“连线图”、“杆图”、“阶梯图”、“条形图”。 2、绘制函数曲线，要求写出程序代码。 (1) 在区间[0:2π]均匀的取50个点，构成向量t (2) 在同一窗口绘制曲线y1=sin(2*t-0.3); y2=3cos(t+0.5)；要求y1曲线为红色点划线，标记点为圆圈；y2为蓝色虚线，标记点为星号。 (3) 分别在靠近相应的曲线处标注其函数表达式。 3.将图形窗口分成两个绘图区域，分别绘制出函数： ???+-=+=1 352221x x y x y 在[0,3]区间上的曲线，并利用axis 调整轴刻度纵坐标刻度，使1y 在[0,12] 区间上，2y 在[-2,1.5]区间上。 4．用mesh 或surf 函数，绘制下面方程所表示的三维空间曲面，x 和y 的

取值范围设为[-3，3]。 10102 2y x z +-= 思考题： 1. 编写一个mcircle(r)函数，调用该函数时，根据给定的半径r ，以原点为圆心,画一个如图所示的红色空心圆。（图例半径r=5）；左图参考 polar 函数的用法，右图绘制圆形的参数方程为x=sin （t ），y=cos （t ）。其中，t 的区间为0~2*pi ，步长为0.1。 2．（1）绘一个圆柱螺旋线（形似弹簧）图。圆柱截面直径为10，高度为5，每圈上升高度为1。如左图所示。（2）利用（1）的结果，对程序做少许修改，得到如右图所示图形。

matlab电力系统潮流计算

m a t l a b电力系统潮流计算 Final approval draft on November 22, 2020

华中科技大学信息工程学院课程设计报告书题目: 电力系统潮流计算专业：电气工程及其自动化班级：学号：学生姓名：指导教师： 2015年 11 月 10 日

信息工程学院课程设计成绩评定表

MATLAB绘图功能大全

MATLAB绘图功能大全-CAL-FENGHAI.-(YICAI)-Company One 1

Matlab绘图强大的绘图功能是Matlab的特点之一，Matlab提供了一系列的绘图函数，用户不需要过多的考虑绘图的细节，只需要给出一些基木参数就能得到所需图形，这类函数称为高层绘图函数。此外，Matlab还提供了直接对图形句柄进行操作的低层绘图操作。这类操作将图形的每个图形元素（如坐标轴、曲线、文字等）看做一个独立的对象，系统给每个对象分配一个句柄，可以通过句柄对该图形元素进行操作，而不影响其他部分。本章介绍绘制二维和三维图形的高层绘图函数以及其他图形控制函数的使用方法，在此基础上，再介绍可以操作和控制各种图形对象的低层绘图操作。一、二维绘图二维图形是将平而坐标上的数据点连接起来的平面图形。可以采用不同的坐标系，如直角坐标、对数坐标、极坐标等。二维图形的绘制是其他绘图操作的基础。（一）绘制二维曲线的基木函数在Matlab中，最基本而且应用最为广泛的绘图函数为plot,利用它可以在二维平面上绘制出不同的曲线。 1.plot函数的基木用法 plot函数用于绘制二维平面上的线性坐标曲线图，要提供一组x 坐标和对应的y坐标，可以绘制分别以x和y为横、纵坐标的二维曲线。plot函数的应用格式

plot(x,y)其中x,y为长度相同的向量，存储x坐标和y坐标。例52在［0,2pi］区间,绘制曲线程序如下：在命令窗口中输入以下命令？ ? x=0:pi/100:2*pi; ? y=2*exp*x).*sin(2*pi*x); ? plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线注意：指数函数和正弦函数之间要用点乘运算，因为二者是向量。例52绘制曲线这是以参数形式给出的曲线方程，只要给定参数向量，再分别求出x,y向量即可输岀曲线： ?t=-pi:pi/100:pi; ? x=t.*cos(3*t); ? y=t.*sin(t).*sin(t); ? plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线 ■ 以上提到plot函数的自变量x,y为长度相同的向量，这是最常见、最基木的用法。实际应用中还有一些变化。 2.含多个输入参数的plot函数

基于MATLAB的电力系统潮流计算

基于MATLAB的电力系统潮流计算 %简单潮流计算的小程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点%编号，将变压器的串联阻抗置于低压侧处理。 %第三列为支路的串列阻抗参数。 %第四列为支路的对地导纳参数。 %第五烈为含变压器支路的变压器的变比 %第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，%“0”为不含有变压器。 %B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点%负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中 %“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。 %X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地%参数。 n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); X=input('节点号和对地参数:X='); Y=zeros(n); Times=1; %置迭代次数为初始值 %创建节点导纳矩阵 for i=1:n1 if B1(i,6)==0 %不含变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %含有变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3);

用matlab电力系统潮流计算

题目：潮流计算与matlab 教学单位电气信息学院姓名学号年级专业电气工程及其自动化指导教师职称副教授

摘要电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLAB仿真。关键词：电力系统潮流计算 MATLAB Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation

MATLAB绘图总结

一、二维数据曲线图 1、MATLAB 最常用的画二维图形的命令是plot ，plot 函数的基本调用格式为：plot(x,y) 其中x 和y 为长度相同的向量，分别用于存储x 坐标和y 坐标数据。例1：在[0,2 ]画)sin(x 。生成的图形如下图1所示： 1 2 3 4 5 6 7 -1-0.8-0.6-0.4-0.200.20.40.60.81 图1 说明：（1）plot 函数的输入参数是矩阵形式时 A 、当x 是向量，y 是有一维与x 同维的矩阵时，则绘制出多根不同颜色的曲线。曲线条数等于y 矩阵的另一维数，x 被作为这些曲线共同的横坐标。 B 、当x,y 是同维矩阵时，则以x,y 对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。 C 、对只包含一个输入参数的plot 函数，当输入参数是实矩阵时，则按列绘制每列元素值相对其下标的曲线，曲线条数等于输入参数矩阵的列数；当输入参数是复数矩阵时，则按列分别以元素实部和虚部为横、纵坐标绘制多条曲线。（2）含多个输入参数的plot 函数调用格式为： plot(x1,y1,x2,y2,…,xn,yn) A 、当输入参数都为向量时，x1和y1，x2和y2，…，xn 和yn 分别组成一组向量对，每一组向量对的长度可以不同。每一向量对可以绘制出一条曲线，这样可以在同一坐标内绘制出多条曲线。 B 、当输入参数有矩阵形式时，配对的x,y 按对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。例2：如下所示的程序： x1=linspace(0,2*pi,100); x2=linspace(0,3*pi,100); x3=linspace(0,4*pi,100);

基于MATLAB的潮流计算源程序代码(优.选)

%*************************电力系统直角坐标系下的牛顿拉夫逊法潮流计算********** clear clc load E:\data\IEEE014_Node.txt Node=IEEE014_Node; weishu=size(Node); nnum=weishu(1,1); %节点总数 load E:\data\IEEE014_Branch.txt branch=IEEE014_Branch; bwei=size(branch); bnum=bwei(1,1); %支路总数 Y=(zeros(nnum)); Sj=100; %********************************节点导纳矩阵******************************* for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2); %末节点 R=branch(m,3); %支路电阻 X=branch(m,4); %支路电抗 B=branch(m,5); %支路对地电纳 k=branch(m,6); if k==0 %无变压器支路情形 Y(s,e)=-1/(R+j*X); %互导纳 Y(e,s)=Y(s,e); end if k~=0 %有变压器支路情形 Y(s,e)=-(1/((R+j*X)*k)); Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/((R+j*X)*k^2); Y(e,e)=-(k-1)/((R+j*X)*k); %对地导纳 end Y(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形 end for t=1:nnum; Y(t,t)=-sum(Y(t,:))+Node(t,12)+j*Node(t,13); %求支路自导纳 end G=real(Y); %电导 B=imag(Y); %电纳 %******************节点分类************************************* * pq=0; pv=0; blancenode=0; pqnode=zeros(1,nnum); pvnode=zeros(1,nnum); for m=1:nnum; if Node(m,2)==3 blancenode=m; %平衡节点编号 else if Node(m,2)==0 pq=pq+1; pqnode(1,pq)=m; %PQ 节点编号 else if Node(m,2)==2 pv=pv+1; pvnode(1,pv)=m; %PV 节点编号 end end end end %*****************************设置电压初值********************************** Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化 for n=1:nnum Uoriginal(1,n)=Node(n,9); %对各点电压赋初值 if Node(n,9)==0;

Matlab作图函数的总结与分析

高等理科教育２００５年第６期（总第６４期）Ｍａｔｌａｂ作图函数的总结与分析＋黄琼湘那斯尔江?吐尔逊（ｔｆｉ疆大学数学与系统科学学院，新疆乌鲁木齐８３００４６）摘要Ｍａｔｌａｂ（ＭＡＴｒｉｘＬＡＢｏｒａｔｏｒｙ的简称）是ＣｌｅｖｅｒＭｏｌｅｒ博士用Ｆｏｒｔｒａｎ语言开发的科学计算工具。它已成为科学研究、工程计算、应用开发的重要工具。国外已将它作为理工科大学的必修课程，国内各大学也开始开设这门课程。Ｍａｔｌａｂ有强大的作图功能，有兴趣的读者可参考文献【卜４’。本文对Ｍａｔｌａｂ的作图函数进行分析和总结，以供教学参考和学生学习之用。关键词Ｍａｔｌａｂ数据可视化作图函数中图分类号Ｇ６４２．０文献标识码Ａ一、Ｍａｔｌ如作图函数的总结Ｍａｔｌａｂ提供了丰富的作图函数，有１００个之多。在教学和学习中显得有点杂乱。我们先对它们进行总结和分类，并提炼出它们的共性和特性。Ｍａｔｌａｂ的作图函数从视角的维数上分有三类：一维作图函数、二维作图函数和三维作图函数。它们的代表分别是ｌｉｎｅ、ｐｌｏｔ和ｐｌｏｔ３等函数。从类型上分大致有四类：通用作图函数（如ｐｌｏｔ函数等）；专业作图函数（如ｃｏｎｔｏｕｒ函数、ｑｕｉｖｅｒ函数等）；动画制作函数（如ｍｏｖｉｅ、ｃｏｍｅｔ３等函数）；图形修饰函数（如ｖｉｅｗ等函数）。Ｍａｔｌａｂ所有的作图函数都可以通过查帮助获得它的功能和用法。这里我们把作图函数按类型分类，列出一些主要和常用的作图函数（见表１），以抓住重点。作图函数虽然功能不同，但它们的调用格式是一致的。我们用ＧｒａｐｈＦ来表示一般的作图函数，它们的调用格式如下：１．ＧｒａｐｈＦ（Ｘ，Ｙ，Ｓ）这是一、二维函数的作图格式。ｘ和Ｙ表示图形的数据点，ｓ表示图形修饰参数组（可以缺省）。当ｘ，Ｙ都是顶点坐标时，ＧｒａｐｈＦ（Ｘ，Ｙ，Ｓ）画出以ｘ，Ｙ为端点，ｓ为参数的线；当ｘ是一组顶点坐标，而Ｙ对应于Ｘ的函数值时，ＧｒａｐｈＦ（Ｘ，Ｙ，Ｓ）画出函数Ｙ的二维图形。２．ＧｒａｐｈＦ（Ｘ，Ｙ，Ｚ，Ｓ）这是三维函数的作图格式。ｚ是ｘ和Ｙ的函数。ｘ，Ｙ以二维坐标形式表示函数值ｚ的作图区域Ｄ，ｓ表示图形修饰参数组（可以缺省）。ＧｒａｐｈＦ（ｘ，Ｙ，Ｚ，Ｓ）画出定义域为Ｄ的函数ｚ的三维图形。值得注意的是，Ｍａｔｌａｂ的作图函数总是描绘数据点（Ｘ，Ｙ）（在平面上）或（ｘ，Ｙ，ｚ）（在空间中）的图形。前者视为Ｙ的函数，而后者视为ｚ的函数。函数ＧｒａｐｈＦ在作图前数据点必须事先给定，在作图时函数ＧｒａｐｈＦ将各数据点用光滑的曲线连接成图形。另外，Ｘ，Ｙ，Ｚ还 ÷收稿日期２００４—０２—１９资助项目新疆大学校基金“应用软件程序设计”重点课程建设项目资助作者简介黄琼湘（１９５８）男，湖南衡阳人，教授，主要从事组合数学与图论、计算机算法研究

基于matlab--psat软件的电力系统潮流计算课程设计

东北电力大学课程设计改革试用任务书: 电力系统潮流计算课程设计任务书设计名称：电力系统潮流计算课程设计设计性质：理论计算，计算机仿真与验证计划学时：两周一、设计目的 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.采用五级评分制：优、良、中、及格、不及格五个等级。

MATLAB下的潮流计算实现-稀疏技术毕业设计

毕业设计（论文）MATLAB下的潮流计算实现-稀疏技术

毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作者签名：日期：指导教师签名：日期：使用授权说明本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名：日期：

摘要电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。因此潮流计算是研究电力系统的一种很重要和很基础的计算。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流问题的方程式阶数越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。本文旨在于研究潮流计算的牛顿—拉夫逊法的基本原理，在Matlab环境中实现牛顿—拉夫逊法潮流计算的数学模型，程序流程以及编制相应程序，并在程序中融合了节点优化编号和稀疏技术，以提高计算效率。最后用IEEE-3O节点标准测试系统验证所编程序。关键词：潮流计算Newtom-Raphson法节点优化稀疏技术Matlab ABSTRACT Power flow calculation is fundanmental of analysis. Network reconfiguration,fault management,state estimator etc also need the data of electrial system power flow.There is important significance to develop power flow calculation in allusion to traits of distribution network. This paper introduces the principle of Newtom-Raphson algorithm, which is developed for calculation of power flow calculation ,where zero sequence network is open.With this algorithm,the three-phase load is resolved into positive/negative sequence power and coupling power,thus,decoupling three phase power flow into sequencet component power flow.The power flow can be obtained by just finding the positive sequence power flow and then finding the negative sequent component from the coupling https://www.360docs.net/doc/0685455.html,pared with the existing methods,the jacobian matrix with the proposed algorithm is of much lower order,thus substantially reducing the computation burden.The proposed algorithm,together with a reference algorithm,has been simulated on an actual IEEE-30 system using statistic load date.And then it will

实验5 Matlab绘图操作实验报告

Tutorial 5 实验报告实验名称：Matlab 绘图操作实验目的： 1、掌握绘制二维图形的常用函数； 2、掌握绘制三维图形的常用函数； 3、掌握绘制图形的辅助操作。实验内容： 1. 设sin .cos x y x x ?? =+ ??+?? 23051，在x=0~2π区间取101点，绘制函数的曲线。 2. 已知： y x =21，cos()y x =22，y y y =?312，完成下列操作：（1）在同一坐标系下用不同的颜色和线性绘制三条曲线；（2）以子图形式绘制三条曲线；（3）分别用条形图、阶梯图、杆图和填充图绘制三条曲线。 3. 已知：ln(x y x x ≤=??+>??0102 ，在x -≤≤55区间绘制函数曲线。 4. 绘制极坐标曲线sin()a b n ρθ=+，并分析参数a 、b 、n 对曲线形状的影响。 5．在xy 平面内选择区域[][],,-?-8888，绘制函数z =的三种三维曲面图。 6. 用plot 函数绘制下面分段函数的曲线。 ,(),,x x f x x x x ?+>? ==??+

8. 在同一坐标轴中绘制下列两条曲线。（1）.y x =-205 （2）sin()cos ,sin()sin x t t t y t t π=?≤≤? =?303 实验结果： 1. 2. （1）

（2）

（3）

Matlab牛拉法潮流计算程序

%本程序的功能是用牛顿——拉夫逊法进行潮流计算 % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点； % 3为PV节点； clear; n=input('请输入节点数:n='); nl=input('请输入支路数:nl='); isb=input('请输入平衡母线节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); % % %--------------------------------------------------- for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧 end %求导纳矩阵 disp('导纳矩阵Y='); disp(Y) %---------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部 e(i)=real(B2(i,3)); f(i)=imag(B2(i,3)); V(i)=B2(i,4); %PV节点电压给定模值 end for i=1:n %给定各节点注入功率 S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量 end %=================================================================== P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率 ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数IT2

MATLAB总结 - 隐函数、符号函数作图

I. 隐函数f(x,y)=0, f(x(t),y(t),z(t))=0; z=f(x,y) ezplot, ezplot3, ezcontour, conctourf, ezpolar, ezmesh, ezmeshc, ezsurf, ezsurfc 1. ezplot：画符号函数图形 ezplot（f）：对于显式函数f=f（x），在默认的范围[-pi> syms x y >> ezplot(2*x^4-y^9) 2. ezplot3：三维曲线图 ezplot3(x,y,z):在默认的范围0> syms t >> ezplot3(t*sin(t),t*cos(t),t,[0,20*pi]) 此外，三维曲线的视角可以通过程序命令来控制，也可以手动设置。 3. ezcontour：画符号函数的等高线图 ezcounter(f):画出二元符号函数f=f(x,y)的等高线图。函数f将被显示在默认的平面区域[-2pi