利用matlab仿真对电力系统谐波治理
摘要：随着国民经济和科学技术的蓬勃发展，冶金、化学等现代化大工业和电气化铁路的发展，电网负荷加大，电力系统中的非线性负荷(硅整流设备、电解设备、电力机车)及冲击性、波动性负荷(电弧炉、轧钢机、电力机车运行)使得电网发生波形畸变(谐波)、电压波动、闪变、三相不平衡，非对称性(负序)和负荷波动性日趋严重。电能质量的下降严重地影响了供用电设备的安全、经济运行，降低了人民的生活质量。所以在世界各国都十分重视电能质量的管理。
引言
新兴负荷的出现对电能质量的要求更高电能质量问题逐渐引起普遍重视，主要原因如下：
（1）大量基于计算机的控制设备和电子装置投入使用，其性能对电压质量非常敏感。
（2）调速电机和无功补偿装置，导致系统谐波水平不断上升，从而对电力系统的容量和安全运行产生影响。
（3）电力用户不断增长的电能质量意识迫使电力公司提高供电质量，设法解决诸如电压中断，电压跌落和开关暂态等电能质量问题。
衡量电能质量的主要指标是电网频率和电压质量。频率质量指标为频率允许偏差；电压质量指标包括允许电压偏差、允许波形畸变率(谐波)、三相电压允许不平衡度以及允许电压波动和闪变。国家技术监督局已公布了上述电能质量的五个国家标准。
电能质量的具体指标。
1．电网频率我国电力系统的标称频率为50Hz，GB／T15945-1995《电能质量一
电力系统频率允许偏差》中规定：电力系统正常频率偏差允许值为±0.2Hz，当系统容量较小时，偏差值可放宽到±0.5Hz，标准中没有说明系统容量大小的界限。在《全国供用电规则》中规定"供电局供电频率的允许偏差：电网容量在300万千瓦及以上者为±0.2HZ；电网容量在300万千瓦以下者，为±0.5HZ。实际运行中，从全国各大电力系统运行看都保持在不大于±0.1HZ范围内。
2．电压偏差GBl2325－90《电能质量一供电电压允许偏差》中规定：35kV及以上供电电压正负偏差的绝对值之和不超过额定电压的10％；10kV及以下三相供电电压允许偏差为额定电压的土7％；220V单相供电电压允许偏差为额定电压的7％～10％。
标准中供电电压为供电部门与用户产权分界处的电压或由供用电协议所规定的电能计量点的电压。
确定允许电压偏差是一个综合的技术经济问题，允许的电压偏差小，有利于用电设备的安全、经济运行，但为此要在电网中增添更多的无功电源和调压设备，需要更多的投入。反过来如果扩大用电设备对电压的适应范围，提高设备在这方面的性能，往往也要增加设备投资。综合国外标准和我国国情制定的供电电压允许偏差的国家标准，能满足绝大部分用电设备的运行要求。
3．三相电压不平衡GB／T15543-1995《电能质量一三相电压允许不平衡度》中规定：电力系统公共连接点正常电压不平衡度允许值为2％，短时不得超过4％。标准还规定对每个用户电压不平衡度的一般限值为1.3％。
但是国标规定的三相电压不平衡度的允许值及计算、测量和取值方法只适用于电力系统正常运行方式下在电网公共连接点由负序分量引起的电压不平衡。因此故障方式引起的不平衡(例如单相接地、两相短路故障等)和零序分量引起的不
平衡均不在考虑之列。由于电网中较严重的不平衡往往是由于单相或三相不平衡负荷所引起的，因此标准衡量点选在电网的公共连接点，以便在保证其它用户正常用电的基础上，给干扰源用户以最大的限值。值得注意的是国标在确定三相电压不平衡度指标时用95％概率作为衡量值。也就是说，标准中规定的"正常电压不平衡度允许值2％"是在测量时间95％内的限值，而剩余5％时间可以超过2％，过大的"非正常值"时间虽短，也会对电网和用电设备造成有害的干扰，特别是对有负序起动元件的快速动作的继电保护和自动装置，容易引起误动。因此标准中对最大的允许值作了"不得大于4％的规定。
4．公用电网谐波GB／T14549--93《电能质量-公用电网谐波》中规定：6～220kV 各级公用电网电压(相电压)总谐波畸变率是0.38kV为5.0％，6～10kV为4.0％，35～66kV为3.0％，110kV为2.0％；用户注入电网的谐波电流允许值应保证各级电网谐波电压在限值范围内，所以国标规定各级电网谐波源产生的电压总谐波畸变率是：0.38kV为2.6％，6～10kV为2.2％，35~66kV为1.9％，110kV为1.5％。对220kV电网及其供电的电力用户参照本标准110kV执行。
5．波动和闪变GBl2326--如《电能质量一电压允许波动和闪变》中规定：在公共供电点的电压波动允许值：10kV及以下为2.5，35～110kV为2％，220kV及以上为1.6％。
电压闪变值主要是表征人眼对灯闪主观感觉的参数。国标推荐的闪变干扰的允许值，对照明要求较高的白炽灯负荷为0.4％，对一般性照明负荷为0.6％
影响电能质量的原因各种各样，大体可以分为：
（1）内因。系统本身接有电弧炉、整流器、单相负荷、大功率电动机等干扰性负荷。这些负荷对电网产生负面影响，如谐波、无功冲击、负序等，而且这些负
面影响可能通过公共连接点（PCC）波及其它终端用户。因此，系统中必须安装相关装置，以及时缓解这些问题，而且还应根据电能质量评估体系，利用经济杠杆约束此类用户对电能质量的影响。
（2）外因。雷电、外力破坏、树枝影响、配电设备故障、电容器投切、线路切换等都可能干扰系统，造成断电或电压变动，甚至影响到相邻线路，导致有害影响蔓延。现在采取的措施，一是减少故障发生的次数和改变排除故障的方式，目前配电系统中的线路主保护是电流保护，该保护最大的缺陷是线路中相当大部分区域上的故障不能无时延地予以切除，此外即使无时延保护，从检测出故障到断路器开断故障，最快也需要3~6 个周波。若是永久性故障，多次重合闸则导致电压的不断波动。二是降低装置对电能质量问题的敏感性，主要是用户侧在敏感负荷或关键负荷处安装补偿装置，这种方法对单个负荷可有直接和明显的效果，但是受限于补偿装置的容量和价格，应用范围也受到限制。
电能质量分析方法
时域仿真方法在电能质量分析中的应用最为广泛，其最主要的用途是利用各种时域仿真程序对电能质量问题中的各种暂态现象进行研究。目前较通用的时域仿真程序有EMTP，EMTDC，NETOMAC等系统暂态仿真程序和SPICE，PSPICE，SABER 等电力电子仿真程序。采用时域仿真计算的缺点是仿真步长的选取决定了可模仿的最大频率范围。因此必须事先知道暂态过程的频率覆盖范围。此外，在模仿开关的开合过程时，还会引起数值振荡。
本次采用的是利用matlab中的sumlink对电力系统中谐波源仿真与治理
以三次谐波为例：
1.各元件的图标及参数：
交流电压源：幅值311 频率50hz
谐波源：幅值100A 频率150hz 产生3次谐波
RL并联支路：充当输电线路
谐波过滤器：L=1.1e-3 H，C=1e-3F
负载：额定电压220V 50hz
测量仪器：示波器电压电流表
2.利用上述器件及所学电路知识搭建如下电路: 治理谐波前:
示波器显示
观察可知此时负载侧电压已经被三次谐波干扰
3.利用LC谐振电路对三次谐波进行消除：
根据公式错误！未找到引用源。求出LC=1.1e-6，所以选取L=1.1e-3，C=1e-3 原理图如下：
示波器显示：
通过前后两次示波器比较可以看出经过LC谐振电路滤波之后电压由原来掺杂三次谐波变为比较完美的正弦波，负载侧电流也变为正弦波，谐波治理成功。
电本0826班
基于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);
电气2013级卓班电力电子技术与电力系统分析 课程实训报告 专业：电气工程及其自动化 班级： 姓名： 学号： 指导教师：
兰州交通大学自动化与电气工程学院 2016 年 1 月日
电力电子技术与电力系统分析课程实训报告 1 电力电子技术实训报告 1.1 实训题目 1.1.1电力电子技术实训题目一 一．单相半波整流 参考电力电子技术指导书中实验三负载，建立MATLAB/Simulink环境下三相半波整流电路和三相半波有源逆变电路的仿真模型。仿真参数设置如下： (1)交流电压源的参数设置和以前实验相关的参数一样。 (2)晶闸管的参数设置如下： R=0.001Ω，L =0H，V f=0.8V，R s=500Ω，C s=250e-9F on (3)负载的参数设置 RLC串联环节中的R对应R d，L对应L d，其负载根据类型不同做不同的调整。 (4)完成以下任务： ①仿真绘出电阻性负载（RLC串联负载环节中的R d= Ω，电感L d=0，C=inf，反电动势为0）下α=30°，60°，90°，120°，150°时整流电压U d，负载电流L 和晶闸管两端电压U vt1的波形。 d ②仿真绘出阻感性负载下（负载R d=Ω，电感L d为，反电动势E=0）α=30°，60°，90°，120°，150°时整流电压U d，负载电流L d和晶闸管两端电压U vt1的波形。 ③仿真绘出阻感性反电动势负载下α=90°，120°，150°时整流电压U d，负载电流L d和晶闸管两端电压U vt1的波形，注意反电动势E的极性。 (5)结合仿真结果回答以下问题： ①该三项半波可控整流电路在β=60°，90°时输出的电压有何差异?
摘要 本论文首先对国内外谐波抑制技术发展现状、有源电力滤波器原理与结构及三相瞬时无功功率理论进行了综述。重点研究了基于瞬时无功功率 理论。检测法及改进的电流移动平均值谐波检测法。在对电流移动平均值原理进行分析的基础上，给出了电流平均值谐波检测方案及实现检测的原理框图。 接着以MATLAB6.1软件包中的SIMULINK仿真环境为平台，构建了平均值谐波检测法的仿真模型;对电流平均值谐波检测方案进行了仿真研 究，并与基于滤波器的。谐波检测法的仿真结果进行了分析对比。结果表明，所采用的仿真方法与所构建的仿真模型不仅有效，而且证实了平均值谐波检测法比滤波器法有良好的动态响应性能。 在仿真基础上，提出了基于LF2407ADSP芯片电流平均值谐波检测法的数字实现方案，进而开发了三相并联型数字有源电力滤波器实验系统。进行了软、硬件设计。搭建的硬件电路包括:过零同步检测、电流和电压检测、PWM输出等几部分。采用模块化设计思路，用DSP汇编语言编写了系统软件，其中包括：ADC及中断处理、捕获及捕获中断处理、三相到两相电流转换、平均值法谐波计算、两相到三相变换、PI调节、PWM输出 控制及主程序等模块，并在软件开发系统下进行了调试。为实现电流同步采样处理，根据LF2407A事件管理器捕获单元特点，提出一种用软件实现锁相环的控制方法。 最后对有源电力滤波器进行了系统调试。实验结果表明，采用电流平均值谐波检测法结合软件锁相环控制方法能有效、准确的检测谐波，用该检测法开发的DSP有源电力滤波器实验系统，能够有效消除由非线性负载产生的谐波。关键词有源电力滤波器，瞬时无功功率，谐波电流检测，电流移动平均值，数字信号处理器 关键词有源电力滤波器，瞬时无功功率，谐波电流检测，电流移动平均值，数字信号处理器
基于Matlab的电力系统故障分析与仿真 摘要：本文介绍了MATLAB软件在电力系统中的应用，以及利用动态仿真工具Simulink和电力系统工具箱PSD进行仿真的基本方法。在仿真平台上，以单机—无穷大系统为建模对象，通过选择模块，参数设置，以及连线，对电力系统的多种故障进行仿真分析。同时，设计一个GUI图形界面，将仿真波形清晰地显示在界面上以便比较和分析。结果表明，仿真波形基本符合理论分析，说明了MATLAB是电力系统仿真研究的有力工具。 关键词：电力系统；仿真；故障；MATLAB；GUI Abstract：This paper introduces the applications of MATLAB in power system analysis, and the basic simulation method of taking use of Simulink and PSD. On MATLAB simulation platform, take a single machine-infinite-bus system as modeling objects, by selecting the module, parameter settings, and connectingmodules to simulate and analysevariousfault of power system. At the same time, in order to facilitate comparison and analysis simulation waveform, design a GUI for showing waveform clearly.The results show that the simulation waveform in line with theoretical analysis, indicates that MATLAB is a powerful tool for researching simulation of power system. Keywords：PowerSystem。 Simulation。 Fault。 Matlab。 GUI 0 前言[1,2] 随着电力工业的发展，电力系统规划、运行和控制的复杂性亦日益增加，电力系统的生产和研究中仿真软件的应用也越来越广泛。现在，我们主要使用的电力系统仿真软件有：EMTP程序，用于电力系统电磁暂态计算，电力系统暂态过电压分析，暂态保护装置的综合选择等。PSCAD／EMTDC程序，典型应用是计算电力系统遭受扰动或参数变化时，参数随时间变化的规律。PSASP，其功能主要有稳态分析、故障分析和机电暂态分析。还有MathWorks公司开发的MATLAB软件。在MATLAB中，电力系统模型可以在Simulink环境下直接搭建，也可以进行封装和自定义模块库，充分显现了其仿真平台的优越性。更重要的是，MATLAB提供了丰富的工具箱资源，以及大量的实用模块，使我们可以更加深入地研究电力系统的行为特性。本篇论文将在熟练掌握MATLAB软件的基础上，对电力系统的故障进行建模、仿真、分析，并且设计一个GUI图形用户界面来反映故障波形。
综合训练① 实验内容：利用matlab绘制频率自定的正弦信号（连续时间和离散时间），复指数信号（连续时间），并举例实际中哪些物理现象可以用正弦信号，复指数信号来表示。绘制成谐波关系的正弦信号（连续时间和离散时间），分析其周期性和频率之间的关系。实验步骤： 一、绘制谐波关系的正弦信号 分析：由于正弦信号可以表示成两个共轭的复指数信号相减，然后再除去两倍的单位虚数得到，故，我们将正弦信号设置为 X=exp(j*pi*n/4)-exp(-j*pi*n/4))/(2*j) 此信号就相当于 x=sin(pi*n/4) 设计程序如下： n=[0:32]; %设置n的取值 x=(exp(j*pi*n/4)-exp(-j*pi*n/4))/(2*j); %限定离散正弦信号 stem(n,x) %绘制该离散正弦信号 通过Matlab所得图形如下：
分析：同样的连续型的正弦信号同样也可以用类似方式绘制. x=sym('(exp(j*pi*t/T)+exp(-j*pi*t/T))/2');%函数表示正弦信号 x5=subs(x,5,'T'); %设置周期大小ezplot(x5,[0,10]) %绘制图形 所得结果如下：
二、绘制复指数信号 分析：由于复指数信号有实数部分和虚数部分，所以绘制其图形，我们采取了分别绘制的方法，将实数和虚数分别画出。 实验程序如下： t=[0:.01:10]; %产生时间轴的等差点 y=exp((1+j*10)*t); %设置复指数信号 subplot(211),plot(t,real(y)); %绘制实数信号图形 grid subplot(212),plot(t,imag(y)); %绘制虚数部分图形 grid 实验所得结果如下：
湖北民族学院 信息工程学院 题目: 基于matlab的电力系统短路电流计算 专业：电气工程及其自动化 班级： 0308407 学号： 030840705 学生姓名： 指导教师： 2011年6 月1 日
信息工程学院课程设计任务书 年月日
信息工程学院课程设计成绩评定表
摘要 随着电力工业的发展，电力系统的规模越来越大，在这种情况下，许多大型的电力科研实验很难进行，尤其是电力系统中对设备和人员等危害最大的事故故障，尤其是短路故障，而在分析解决事故故障时要不断的实验，在现实设备中很难实现，一是实际的条件难以满足；二是从系统的安全角度来讲也是不允许进行实验的。考虑这两种情况，寻求一种最接近于电力系统实际运行状况的数字仿真工具十分重要，而MATLAB软件中的SIMULINK是用来对动态系统进行建模、仿真和分析的集成开发环境，是结合了框图界面和交互仿真能力的非线性动态系统仿真工具，为解决具体的工程问题提供了更为快速、准确和简洁的途径。 关键词：短路电流计算，MATLAB，仿真 Abstract Along with the development of the electric power industry, the scale of the power system is more and more big, in this case, many large power research is difficult to, especially in the power system, equipment and personnel to the harm such as the biggest accident, especially fault fault location, and on the analysis of the accident to solve the fault of the experiment, in the reality constantly in equipment, it is difficult to accomplish a is practical conditions to meet; The security of the system from the perspective is not allowed in the experiment. Consider the two kinds of circumstances, for one of the most close to power system actual the operation condition of digital simulation tool is very important, and MATLAB software SIMULINK is used for the dynamic system modeling, simulation and analysis of the integrated development environment, is combined with the block diagram interface and interactive simulation of nonlinear dynamic system ability of simulation tools, to solve the specific engineering problem, provides a more rapid, accurate and simple way. Keywords: short-circuit current calculation, MATLAB, the simulatio
为了便于进行周期信号的分析与处理，常要把复杂的周期信号进行分解，即将周期信号分解为正余弦等此类基本信号的线性组合，通过对这些基本信号单元在时域和频域特性的分析来达到了解信号特性的目的。本文主要阐述了傅立叶级数的推演过程，从而得出周期信号的分解与合成的基本原理。 1 绪论 研究信号是为了对信号进行处理和分析，信号处理是对信号进行某些加工或变换，目的是提取有用的部分，去掉多余的部分，滤除各种干扰和噪声，或 将信号进行转化，便于分析和识别。信号的特性可以从时间特性和频率特性两 方面进行描述，并且信号可以用函数解析式表示（有时域的，频域的及变化域 的），也可用波形或频谱表示。 系统分析的主要任务是分析系统对指定激励所产生的影响。其分析过程主要包括建立系统模型，根据模型建立系统的方程，求解出系统的响应，必要时 对解得的结果给出物理解释。系统分析是系统综合与系统诊断的基础。 任何满足狄里赫利条件的周期信号都是由各种不同频率、幅度和初相的正弦波叠加而成的。对周期信号由它的傅立叶级数展开式可知，各次谐波为基波 频率的整数倍。而非周期信号包括了从零到无穷大的所有频率成分，每一个频 率成分的幅度均趋向无穷小，但其相对大小式不同的。 信号的分解与合成 周期信号的信号分解与合成 设有周期信号，它的周期为T，角频率，则的三角傅里 叶级数表示的一般形式为 （2.2-1） 其中
可以写成更紧凑的和式为： 式（2.2-1）中的系数、称为傅里叶系数，为在函数中的分量（相对大小）；为在函数中的分量，它可由式（2.1-7）求 得。为简便，式（2.1-7）的积分区间取为或。考虑到正、余弦函数的正交条件（2.1-3），由式（2.1-7），可得傅立叶系数 （2.2-2） 周期信号也可分解为一系列余弦信号，即： 其中 方波信号的分解与合成 现以周期为T、幅值为1的方波信号为例 方波信号的分解与合成【12】 由式（2）可得
《电力系统设计》报告题目: 基于MATLAB的电力系统仿 学院：电子信息与电气工程学院 班级: 13级电气 1 班 姓名:田震 学号: 日期：2015年12月6日 基于MATLAB的电力系统仿真 摘要：目前，随着科学技术的发展和电能需求量的日益增长，电力系统规模越来 越庞大，超高压远距离输电、大容量发电机组、各种新型控制装置得到了广泛的应用，这对于合理利用能源，充分挖掘现有的输电潜力和保护环境都有重要意义。另一方面，随着国民经济的高速发展，以城市为中心的区域性用电增长越来越快，大电网负荷中心的用电容量越来越大，长距离重负荷输电的情况日益普遍，电力系统在人们的生活和工作中担任重要角色，电力系统的稳定运行直接影响着人们的日常生活。从技术和安全上考虑直接进行电力试验可能性很小，因此迫切要求运用电力仿真来解决这些问题。 电力系统仿真是将电力系统的模型化、数学化来模拟实际的电力系统的运行，可以帮助人们通过计算机手段分析实际电力系统的各种运行情况，从而有效的了解电力系统概况。本文根据电力系统的特点，利用MATLAB的动态仿真软件Simulink搭建了无穷大电源的系统仿真模型，得到了在该系统主供电线路电源端发生三相短路接地故障并由故障器自动跳闸隔离故障的仿真结果，并分析了这一暂态过程。通过仿真结果说明MATLAB 电力系统工具箱是分析电力系统的有效工具。 关键词：电力系统；三相短路；故障分析；MATLAB仿真 目录 一．前言.............................................. 二．无穷大功率电源供电系统仿真模型构建............... 1.总电路图的设计......................................
基于matlab的电路仿真 杨泽辉51130215 %基于matlab的电路仿真 %关键词: RC电路仿真, matlab, GUI设计 % 基于matlab的电路仿真 %功能：产生根据输入波形与电路的选择产生输出波形 close all;clear;clc; %清空 figure('position',[189 89 714 485]); %创建图形窗口，坐标（189，89），宽714，高485；Na=['输入波形[请选择]|输入波形:正弦波|',... '输入波形:方形波|输入波形:脉冲波'];%波形选择名称数组； Ns={'sin','square','pulse'}; %波形选择名称数组； R=2; % default parameters: resistance 电阻值 C=2; % default parameters: capacitance电容值 f=10; % default parameters: frequency 波形频率 TAU=R*C; tff=10; % length of time ts=1/f; % sampling length sys1=tf([1],[1,1]); % systems for integral circuit %传递函数； sys2=tf([1,0],[1,1]); % systems for differential circuit a1=axes('position',[0.1,0.6,0.3,0.3]); %创建坐标轴并获得句柄； po1=uicontrol(gcf,'style','popupmenu',... %在第一个界面的上方创建一个下拉菜单'unit','normalized','position',[0.15,0.9,0.2,0.08],... %位置 'string',Na,'fontsize',12,'callback',[]); %弹出菜单上的字符为数组Na，字体大小为12， set(po1,'callback',['KK=get(po1,''Value'');if KK>1;',... 'st=char(Ns(KK-1));[U,T]=gensig(st,R*C,tff,1/f);',... 'axes(a1);plot(T,U);ylim([min(U)-0.5,max(U)+0.5]);',... 'end;']); %pol触发事件：KK获取激发位置，st为当前触发位置的字符串，即所选择的波形类型； %[U,T]，gensing,产生信号，类型为st的值，周期为R*C，持续时间为tff， %采样周期为1/f，U为所产生的信号，T为时间； %创建坐标轴al；以T为x轴，U为y轴画波形，y轴范围。。。 Ma=['电路类型[请选择]|电路类型:积分型|电路类型:微分型']; %窗口2电路类型的选择数组； a2=axes('position',[0.5,0.6,0.3,0.3]);box on; %创建坐标轴2； set(gca,'xtick',[]);set(gca,'ytick',[]); %去掉坐标轴的刻度 po2=uicontrol(gcf,'style','popupmenu',... %在第二个窗口的位置创建一个下拉菜单，同1 'unit','normalized','position',[0.55,0.9,0.2,0.08],... 'string',Ma,'fontsize',12,'callback',[]); set(po2,'callback',['KQ=get(po2,''Value'');axes(a2);',... %po2属性设置，KQ为选择的电路类型，'if KQ==1;cla;elseif KQ==2;',... %1则清除坐标轴，2画积分电路，3画微分电路 'plot(0.14+0.8i+0.02*exp(i*[0:.02:8]),''k'');hold on;',... 'plot(0.14+0.2i+0.02*exp(i*[0:.02:8]),''k'');',... 'plot(0.84+0.2i+0.02*exp(i*[0:.02:8]),''k'');',... 'plot(0.84+0.8i+0.02*exp(i*[0:.02:8]),''k'');',... 'plot([0.16,0.82],[0.2,0.2],''k'');',... 'plot([0.16,0.3],[0.8,0.8],''k'');',... 'plot([3,4,4,3,3]/10,[76,76,84,84,76]/100,''k'');',... 'plot([0.4,0.82],[0.8,0.8],''k'');',... 'plot([0.6,0.6],[0.8,0.53],''k'');',... 'plot([0.6,0.6],[0.2,0.48],''k'');',... 'plot([0.55,0.65],[0.53,0.53],''k'');',... 'plot([0.55,0.65],[0.48,0.48],''k'');',... 'text(0.33,0.7,''R'');',...
信号与线性系统 课程设计 题目 学号姓名学号姓名学号姓名学号姓名 院系 年级 专业 日期
摘要 的方波信号进行傅里叶级数展开，并绘制离散幅度利用MATLAB对周期为T 谱和不同次谐波叠加后的图形。通过观察绘制的各个图像，加深对傅立叶变换和信号的分解与合成的理解。 Abstract Expanded the square wave signal with periodic T0 to Fourier series by MATLAB , and drew the discrete spectrum and plot the patterns after different sub harmonics are superimposed. Through the observation of each image, deepen the understanding of the Fourier transform and signal decomposition and synthesis. 关键词：矩形信号傅里叶级数谐波叠加分解与合成 Keywords: Squarewave signal.Fourier series. Harmonic superposition. Decomposition and synthesis 一、设计目的和要求 本设计主要利用MATLAB绘制信号的离散幅度谱和各次谐波叠加后的波形，通过观察谐波展开次数增加后的波形，进一步掌握信号分解与合成的原理。 培养运用所学知识分析解决问题的能力。 掌握用MATLAB实现通信系统仿真实验的能力。 这里要做一个信号的分解与合成的仿真系统，利用matlab软件的仿真模拟能力来体现信号的分解与合成过程中出现的情况。
收稿日期:2006-08-25 作者简介:曾江华, 女, 长江水利委员会设计院机电处, 工程师, 硕士。 文章编号:1001-4179(2006 11-0041-02 MAT LAB 在电力系统仿真中的运用 曾江华陈晓明金伟江万里李远青 (长江水利委员会设计院, 湖北武汉 摘要:MAT LAB 是将计算、可视化、真中运用很广泛。, 由于电力系统是个复杂的系统, , 也不直观。M AT LAB 的M LAB 的POWERSY STE M BLOCK 对避雷器在有电抗器补, 。关键; ; 仿真; 运用T 文献标识码:A 1概述 M AT LAB 是由美国Mathw orks 公司开发的大型软件, 它是以 矩阵运算为基础, 把计算、可视化、程序设计融合在一个交互的工作环境中, 在此环境中可以实现工程计算、算法研究、建模和仿真、应用程序开发等。在M AT LAB 中包括了两大部分, 数学计算和工程仿真, 其中在工程仿真方面,M AT LAB 提供的软件支持涉及到各个工程领域, 并且在不断完善。M AT LAB 所具有的程序设计灵活, 直观, 图形功能强大的优点使其已经发展成为多学科, 多平台的强大的大型软件。M AT LAB 提供的S imulink 工具箱是一个在M AT LAB 环境下用于对动态系统进行建模、仿真和分析的软件包, 它提供了用方框图进行建模的接口, 与传统的仿真建模相比, 更加直观、灵活。S imulink 的作用是在程序块间的互联基础上建立
起一个系统。每个程序块由输入向量, 输出向量以及表示状态变量的向量等3个要素组成。在计算前, 需要初始化并赋初值, 程序块按照需要更新的次序分类, 然后用ODE 计算程序通过数值积分来模拟系统。M AT LAN 含有大量的ODE 计算程序, 有固定步长的, 有可变步长的, 为求解复杂的系统提供了方便。 M AT LAB 在电力系统建模和仿真的应用主要由电力系统仿真模块(P ower System Blockset 来完成, P ower System Block 是由TE QSI M 公司和魁北克水电站开发的。PS B 是在S imulink 环境下 使用的模块, 采用变步长积分法, 可以对非线性、刚性和非连续系统进行精确的仿真, 并精确地检测出断点和开关发生时刻, PS B 程序库含有代表电力网络中一般部件和设备的S imulink 程 序块, 通过PS B 可以迅速建立模型, 并立即仿真。PS B 程序块程序库中的测量程序和控制源起到电信号与S imulink 程序之间连 接作用。2电力系统元件库 电力系统元件库包括了电路、电力电子、电机和电力系统等 常用的基本元件和系统的仿真模型。其包含以下库元件: (1 电源元件。包括了交流电压源和电流源、直流电压源、可控电源及三相电源等产生电信号的元件。 (2 线路元件。包括各种线性网络电路元件和非线性网络电路元件。 (3 电力电子元件。包括如二级管、晶闸管等各种电力电子元件。 (4 电机元件。包括各种电机模型元件。 (5 连接器元件。包含有在各种不同情况下用于相互连接的元件。
《电力系统设计》报告 题目: 基于MATLAB的电力系统仿学院：电子信息与电气工程学院 班级: 13级电气 1 班 姓名:田震 学号: 20131090124 日期：2015年12月6日
基于MATLAB的电力系统仿真 摘要：目前，随着科学技术的发展和电能需求量的日益增长，电力系统规模越来越庞大，超高压远距离输电、大容量发电机组、各种新型控制装置得到了广泛的应用，这对于合理利用能源，充分挖掘现有的输电潜力和保护环境都有重要意义。另一方面，随着国民经济的高速发展，以城市为中心的区域性用电增长越来越快，大电网负荷中心的用电容量越来越大，长距离重负荷输电的情况日益普遍，电力系统在人们的生活和工作中担任重要角色，电力系统的稳定运行直接影响着人们的日常生活。从技术和安全上考虑直接进行电力试验可能性很小，因此迫切要求运用电力仿真来解决这些问题。 电力系统仿真是将电力系统的模型化、数学化来模拟实际的电力系统的运行，可以帮助人们通过计算机手段分析实际电力系统的各种运行情况，从而有效的了解电力系统概况。本文根据电力系统的特点，利用MATLAB的动态仿真软件Simulink搭建了无穷大电源的系统仿真模型，得到了在该系统主供电线路电源端发生三相短路接地故障并由故障器自动跳闸隔离故障的仿真结果，并分析了这一暂态过程。通过仿真结果说明MATLAB电力系统工具箱是分析电力系统的有效工具。 关键词：电力系统；三相短路；故障分析；MATLAB仿真
目录 一．前言 (4) 二．无穷大功率电源供电系统仿真模型构建 (5) 1.总电路图的设计 (5) 2.各个元件的参数设定 (6) 2.1供电模块的参数设定 (6) 2.2变压器模块的参数设置 (6) 2.3输电线路模块的参数设置 (7) 2.4三相电压电流测量模块 (8) 2.5三相线路故障模块参数设置 (8) 2.6三相并联RLC负荷模块参数设置 (9) 3.仿真结果 (9)
1.目前常用的电力系统仿真软件有：BPA 程序和EMTP(程序；PSCAD /EMTDC; NETOMAC; PSASP;MATLAB 2.SimPowerSystems库产品 SimPowerSystems 4.0中含有130 多个模块，分布在7个可用子库中。这7个子库分别为“应用子库(Application Libraries)”、“电源子库(Electrical Sources)”、“元件子库(Elements)”、“附加子库(Extra Library)”、“电机子库(Machines)”、“测量子库(Measure-ments)”和“电力电子子库(Power Electronics)”。此外，SimPowerSystems 4.0中还含有一个功能强大的图形用户分析工具Powergui和一个废弃的“相量子库”(Phasor Elements) 3.MATLAB的特点:(1) 提供了便利的开发环境。(2) 提供了强大的数学应用功能。(3) 编 程语言简易高效。(4) 图形功能强大(5) 提供了功能强大的工具箱。(6) 应用程序接口功能强大。(7) MATLAB的缺点。和其它高级程序相比，MATLAB程序的执行速度较慢。 4.SIMULINK的特点:(1) 建立动态系统的模型并进行仿真。(2) 以直观的方式建模。(3) 增 添定制模块元件和用户代码。(4) 快速、准确地进行设计模拟。(5) 分层次地表达复杂系统。(6) 交互式的仿真分析。 5.SimPowerSystems库的特点:(1) 使用标准电气符号进行电力系统的拓扑图形建模和仿 真。(2) 标准的AC和DC电机模型模块、变压器、输电线路、信号和脉冲发生器、HVDC 控制、IGBT 模块和大量设备模型。(3) 使用SIMULINK强有力的变步长积分器和零点穿越检测功能，给出高度精确的电力系统仿真计算结果。(4) 利用定步长梯形积分算法进行离散仿真计算，为快速仿真和实时仿真提供模型离散化方法(5) 利用Powergui交互式工具模块可以修改模型的初始状态，从任何起始条件开始进行仿真分析(6) 提供了扩展的电力系统设备模块，如电力机械、功率电子元件、控制测量模块和三相元器件。 (7) 提供大量功能演示模型，可直接运行仿真或进行案例学习。 6.默认的程序主界面主要包括下列区域：①菜单；②工具栏；③命令窗口；④当前 路径浏览器；⑤工作空间浏览器；⑥命令历史浏览器。 7.菜单功能:(1) [File>New>M-File]：进入文本编辑窗界面，建立一个文本文件，实现 MATLAB命令文件的输入、编辑、调试、保存等处理功能，保存时文件后缀名为.m。(2) [File>New>Figure]：进入图形窗界面，建立一个图形文件，实现MATLAB图形文件的显示、编辑、保存等处理功能，保存时文件名后缀为.fig.(3) [File>New>Model]：建立一个SIMULINK模型文件，实现SIMULINK仿真模型的建模、仿真、调试、保存等处理功能，保存时文件名后缀为.mdl。 8.：进入SIMULINK仿真环境界面，作用相当于在MATLAB的命令窗口中输入 simulink命令并按回车键。 9.Matlab默认工作路径：安装路径\Matlab\work .修改路径（1）利用图标 （2）利用菜单项[File>Set Path>Add Folder]将用户拟采用的 目录添加到Matlab 搜索路径中。 10.MATLAB编程有两种工作方式：一种称为行命令方式，就是在工作窗口中一行一行地输 入程序，计算机每次对一行命令做出反应，因此也称为交互式的指令行操作方式；另一种工作方式为M文件编程工作方式。 11.变量是保存数据信息的一种最基本的数据类型。变量的命名应遵循如下规则：(1) 变量 名必须以字母开头；(2) 变量名可以由字母、数字和下划线混合组成；(3) 变量名区分字母大小写；(4) MATLAB保留了一些具有特定意义的默认变量(见表2-3)，用户编程时可以直接使用，并尽量避免另外自定义例如，Long和My_long1均是有效的变量名，Long
一. 生成一个基频为20Hz ，幅度为±1的周期方波信号； （1）试确定该信号的周期？ （2）画出该信号5个周期内的波形；画出该信号的幅度频谱图； （3）该信号谱线间隔和频带宽度是多少，对它进行时域采样，采样频率Fs 至少为多少？ （4）设计一个IIR/FIR 低通滤波器，滤去该周期信号中80Hz 以后的频率成分，要求通带截止频率80Hz ，阻带截止频率90Hz ，通带纹波小于3dB ，阻带衰减大于50dB ,用matlab 命令画出滤波器的频率响应，以及滤波后的信号波形和频谱。 答：1.周期为0.05秒； 3.谱线间隔为20Hz, 带宽为40Hz ，采样频率至少为80Hz 。 Fs=1000; N=256; % 时域/频率采样点 t=0:1/Fs:(N-1)/Fs; x=square(2*pi*20*t); %产生幅度为[-1,1]，频率为20Hz 的方波 subplot(2,1,1) %方波绘图 plot(t,x); title('5个周期的矩形波') axis([0,0.1,-1.5,1.5]) T=0.05; %求方波的频谱 t1=0:1/Fs:(T-1/Fs); xn=square(2*pi*20*t1); Xk=fft(xn,N); Xk=fftshift(abs(Xk)); Xkl=Xk/N; subplot(2,1,2) plot(Fs/N*[0:1:N-1],Xkl) title('方波信号的频谱') [n,Wn]=buttord(80/500,90/500,3,50); %确定butterworth 滤波器的阶次 butterworth 低通滤波 %器：通带截止频率80Hz ，阻带截止频率90Hz ，通带纹波小于3dB ，阻带衰减大于50dB , %以上参数可以变化 [b,a]=butter(n,Wn); %确定滤波器传递函数 y=filter(b,a,x); %滤波 figure %滤波器频率特性 freqs(b,a); title('低通滤波器频率特性') ; figure subplot(2,1,1) plot(t,y) ;title('矩形波通过低通滤波器后的波形') ;axis([0,0.1,-1.5,1.5]) subplot(2,1,2) N=512; %可适当增减 fy=fftshift(fft(y,N)); %矩形波通过低通滤波器后的频谱 f=(-N/2:(N/2-1))*2; f=f*Fs/N; plot(f,abs(fy)/N) ;title('矩形波通过低通滤波器后的频谱') 二.设有连续信号()cos(240)f t t π=??， （1）利用DFT 计算该连续时间信号的频谱，采样频率至少应取多少？80HZ （2）以采样频率Hz F s 1000=对信号)(t f 采样得)(t f a ，请画出两个周期内的离散 信号)(t f a 的波形；[要求横
考核人数______ 考核班次_______________ 任课教员_________ 出题教员签名________ 任课教研室主任签名_______日期_______ 队别__________ 教学班次___________ 学号___________ 姓名____________ …………………………密………………………………封………………………………线……………………………………… 通信系统仿真题目 1．学习电路时已知LC 谐振电路具有选择频率的作用，当输入正弦信号频率与LC 电路的谐 振频率一致时，将产生较强的输出响应，而当输入信号频率适当偏离时，输出响应相对值很弱，几乎为零（相当于窄带通滤波器）。利用这一原理可以从非正弦周期信号中选择所需的正弦频率成分。题图所示RLC 并联电路和电流1()i t 都是理想模型。已知电路的谐振频率为 0100f kHz = =,100R k =Ω谐振电路品质因素Q 足够高（可滤除邻近频率成分） 1()i t 为周期矩形波，幅度为1 mA 当1()i t 的参数(,)T τ为下列情况时，粗略地画出输出电压 2()t υ的波形，并注明幅度值。 （1）510s T s τμμ== （2）1020s T s τμμ== （3）1530s T s τμμ== 2．设()x n 为一限长序列，当0n <和n N ≥时，()0x n =，且N 等于偶数。已知[()]DFT x n = ()X k ，试用()X k 表示以下各序列的DFT 。 (1)1()(1)x n x N n =-- (2)2()(1)()n x n x n =- (3) 3() (01)()()(21)0()x n n N x n x n N N n N n ≤≤-?? =-≤≤-??? 为其他值 (4) 4()()(01) ()2 2 () N N x n x n n x n n ?≠+≤≤ -? =???为其他值 (5) 5()(01)()0 (21)0() x n n N x n N n N n ≤≤-?? =≤≤-??? 为其他值 (6) 6() ()20()n x n x n n ??? ? ?=????? 为偶数为奇数 (DFT 有限长度取2N ，k 取偶数。) (7) 7()(2)x n x n =(DFT 有限长度取 2 N )。 3．已知三角脉冲1() f t 的傅里叶变换为21()24E F Sa τωτω??= ??? 试利用有关定理求210()cos()2f t f t t τω?? =- ??? 的傅里叶变换2()F ω。1()f t 、2()f t 的波形如下图所示。 4．求下图所示半波余弦信号的傅里叶级数。若E=10V ，f=10kHz ，大致画出幅度谱。 5．求下图所示()F ω的傅里叶逆变换()f t 。
课程设计 课程名称：电力系统分析 设计题目：基于Matlab计算程序的电力系统运行分析学院：电力工程学院 专业：电气工程自动化 年级： 学生姓名： 指导教师： 日期： 教务处制
目录 前言 (1) 第一章参数计算 (2) 一、目标电网接线图 (2) 二、电网模型的建立 (3) 第二章潮流计算 (6) 一．系统参数的设置 (6) 二．程序的调试 (7) 三、对运行结果的分析 (13) 第三章短路故障的分析计算 (15) 一、三相短路 (15) 二、不对称短路 (16) 三、由上面表对运行结果的分析及在短路中的一些问题 (21) 心得体会 (26) 参考文献 (27)
前言 电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。 在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。产生短路故障的主要原因是电力设备绝缘损坏。短路故障分为三相短路、两相短路、单相接地短路及两相接地短路。其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。短路故障大多数发生在架空输电线路。电力系统设计与运行时，要采取适当的措施降低短路故障的发生概率。短路计算可以为设备的选择提供原始数据。
MATLAB谐波信号双谱估计 clc clear close all N=128*64; n=[0:N-1]; Nfft=128; Fs=1; t=n/Fs; f1=0.6381/2/pi; f2=0.8345/2/pi; f3=f1+f2; f4=0.4909/2/pi; f5=1.7671/2/pi; f6=f4+f5; fai1=pi/6; fai2=pi/3; fai3=fai1+fai2; fai4=5*pi/4; fai5=6*pi/4; fai6=pi; s1=cos(2*pi*f1*t+fai1); s2=cos(2*pi*f2*t+fai2); s3=cos(2*pi*f3*t+fai3); s4=cos(2*pi*f4*t+fai4); s5=cos(2*pi*f5*t+fai5); s6=cos(2*pi*f6*t+fai6); s=s1+s2+s3+s4+s5+s6; %s=awgn(s,0,'measured'); noi=zeros(1,N); %高斯白噪声 noi=randn(1,N); s=s+noi; %有色噪声 nys=0; if(nys==1) Ln=length(s); Nd=[1-0.50.70.1]; Nc=[10.50.2]; nd=length(Nd)-1;
nc=length(Nc)-1; xik=zeros(nc,1); ek=zeros(nd,1); xi=randn(Ln,1); for kn=1:Ln e(kn)=-Nd(2:nd+1)*ek+Nc*[xi(kn);xik]; for nn=nd:-1:2 ek(nn)=ek(nn-1); end ek(1)=e(kn); for nn=nc:-1:2 xik(nn)=xik(nn-1); end xik(1)=xi(kn); end s=s+0.25*e; noi=0.25*e; end %信噪比 S_N=10*log10(sum(abs(s).^2)/sum(abs(noi).^2)); [bspec,waxis]=bispecd(s,Nfft,1,64,32); close all %双谱信号（包含相干信号的幅度和相位信息，幅度受信号的功率谱影响）SS=Nfft/2+1; waxis1=waxis(SS:end)*Fs; figure subplot(121) bspec1=bspec(SS:end,SS:end); contour(waxis1,waxis1,abs(bspec1),4); grid on title('双谱估计二维等高线图') xlabel('f1/Hz'),ylabel('f2/Hz') [X1,Y1]=meshgrid(waxis1,waxis1); subplot(122) mesh(X1,Y1,abs(bspec1)) axis tight xlabel('f1/Hz') ylabel('f2/Hz') zlabel('幅度') title('双谱估计三维图');