MATLAB应用技术5MATLAB与电力系统仿真5.1电力系统的数学模型电力系统一般由发电机、变压器、电力线路和电力负荷构成。电力系统的数学模型一般是由电力系统元件的数学模型组合构成。MATLAB为电力系统的建模提供了简洁的工具,通过电力系统的电路图绘制,可以自动生成数学模型。电路图模型的主要特点是具有良好的人机界面,便于进行简单的操作,省去了利用程序建立电力系统模型的反覆步骤。利用这种方式构成的数学模型相对于控制系统中的微分方程模型、状态方程模型、传递函数模型有着更直观和实用的优点。另外,在电路图模型建立以后,在MATLAB软件中,提供了power2sys函数作为短路模型的结构分析函数,可以利用power2sys函数将电力系统的电路图模型向状态方程模型和传递函数模型进行转换。MATLAB应用技术5.1.1电力系统元件库1.启动和退出电力系统元件库启动电力系统元件库的方法有几种,下面介绍两种最简单的方法。(1)利用指令窗口(CommandWindows)启动:在指令窗口中键入powerlib单击回车,则MATLAB软件中弹出电力系统元件对话框(powerlib)(2)利用开始(Start)导航区启动:单击开始按钮,选择仿真(Simulink)命令,再选择电力系统仿真命令(SimPowerSystem),在弹出的对话框中选择电力系统元件库(BlockLibrary)命令即可MATLAB应用技术MATLAB应用技术2.电力系统元件库简介在电力系统元件库对话框中包含了10类库元件,分别是电源元件(ElectricalSources)、线路元件(Elements)、电力电子元件(PowerElectronics)、电机元件(Machines)、连接器元件(Connectors)、电路测量仪器(Measurements)、附加元件(Extras)、演示教程(Demos)、电力图形用户接口(Powergui)、电力系统元件库模型(Powelib_models)。1)电源元件电源元件库包括7类元件,分别为:直流电压源元件(DCVoltageSource),交流电压源元件(ACVoltageSource),交流电流源元件(ACCurrentSource),受控电压源元件(ControlledVoltageSource),受控电流源元件(ControlledCurrentSource),三相电源元件(3-PhaseSource),三相可编程电压源元件(3-PhaseProgrammableVoltageSource)。MATLAB应用技术(1)直流电压源元件(DCVoltageSource)直流电压源元件在电力系统中可以用来实现一个直流的电压源,如操作电源等。MATLAB软件提供的直流电源为理想的直流电压源。(2)交流电压源元件(ACVoltageSource)交流电压源可以用来实现理想的单相正弦交流电压。(3)交流电流源元件(ACCurrentSource)MATLAB软件提供的交流电流源为一理想电流源(4)受控电压源元件(ControlledVoltageSource)MATLAB软件提供的受控电压源是由激励信号源控制的,激励源可以是交流激励源也可以是直流激励源。(5)受控电流源元件(ControlledCurrentSource)MATLAB应用技术(6)三相电源元件(3-PhaseSource)三相电源元件是电力系统设计中最常见的电路元件,也是最重要的元件,其运行特性对电力系统的运行状态起到决定性的作用。三相电源元件提供了带有串联RL支路的三相电源。(7)三相可编程电压源元件(3-PhaseProgrammableVoltageSource)三相可编程电压源是可以对其进行编程的三相电压源,它的幅值、相位、频率、谐波均可随时间进行变化,应用非常灵活。其主要作用是提供一个幅值、相位、频率、基频分量进行实时变性编程的三相电压源;此外,还可以提供两个谐波分量,作用于基频信号。MATLAB应用技术2)线路元件线路元件库包括各种线性网络电路元件和非线性网络电路元件,线路元件共有4类,分别是:支路元件(Elements)、输配电线路元件(Lines)、断路器元件(CircuitBreakers)、变压器元件(Transformers)。(1)支路元件(Elements)支路元件用来实现各种串并联支路或者负载元件,它包括12种元件(2)输配电线路元件(Lines)在电力系统设计和分析中,输配电线路一般用各种类型的等值电路来进行简化以便于简化分析。输配电线路元件的作用就是构成各种线路的等值电路,在输配电线路元件中包括3种元件MATLAB应用技术(3)断路器元件(CircuitBreakers)在电力系统中,断路器的作用是通断高压电力线路,可靠地接通或切断有载电路和故障电路。断路器元件就是用来实现各种电路中的高压断路器。在断路器元件中包括3种元件(4)变压器元件(Transformers)在电力系统中,电力变压器是最重要的电气设备,其作用是进行能量的传输并改变电压的等级。变压器的种类有很多种,变压器元件就是用来设计实现各种类型的变压器。在变压器元件种包括6种元件3)其他元件在电力系统元件库中还有其他元件:电力电子元件(PowerElectronics)、电机元件(Machines)、连接器元件(Connectors)、电路测量仪器(Measurements)、附加元件(Extras),这些元件都具有特定的功能MATLAB应用技术3.示例下面以几个简单的例子来介绍如何使用这些电气元件。例1:交流电压源的叠加设计的交流电路如图5-20所示,在此电路图中,交流电压源的幅值、频率、相位均不相同,可以通过仿真结果直接对各自电压源的输出和他们的叠加结果进行分析,这种分析方法简单、直接。图5-20交流电压源的叠加电路图电路图设计步骤(1)从电源元件库选择交流电压源元件,复制后粘贴在电路图中。步骤1:将电压源元件改名为U1步骤2:双击交流电压源元件,对交流电压源元件的参数进行如下设置:图5-20交流电压源的叠加电路图MATLAB应用技术峰值振幅(PeakAmplitude):100初始相位(Phase):30频率(Frequency):60采样时间(Sampletime):0测量选项(Measurements):选择不测量电气量步骤1:复制交流电压源元件并改名为U2步骤2:双击交流电压源元件,对交流电压源元件的参数进行如下设置:峰值振幅(PeakAmplitude):75初始相位(Phase):60频率(Frequency):50采样时间(Sampletime):0测量选项(Measurements):选择不测量电气量)6120sin(1001tu)3100sin(752tuMATLAB应用技术(2)从线路元件库中选择串联RLC支路对串联RLC支路元件的参数进行如下设置:电阻(ResistanceR)200电感(InductanceL)100e-3电容(CapacitanceC)150e-6测量选项(Measurements):选择不测量电气量(3)从电路测量仪器中选择电压计元件,复制后粘贴于电路图中。(4)在仿真元件库(SimulinkLiberary)中选择示波器,复制示波器并改变其名称为1。(5)从连接元件库(Connectors)中选择接地及相应的元件进行合理的放置。对该电路图进行接线,就可以完成电路图的绘制。在接线时,如果提示颜色为红色,则表示在接线时出现了错误。MATLAB应用技术仿真参数设置在电路图菜单选项中,选择仿真(Simulation)菜单,激活仿真参数(SimulationParameters)命令,即可弹出仿真参数对话框,根据相应选项对其进行设置:开始时间(Starttime):0停止时间(Stoptime):0.4求解程序类型(Type)选项:可变步长(Variable),Ode45(Domand-Price)最大步长(Maxstepsize)选项:自动(Auto)最小步长(Minstepsize)选项:自动(Auto)初始步长(Initialstepsize)选项:自动(Auto)相对容差(Relativetolerance)选项:1e-3绝对容差(Absolutetolerance)选项:1e-6MATLAB应用技术仿真结果及分析合理设置示波器参数后,激活仿真按钮,得到仿真结果如图所示。示波器1输出的电压波形为交流电压源U1和U2的叠加,横轴为时间轴,纵轴为电压幅值。从仿真结果可见,在交流电路中,多个交流电压源共同作用的结果等效于一个非线性电压源。图例1仿真结果MATLAB应用技术5.1.2电力系统电路图模型结构分析利用电力系统工具箱建立电路图模型,操作简单,熟悉电路元件的人员可以很容易地掌握建立电力系统数学模型的方法,避免了利用程序建模的复杂步骤。根据上一小节示例的建模方法可以很容易建立起电力系统的电路图模型。在MATLAB软件中,提供了一种对电路图进行分析的方法,这就是power2sys函数。利用该函数,可以对电路图的结构特征、状态方程等进行较为全面的分析。power2sys函数的表达式如下:psb=power2sys('sys','structure')用来显示电路图的结构psb=power2sys(‘sys’,‘sort’)用来显示电路图中元件和支路的信息psb=power2sys(‘sys’,‘ss’)将电路图模型转换为状态方程[A,B,C,D,x0,states,inputs,outputs,uss,xss,yss,freqyss,Hlin]=psb=power2sys('sys')用来显示电力系统模型的结构信息psb=power2sys('sys','net')用来显示电力系统的网络结构MATLAB应用技术5.1.3Park变换同步电机是电力系统中的重要元件,它实质上是由定子和转子两个部件组成。在研究同步机的数学模型时,我们都假设定子三相绕组的结构完全相同,空间位置彼此相差120度,转子铁芯及绕组对极中心轴和极间轴完全对称。一般,在推导同步机的数学模型时应用的是用abc坐标系统表示的电压和磁链方程。abc三轴就是定子三相绕组的中心轴线。定子三相绕组中的电流分别表示如下:)120cos()120cos(cos)cos(0mcmbmsmaIiIiItIiMATLAB应用技术利用该坐标系统建立同步机的电压和磁链方程时非常容易理解,但是所建立的方程为变系数的微分方程,它们的求解非常的困难。为了克服这个困难,最简单有效的方法时将定子abc三相绕组的磁链和电压方程用一组新的变量替换,这样使方程更易于求解。变量变换又称为坐标变换,最常用的坐标变换,即Park变换。这样变换后电流的表示方式如下:cqcdcbqbdbaqadaqdiiiiiiiiiiiiMATLAB应用技术(1)abc坐标系统变换为dq0坐标系统的变换公式如下:在MATLAB中,使用abc坐标系统转换为dq0坐标系统(abc_to_dq0Transformation)元件可以实现这种变换。abc_to_dq0Transformation在电力系统元件库(PowerLib)中的附加元件(Extras)下的测量元件(Measurements)中。cbaqdiiiiii212121)32sin()32sin(sin)32cos()32cos(cos320MATLAB应用技术MATLAB应用技术(2)dq0坐标系统变换为abc坐标系统的变换公式如下:01)32sin()32cos(1)32si