电力系统的仿真与应用教材.pptx

1、第8章 定 制 模 块 第8章 定 制 模 块 8.1 定制非线性模块定制非线性模块 8.2 S函数编写及应用函数编写及应用 电力系统的仿真与应用教材第1页第8章 定 制 模 块 8.1 定制非线性模块定制非线性模块8.1.1 非线性电感模块非线性电感模块【例【例8.1】定制一个非线性电感元件，当电压在0120 V时，电感恒定为2 H；当电压超出120 V时，电感元件饱和，电感降低到0.5 H。图8-1所表示为该非线性电感对应磁通电流特征曲线，单位为p.u.。其中VB=120/sqrt(2)V，fB=50 Hz。电力系统的仿真与应用教材第2页第8章 定 制 模 块 图8-1 例8.1磁通电流特

2、征曲线电力系统的仿真与应用教材第3页第8章 定 制 模 块 解：解：(1)理论分析。显然电感元件上电压v和电流i含有以下关系：其中，为电感元件上自感磁链。由式(8-1)能够得到磁链为所以电感上电流i为(8-3)(8-2)(8-1)电力系统的仿真与应用教材第4页第8章 定 制 模 块 可见，能够用受控电流源表示该非线性电感元件，该电流源受控于电流源两端电压。(2)按图8-2搭建非线性电感模型。该模型包含一个电压表模块、一个可控电流源模块(电流源电流方向为箭头所表示方向)、一个积分模块和一个用于描述磁通电流饱和特征查表模块。选取各模块名称及提取路径见表8-1。图中有一个信号输出口m，输出非线性电感

3、模块上磁通和该模块两端电压值。电力系统的仿真与应用教材第5页第8章 定 制 模 块 图8-2 非线性电感模型电力系统的仿真与应用教材第6页第8章 定 制 模 块 表表8-1 例例8.1非线性电感模型中包含模块非线性电感模型中包含模块电力系统的仿真与应用教材第7页第8章 定 制 模 块 图8-3 例8.1查表模块参数设置打开查表模块参数对话框，按图8-3设置参数。该图中参数实际上就是图8-1磁通电流特征。其余模块参数采取默认设置。电力系统的仿真与应用教材第8页第8章 定 制 模 块(3)将搭建好非线性电感模型组合为一个子系统并命名为Nonlinear Inductance后，按图8-4所表示搭建

4、仿真系统。选取各模块名称及提取路径见表8-2。电力系统的仿真与应用教材第9页第8章 定 制 模 块 图8-4 例8.1仿真系统图电力系统的仿真与应用教材第10页第8章 定 制 模 块 表表8-2 例例8.1仿真电路模块名称及提取路径仿真电路模块名称及提取路径电力系统的仿真与应用教材第11页第8章 定 制 模 块 该系统中含有两个电压源，一个峰值为120 V、50 Hz、相角为90交流电压源Vs和一个幅值为0直流电压源VDC。串联RLC支路为纯电阻电路，其中电阻元件R=5 W。在仿真参数对话框中设置变步长ode23tb算法，仿真结束时间为1.5 s。(4)开始仿真。设置直流电压源幅值为0，开始仿

5、真。图8-5所表示为仿真最终5个周期波形，图中波形从上到下依次为非线性电感元件上磁通、电流和电压。此时，电压未超出极限值120 V，电感为2 H。对应电流幅值Im为(8-4)电力系统的仿真与应用教材第12页第8章 定 制 模 块 磁链为与观察到波形一致。(8-5)电力系统的仿真与应用教材第13页第8章 定 制 模 块 图8-5 仿真波形图(VDC=0 V)电力系统的仿真与应用教材第14页第8章 定 制 模 块 将直流电压源幅值改为1 V，再次仿真。按理论分析，直流电压源单独作用时，电路中电感相当于短路，所以观察仿真电流中应该含有一个直流电流，该电流值为1 V/5=0.2 A。观察仿真最终5个周

6、期波形如图8-6所表示，图中波形从上到下依次为非线性电感元件上磁通、电流和电压。因为对1 V直流电压源分量进行了积分运算，造成磁通饱和，所以电流波形发生畸变。由图可见，此时电流增大到0.575 A。电力系统的仿真与应用教材第15页第8章 定 制 模 块 图8-6 仿真波形图(VDC=1 V)电力系统的仿真与应用教材第16页第8章 定 制 模 块 经过FFT模块提取电流信号中基频和直流分量，如图8-7所表示，图中波形为非线性电感元件电流基频分量和直流分量。可见，基频电流增大到0.27 A，直流电流分量为0.2 A，与理论分析值一致。电力系统的仿真与应用教材第17页第8章 定 制 模 块 图8-7

