《MATLAB变压器仿真设计》由会员分享，可在线阅读，更多相关《MATLAB变压器仿真设计（26页珍藏版）》请在金锄头文库上搜索。
1、大学专业软件应用综合设计报告 水能 学院 13 级 电气 专业题目 变压器综合仿真设计二 学生 某某某 学号 131504207 指导教师 建华 2015年 12 月 30 日目录1、 设计题目.22、 正文.2 1、引言.2 2、设计依据及框图.3 2.1 设计平台.3 2.2 设计思想.4 2.3 设计结构框图或流程图.6 2.4各模块功能简介.6 3、软件调试分析.10 4、结语.23 5、参考文献.25 6、致.25.变压器综合仿真设计二摘要：随着变压器技术的进步，传统仿真已经受到了很大的限制。并且当下要推动变压 器技术的发展，已经不能再依靠传统仿真。因此，对于变压器的计算机仿真技术势
2、在必行。本为通过MATLAB软件，对变压器的运行特性进行了仿真。主要仿真的容包括：变压器磁路电流畸变以及变压器负载运行特性曲线研究。仿真所用到的方法为数值计算方法，通过插值的方法实现了对曲线的拟合。仿真时，结合实际情况可输入不同参数便于研究。文中给出了各种运行特性的仿真结果图，并且结合理论对其做了简单的分析，验证了仿真方法的准确性和可行性。关键字：变压器；MATLAB仿真分析；曲线拟合1 引言设随着科学技术进步，电工电子新技术的不断发展，新型电气备不断涌现，人们使用电的频率越来越高，人与电的关系也日益紧密，对于电性能和电气产品的了解，已成为人们必需的生活常识。变压器是一种静止的电气设备，它是利
3、用电磁感应原理把一种电压的交流电能转变成同频率的另一种电压的交流电能，以满足不同负载的需要。在电力系统中，变压器是一个重要的电气设备，它对电能的经济传输，灵活分配和安全使用具有重要的作用，此外，也使人们能够方便地解决输电和用电这一矛盾。由于计算机仿真技术的出现，传统的物理仿真系统逐渐的被计算机仿真系统代替。计算机仿真系统所具有的效率高、精度高、重复性和通用性好、容易改变仿真参数等优点，还可以实现物理仿真无法实现的有危险性的或者是成本昂贵的仿真。在我国电力行业发展迅速的今天，变压器的仿真技术不能够再依托于传统的物理仿真系统，而是需要能够采用能够促进变压器技术发展的仿真技术。对变压器特性的仿真涉及
4、到很多方面，比如变压器空载励磁电流在饱和和磁滞影响时的特性、变压器磁滞回环在不同电压等级下的数据仿真、变压器空载合闸时的过电流现象、变压器在突发短路时的过电流现象，还有基本的比如效率特性、外特性、短路试验、空载试验等。在学习完本课程后，运用MATLAB相关仿真技术对变压器进行仿真研究，本文的仿真主要以变压器磁路电流畸变以及变压器负载运行特性曲线为主要研究对象，通过结合实际进行曲线拟合、波形分析，得出相应结论。2 设计依据及框图2.1 设计平台计算机技术的发展使得大量的数据计算变得方便快捷，一些因为需要不断的迭代而数据量庞大的数学算法也可以在实际中得到应用。不仅将工作者从繁忙的数据计算中解脱出来
5、，而且还可以做到不同精度的计算。MATLAB软件在数值计算方面独占鳌头，由于其提供了数据视图，文字处理的同一环境而受到欢迎。MATLAB的中文意思为矩阵实验室，起源于20世纪80年代，由其开创者Cleve Moler开发。经过后期的不断完善，MATLAB最终走向正轨，并且由MATH WORKS公司以商品形式发布。从MATLAB的开创到现今，随着其版本的更替，功能也变得愈加强大。其核心编写所采用的语言最终也从FORTRAN语言变为了C语言。MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式
6、程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、
7、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+，JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。其主要功能有：数值分析、数值和符号计算、工程与科学绘图、控制系统的设计
8、与仿真、数字图像处理、数字信号处理、通讯系统设计与仿真、财务与金融工程等MATLAB自产生之日起就具有方便的数据可视化功能，以将向量和矩阵用图形表现出来，并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图。可用于科学计算和工程绘图。新版本的MATLAB对整个图形处理功能作了很大的改进和完善，使它不仅在一般数据可视化软件都具有的功能（例如二维曲线和三维曲面的绘制和处理等）方面更加完善，而且对于一些其他软件所没有的功能（例如图形的光照处理、色度处理以及四维数据的表现等），MATLAB同样表现了出色的处理能力。同时对一些特殊的可视化要求，例如图形对话等，MA
9、TLAB也有相应的功能函数，保证了用户不同层次的要求。另外新版本的MATLAB还着重在图形用户界面（GUI）的制作上作了很大的改善，对这方面有特殊要求的用户也可以得到满足。MATLAB对许多专门的领域都开发了功能强大的模块集和工具箱。一般来说，它们都是由特定领域的专家开发的，用户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码。领域，诸如数据采集、数据库接口、概率统计、样条拟合、优化算法、偏微分方程求解、神经网络、小波分析、信号处理、图像处理、系统辨识、控制系统设计、LMI控制、鲁棒控制、模型预测、模糊逻辑、金融分析、地图工具、非线性控制设计、实时快速原型及半物理仿真、嵌入式系
10、统开发、定点仿真、DSP与通讯、电力系统仿真等，都在工具箱（Toolbox）家族中有了自己的一席之地。MATLAB包括拥有数百个部函数的主包和三十几种工具包。工具包又可以分为功能性工具包和学科工具包。功能工具包用来扩充MATLAB的符号计算，可视化建模仿真，文字处理及实时控制等功能。学科工具包是专业性比较强的工具包，控制工具包，信号处理工具包，通信工具包等都属于此类。开放性使MATLAB广受用户欢迎。除部函数外，所有MATLAB主包文件和各种工具包都是可读可修改的文件，用户通过对源程序的修改或加入自己编写程序构造新的专用工具包。现在的MATLAB，已经发生了质的飞跃。完善的数值计算系统和简单的
11、程序编写环境，使得MATLAB软件不需要太过专业的程序编写技能基础就可以轻松的使用。MATLAB以矩阵作为运算的基本单位，使得矩阵运算变得方便快捷。此外，MATLAB所提供的丰富的函数可以很容易的实现各种数值算法。MATLAB最为突出的特点就是编程语言的简洁、直观。而且其对语法要求不是特别严格，不像其他编程语言，比如C语言等。MATLAB语言程序文件为文本文件，后缀为.m，一般称为M文件。M文件具有保存和容易修改命令的优点。MATLAB提供了专门的M文件编辑器。通过M文件还可以自行的定义具有具体功能的函数，使的程序的编写得到简化。MATLAB提供的数据可视功能为实现仿真的结果分析提供了方便，通
12、过图像可以很容易的对大量数据的变化实现分析。2.2 设计思想对变压器特性的仿真涉及到很多方面，比如变压器空载励磁电流在饱和和磁滞影响时的特性、变压器磁滞回环在不同电压等级下的数据仿真、变压器空载合闸时的过电流现象、变压器在突发短路时的过电流现象，还有基本的比如效率特性、外特性、短路试验、空载试验等。变压器动态特性的分析主要是为计算和分析变压器在空载合闸和突然短路或者是其他故障时出现的暂态过电流和过电压。在电力系统中，变压器从发电厂到输配电网中都充当着重要的角色。变压器的运行特性直接影响到电力系统的正常运行，所以对变压器的运行特性进行研究是非常必要的。电力变压器的保护受到变压器励磁涌流的影响，单
13、相变压器空载合闸产生励磁涌流的大小与变压器合闸角有关，通过对励磁电流特点的研究可以制定出防止励磁涌路引起误动的方法。随着人们对变压器技术的不断探索，变压器技术已从刚开始的基本感应定律发展到能够对变压器的结构进行合理的设计、在理论上对变压器运行时的各种现象进行详细分析、使变压器的理论模型达到更高的精确程度。在20世纪80年代末期，国外学者引入了磁化曲线和各频率下的等效激磁并联电导，建立了变压器的频域模型，但是在计算中存在了一定的误差，并且数据参数的获取非常复杂。另一种方法为引入磁滞回线，建立时域模型，这种方法较为精确，但是在拟合磁化曲线建立较为准确的曲线模型存在着困难。有的学者通过传统的T型等效
14、电路的励磁支路来实现变压器非线性的特性，但是在仿真中要求对铁芯损耗进行估计，而且采用的计算方法相当的复杂，难以实现。另外有学者对变压器的非线性特性进行了详细的分析，但是并没有给出具体的模型和仿真结果。很多的文献中提到了分段磁化曲线的简单模拟方法，但是通过这种方法的仿真结果不够精确，模型的准确程度不够。后来出现的通过拟合曲线的方法来准确的反应铁芯磁化的非线性。变压器的Jiles.Atherton模型详细的说明了铁芯的磁化过程，并且通过磁学理论明了铁芯磁滞现象的原理。在这个基础上学者通过MATLAB软件进行了仿真，仿真所涉到的五个模型参数包括饱和磁化强度、表明滞后磁化曲线形状的参数、可逆磁化系数、
15、磁畴对运动阻碍作用的参数、磁畴间相互作用的参数，通过实验的方法可以获得。为了能够使计算速度更快，没有采用经典的四阶龙格库塔算法来求解，而是采用了较为简单的欧拉算法。后来较为常见的描述变压器磁滞回线的方法是将磁滞回线分为了主磁滞环和饱和区，这种描述可已较为精确的考虑到磁滞特性。通过拟合曲线的方法，将实验数据拟合成为接近实际的磁化曲线。主磁滞环采用修正的反正切函数加以拟合，饱和区特性将其视为线性可逆，即磁化曲线用直线表示。在拟合出最大的磁滞回环后，通过对曲线按照比例的压缩，就可以得到一系列的反应不同程度的磁滞回线。通过龙格库塔算法求解非常系数的微分方程，考虑到了励磁阻抗随着饱和程度而变化地情况。随
16、着仿真软件所提供的功能越来越强大，变压器的仿真可以省去编写程序的负担，而是直接用仿真软件中的功能块来实现，通过对模块的合理的联结和参数设置就可以完成仿真。在MATLAB环境下，基于单相变压器数学模型的基础，建立单独的功能模块，通过有机的整合功能模块搭建单相变压器的仿真模型。变压器以磁链为状态变量的数学模型与以电流作为状态变量的数学模型比较，前者在物理概念上更加的清晰。而且在采用数值积分的方法来计算时，此模型要更加稳定。上述所提到的利用数值计算来实现仿真的方法虽然都较为的成熟了，但是数值计算方法在仿真时对于初值的依赖较强。如果初值的设定不合理，仿真程序很容易出现错误。甚至会得到发散的仿真结果，使
17、得程序无止境的循环下去。本文将通过MATAB软件对变压器的稳态和动态特性进行仿真。稳态特性仿真主要是考虑变压器磁路电流畸变以，动态特性仿真为变压器负载运行特性曲线研究。方法依托于MATLAB仿真软件的数值计算能力。在对变压器磁路电流畸变数据处理时将用到三次样条插值拟合曲线的方法。动态特性仿真同样是利用了MATLAB的数值计算功能。通过对仿真结果的分析，验证数值算法仿真的正确性和可行性。主要针对单相变压器负载运行时的参数分析，关注负载运行效率问题，研究负载运行效率与不同性质负载及不同负载系数之间的关系，通过仿真得出相关波形再进行讨论研究。2.3 设计结构框图或流程图设计仿真模型变压器综合仿真编写
18、仿真程序不同负载时的外特性变压器电流畸变研究负载特性与效率曲线设置参数进行调试波形、数据情况观察分析数据、波形图1 变压器综合仿真框图 图2变压器综合仿真流程图2.4各模块功能简介SIMULINK是一个对动态系统（包括连续系统、离散系统和混合系统）进行建模、仿真和综合分析的集成软件包，是MATLAB的一个附加组件，其特点是模块化操作、易学易用，而且能够使用MATLAB提供的丰富的仿真资源。在SIMULINK环境中，用户不仅可以观察现实世界中非线性因素和各种随机因素对系统行为的影响，而且也可以在仿真进程中改变感兴趣的参数，实时地观察系统行为的变化，因此已在许多领域，如通信、信号处理、DSP、电力
19、、金融、生物系统等获得重要应用。对电类专业的学生来说，无论是学习专业课程或相关课程设计，还是在今后的工作中，SIMULINK都是一个重要的仿真建模工具。SimPowerSystems是在Simulink环境下进行电力电子系统建模和仿真先进工具。SimPowerSystems是Simulink下面的一个专用模块库，包含电气网络中常见的元器件和设备，以直观易用的图形方式对电气系统进行模型描述。模型可与其它Simulink模块的相连接，进行一体化的系统级动态分析。一、SimPowerSystems专用模块库的特点：1.使用标准电气符号进行电力系统的拓扑图形建模和仿真；2.标准的AC和DC电机模型模块
20、；变压器；传输线；信号和脉冲发生器；HVDC控制；IGBT模块和大量设备模型，有断路器，二极管，IGBT，GTO，MOSFET和晶闸管；3.使用Simulink强有力的变步长积分器和零点穿越检测功能，给出高度精确的电力系统仿真计算结果4.为快速仿真和实时仿真提供了模型离散化方法；5.提供多种分析方法，可以计算电路的状态空间表达、计算电流和电压的稳态解、设定或恢复初始电流/电压状态、电力机械的潮流计算；6.提供了扩展的电气系统网络设备模块，如电力机械，功率电子元件，控制测量模块和3相元器件；7.提供36个功能演示模型，可直接运行仿真；二、SimPowerSystems专用模块库的强大功能：Sim
21、PowerSystems中的模块SimPowerSystems中模块的数学模型基于成熟的电磁和机电方程，用标准的电气符号表示。它们可以同标准的Simulink模块一起使用建立包含电气系统和控制回路的模型。连接通过与SimPowerSystems提供的测量模块实现。SimPowerSystems拥有近100个模块，分别位于7个子模块库中。这些库模块涵盖了以下应用围：1.电气网络（ElectricalSources&Elements）RLC支路和负载，型传输线，线性和饱和变压器，浪涌保护，电路分离器，互感，分布参数传输线，3相变压器（2个和3个绕组），AC和DC电压源，受控电压源和受控电流源。2.
22、电力机械（Machines）完整或是简化形式的异步电动机，同步电动机，永磁同步电动机，直流电动机，激磁系统和水轮机涡轮机/调速系统模型。3.电力电子（PowerElectronics）.二极管，简化/复杂晶闸管，GTO，开关，MOSFET，IGBT和通用型桥接管模型。4.控制和测量模块（Measurements）电压、电流和电抗测量，RMS测量，有功和无功功率计算，计时器，万用表，傅立叶分析，HVDC控制，总谐波失真，abc到dq0和dq0到abc轴系变换，3相V-I测量，3相脉冲和信号发生，3相序列分析，3相PLL和连续/离散同步6-，12-脉冲发生器。5.三相网络元器件（Electrica
23、lSources&Elements）3相RLC负载和支路，3相断路器，3相，3相电抗，型传输线，AC电压源，6-脉冲二极管和晶闸桥管，整流二极管，Y-/Y-/Y-/Y-Y-可配置3相变压器。3、 SimPowerSystems常用模块库简介1. Continuous:（连接器元件库有10种模块）2. ElectricalSources:（电源元件库有7种电源功能模块） 1）DCVoltageSource:直流电压源 2）ACVoltageSource:交流电压源 3）ACCurrentSource:交流电流源 4）ControlledVoltageSource:受控电压源 5）Controll
24、edCurrentSource:受控电流源 6）3-PhaseSource:三相电源 7）3-PhaseProgrammableVoltageSource:三相可编程电压源 3.Elements:（线路元件库有24种模块） 4.Machines:（电机元件库有16种模块） 5.Measurements:（电路测量模块元件库有5种模块）而本文所采用的仿真方法大多是是利用m文件进行仿真研究。下面就这种方式及相关知识进行介绍。Matlab输入命令的常用方式有两种：一种是直接在Matlab的命令窗门中逐条输入Matlab命令；二是m文件工作方式。当命令行很简单时，使用逐条输入方式还是比较方便的。但当命
25、令行很多时(比如说几十行乃至全成百上千行命令)，显然再使用这种方式输入MATLAB命令，就会显得杂乱无章，不易于把握程序的具体走向，并且给程序的修改和维护带来了很大的麻烦。这时，建议采用Matlab命令的第二种输入形式m文件工作方式。m文件工作方式，指的是将要执行的命令全部写在一个文本文件中，这样既能使程序显得简洁明了，又便于对程序的修改与维护。m文件直接采用Matlab命令编写，就像在Matlab的命令窗口直接输入命令一样，因此调试起来也十分方便，并且增强了程序的交互性。m文件与其他文本文件一样，可以在任何文本编辑器中进打编辑、存储、修改和读取。利用m文件还可以根据白己的需要编写一些函数，这
26、些函数也可以橡Matlab提供的函数一样进行调用。从某种意义上说，这也是对MATLAB的二次开发。m文件有两种形式：一种是命令方式或称脚本方式；另一种就是函数文件形式。两种形式的文件扩展名均是.m。1、 M文件 当遇到输入命令较多以及要重复输入命令的情况时，利用命令文件就显得很方便了。将所有要执行的命令按顺序放到一个扩展名为.m的文本文件中，每次运行时只需在命令窗口输入m文件的文件名就可以了。需要注意的是，m文件最好直接放在Matlab的默认搜索路径下(一般是Matlab安装目录的子目录work中)，这样就不用设置m文件的路径了，否则应当用路径操作指令path重新设置路径。另外，m文件名不应该
27、与Matlab的置函数名以及工具箱中的函数重名，以免发生执行错误命令的现象。Matlab对命令文件的执行等价于从命令窗口中顺序执行文件中的所有指令。命令文件可以访问Matlab工作空间里的任何变量及数据。命令文件运行过程中产生的所有变量都等价于从Matlab工作空间中创建这些变量。因此，任何其他命令文件和函数都可以自由地访问这些变量。这些变量一旦产生就一直保存在存中，只有对它们重新赋值，它们的原有值才会变化。关机后，这里变量也就全部消失了。另外，在命令窗口中运行clear命令，也可以把这些变量从工作空间中删去。当然，在Matlab的工作空间窗口中也可以用鼠标选择想要删除的变量，从而将这些变量从
28、工作空间中删除。 2、M函数m函数文件是一个特殊的m文件，其常见格式如下：funcdon返问变量列表函数名(输入变量列表)注释说明语句段函数体语句。需要说明的是，这里输入变量的个数以及输出变量的个数是由MATLAB本身提供的两个保留变量nargin和nargout来给出的，它们分别是Numberoffunctioninputarguments和Numberoffunctionoutputarguments的缩写形式。输入变量要用逗号隔开，输出变量多于1个时，要用方括号括起来。用户可以借助于help命令显示其中的注释说明语句段。通过这样的方法就可以建立函数文件或者称m函数，其调用方法与一放的Ma
29、tlab函数的调用方法相同。函数文件相当于对Matlab进行了二次开发。其作用与其他高级语言子函数的作用基本相同，都是为了实现特定目的而由用户自己编写的子函数。函数文件与命令文件有着鲜明的区别：函数文件的第一行必须包含function字符；命令文件无此要求。函数文件的第一行必须指定函数名、输入参数及输出参数，命令文件无此要求。一个函数文件可以合0个、1个或多个输入参数和返回值。函数文件要在文件的开头定义函数名，如functiony1,y2=func(x,a,b,c)，则该函数文件名必须存为func.m，而命令文件无此要求。命令文件的变量在文件执行结束以后仍然保存在存中而不会丢失，而函数文件的变
30、量仅在函数运行期间有效(除非用global把变量说明成全局变量，否则函数文件中的变量均为局部变量)，当函数运行完毕后，这些变量也就消失了。需要说明的是，调用函数时所用的输入输出变量名并不要求与编写函数文件时所用的输入输出变量名相同。 3 软件调试分析3.1励磁电流的研究图3磁通、磁化曲线、励磁电流波形图 铁芯在反复磁化时，磁化特性曲线就会偏离基本磁化曲线，励磁电流会滞后于磁通的变化。磁滞损耗使得励磁电流的波形超前于磁通波形。实际上铁芯中涡流产生的损耗也会使得励磁电流超前于磁通波形。 空载励磁电流波形除了受到铁芯磁化的非线性影响外，还受到铁芯磁滞现象的影响。如果铁芯在剩磁为零时开始磁化，那么铁芯
31、磁化特性将沿着基本磁化曲线进行，但是当磁通开始发生转折后，磁化特性将不再按照基本磁化曲线进行，而是磁通变化滞后于励磁电流的变化，这种现象为铁芯的磁滞现象。由于铁芯磁滞想象的存在，铁芯的磁化呈现为基本磁化曲线两侧的回线。在这种情况下，同一个数值的磁通会对应两个励磁电流。下图以为MATLAB程序绘制的磁磁滞回线。有一实际例子磁滞回线磁通和电流数据为分别为：1=-1.0，-0.90，-0.79，-0.50，0.30，0.70，0.85，0.91，0.95，0.97，1.0，2=-1.0，-0.97，-0.91，-0.85，-0.70，-0.30，0.50，0.79，0.90，1.0，i= -0.05
32、：0.01：0.05，单位为.p.u。其中1 与i对应磁滞回线中的下降曲线，2 与i对应磁滞回线中的上升曲线。当输入磁通信号时，依旧用磁通信号的数据对磁滞回线进行插值计算。将插值计算的结果可视化输出，即得到受到铁芯磁滞影响的励磁电流波形图。从图中知道励磁电流在铁芯磁滞现象的影响下，励磁电流波形超前于磁通波形。这说明铁芯的磁滞现象在铁芯中造成了损耗，称之为磁滞损耗。当改变磁通信号幅值时，就会得到受不同程度磁滞影响的励磁电流波形。图4 磁滞回线图 图5 考虑磁滞后的电流波形图3.2 变压器磁路电流畸变程序设计变压器铁芯磁化数据如表所示，设磁路中磁通按正弦规律变化，分析励磁电流的变化规律。铁磁材料磁
33、化数据1234567891011121314151617I1.381.541.711.912.162.502.933.494.225.216.528.3610.914.822.437.864.00.400.480.560.640.720.800.880.961.041.121.201.281.361.441.521.601.68(1)研究磁路中按正弦规律变化的磁通, 研究磁通与时间的关系，即=f(t)。M文件具体如下：idata=1.38,1.54,1.71,1.91,2.16,2.50,2.93,3.49,4.22,5.21,6.52,8.36,10.91,14.8,22.4,37.8,64
34、.0;phidata=0.40:0.08:1.68;t=0:0.001:pi;phi=2*sin(t);i=spline(phidata,idata,phi);plot(t,phi)xlabel(t/s);ylabel(/wb)title(=f(t)可得仿真波形，从仿真波形可以看出磁路中磁通按正弦规律变化。图6 磁路磁通图(2)研究磁通与电流的关系，即=f(I)。M文件具体如下：idata=1.38,1.54,1.71,1.91,2.16,2.50,2.93,3.49,4.22,5.21,6.52,8.36,10.91,14.8,22.4,37.8,64.0;phidata=0.40:0.08
35、:1.68;t=0:0.001:pi;phi=2*sin(t);i=spline(phidata,idata,phi);plot(i,phi)xlabel(I/A)ylabel(/Wb)title(=f(I)axis(0 400 0 3)图7 磁通电流关系图(3)研究电流的畸变规律，即I=f(t)idata=1.38,1.54,1.71,1.91,2.16,2.50,2.93,3.49,4.22,5.21,6.52,8.36,10.91,14.8,22.4,37.8,64.0;phidata=0.40:0.08:1.68;t=0:0.001:pi;phi=2*sin(t);i=spline(p
36、hidata,idata,phi);plot(t,i)xlabel(t/s)ylabel(I/A)title(I=f(t)axis(0,3.5,0,400)图8 畸变电流图小结： 由于变压器的铁芯存在着饱和现象，即铁芯磁化的非线性现象，使得磁通与与励磁电流之间的关系= f(i0)，呈现非线性。磁通和电压按照同样的规律变化，即正弦变化，只是要滞后电压0.5*pi的电角度。根据磁化曲线可以求得对应每一时刻的励磁电流，从而绘制励磁电流波形图。上面的波形图为正弦变化地磁通曲线、基本磁化曲线和励磁电流曲线。图中励磁电流波形偏离了正弦波，呈现为尖顶波。3、变压器负载运行 变压器的负载运行是指原绕组接入电源
37、电压，副绕组接负载时的工作状况。这时，变压器的副边也有电流流通，原边的接入电路与空载相比相应增大，副边端电压将受到负载的影响而发生变化，这是负载运行与空载运行的主要区别。 在变压器负载运行时，由于变压器的一次侧漏阻抗很小，即使负载时的一次侧电流变为了空载电流的很多倍，依旧存在了I1nZ1U1的关系。所以忽略一次侧电压在漏阻抗上的电压降落，认为U1E1=4.44fN1m，所以负载运行时的磁通与空载运行时的磁通相差无几。二次侧电流产生的磁动势势必会影响到励磁磁动势，要使得缠身磁通的励磁磁动势不发生变化，一次侧必须有中和二次侧电流产生的磁动势的电流。所以负载运行时，一次侧电流可以将其视为励磁分量和负
38、载分量之和。这里直接给出变压器的基本方程，如下： (3-1)将二次侧量折合到一次侧，基本方程变为下式： (3-2)图9 T形等效电路根据式(2-10)，可以得到变压器的T形等效电路，如图9。负载时的一次侧电流要比励磁电流大很多，如果忽略励磁电流，便可简化T形等效电路为变压器的简化等效电路。如图10。另，电路图变为如图11。图11 短路阻抗表示的等效电路图10 简化电路图变压器在空载时的二次侧电压为其额定电压。而当变压器二次侧接负载时，变压器的二次侧电压会发生变化，用表示：（3-3）根据给出的等效电路图和简化向量图，给出另一个的计算公式： （3-4）其中为负载因数，Rk和Xk为短路电阻和短路电抗
39、，为短路阻抗的阻抗角。变压器带感性负载时电压降低，带容性负载时电压升高。变压器效率受到铁损和铜损的影响。由于变压器空载和负载时铁芯中主磁通基本不变，所以我们认为变压器铁损基本不变。而变压器铜损随着变压器负载的变化会发生变化。设额定电流时，变压器的空载损耗为，那么我们可以将铜损表达如下：（3-5）这里我们给出变压器效率的计算公式：（3-6）本次课程设计中，针对上述理论结合实际利用m文件编程，所编写的程序应用性强，可适用于不同规格型号的变压器，只需提供输入相关参数即可。下面是各程序的主体：(1)求解变压器相关参数数值的程序%输入基本参数SN=input(请输入额定容量SN=);U1N=input(
40、请输入一次侧额定电压U1N=);U2N=input(请输入一次侧额定电压U2N=);r1=input(请输入一次绕组漏阻值r1=);r2=input(请输入二次绕组漏阻值r2=);x1=input(请输入一次绕组漏抗值x1=);x2=input(请输入二次绕组漏抗值x2=);rm=input(请输入励磁电阻rm=);xm=input(请输入励磁电抗xm=);ZL=input(请输入负载阻抗ZL=);%计算额定电流和变比I1N=SN/U1N;I2N=SN/U2N;k=U1N/U2N;%计算T型等效电路中的未知参数，其中Z1为原边漏阻抗，rr2为副边漏阻抗折算值，xx2*rr2为副边漏电抗折算值%
41、ZZ2为副边漏阻抗折算值，ZZL为负载阻抗折算值，Zm为励磁阻抗，Zd为输入阻抗Z1=r1+j*x1;rr2=k2*r2;xx2=k2*x2;ZZ2=rr2+j*xx2;ZZL=k2*ZL;Zm=rm+j*xm;Zd=Z1+1/(1/Zm+1/(ZZ2+ZZL);%U1I为原边电压，I1I原边电流，E1I励磁电动势，I22I折算后副边电流，I2I副边电流，U22I折算后副边电压%U2I副边电压U1I=U1N;I1I=U1I/Zd;E1I=-(U1I-I1I*Z1);I22I=E1I/(ZZ2+ZZL);I2I=k*I22I;U22I=I22I*ZZL;U2I=U22I/k;%cospsil输入
42、侧功率因数，cospsi2负载侧功率因数，P1输入有功功率，P2输出有功功率cospsi1=cos(angle(Zd);cospsi2=cos(angle(ZL);P1=abs(U1I)*abs(I1I)*cospsi1;P2=abs(U2I)*abs(I2I)*cospsi2;eta=P2/P1;%ImI励磁电流，pFe铁损耗，pcu1原边铜损耗，pcu2副边铜损耗ImI=E1I/Zm;pFe=abs(ImI)2*rm;pcu1=abs(I1I)2*r1;pcu2=abs(I2I)2*r2;%输出数据disp(原边电流=),disp(abs(I1I)disp(副边电流=),disp(abs(
43、I2I)disp(副边电压=),disp(abs(U2I)disp(原边功率因数=),disp(cospsi1)disp(原边功率=),disp(P1)disp(副边功率因数=),disp(cospsi2)disp(副边功率=),disp(P2)disp(效率=),disp(eta)disp(励磁电流=),disp(abs(ImI)disp(铁损耗=),disp(pFe)disp(原边铜损耗=),disp(pcu1)disp(副边铜损耗=),disp(pcu2)下面通过一个具体事例进行分析，有一单相变压器参数：SN=10kVA，U1N/ U2N=380/220V，r1=0.14，r2=0.03
44、5，x1=0.22，x2=0.055，rm=30，xm=310。通过输入相关参数可以得到仿真结果为：请输入额定容量SN=10e3请输入一次侧额定电压U1N=380请输入二次侧额定电压U2N=220请输入一次绕组漏阻值r1=0.14请输入二次绕组漏阻值r2=0.035请输入一次绕组漏抗值x1=0.22请输入二次绕组漏抗值x2=0.055请输入励磁电阻rm=30请输入励磁电抗xm=310请输入负载阻抗ZL=4+j*3原边电流= 25.5752副边电流= 42.7447副边电压= 213.7237原边功率因数= 0.7725原边功率= 7.5072e+003副边功率因数= 0.8000副边功率= 7
45、.3085e+003效率= 0.9735励磁电流= 1.1998铁损耗= 43.1836原边铜损耗= 91.5725副边铜损耗= 63.9489可以看出，上面的程序对于求解变压器参数很便捷，只要输入一些已知的基本参数便可求出另外的一些重要参数，另外程序在编写的过程中也体现出“自动”的概念，及实现了输入参数与程序相结合的方式，而不仅仅简单的把数据直接编写在程序中，实现了一定的“智能化”。（2）变压器不同负载性质运行时的效率的程序M文件程序主体：c0=input(请输入阻性功率因数c0=);s0=input(请输入阻性因数s0=);c1=input(请输入感性功率因数c1=);s1=input(请
46、输入感性因数s1=);c2=input(请输入容性功率因数c2=);s2=input(请输入容性因数s2=);xk=input(请输入短路电抗xk=);rk=input(请输入短路阻抗rk=);x=0:0.1:1.5;du2=x*(rk*c2+xk*s2);du1=x*(rk*c1+xk*s1);du0=x*(rk*c0+xk*s0);u0=1-du0;u1=1-du1;u2=1-du2;plot(x,u2,x,x,u1,o,x,u0,x,1,-.);xlabel(负载系数);ylabel(u2/.p.u);title(变压器的外特性);legend(感性负载,容性负载,纯电阻负载,3)下面以
47、一个例子加以说明，若有一变压器不同负载参数如下：c2=0.8,s2=0.6,xk=0.0504,rk=0.024;c1=0.8,s1=-0.6;c0=1,s0=0;对其进行仿真可以得到下面的波形。图12外特性图（3）变压器运行效率与负载系数的程序M文件程序如下：f=inline(x*sn*c/(x*sn*c+p0+x.2*pkn),sn,c,p0,pkn,x);c=input(请输入功率因数值c：);sn=input(请输入额定容量值sn：);p0=input(请输入空载损耗值p0：);pkn=input(请输入额定损耗值pkn：);n=1;for x=0:0.01:1.5h(n)=f(sn,
48、c,p0,pkn,x);n=n+1;endx=0:0.01:1.5plot(x,h)axis(0 1.5 0 1)ymax=max(h)xm=find(diff(sign(diff(h)=-2)+1;xmax=xm(1)*0.01plot(x,h,xm(1)*0.01,max(h),r*)xlabel(负载系数)ylabel(变压器效率)title(效率特性曲线)text(xm(1)*0.01,max(h)+0.04,最高效率点)下面以本文前面所引出的一个例子进行参数设置仿真：请输入功率因数值c：0.8 请输入额定容量值sn：10e3请输入空载损耗值p0：0.068e3请输入额定损耗值pkn：
49、0.306e3图12 效率特性图结果：最高效率=0.9652负载系数=0.4800小结：（1） 从曲线上可以看出，空载时输出功率为0，所以效率=0。（2） 负载较小时，由于变压器铁损耗相对较大，故较小。（3）随着负载的增加，效率也随之增加。但超过某一负载时，因铜耗占的比例增加，效 率反而降低。4 结语4.1 结论与讨论课程设计诚然是一门专业课，给我很多专业知识以及专业技能上的提升，同时又是一门讲道课，一门辩思课，给了我许多道，给了我很多思，给了我莫大的空间。同时，设计让我感触很深。使我对抽象的理论有了具体的认识。通过这次课程设计，不仅培养了独立思考、动手操作的能力，在各种其它能力上也都有了提高
50、。更重要的是，在实验课上，我们学会了很多学习的方法。而这是日后最实用的，真的是受益匪浅。要面对社会的挑战，只有不断的学习、实践，再学习、再实践。这对于我们的将来也有很大的帮助。以后，不管有多苦，我想我们都能变苦为乐，找寻有趣的事情，发现其中珍贵的事情。就像中国提倡的艰苦奋斗一样，我们都可以在实验结束之后变的更加成熟，会面对需要面对的事情。回顾起此课程设计，至今我仍感慨颇多，从理论到实践，在这段日子里，可以说得是苦多于甜，但是可以学到很多很多的东西，同时不仅可以巩固了以前所学过的知识，而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远
51、远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务，从而提高自己的实际动手能力和独立思考的能力。在设计的过程中遇到问题，可以说得是困难重重，但可喜的是最终都得到了解决。这次课程设计虽然只持续了不到一周的时间，中间还夹着上课，真正用在实验上的时间并不是很多，但是同学们却一致认为我们从中学到了很多东西。首先，之前我们在学习MATLAB时，注重它的基本使用方法和语句，而且也还是停留在很浅显的层面上的，包括之前的考试也是注重的理论知识。而我们知道MATLAB是一门实用性的语言，它有着强大的图形图像功能，很多东西都需要我们在实践中，在日常使用中去发现、去摸索、去掌握。而
52、这次课程设计则给我们提供了这样的一个很好的平台。在使用过程中我们深刻体会到其功能的丰富与强大。作为我们这个专业的学生，除了学习好规定的理论知识外，更重要的就是努力将自己培养成实用性的人才，努力学习使用专业相关的软件知识。同时，通过这次课程设计我们对变压器的特性有了更进一步的认识。这次课程设计我主要通过MATLAB以m文件方式对变压器磁路畸变电流以及变压器负载运行时的效率问题进行了仿真。另外通过查阅相关资料，本次课程设计也做了一些关于励磁电流的仿真。通过这次实验我对于变压器的认识有了更进一步的认识，由电机书上的理论的平面知识，逐渐开始了解其运行的工作特性和运作方式。参考文献1亮.MATLAB语言
53、与控制系统仿真.工业大学.2胡虔生.胡敏强.电机学.中国电力.3薛定宇.控制系统计算机辅助设计MATLAB语言及应用M.:清华大学,1996.社,4垛生.电力拖动自动控制系统M.:高等教育,1998.致此次课程设计虽然说持续的时间，但为了完成此次课程设计，我们都花费了很大的心思和时间，所以我首先相感的是他们的支持和努力！同时，在利用m文件编程过程中，由于之前并没有好好的学习去了解，在实验的过程中遇到了很多这样那样的问题，但每一次都得到了建华老师的悉心讲解，所以在这里同样需要感他在此次课课程设计中给予我们小组的指导！另外，我还想感其他同学在实践中给我们提供的建议与指导！我相信，没有以上这些人的帮助我是不会如此顺利的完成此次课程设计的。