7、 基频和直流电流分量电力系统的仿真与应用教材第18页第8章 定 制 模 块 尽管在例8.1中已经建立了一个非线性电感模块，不过该模块外观并不完美。下面将对该非线性电感模块进行封装，创建对话框、图标和说明文档，并经过对话框来设定非线性电感磁通电流特征，使它看起来和SIMULINK库提供其它模块一样完整。单击图8-4中“非线性电感元件”图标，打开菜单Edit Mask subsystem，弹出封装子系统编辑窗口如图8-8所表示。选择“文档”(Documentation)标签页，在“封装类型”(Mask type)文本框中输入文字“非线性电感”，在“封装说明”(Mask description)多行

8、文本框中输入该模块简单说明和注意事项，在“封装帮助”(Mask help)多行文本框中输入该模块帮助文件。电力系统的仿真与应用教材第19页第8章 定 制 模 块 图8-8 封装子系统编辑窗口(文档标签页)电力系统的仿真与应用教材第20页第8章 定 制 模 块 选择“参数”(Parameters)标签页如图8-9所表示，经过点击 按键添加、删除、移动项目。按图8-9分别添加额定电压、额定频率、线性电感和饱和特征。其中，在“变量说明”(Prompt)列中输入各变量简单说明，在“变量名”(Variable)列中输入各变量名称。注意这些变量是能够被封装子系统作为已知参数调用，所以，这些变量名应该是轻易

9、记忆，同时SIMULINK不区分大小写。在“类型”(Type)列中选择参数类型，可选类型有“文本框”(edit)、“列表框”(checkbox)和“下拉框”(popup)。“可计算”(Evaluate)和“可调用”(Tunable)列为可选项。电力系统的仿真与应用教材第21页第8章 定 制 模 块 选中“可计算”(Evaluate)列后，SIMULINK首先对用户输入表示式进行计算，然后再将计算结果赋值给变量，不然SIMULINK直接把用户输入表示式作为一个字符串赋值给变量。选中“可调用”(Tunable)列将允许该参数在仿真过程中被修改。电力系统的仿真与应用教材第22页第8章 定 制 模 块

10、 图8-9 封装子系统编辑窗口(参数标签页)电力系统的仿真与应用教材第23页第8章 定 制 模 块 选择“初始化”(Initialization)标签如图8-10所表示，在“初始化命令”(Initialization commands)窗口中输入以下命令并提取电流变量Current_vect和磁通变量Flux_vect。电力系统的仿真与应用教材第24页第8章 定 制 模 块 图8-10 封装子系统编辑窗口(初始化标签页)电力系统的仿真与应用教材第25页第8章 定 制 模 块 电力系统的仿真与应用教材第26页第8章 定 制 模 块 电力系统的仿真与应用教材第27页第8章 定 制 模 块 单击“确

11、定”(OK)按键，关闭封装子系统编辑窗口。接下来，将定义电流变量Current_vect和磁通变量Flux_vect传递到非线性电感元件查表模块中。经过菜单 Editlook under Mask进入图8-2所表示“非线性电感元件”窗口，打开查表模块对话框，设置输入、输出参数如图8-11所表示。确认后退出“非线性电感元件”子系统，回到主窗口中，双击“非线性电感元件”图标，出现图8-12所表示参数对话框，输入额定电压、额定频率、线性电感值和磁通电流饱和特征。现在能够开始仿真了，仿真波形和图8-6完全相同。电力系统的仿真与应用教材第28页第8章 定 制 模 块 图8-11 查表模块中参数设置电力系

12、统的仿真与应用教材第29页第8章 定 制 模 块 图8-12 非线性电感元件参数设置电力系统的仿真与应用教材第30页第8章 定 制 模 块 回到仿真主窗口，继续为该模块定制一个图标，使得该模块像SIMULINK任何一个模块一样漂亮。选中非线性电感模块，经过菜单EditEdit Mask打开封装子系统编辑窗口，选择“图标”(Icon)标签如图8-13所表示。电力系统的仿真与应用教材第31页第8章 定 制 模 块 图8-13 封装子系统编辑窗口(图标标签页)电力系统的仿真与应用教材第32页第8章 定 制 模 块 在“画图命令”(Drawing commands)窗口中输入命令plot(Curren

13、t_vect,Flux_vect);在“透明度”(Transparency)下拉框中选择“透明”(Transparent)。单击“确定”(OK)按键或者“应用”(Apply)按键后退出子系统编辑窗口，能够看见磁通电流饱和特征曲线出现在非线性电感模块上，因为选择透明处理，输入、输出端口名称也一并显示在该模块上。封装后非线性电感模块图标如图8-14所表示。当然，也能够选择“不透明”(Opaque),这么，输入、输出端口名称被图形覆盖了。电力系统的仿真与应用教材第33页第8章 定 制 模 块 图8-14 非线性电感模块图标电力系统的仿真与应用教材第34页第8章 定 制 模 块 8.1.2 定制非线性

14、电阻元件定制非线性电阻元件非线性电阻元件建模方法和非线性电感元件建模方法类似。本节不再利用查表方式建立电压电流关系，而是直接由电压电流数学关系建立表示式。【例【例8.2】定制一个非线性MOV电阻元件，搭建电路，观察效果。解：解：(1)理论分析。金属氧化物压敏电阻MOV电压电流含有以下关系：(8-6)电力系统的仿真与应用教材第35页第8章 定 制 模 块 其中，v、i为瞬时电压和电流；V0为钳制电压；I0为钳制电压对应参考电流；用来定义非线性特征，通常在10，50间取值。所以，本例能够用受控电流源来表示该非线性电阻元件，受控电流源受控于该电源两端电压。(2)按图8-15搭建非线性MOV电阻模型。

15、该模型包含一个电压表模块、一个可控电流源模块(电流源电流方向为箭头所表示方向)、一个传递函数模块和一个自定义函数模块。选取各模块名称及提取路径见表8-3。图中有一个信号输出口m，输出非线性电阻模块电压值。电力系统的仿真与应用教材第36页第8章 定 制 模 块 图8-15 非线性电阻模型电力系统的仿真与应用教材第37页第8章 定 制 模 块 表表8-3 例例8.2非线性电阻模型中包含模块名称及提取路径非线性电阻模型中包含模块名称及提取路径电力系统的仿真与应用教材第38页第8章 定 制 模 块 本例利用数学函数模块直接建立非线性MOV电压电流特征。因为纯电阻模块不含状态变量，这么在SIMULINK

16、内部运算时将产生一个代数循环，造成运算速度降低。为了解开代数环，在电压测量模块输出口加入了一个1阶滞后传递函数。按图8-16封装该非线性电阻模块，设置钳位电压为2倍额定电压，即2120e3sqrt(2)/sqrt(3)=195.96 kV。参考电流为500 A，为25。电力系统的仿真与应用教材第39页第8章 定 制 模 块 图8-16 封装非线性电阻模块 电力系统的仿真与应用教材第40页第8章 定 制 模 块 继续给该非线性电阻加一个图标。打开封装子系统编辑窗口，选择“初始化”标签，在“初始化命令”(Initialization commands)窗口中输入以下命令：t=0:0.0001:0.

17、04;x=sin(100*pi*t);y=500*x.25;选择“图标”(Icon)标签，在“画图命令”(Drawing commands)窗口中输入命令：plot(y,x);确定后退出该编辑窗口。非线性电阻模块图标如图8-17所表示。电力系统的仿真与应用教材第41页第8章 定 制 模 块 图8-17 非线性电阻模块图标电力系统的仿真与应用教材第42页第8章 定 制 模 块 图8-18 例8.2仿真系统图(3)将定制非线性电阻模块用于保护线电压为120 kV、按图8-18所表示搭建仿真系统。选取各模块名称及提取路径见表8-4。电力系统的仿真与应用教材第43页第8章 定 制 模 块 表表8-4

18、例例8.2仿真电路模块名称及提取路径仿真电路模块名称及提取路径电力系统的仿真与应用教材第44页第8章 定 制 模 块 将交流电压源峰值设为2.3120sqrt(2)/sqrt(3)kV、频率为50 Hz、相角为0。串联RLC支路RL中，电阻元件R=1.92，电感元件L=26 mH。串联RLC负荷额定电压为120/sqrt(3)kV，频率为50 Hz，有功功率为10 MW，感性和容性无功功率为0。在仿真参数对话框中设置变步长ode23tb算法，仿真结束时间为0.1 s。开始仿真。图8-19所表示为仿真波形图，图中波形为非线性电阻元件上电流和电压。电力系统的仿真与应用教材第45页第8章 定 制 模

19、 块 图8-19 仿真波形图电力系统的仿真与应用教材第46页第8章 定 制 模 块 由波形图可见，在正常电压条件下，MOV相当于一个阻值极大线性电阻。当电压大于200 kV时，MOV内阻急剧下降并快速导通，其上电流增加几个数量级，而电压被钳位在200 kV，从而有效地保护了电路中其它元器件不至于过压而损坏。对应电压电流特征如图8-20所表示。电力系统的仿真与应用教材第47页第8章 定 制 模 块 图8-20 实测V-I特征电力系统的仿真与应用教材第48页第8章 定 制 模 块 8.1.3 定制模块库定制模块库为了方便地对定制模块进行管理，用户还能够像SIMULINK一样创建自己模块库，在自己库

20、中能够添加、删除模块，或者在这个库中建立子库。当对库中某模块进行修改时，全部模型文件中该模块都被自动修改，这么，就不需要一个个打开模型文件去进行模块修改了。在SIMULINK浏览器窗口中，点击菜单FileNew Library，将出现一个新名为Library:untitled窗口。将定制非线性电阻模块和非线性电感模块复制到新窗口中，保留新窗口，比如取名为mypsb_library。新模块库窗口如图8-21。电力系统的仿真与应用教材第49页第8章 定 制 模 块 图8-21 定制模块库电力系统的仿真与应用教材第50页第8章 定 制 模 块 8.2 S函数编写及应用函数编写及应用S函数是Syste

21、m Function简称，即系统函数，它是扩展SIMULINK功效强有力工具，在很多情况下非常有用。用户能够利用MATLAB、C、C+以及FORTRAN等语言来编制程序，组成S函数模块，并像标准SIMULINK模块一样直接调用。S函数使用一个特殊调用规则，使得用户能够与SIMULINK内部算法进行交互。这种交互和SIMULINK内部算法与内置模块之间交互非常相同，而且能够适合用于不一样性质系统。电力系统的仿真与应用教材第51页第8章 定 制 模 块 图8-22 S函数模块图标本节主要介绍利用MATLAB语言设计S函数，并经过例子介绍S函数应用技巧。8.2.1 S函数模块函数模块S函数模块在Si

22、mulink/User-Defined Functions库中，用此模块能够创建包含S函数SIMULINK模块。S函数模块图标如图8-22所表示。电力系统的仿真与应用教材第52页第8章 定 制 模 块 S函数模块只有一个信号输入口和一个信号输出口，分别对应输入变量和输出变量。S函数模块输入端口只能接收一维向量信号。假如S函数模块含有多个输入变量和输出变量，需要使用“信号合成”(Mux)模块和“信号分离”(Demux)模块，同时需要指定sizes结构适当区域值。注意输入变量个数和输出变量个数必须和S函数内部定义输入变量和输出变量个数相同。S函数模块自动将输入列向量第一个值和S函数内部第一个输入变

23、量对应，第二个值和S函数内部第二个输入变量对应，等等。输出变量也存在一样关系。电力系统的仿真与应用教材第53页第8章 定 制 模 块 双击S函数模块，弹出该模块参数对话框，如图8-23所表示。该对话框中包含以下参数：(1)“S函数文件名”(S-Function Name)文本框：填写S函数文件名，确定后图标上将显示该文件名。注意，不需要扩展名，但该文本框不能为空。(2)“S函数参数”(S-Function Parameters)文本框：填写S函数需要外部参数。S函数模块允许使用外部参数。函数参数能够是MATLAB表示式，也能够是变量。参数并列给出，各参数间以逗号分隔，但参数列不需要用小括号括住

24、表示。若无外部变量，该文本框为空。电力系统的仿真与应用教材第54页第8章 定 制 模 块 图8-23 S函数模块参数对话框电力系统的仿真与应用教材第55页第8章 定 制 模 块(3)“编辑”(Edit)按键：点击该按键将打开指定文件名S函数文件编辑窗口，在该窗口中能够进行编辑工作。【例【例8.3】假设S函数需要外部参数有三种，分别为2、矩阵1,2;3,4和字符串miles。试将这三种参数输入“S函数参数”参数文本框。解：解：在“S函数参数”参数文本框中输入2,1,2;3,4;,miles即可。电力系统的仿真与应用教材第56页第8章 定 制 模 块 8.2.2 S函数编写函数编写S函数模块M文件

25、含有一套固定调用变量规则。能够经过创建一个新M文件并在其中按规则编写S函数，或者直接用SIMULINK提供模板文件sfuntmpl.m，并从这个模板出发构建需要S函数。这类M文件中第一行程序是一个函数语句：function sys,x0,str,ts=f(t,x,u,flag,p1,p2,pn)其中，f为S函数函数名，比如S函数以函数名sfun.m保留，则f应为sfun。输入变量和输出变量说明分别见表8-5和表8-6。电力系统的仿真与应用教材第57页第8章 定 制 模 块 表表8-5 S函数输入变量表函数输入变量表电力系统的仿真与应用教材第58页第8章 定 制 模 块 表表8-6 S函数输出变

26、量表函数输出变量表 电力系统的仿真与应用教材第59页第8章 定 制 模 块 S函数M文件在运行过程中不停检测输入变量flag值，并按照表8-7所表示规则调用不一样内部函数。表表8-7 S函数函数flag参数表参数表电力系统的仿真与应用教材第60页第8章 定 制 模 块 由表8-7可知，变量flag共有六种可能值，下面将分别讨论这六种值所对应操作。1.初始化设置初始化设置flag=0时，调用函数mdlInitializeSizes进行初始化设置。在mdlInitializeSizes函数中，需要对四种参数进行初始值设置。1)sizes首先经过语句sizes=simsizes取得默认系统参数变量s

27、izes。该变量是一个没有初始化结构体变量，共包含有六个字段，各字段说明如表8-8所表示。电力系统的仿真与应用教材第61页第8章 定 制 模 块 表表8-8 结构体变量结构体变量sizes字段字段电力系统的仿真与应用教材第62页第8章 定 制 模 块 对sizes结构体变量中每个元素赋值。sizes中每个元素都必须有值，即使这个值为0。设置好后，经过语句sys=simsizes(sizes)将变量sizes赋值给输出变量sys。【例【例8.4】某S函数不含连续状态变量和离散状态变量，只有一个输入变量和一个输出变量，且输入变量直接控制输出变量，采样周期唯一。试对该S函数结构体变量sizes进行初

28、始化设置。解：解：输入以下语句进行初始化设置：电力系统的仿真与应用教材第63页第8章 定 制 模 块 电力系统的仿真与应用教材第64页第8章 定 制 模 块 2)x0设置状态变量初始值。【例【例8.5】某S函数含有两个连续状态变量和两个离散状态变量。连续状态变量初始值为1.0，离散状态变量初始值为0.0。试对该S函数状态变量初始值进行设置。解：解：x0=1.0,1.0,0.0,0.0;3)str变量str总为空：str=;电力系统的仿真与应用教材第65页第8章 定 制 模 块 4)ts经过设置采样时间和采样延迟时间矩阵ts来定义SIMULINK调用函数时间。SIMULINK提供了各种设置ts方

29、式，如表8-9所表示。表表8-9 采样时间设置采样时间设置电力系统的仿真与应用教材第66页第8章 定 制 模 块 若系统含有多个采样周期，能够用一样方法进行设置。【例【例8.6】某S函数含有两个离散采样周期，第一个采样周期为0.25 s，无延迟；第二个采样周期为1.0 s，采样延迟0.1 s。试对该S函数ts进行初始化设置。解：解：ts=0.25,0;1.0,0.1;电力系统的仿真与应用教材第67页第8章 定 制 模 块 2.连续状态变量微分计算连续状态变量微分计算flag=1时，调用函数mdlDerivatives将连续状态变量微分值赋给sys。【例【例8.7】已知非线性系统当flag=1时

30、，试设置sys属性。解：解：sys(1)=x(2);sys(2)=x(1)-3*x(2)2+u(1);电力系统的仿真与应用教材第68页第8章 定 制 模 块 3.离散状态变量更新离散状态变量更新flag=2时，调用函数mdlUpdates将离散状态变量更新后值赋给sys。【例【例8.8】已知一个一阶单变量离散子系统当flag=2时，试设置sys属性。解：解：sys=x(1)+u(1);4.输出变量计算输出变量计算flag=3时，调用函数mdlOutputs对输出变量进行计算，并把输出变量赋值给sys。电力系统的仿真与应用教材第69页第8章 定 制 模 块 5.下次采样时间计算下次采样时间计算f

31、lag=4时，调用函数mdlGetTimeOfNextVarHit对下一次采样时间进行计算，并把下一次采样时间赋值给sys。只有在初始值设置中已将采样时间设置为可变后，该函数才会被调用。6.仿真结束仿真结束当仿真以某种原因结束时，flag=9，调用函数mdlTerminate执行必要任务终止仿真，不对sys赋值。电力系统的仿真与应用教材第70页第8章 定 制 模 块【例【例8.9】按文件14定制一个非线性电弧炉元件以反应电弧炉电压电流特征，搭建电路，观察效果。解：解：(1)理论分析。文件14所述模型电弧电压v和电流i关系为其中，g为电弧电导，定义为(8-8)(8-7)电力系统的仿真与应用教材第

32、71页第8章 定 制 模 块 r为电弧半径，定义为其中，C对应电弧半径初值；Im为电流幅值；k1、k2、k3是常数，分别为3000、1、12.5。考虑414 Hz频率范围白噪声后，电弧半径为(8-9)(8-10)电力系统的仿真与应用教材第72页第8章 定 制 模 块(2)按图8-24搭建电弧炉电气特征模型。该模型中，电弧炉电气特征用一个可控电压源等效替换，可控电压源大小主要由电弧炉电弧半径决定，而电弧半径由电流幅值决定，同时电弧半径上还叠加了一个噪声。该模型中选取各模块名称及提取路径见表8-10。电力系统的仿真与应用教材第73页第8章 定 制 模 块 表表8-10 例例8.9电弧炉电气特征模型

33、中选取模块电弧炉电气特征模型中选取模块名称及提取路径名称及提取路径电力系统的仿真与应用教材第74页第8章 定 制 模 块 图8-24 电弧炉电气特征模型电力系统的仿真与应用教材第75页第8章 定 制 模 块 设置正弦波发生器模块中电流幅值为3 A、频率为100pi、相角为0。FFT模块基频为50 Hz，提取基频分量(1次分量)，滤波器为带通滤波器，截断频率为14 Hz。(3)编写s函数。在MATLAB/toolbox/Simulink/Blocks目录下复制模块文件sfuntmpl.m到MATLAB/work工作目录下，打开该模板，另存为myarc_r.m和myarc_Ron.m，其中myar

34、c_r为电弧半径计算模块，myarc_Ron为电弧炉等效电阻计算模块。myarc_r程序代码以下：电力系统的仿真与应用教材第76页第8章 定 制 模 块 电力系统的仿真与应用教材第77页第8章 定 制 模 块 电力系统的仿真与应用教材第78页第8章 定 制 模 块 电力系统的仿真与应用教材第79页第8章 定 制 模 块 电力系统的仿真与应用教材第80页第8章 定 制 模 块 电力系统的仿真与应用教材第81页第8章 定 制 模 块 电力系统的仿真与应用教材第82页第8章 定 制 模 块(4)封装S函数模块。以上模型还不能运行，因为其中外部参数未设置。所以，对图8-24所表示模型按8.1节所述方法

35、进行封装，并设置参数对话窗口，如图8-25所表示。电力系统的仿真与应用教材第83页第8章 定 制 模 块 图8-25 电弧炉模型参数对话框电力系统的仿真与应用教材第84页第8章 定 制 模 块 图8-26 电弧炉模型图标封装完成后电弧炉模型图标如图8-26所表示。电力系统的仿真与应用教材第85页第8章 定 制 模 块(5)开始仿真，仿真结果如图8-27所表示。图中从上到下依次为电弧炉上电流、电弧炉端口电压和电弧炉电弧半径。图8-28为实测电弧炉电压电流特征。从电压电流特征曲线上能够看出，仿真效果良好。从上例能够看出，经过编写S函数，并与MATLAB/SIMULINK结合，定义、搭建新模型是一个

36、非常简单、直观过程，它极大地扩展了SIMULINK功效。电力系统的仿真与应用教材第86页第8章 定 制 模 块 S函数还允许采取C、C+、FORTRAN等语言进行编写。不过，C语言编写和调试过程比MATLAB语言编写S函数过程要复杂得多，所以在纯仿真中最好使用MATLAB语言去编写。当然，在有些应用中，因为MATLAB语言S函数不能转换成C语言程序，生成独立文件，所以应该采取C语言去编写S函数。关于采取C、C+、FORTRAN等语言编写S函数方法，读者能够参考MATLAB相关帮助。电力系统的仿真与应用教材第87页第8章 定 制 模 块 图8-27 电弧炉模型图标电力系统的仿真与应用教材第88页第8章 定 制 模 块 图8-28 电弧炉电压电流特征电力系统的仿真与应用教材第89页