APF matlab仿真建模要点

APF matlab仿真建模要点
APF matlab仿真建模要点

电力电子系统建模与仿真

学院:电气工程学院

年级:2012级

学号:12031236

姓名:周琪俊

指导老师:舒泽亮

二极管钳位多电平APF电压平衡SPWM仿真报告

1 有源电力滤波器的发展及现状

有源电力滤波器的发展最早可以追溯到20 世纪60 年代末,1969 年B.M.Bird 和J.F.Marsh发表的论文中,描述了通过向电网注入三次谐波电流来减少电源电流中的谐波成分,从而改善电源电流波形的新方法,这种方法是APF 基本思想的萌芽。1971年日本的H.Sasaki 和T.Machida 首先提出APF 的原始模型。1976 年美国西屋电气公司的L.Gyugyi 等提出了用PWM 变流器构成的APF 并确立了APF 的概念。这些以PWM 变流器构成的APF 已成为当今APF 的基本结构。但在70 年代由于缺少大功率的快速器件,因此对APF 的研究几乎没有超出实验室的范围。80 年代以来,随着新型电力半导体器件的出现,脉宽调制的发展,以及H.Akagi 的基于瞬时无功功率理论的谐波电流瞬时检测方法的提出,APF有了迅速发展。

现在日本、美国、德国等工业发达国家APF已得到了高度重视和日益广泛的应用。由于理论研究起步较早,目前国外有源电力滤波器的研究已步入工业化应用阶段。随着容量的逐步提高,其应用范围也从补偿用户自身的谐波向改善整个电网供电质量的方向发展。有源电力滤波器的工业化应用对理论研究起了非常大的推动作用,新的理论研究成果不断出现。1976 年美国西屋公司的L.Gyugyi 率先研制出800kV A的有源电力滤波器。在此以后的几十年里,有源电力滤波器的实践应用得到快速发展。在一些国家,已经投入工业应用的有源电力滤波器容量已增加到50MV A。目前大部分国际知名的电气公司如西屋电气、三菱电机、西门子和梅兰日兰等都有相关的部门都已有相关的产品。

我国在有源电力滤波器的研究方面起步较晚,直到20 世纪80 年代末才有论文发表。90 年代以来一些高等院校和科研机构开始进行有源电力滤波器的研究。1991 年12 月由华北电科院、北京供电局和冶金部自动化研究所研制的国内第一台400V/50kV A 的有源电力滤波器在北京某中心变电站投运,2001 年华北电科院又将有源电力滤波器的容量提高到了10kV/480kV A。由中南大学和湖南大学研制的容量为500kV A 并联混合型有源电力滤波器已在湖南娄底早元220kV 变电站挂网运行。在近几年国内的有源电力滤波器产品已有很多应用,本文研制的两种APF都已应用于工业现场。

2 二极管箝位式多电平逆变器

自从日本学者南波江章于1980 年提出三电平中性点箝位逆变器以来,多电平逆变器的拓扑结构就受到人们的普遍关注,很多学者相继提出了一些实际应用

性强的多电平电路结构,主要有箝位式、级联式、层叠式等多电平逆变器,其中箝位式又包含二极管箝位式、飞跨电容箝位式和混合箝位式等结构,本文研究的对象为二极管箝位式多电平逆变器。二极管箝位式多电平逆变器的显著特点是采用二极管对相应的开关管进行箝位,利用不同的开关状态组合得到不同的输出电压电平数。假定输出电压的电平数为m,相位数为n,则直流支撑电容的个数为m-1,开关管个数为 2 n (m-1),箝位二极管的个数为 2 n (m-2)。以三相五电平为例,直流分压电容的个数为4,开关管的个数为24,箝位二极管的个数为18 个。三相二极管箝位式五电平逆变器的主电路如图1-1 所示。

图1-1 三相二极管箝位式五电平逆变器主电路

对于多电平变换器来说,若其电平数为M,则它的直流侧需要(M-1)个分压电容,输出相电压的电平数为M,输出的线电压电平数为(2M-1)。二极管箝位位式三相五电平逆变器结构如图1-1所示。

图1-1中E为直流侧电源,C1,C2,C3,C4为个直流侧箝位电容,把直流侧电压分为五个部分。图中S41、S42…S47、S48共8个IGBT串联组成一个桥臂,S51、S52…S57、S58这8个IGBT组成一个桥臂,S61、S62…S67、S68组成一个桥臂,这三个桥臂的中点引出变换器交流侧的三相电压。

从图中可以看出三相五电平变换器电路的每一个桥臂有8个IGBT组成,这8个IGBT又可以分为4对对管,例如第一个桥臂中S41和S45,S42和S46,

S43和S47,S44和S48。每一个对管中的两个IGBT不能同时导通,否则会造成短路,正常工作时对管的开关状态互补。下面以第一个桥臂为例研究变换器多电平输出时各个开关的状态。

以表1-1为例,每一个桥臂输出分为V0,V1,V2,V3,V4时,对应的8个IGBT开关状态。表中“1”表示导通状态,“0”表示关断状态。

表1-1 输出电平电压和开关管的状态

开关状

SV1 SV2 SV3 SV4 SV5 SV6 SV7 SV8 态

V0 0 0 0 0 1 1 1 1 V1 0 0 0 1 1 1 1 0 V2 0 0 1 1 1 1 0 0 V3 0 1 1 1 1 0 0 0 V4 1 1 1 1 0 0 0 0

3 多电平变换器的PWM控制方法

多电平变换器脉宽控制技术(Pulse width Modulation,PWM)是用一种以正弦波参考波作为调制波,以N倍于调制波频率的三角波为载波,由于三角波的上下限是线性变化的,所以通过调制波与载波进行比较,调制波大于载波的部分可以得到一组幅值相等,宽度正比于调制波的矩形序列脉冲,用开关量取代模拟量,通过开关管的通断,把直流电能变换成交流电能。我们通常把这种控制技术称为PWM控制技术。

多电平变换器的PWM控制方法主要分为三个大类:载波调制PWM控制法,空间电压相量调制(SVPWM),消除特定谐波PWM控制法。载波调制法又分为载波移相、开关频率优化、阶梯波EPWM、载波层叠法和分段载波层叠法等五种。不同的电路结构和要求,就需要不同的PWM控制法。空间相量控制法不适合于五电平以上的多电平逆变器,以为此时电路会非常复杂。二极管箝位式电路多采用载波层叠法和开关频率优化法,下面我们以载波层叠法来完成二极管箝位式的PWM控制。

对M电平变换器来说,利用(M-1)个频率相同,幅值相等的三角波与一个正弦波进行比较,(M-1)个三角波对称分布在参考量的正负两侧根据正弦调制波与各个三角波的比较结果输出不同的电平,并决定相应的开关管的开关状态。以五电平为例,就需要4个频率、幅值相同的三角载波和一个正弦波进行比较,4个三角波的总幅值要大于等于正弦波的幅值,否则载波层叠PWM法就没有意义。

下图3-1所示,以五电平为例,正弦波与各个三角波进行比较时,当正弦波的幅值大于某个三角波的幅值时,就令相应的IGBT管导通,否则关断。V0,V1,V3,V4为4个频率,幅值完全相同的三角载波,电压依次升高,但中间不间断,V2为正弦调制波。

图3-1 五电平层叠PWM

对IGBT编号为1和5的开关管进行PWM控制,三角波V4与调制正弦波相比较,当正弦波电压高于三角波时,PWM1输出高电平,IGBT管S41导通,相应的S41的对管S45截止。反之,当三角波高于正弦波时,PWM1输出低电平,开关管S41截止,对管S45导通。从而完成PWM控制。下图3-2即为开关管S41和S45的PWM控制。

图3-2 正弦波与三角波进行比较

对IGBT编号为2的开关管进行PWM控制,三角波V3与正弦波进行比较,当正弦波电压高于三角波时,PWM2输出高电平,IGBT管S42导通,相应的S42的对管啥S46截止。反之,当三角波高于正弦波时,PWM2输出低电平,管S42截止,对管S46导通。

对IGBT编号为3和4的开关管进行PWM控制,三角波V1和V0分别与

正弦波进行比较,当正弦波电压高于三角波时,输出高电平,IGBT 管导通,相应它们的对管截止。反之,当三角波高于正弦波时,输出低电平。通过上述过程,从而完成一个桥臂的PWM 控制,图3-3为PWM 控制状态和对应的输出波形图。

图3-3 一个周期的开关状态及输出的电压波形

下图3-4为STATCOM 仿真时的A 相层叠PWM 法控制波形图。

0.1450.150.155

0.160.1650.17-1-0.8

-0.6

-0.4

-0.2

00.20.4

0.6

0.8

1

Time(s)层叠P W M 调制(V )

图3-4 五电平STA TCOM 仿真单相层叠PWM 调制波形

在层叠PWM 控制法中,三角波的频率与幅值通常是固定的,正弦波的幅值与三角波信号的幅值之比称为幅值调制比,用M 表示,通称为调制度。M 的值在0-1之间时,逆变器输出电压的基波分量与M 成线性关系,当M 大于1时,脉冲宽度就不会按照正弦规律进行变化,这种情况称为过调制。

4 直流侧电容电压不平衡的原因及危害

二极管箝位式STATCOM 的主要组成部分为三相多电平变换器,变换器的交流侧通过连接电抗器或者变压器与电网相连,直流侧根据电平数M 的数值,连接M-1个电容,它们起到承载逆变器输入的纹波电流,支撑直流侧电压的作用,多电平的输出就是通过连通不同的电容,达到输出的电压不同。直流侧电容的另一个作用是为电压型变换器提供一个稳定的直流电压源,为系统提供稳定的直流电压和少量的有功功率,以补偿系统的有功损耗。

4.1 直流侧电容电压不平衡的原因

当二极管箝位式多电平变换器在传递有功功率时,每个直流侧电容充放电的时间会有所差异,这就造成它们之间的电压不平衡,从而导致输出的电压产生畸变,甚至得不到期望的输出电平数。所以,怎么解决直流侧电容电压不平衡问题,是多电平STATCOM 的广大研究者必须面对的。

下图2-12即为二极管箝位式五电平交直交电路,a u 、b u 、c u 为三相交流电

源, C1,C2,C3,C4为直流侧电容,其左侧为整流电路的三个桥臂,电容的右侧为逆变电路,电流in i 和out i 为流进流出电容的电流,通过这几个电流的流向

讨论直流侧电容的电压值。 S41S42

S43S44S45S46S47S48D41D42D43D45D44D46S51S52S53S54S55S56S57S58D51D52D53D55D54D56S61S62S63S64S65S66S67S68

D61D62D63

D65D64D66

L1整流电路C1

C2C3C4S11S12S13S14S15S16S17S18D11D12D13D15D14D16S21S22S23S24S25

S26

S27S18D11D22D13D25D24D26S31

S32S33S34S35S36S37S38D31D32D33D35D34D36L2L3

直流侧电容M Load i in5

i in4 i in3 i in2 i in1 i out5 i out4 i out3 i out2

i out1 逆变电路

u a u b

u c u a

V5V4V3V2V1图4.1 二极管箝位式五电平交直交电路

有功情况下电压和电流同相位,电流5in i 工作时图4-2,电压和电流都为正,

电压V5升高,无功时电流和电压相差90度,在输出电压为V5的时间内,前半部份电流为正,后半部分电压为负,相互抵消,所以V5的值不变。电流4in i 工作,有功时,电压和电流也同时为正,电压V4升高。无功时V4电压不变。电流3in i 工作,有功时出于电压和电流也是一半为正,一半为负,所以V3不变。同此,电流2in i 、1in i 工作,有功时V2,V1降低,无功时V2,V1不变。

图4-2 电流5in i 和电压的关系

纯有功情况下分析,电流和电压同相位。

图4-3 有功情况下电压和电流的关系

纯无功情况下:由图4-4可以看出无功情况下,电流和电压相位一直相差90度。

图4-4 无功状态下电流和电压的关系

5 APF控制策略基本思想

对并联型APF 的控制策略进行讨论。APF的控制关键就是对输出电流及直流侧电压进行控制。电压外环和电流内环的双闭环串级控制结构是一种常用的控制策略。双闭环串级控制的主要特点是物理意义清晰,控制结构简单,控制性能优良,在同步旋转坐标系下可以很方便地引入电流状态反馈实现d、q 轴电流的解耦控制。解耦后的控制系统和直流电机双闭环调速系统非常相似。由于双闭环控制在电力电子及其它工业领域中都已得到广泛应用,其控制器参数的工程化整定方法已趋成熟,所以双闭环串级系统的控制器设计几乎可以完全借用这种工程化设计方法。

由于APF 的主要应用是补偿非线性负载的谐波电流,因此其输出的电流为高次谐波电流,在双闭环控制策略中,决定APF输出电流精度的主要因素是谐波检测算法、电流调节器和PWM策略。

5.1 谐波检测算法及其对直流侧电压的影响

并联型APF的一个关键技术就是补偿指令电流的快速精确的获取。谐波检测算法的有效性将直接影响到APF 的补偿性能,根据补偿的目标不同,并联型APF 可以进行单次,多次或全部谐波电流补偿、无功功率补偿以及不平衡电流补偿等。因此要求谐波检测算法从采集到的电网侧或负载侧电流中,能有效迅速的分离出谐波电流分量、无功电流分量或不平衡电流分量。APF 发展到现在,常用的谐波检测算法大体上有以下几种:

①、基于频域分析模拟实现的带通、带阻滤波器或陷波器检测法带通、带阻滤波器或陷波器从采集到的电流信号自己中分离出某一预定频率分量,是用模

拟的方法实现频域分析的一种方法,也是最早被采用的谐波检测算法。通常,通过带通滤波器获得基波电流分量,再与被检测电流相减得到谐波分量,或采用陷波器滤除基波电流分量,得到谐波电流分量。

该检测方法的优点在于电路结构简单、造价低廉、品质因数易于控制。但缺点也较多,比如滤波器的中心频率对元件参数十分敏感等,难设计,误差大,目前已极少采用。

②、基于频域分析的离散傅立叶变换检测法。利用离散傅立叶变换(Discrete Fourier Transform, DFT),对一系列的采样值进行处理,则可以计算出每一个指定的谐波分量。该方法建立在傅立叶分析的基础上,因此要求被补偿的波形是周期变化的,否则会带来较大的误差。通过DFT 将检测到的一个周期的谐波信号进行分解,得到各次谐波的幅值和相位系数,将拟抵消的谐波分量通过带通滤波器或傅立叶变换得到所需的误差信号,再将该误差信号进行DFT 反变换,即可得到补偿信号。

该算法的优点是可以选择拟消除的谐波次数,受电网的影响也比较小。但是这种方法需要进行FFT变换及其反变换,计算量非常大,有较大的时间延迟。而且当电网波形畸变严重或者频率波动时,将引入较大的非同步采样误差,对谐波电流的检测精度影响很大。

③、基于Fryze 时域分析的有功电流检测方法。该方法的基本原理是将负载电流分解为两个正交分量:一个是与电网电压波形完全一致的电流分量,称为有功电流分量;另外一个分量为负载电流与有功电流的差值,包含基波无功和谐波,称为广义无功电流分量。

这种方法的主要缺点是必须计算负载的有功功率和电网电压的有效值,用该方法计算出广义无功电流瞬时值至少有一个周期以上的时间延迟,故不适用于频繁变化负载的补偿。

④、基于采样保持原理的谐波检测算法。此方法将负载电流经过一个带通滤波器得到相电流基波瞬时值,经过输入采样和保持电路得到相应的数字信号。采样保持电路与相电压峰值同步,获得的直流信号正比于电流有功分量幅值。用同样的方法得到另外两相的直流信号。根据负载的有功功率,对这三相直流信号进行平均,平均后的直流信号与三个对称的正弦参考波形相乘,可以获得各相对称基波有功电流,然后它们与实际负载电流相减就得到了所需的补偿电流值。

该检测方法得到的检测电路既可以用于谐波抑制、补偿功率因数;又能平衡系统三相功率

⑤、基于瞬时无功功率理论的检测方法,三相电路瞬时无功功率理论,在谐波和无功电流的实时检测方面得到了成功的应用。目前APF中,基于瞬时无功

功率理论的谐波和无功电流检测方法应用最多。三相瞬时无功功率理论首先由日本学者H.Akagi 于1984 年提出,此后经不断研究逐渐得到完善。现已包括p-q 法、ip- iq法和d-q 法。其中p-q 法最早应用,仅适用于对称三相且无畸变的市电电网;ip- iq 法不仅对电源电压畸变有效,而且在不对称三相市电电网的检测中,相对于p-q法来说,检测误差要小一些;基于Park 变换的d-q法,不仅简化了电网对称无畸变的电流检测,而且也适用于不对称、有畸变的市电电网检测。瞬时无功功率理论相对于传统理论中的有功功率、无功功率在平均值的基础上定义,只适用于电压、电流的波形都是正弦形式的情况。在瞬时无功功率理论中,有功功率、无功功率都是在瞬时值的基础上定义的。

因此,它不仅适用于正弦波,也适用于非正弦波和任何过渡过程。此外检测谐波电流时,因被检测对象电流中谐波构成和所采用的低通滤波器的不同,会有不同的延时,同时谐波电流检测的有效性很大程度上受低通滤波器性能的好坏的影响。

⑥、自适应检测法

该方法是根据信号处理技术中的自适应干扰对消的原理发展起来的。自适应检测法对于市电电压畸变,频率偏移以及市电参数变化有较好的自适应调整能力,可以应用于三相或者单相电路。

⑦、基于神经网络控制法

该方法是随着神经网络理论在系统中的应用而发展起来的一种新型智能控制检测手段。人工神经网络自学功能性强,将进化算法和反向传播用于神经网络的训练,避免了对于给定补偿电流的复杂计算。分析上述这些常用的谐波检测算法,可以将其分为下列两类算法:基波提取法(间接提取法):从负载电流il中提取出基波正序有功分量ilb,使补偿电流if等于il-ilb。常用的算法有基于频域分析模拟实现的带通、带阻滤波器或陷波器检测法,基于瞬时无功功率理论的检测方法,同步旋转坐标变换检测算法等。谐波直接提取法:从负载电流il中直接提取出需要补偿的谐波、无功或不平衡分量。常用的算法有基于频域分析的离散傅立叶变换检测法以及广义同步旋转坐标变换算法等。

5.2 同步旋转坐标变换检测算法

同步旋转坐标变换检测算法的特点是不仅简化了电网电压对称、无畸变下的电流检测,而且也适用于不对称、有畸变的电网电压情况。设三相负载电流为:

式中ilk+、ilk-、ilk0 分别为第k 次谐波的正序、负序和零序分量的幅值。经过式(2-23)的恒模长同步旋转d-q-0坐标变换,可得式(5-2)。由式(5-2)可以看出,三相负载电流经同步旋转d-q-0 坐标后:正序基波分量为直流;其它正序分量则降低一阶次,负序分量则升高一阶次。例如正序7次谐波在同步旋转坐标系下为6次谐波,而负序基波为2次谐波。

基于同步旋转d-q-0坐标系下的检测方法如图5.1所示,负载电流ila、ilb、ilc经过变换后得到ild、ilq、il0。其中,ild 经过低通滤波器LPF,可以获得d 轴直流分量(即基波有功分量);同样,ilq 经过低通滤波器LPF,可以获得q 轴直流分量(即基波无功分量)。再经过d-q-0到a-b-c的反变换即可得到谐波指令电流ira、irb、irc。当开关S分开时,谐波指令电流包含有基波无功分量,APF系统将补偿负载电流中的无功分量;当开关S合上时,APF 系统将不补偿负载电流中的无功分量。

6 APF系统仿真

一个带有三电容辅助稳压电路的APF仿真电路,用于验证辅助稳压电路对多电平变换器直流侧电容电压的均压作用。

APF仿真主要参数如下:

1) 容量为200KV A;

2) 电源额定电压1000V,频率50Hz;

3) 额定电流为300A;

4) 直流侧电压为2200V;

5) 层叠PWM控制频率5KHz。

其二极管钳位多电平APF主电路系统仿真模型如下图6-1所示:

图6-1 二极管钳位多电平APF主电路系统仿真模型

其中,五电平变换器的直流侧支撑电容值均为C=2e-2F。

系统设定直流侧电压值为2200V,0.08s是负载电路负载参数发生变化,但是在很短的时间内,电路迅速恢复到正常值2200V。辅助平衡电路对直流侧电压没有影响,如果直流侧电压发生变化或者出现大的波动,就会影响整个APF 的性能。直流侧电容电压值如图6-2所示:

图6-2 直流侧电容电压值如

其五电平变换器稳压电容仿真电压平衡结果如下图6-3所示:

图6-3 五电平变换器稳压电容仿真电压平衡

图6-3中给出四个均压电容的实际仿真值,基本都维持在Udc/4=550V左右,由此可知,辅助稳压电路实现了电容电压平衡。

二极管钳位多电平APF主电路系统网测电压Vabc,电流Iabc,可控逆变侧电流Ic,及不控整流负载端电流Iioad,如下图6-4所示:

图6-4 系统网测电压Vabc,电流Iabc,可控逆变侧电流Ic,

及不控整流侧负载侧电流Iload

分别对主电路系统网测电压Vabc,电流Iabc,可控逆变侧电流Ic,及不控整流侧负载侧电流Iload,进行THD值分析,分析结果如下图6-5所示:

图6-5a 网侧电压V abc的THD值

图6-5b 网侧电流Iabc的THD值

图6-5c 可控逆变侧电流Ic的THD值

图6-5d 不控整流负载端电流Iload的THD值

通过对二极管钳位多电平APF主电路系统仿真,分别得到网测电压Vabc的THD值为0.54%,网侧电流Iabc的THD值为1.50%,可控逆变侧电流Ic的THD 值为1.85%,及不控整流负载端电流Iioad的THD值为1.17%,主要存在3次,11次,15次谐波畸变比较严重,其系统网测电压Vabc,电流Iabc,可控逆变侧电流Ic,及不控整流侧负载侧电流Iload四个主要THD指标均在标准范围之内,验证此二极管钳位多电平APF控制系统是可行的。

总结和致谢:

针对逆变器直流侧电容电压的不平衡现象,本文提出了一种基于电容的辅助稳压电路来实现电压平衡的措施,给出了一种辅助电路的结构和控制方法,仿真和实验结果显示,这种辅助稳压电路完全能够使直流侧电容电压保持均等。在此基础上搭建了五电平APF仿真系统,并完成了具有辅助稳压电路的五电平APF 的仿真。仿真结果都表明了:辅助稳压电路能够有效实现直流侧电容的均衡,五电平APF能够正常工作。

借此机会,感谢我的导师舒泽亮教授在这一学期的谆谆教导,使我作为一个Matlab仿真软件零基础的学生,在18周的课程学习当中,不仅学习了专业的电力电子仿真技术,还学会熟练的运用simulink进行学习,和其他科目的仿真应用。课堂上丰富的课程内容,便于同学吸收应用的专业知识,以及很有专业针对性的模块仿真。舒老师对我们的难题解答,也是耐心之至,直到听懂为止,使我对电力电子建模与仿真这门课程有了更深层次的认识,同时也对这门课程产生了浓厚的兴趣。

参考文献

[1] 王兆安,杨君,刘进军.谐波抑制和无功功率补偿[M]. 机械工业出版社,2002.

[2] 林中岳. 现代电力电子应用技术[M]. 科学出版社,2007.

[3] 刘风君. 多电平逆变技术及其应用[M]. 机械工业出版社,2007.

[4] 李永东肖曦高跃大容量多电平变化器[M].科学出版社,2005.

[5] 游小杰,李永东,Victor Valouch,等.并联型有源电力滤波器在非理想电源

电压下的控制[J].中国电机工程学报,2004,24(2):55-60.

[6] 廖志凌,刘国海,贾洪平.应用HPF同时检测高次谐波和无功电流方法的讨

论[J].电力系统自动化,2002,26(23):67-68.

[7] 王广柱. 二极管箝位式多电平逆变器直流侧电容电压不平衡机理的研究[J].

中国电机工程学报,2002,12.

[8] 王广柱,洪春梅多电平逆变器直流侧电容电压的平衡与控制[J].电力系统自

动化,2002,11.

[9] 舒泽亮,电力电子技术建模与仿真.课件

现代信号处理Matlab仿真——例611

例6.11 利用卡尔曼滤波估计一个未知常数 题目: 设已知一个未知常数x 的噪声观测集合,已知噪声v(n)的均值为零, 方差为 ,v(n)与x 不相关,试用卡尔曼滤波估计该常数 题目分析: 回忆Kalman 递推估计公式 由于已知x 为一常数,即不随时间n 变化,因此可以得到: 状态方程: x(n)=x(n-1) 观测方程: y(n)=x(n)+v(n) 得到A(n)=1,C(n)=1, , 将A(n)=1,代入迭代公式 得到:P(n|n-1)=P(n-1|n-1) 用P(n-1)来表示P(n|n-1)和P(n-1|n-1),这是卡尔曼增益表达式变为 从而 2v σ1??(|1)(1)(1|1)(|1)(1)(1|1)(1)()()(|1)()[()(|1)()()]???(|)(|1)()[()()(|1)](|)[()()](|1)H w H H v x n n A n x n n P n n A n P n n A n Q n K n P n n C n C n P n n C n Q n x n n x n n K n y n C n x n n P n n I K n C n P n n --=----=----+=--+=-+--=--2()v v Q n σ=()0w Q n =(|1)(1)(1|1)(1)()H w P n n A n P n n A n Q n -=----+21 ()(|1)[(|1)]v K n P n n P n n σ-=--+22(1)()[1()](1)(1)v v P n P n K n P n P n σσ-=--=-+

APF matlab仿真建模要点

电力电子系统建模与仿真 学院:电气工程学院 年级:2012级 学号:12031236 姓名:周琪俊 指导老师:舒泽亮

二极管钳位多电平APF电压平衡SPWM仿真报告 1 有源电力滤波器的发展及现状 有源电力滤波器的发展最早可以追溯到20 世纪60 年代末,1969 年B.M.Bird 和J.F.Marsh发表的论文中,描述了通过向电网注入三次谐波电流来减少电源电流中的谐波成分,从而改善电源电流波形的新方法,这种方法是APF 基本思想的萌芽。1971年日本的H.Sasaki 和T.Machida 首先提出APF 的原始模型。1976 年美国西屋电气公司的L.Gyugyi 等提出了用PWM 变流器构成的APF 并确立了APF 的概念。这些以PWM 变流器构成的APF 已成为当今APF 的基本结构。但在70 年代由于缺少大功率的快速器件,因此对APF 的研究几乎没有超出实验室的范围。80 年代以来,随着新型电力半导体器件的出现,脉宽调制的发展,以及H.Akagi 的基于瞬时无功功率理论的谐波电流瞬时检测方法的提出,APF有了迅速发展。 现在日本、美国、德国等工业发达国家APF已得到了高度重视和日益广泛的应用。由于理论研究起步较早,目前国外有源电力滤波器的研究已步入工业化应用阶段。随着容量的逐步提高,其应用范围也从补偿用户自身的谐波向改善整个电网供电质量的方向发展。有源电力滤波器的工业化应用对理论研究起了非常大的推动作用,新的理论研究成果不断出现。1976 年美国西屋公司的L.Gyugyi 率先研制出800kV A的有源电力滤波器。在此以后的几十年里,有源电力滤波器的实践应用得到快速发展。在一些国家,已经投入工业应用的有源电力滤波器容量已增加到50MV A。目前大部分国际知名的电气公司如西屋电气、三菱电机、西门子和梅兰日兰等都有相关的部门都已有相关的产品。 我国在有源电力滤波器的研究方面起步较晚,直到20 世纪80 年代末才有论文发表。90 年代以来一些高等院校和科研机构开始进行有源电力滤波器的研究。1991 年12 月由华北电科院、北京供电局和冶金部自动化研究所研制的国内第一台400V/50kV A 的有源电力滤波器在北京某中心变电站投运,2001 年华北电科院又将有源电力滤波器的容量提高到了10kV/480kV A。由中南大学和湖南大学研制的容量为500kV A 并联混合型有源电力滤波器已在湖南娄底早元220kV 变电站挂网运行。在近几年国内的有源电力滤波器产品已有很多应用,本文研制的两种APF都已应用于工业现场。 2 二极管箝位式多电平逆变器 自从日本学者南波江章于1980 年提出三电平中性点箝位逆变器以来,多电平逆变器的拓扑结构就受到人们的普遍关注,很多学者相继提出了一些实际应用

MATLAB Simulink系统建模与仿真 实验报告

MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名:****** 专业:电气工程及其自动化 班级:******************* 学号:*******************

实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件,点击Simulink模型构建,根据电路原理图,添加下列模块: (1)无穷大功率电源模块(Three-phase source) (2)三相并联RLC负荷模块(Three-Phase Parallel RLC Load) (3)三相串联RLC支路模块(Three-Phase Series RLC Branch) (4)三相双绕组变压器模块(Three-Phase Transformer (Two Windings)) (5)三相电压电流测量模块(Three-Phase V-I Measurement) (6)三相故障设置模块(Three-Phase Fault) (7)示波器模块(Scope) (8)电力系统图形用户界面(Powergui) 按电路原理图连接线路得到仿真图如下: 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV,相角0°,频率50Hz,接线方式为中性点接地的Y形接法,电源电阻0.00529Ω,电源电感0.000140H,参数设置如下图:

1.2.2 变压器模块 变压器模块参数采用标幺值设置,功率20MVA,频率50Hz,一次测采用Y型连接,一次测电压110kV,二次侧采用Y型连接,二次侧电压11kV,经过标幺值折算后的绕组电阻为0.0033,绕组漏感为0.052,励磁电阻为909.09,励磁电感为106.3,参数设置如下图: 1.2.3 输电线路模块 根据给定参数计算输电线路参数为:电阻8.5Ω,电感0.064L,参数设置如下图: 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号,相当于电压、电流互感器的作用,勾选“使用标签(Use a label)”以便于示波器观察波形,设置电压标签“Vabc”,电流标签“Iabc”,参数设置如下图:

Matlab仿真实例-卫星轨迹

卫星轨迹 一.问题提出 设卫星在空中运行的运动方程为: 其中是k 重力系数(k=401408km3/s)。卫星轨道采用极坐标表示,通过仿真,研究发射速度对卫星轨道的影响。实验将作出卫星在地球表面(r=6400KM ,θ=0)分别以v=8KM/s,v=10KM/s,v=12KM/s 发射时,卫星绕地球运行的轨迹。 二.问题分析 1.卫星运动方程一个二阶微分方程组,应用Matlab 的常微分方程求解命令ode45求解时,首先需要将二阶微分方程组转换成一阶微分方程组。若设,则有: 2.建立极坐标如上图所示,初值分别为:卫星径向初始位置,即地球半径:y(1,1)=6400;卫星初始角度位置:y(2,1)=0;卫星初始径向线速度:y(3,1)=0;卫星初始周向角速度:y(4,1)=v/6400。 3.将上述一阶微分方程及其初值带入常微分方程求解命令ode45求解,可得到一定时间间隔的卫星的径向坐标值y(1)向量;周向角度坐标值y(2)向量;径向线速度y(3)向量;周向角速度y(4)向量。 4.通过以上步骤所求得的是极坐标下的解,若需要在直角坐标系下绘制卫星的运动轨迹,还需要进行坐标变换,将径向坐标值y(1)向量;周向角度坐标值y(2)向量通过以下方程转换为直角坐标下的横纵坐标值X,Y 。 5.卫星发射速度速度的不同将导致卫星的运动轨迹不同,实验将绘制卫星分别以v=8KM/s ,v=10KM/s ,v=12KM/s 的初速度发射的运动轨迹。 三.Matlab 程序及注释 1.主程序 v=input('请输入卫星发射速度单位Km/s :\nv=');%卫星发射速度输入。 axis([-264007000-1000042400]);%定制图形输出坐标范围。 %为了直观表达卫星轨迹,以下语句将绘制三维地球。 [x1,y1,z1]=sphere(15);%绘制单位球。 x1=x1*6400;y1=y1*6400;???????-=+-=dt d dt dr r dt d dt d r r k dt r d θ θθ2)(2 22222θ==)2(,)1(y r y ?????????????**-=**+*-===)1(/)4()3(2)4()4()4()1()1()1()3()4()2() 3()1(y y y dt dy y y y y y k dt dy y dt dy y dt dy ???*=*=)] 2(sin[)1(Y )]2(cos[)1(X y y y y

基于MATLAB的变压器仿真 与分析

于MATLAB_Simulink的牵引变压器建模与仿真 基于MATLAB/Simulink的牵引变压器建模与仿真徐(西安铁路局安康供电段新陕西汉中 723000)摘要:针对多种牵引变压器接线方式,建立数学模型,基于Matlab/Simulink仿真软件,建立牵引变压器的仿真模型,并验证数学模型和仿真模型的一致性。利用所建立仿真模型对不同接线形式牵引变压器在不同条件下对公用电网产生的谐波和负序影响进行仿真试验,对研究各种类型的牵引变压器特性在我国电气化铁路的应用提供条件。关键词:牵引变压器;数学模型;仿真模型;Matlab/Simulink 中图分类号:U223.6 文献标识码:A 文章编号:1671-7597(2011)0610061-03 牵引变压器按其特性可分为平衡接线和不平衡接线。其中不平衡接线有单相接线、Vv接线和YNd11接线;平衡接线是试图实现三相两相对称变换而提出的,主要代表方式有Scott,Leblanc、Kubler、Wood-bridge、阻抗匹配接线等。本次主要总结了常用牵引变压器的特点并建立数学模型,包括每种牵引变压器的原理结构、原次边电气量关系等,基于Matlab/Simulink软件建立牵引变压器仿真模型,并对牵引变压器在不同条件下的负序、谐波特性的进行了研究. 1 牵引变压器数学模型研究 1.1 YNd11接线 YNd11变压器接线原理如下图所示,如果忽略激磁电流及其漏阻抗压降,二次侧绕组ac相与一次侧绕组A相同相,cb相与C相同相。由于变压器一次侧绕组A,B,C相与电力系统的相序一致,A相滞后C相,对应的二次侧ac也滞后cb相[2]。其中Z为牵引端口对应变压器漏抗,和β相的端口电压。 1.2 Vv接线 Vv接线牵引变压器接线原理如图2所示。为二次侧空载相即α相图2 Vv接线牵引变压器设Vv接线变压器一次侧、二次侧绕组匝数分别为可得电流输入输出关系[3]:和,电压输入输出关系如下:图1 YNd11接线牵引变压器设YNd11接线变压器一次侧、二次侧绕组匝数分别为和假设变压器原边中性点接地,可以得出一次侧三相电流。,其中为牵引端口对应变压器漏抗,为二次侧空载相即α相和β相的端口电压。 1.3 Scott接线 Scott接线变压器(又称T形接法变压器)属于能完成三相-两相变换的平衡变压器,Scott接线牵引变压器接线原理如图3所示。图3 Scott牵引变压器接线原理图 1 61 设一次侧绕组BC的匝数为次侧绕组AD的匝数为,记,二次的绕组ad、bc的匝数为,则一。可得电流输入输出关系[4]:把一次侧绕组电流用相电流替换,即为:式中,为从三相端子流进变压器的电流。输出端口电压方程为:图6 YNd11接线牵引变压器两供电臂输出电压波形从电压输出波形中可以得到α供电臂电压波形超前β供电臂电压波形120°,在对称阻性负载下,两臂电流输出波形幅值相同,相位相差120°,满足理论值。 2.2 Vv接线牵引变压器 Vv 接线牵引变压器是由两个单相牵引变压器并联而成,仿真模型如图7所示.在仿真模型中牵引变压器T1和T2的原、次边变比设置为110kV/27.5kV。对,于

倒立摆系统的建模及Matlab仿真资料

第1 页共11 页 倒立摆系统的建模及Matlab仿真 1.系统的物理模型 考虑如图(1)所示的倒立摆系统。图中,倒立摆安装在一个小车上。这里仅考虑倒立摆在图面内运动的二维问题。 图(1)倒立摆系统 假定倒立摆系统的参数如下。 摆杆的质量:m=0.1g l=1m小车的质量:摆杆的长度:2重力加速度:g=9.8m/M=1kg s摆杆的质量在摆杆的中心。 设计一个控制系统,使得当给定任意初始条件(由干扰引起)时,最大超调量?≤10%,调节时间ts ≤4s ,通过小车的水平运动使倒立摆保持在垂直位置。 2.系统的数学模型 2.1建立倒置摆的运动方程并将其线性化。 为简化问题,在数学模型中首先假设:1)摆杆为刚体;2)忽略摆杆与支点之间的摩擦;3)忽略小车与接触面间的摩擦。 ?),在u设小车瞬时位置为z,摆心瞬时位置为(作用下,小车及摆均产生加速远 动,sin?lz根据牛顿第二定律,在水平直线远动方向的惯性力应与u平衡,于是有 22dzd?)?sinu?M?m(zl22dtdt???2????z(M?mml?)cos?mlusin? 即:??①

绕摆轴转动的惯性力矩与重力矩平衡,因而有. 第2 页共11 页 2??d??? sin??lcosm(z?lsinmgl)??2dt?????22???????即: nis?l?ocgcosincoszs?ls??② 以上两个方程都是非线性方程,为求得解析解,需作线性化处理。由于控制的目的是保持倒立摆直?2?????且可忽略则,立,在试驾合适的外力条件下,假定θ很小,接近于零时合理的,1sincos??,项。于是有 ???M?zm?u?ml??)(③ ????g?z?l??④联立求解可得1mg?u?z????MM 1)?m(M????u??MlMl 列写系统的状态空间表达式。2.2??T xx,x,x,,选取系统变量则 xx,x,xx?,42134123xx??211mgux???x?32MM x?x?431)(M?mu?x?x? 34MlMl 即00100????z??1mg??????000?z?????d MM??Bu?Ax?xux????????00001???dt????1gm?(M)????000??????? MlMl??????Cx?0?y?xx1001代入数据计算得到:0100????000?1??????T0D,?0??1BA?,?001,C100??1000??00011?? 11 页3 页共第 3.设计控制器3.1判断系统的能控性和稳定性 1100????0011????23BBAABAB?Q?故被控对象完全可控, rank()=4,Q kk??11?0?10??011?10???22???11?。出现大于零的特征值,故被,,0 解得特征值为 0由特征方程0??11I?A?)(控对象不稳定3.2确定希望的极点, 另一对为远极点,认为系统性能主要由主导,选其中一对为主导极点和希望的极点n=4ss21极点决定,远极点只有微小影响。根据二阶系统的关系式,先确定主导极点???42??1????10.?e??t1.67?有,闭环可得;取误差带,于是取,则6.?059?0.02.?0? pns??n2????1?js??=-10.8j,远极点选择使它和原点的距离大于主导极点与原点 距离主导极点为?n,21s??15倍,取的54,33.3采用状态反馈方法使系统稳定并配置极点 ??kkkk?k;状态反馈系统的状态方程,馈状态反的控制规律为为kxu??3102?,其

MATLAB的建模和仿真

课程设计说明书 题目:基于Matlab的IIR滤波器设计与仿真班级:2012 级电气五班 姓名:王璐 学号:201295014178 指导教师:张小娟 日期:2015年 1 月12日

课程设计任务书

基于MATLAB的IIR滤波器设计与仿真 前言 数字信号处理(digital signal processing,DSP)是从20世纪60年代以来,随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科。数字信号处理是把信号用数字或符号表示的序列,通过计算机或通用(专用)信号处理设备,用数字的数值计算方法处理(例如滤波、变换、压缩、增强、估计、识别等),以达到提取有用信息便于应用处理的目的。数字信号处理系统有精度高、灵活性高、可靠性高、容易大规模集成、时分复用、可获得高性能指标、二维与多维处理等特点。正是由于这些突出的特点,使得它在通信、语音、雷达、地震测报、声呐、遥感、生物医学、电视、仪器中得到愈来愈广泛的应用。在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器(DF,Digital Filter),根据其单位冲激响应函数的时域特性可分为两类:无限冲激响应IIR(Infinite Impulse Response)滤波器和有限冲激响应FIR(Finite Impulse Response)滤波器。MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来结算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的有点,使MATLAB成为一个强大的数学软件,在新的版本中也加入了对C,FORTRAN,C++,JA V A的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。 1 数字滤波器概述 数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波实质上是一种运算过程,实现对信号的运算处理。输入数字信号(数字序列)通过特定的运算转变为输出的数字序列,因此,数字滤波器本质上是一个完成特定运算的数字计算过程,也可以理解为一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则,使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性:Y(eωj)=X(eωj)H(eωj) 其中Y(eωj)、X(eωj)分别是数字滤波器的输出序列和输入序列的频域特性(或称为

MATLAB实现通信系统仿真实例

补充内容:模拟调制系统的MATLAB 仿真 1.抽样定理 为了用实验的手段对连续信号分析,需要先对信号进行抽样(时间上的离散化),把连续数据转变为离散数据分析。抽样(时间离散化)是模拟信号数字化的第一步。 Nyquist 抽样定律:要无失真地恢复出抽样前的信号,要求抽样频率要大于等于两倍基带信号带宽。 抽样定理建立了模拟信号和离散信号之间的关系,在Matlab 中对模拟信号的实验仿真都是通过先抽样,转变成离散信号,然后用该离散信号近似替代原来的模拟信号进行分析的。 【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率,即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样,并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小,要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例 【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小,要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样

三相变压器建模及仿真及MATLAB仿真

XXXXXXX学院课程设计报告 课程名称: 系部: 专业班级: 学生姓名: 指导教师: 完成时间: 报告成绩: 学院教学工作部制

目录 摘要 (3) 第一章变压器介绍 (4) 1.1 变压器的磁化特性 (4) 1.2 变压器保护 (4) 1.3 励磁涌流 (7) 第二章变压器基本原理 (9) 2.1 变压器工作原理 (9) 2.2 三相变压器的等效电路及联结组 (10) 第三章变压器仿真的方法 (11) 3.1 基于基本励磁曲线的静态模型 (11) 3.2基于暂态磁化特性曲线的动态模型 (13) 3.3非线性时域等效电路模型 (14) 第四章三相变压器的仿真 (16) 4. 1 三相变压器仿真的数学模型 (16) 4.2电源电压的描述 (20) 4.3铁心动态磁化过程简述 (21) 第五章变压器MATLAB仿真研究 (25) 5.1 仿真长线路末端电压升高 (25) 5.2 仿真三相变压器 T2 的励磁涌流 (28) 5.3三相变压器仿真模型图 (34) 5.4 变压器仿真波形分析 (36) 结论 (40) 参考文献 (41)

摘要 在电力变压器差动保护中,励磁涌流和内部故障电流的判别一直是一个关键问题。文章阐述了励磁涌流的产生及其特性,利用 MATLAB 对变压器的励磁涌流、内部故障和外部故障进行仿真,对实验的数据波形分析,以此来区分故障和涌流,目的是减少空载合闸产生的励磁涌流对变压器差动保护的影响,提高保护的灵敏性。 本文在Matlab的编程环境下,分析了当前的变压器仿真的方法。在单相情况下,分析了在饱和和不饱和的励磁涌流现象,和单相励磁涌流的特征。在三相情况下,在用分段拟和加曲线压缩法的基础上,分别用两条修正的反正切函数,和两条修正的反正切函数加上两段模拟饱和情况的直线两种方法建立了Yd11、Ynd11、Yny0和Yy0四种最常用接线方式下三相变压器的数学仿真模型,并在Matlab下仿真实现。通过对三相励磁涌流和磁滞回环波形分析,三相励磁涌流的特征分析,总结出影响三相变压器励磁涌流地主要因素。最后,分析了两种方法的优劣,建立比较完善的变压器仿真模型。 关键字: 变压器;差动保护;励磁涌流;内部故障;外部故障;波形分析;仿真;数学模型

基于Matlab、Simulink 的AM通信系统仿真设计与研究

天津理工大学计算机与通信工程学院通信工程专业设计说明书 基于Matlab/Simulink 的AM通信系统仿真设计与研究 姓名杜艳玮 学号 20092177 班级 09通信-2 指导老师赵健 日期2012/12/16

目录 摘要 (3) 第一章前言 (4) 1.1专业设计任务及要求 (4) 1.2 Matlab简介 (4) 1.4 通信系统模型 (6) 第二章 AM调制原理及仿真 (7) 2.1 AM调制原理 (7) 2.1.1 AM介绍 (7) 2.1.2 AM调制原理框图 (8) 2.2 AM调制方式的Matlab仿真 (8) 2.2.1 载波信号分析 (8) 2.2.2 AM调制 (9) 2.3 AM调制方式Matlab-simulink仿真 (10) 2.3.1 仿真框图 (10) 2.3.2 仿真结果 (11) 第三章 AM解调 (13) 3.1 AM解调原理 (13) 3.2 AM解调方式Matlab仿真 (13) 3.2.1 滤波前AM解调信号波形 (13) 3.2.2 AM调制信号解调 (15) 3.3 AM解调方式的Matlab-simulink仿真 (17) 3.3.1 仿真框图 (17) 3.3.2 仿真结果 (18) 第四章结论 (19) 参考文献 (20)

摘要 学习AM调制原理,AM调制就是由调制信号去控制高频载波的幅度,使之随调制信号作线性变化的过程。在波形上,幅度已调信号的幅度随基带信号的规律而呈正比地变化。解调方法利用相干解调。解调就是实现频谱搬移,通过相乘器与载波相乘来实现。通过相干解调,通过低通滤波器得到解调信号。相干解调时,接收端必须提供一个与接受的已调载波严格同步的本地载波,它与接受的已调信号相乘后,经低通滤波器取出低频分量,得到原始的基带调制信号。通过信号的功率谱密度的公式,得到功率谱密度。利用Matlab和Matlab-Simulink仿真建立AM调制的通信系统模型,用Matlab仿真程序画出调制信号、载波、已调信号、相干解调之后信号的波形以及功率频谱密度,分析所设计系统性能。用Matlab-Simulink仿真建立基于相干解调的AM仿真模型,详细叙述模块参数的设置,分析仿真结果。 关键字:AM调制相干解调 Matlab仿真 Matlab-Simulink仿真

倒立摆系统的建模及Matlab仿真

倒立摆系统的建模及Matlab 仿真 1.系统的物理模型 考虑如图(1)面内运动的二维问题。 图(1)倒立摆系统 假定倒立摆系统的参数如下。 摆杆的质量:m=0.1g 摆杆的长度:l =1m 小车的质量: M=1kg 重力加速度:g=9.8m/2s 摆杆的质量在摆杆的中心。 设计一个控制系统,使得当给定任意初始条件(由干扰引起)时,最大超调量δ ≤10%,调节时 间ts ≤4s ,通过小车的水平运动使倒立摆保持在垂直位置。 2.系统的数学模型 2.1建立倒置摆的运动方程并将其线性化。 为简化问题,在数学模型中首先假设:1)摆杆为刚体;2)忽略摆杆与支点之间的摩擦;3)忽略小车与接触面间的摩擦。 设小车瞬时位置为z,摆心瞬时位置为(θsin l z +),在u 作用下,小车及摆均产生加速远动,根据牛顿第二定律,在水平直线远动方向的惯性力应与u 平衡,于是有 u l z dt d m dt z d M =++)sin (22 22θ 即: u ml ml z m M =-++θθθθsin cos )(2&&&&& ① 绕摆轴转动的惯性力矩与重力矩平衡,因而有

θθθsin cos )sin (22mgl l l z dt d m =??? ????+ 即: θθθθθθθsin cos sin cos cos 22g l l z =-+&&&&& ② 以上两个方程都是非线性方程,为求得解析解,需作线性化处理。由于控制的目的是保持倒立摆直 立,在试驾合适的外力条件下,假定θ很小,接近于零时合理的,则1cos ,sin ≈≈θθθ,且可忽略θ θ2&项。于是有 u ml z m M =++θ&&&& )( ③ θθg l z =+&&&& ④ 联立求解可得 u Ml Ml m M u M M mg z 1)(1 -+=+- =θθθ&&&& 2.2列写系统的状态空间表达式。 选取系统变量4321,,,x x x x , []T x x x x x 4321,,,=则 u Ml x Ml m M x x x u M x M mg x x x 1 )(134433221-+= =+-==&&&& 即 []Cx x x y Bu Ax u Ml M x Ml g m M M mg z z dt d x ===+=?????? ? ???????-+?????????? ??? ? +- =???? ????????=000110100)(0 010 0000000 1 1θθ&&& 代入数据计算得到: [][]0,0001,1010,01100 1000010000 1 0==-=? ? ??? ? ??? ???-=D C B A T

四旋翼飞行器建模与仿真Matlab

四轴飞行器的建模与仿真 摘要 四旋翼飞行器是一种能够垂直起降的多旋翼飞行器,它非常适合近地侦察、监视的任务,具有广泛的军事和民事应用前景。本文根据对四旋翼飞行器的机架结构和动力学特性做详尽的分析和研究,在此基础上建立四旋翼飞行器的动力学模型。四旋翼飞行器有各种的运行状态,比如:爬升、下降、悬停、滚转运动、俯仰运动、偏航运动等。本文采用动力学模型来描述四旋翼飞行器的飞行姿态。在上述研究和分析的基础上,进行飞行器的建模。动力学建模是通过对飞行器的飞行原理和各种运动状态下的受力关系以及参考牛顿-欧拉模型建立的仿真模型,模型建立后在Matlab/simulink软件中进行仿真。 关键字:四旋翼飞行器,动力学模型,Matlab/simulink Modeling and Simulating for a quad-rotor aircraft ABSTRACT The quad-rotor is a VTOL multi-rotor aircraft. It is very fit for the kind of reconnaissance mission and monitoring task of near-Earth, so it can be used in a wide range of military and civilian applications. In the dissertation, the detailed analysis and research on the rack structure and dynamic characteristics of the laboratory four-rotor aircraft is showed in the dissertation. The dynamic model of the four-rotor aircraft areestablished. It also studies on the force in the four-rotor aircraft flight principles and course of the campaign to make the research and analysis. The four-rotor aircraft has many operating status, such as climbing, downing, hovering and rolling movement, pitching movement and yawing movement. The dynamic model is used to describe the four-rotor aircraft in flight in the dissertation. On the basis of the above analysis, modeling of the aircraft can be made. Dynamics modeling is to build models under the principles of flight of the aircraft and a variety of state of motion, and Newton - Euler model with reference to the four-rotor aircraft.Then the simulation is done in the software of Matlab/simulink. Keywords: Quad-rotor,The dynamic mode, Matlab/simulink

基于Matlab的电力系统自动重合闸建模与仿真讲解

实践课程设计报告 课程名称:Matlab上机 题目:基于MATLAB的电力系统自动重合闸 所在学院: 学科专业: 学号: 学生姓名: 指导教师: 二零一五年四

摘要 分析了单相自动重合闸的工作特性,并利用MATLAB软件搭建了220kv电力系统的自动重合闸的仿真模型,模拟系统发生单相接地、三相相间短路故障,断路器跳闸后自动重合闸的工作过程。 关键词:电力系统自动重合闸MATLAB 短路故障

目录 1 引言 (1) 2 模型中主要模块的选择和参数 (2) 2.1同步发电机模块 (2) 2.2 变压器模块 (2) 2.3 输电线路模块 (3) 2.3.1 150km线路模块 (3) 2.3.2 100km线路模块 (4) 2.1 电源模块 (5) 2.3 负载模块 (6) 2.3.1 三相串联RLC负载Load1 (6) 2.3.2 三相串联RLC负载Load4 (7) 2.4 断路器模块 (8) 2.5 测量模块 (9) 2.6 显示模块 (9) 2.7 其他模块 (9) 2.8 仿真参数设置 (10) 3 仿真结果及波形分析 (10) 3.1 线路单相重合闸 (10) 3.2 线路三相重合闸 (12) 总结 (13) 参考文献 (14)

基于Matlab的电力系统自动重合闸 1 引言 随着技术的发展,电力系统的规模越来越复杂。从实际条件与安全角度考虑,不太可能进行电力系统科研实验,因而电力系统数字仿真成为了电力系统研究、规划和设计的重要手段。电力系统仿真软件如BPA,EMTP,PSCAD/ EMTDC ,NETOMAC,PSASP,MATLAB等,正向着多功能,具有更高的可移植性方向发展。其中在MATLAB 中,电力系统模型可以在Simulink环境下直接搭建,Simulink电力系统元件库中有多种多样的电气模块,电力系统大多数元件都包含。其中,可以直接调用。电力系统大部分故障是瞬时性故障,因此采用自动重合闸后,电力系统发生瞬时性故障时供电的连续性、系统的稳定性得到很大的提高。此外,自动重合闸有效纠正由于断路器或继电保护误动作引起的误跳闸。 本文以MATLAB为工具,对简单系统的线路单相重合闸和线路三相重合闸进行分析与研究。 1.1 仿真模型的设计和实现 电力系统正常运行时可以认为是三相对称的,即电压、电流对称,且具有正弦波形。下图为理想情况下220kv电力系统的模型。 图 1 220kv电力系统模型

Matlab通信系统建模与仿真例题源代码-第三章

% ch3example1A.m clear; f_p=2400; f_s=5000; R_p=3; R_s=25; % 设计要求指标 [n, fn]=buttord(f_p,f_s,R_p,R_s, 's'); % 计算阶数和截止频率 Wn=2*pi*fn; % 转换为角频率 [b,a]=butter(n, Wn, 's'); % 计算H(s) f=0:100:10000; % 计算频率点和频率范围 s=j*2*pi*f; % s=jw=j*2*pi*f H_s=polyval(b,s)./polyval(a,s); % 计算相应频率点处H(s)的值 figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel('频率Hz');ylabel('幅度dB'); subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabel('频率Hz');ylabel('相角rad'); figure(2); freqs(b,a); % 也可用指令freqs直接画出H(s)的频率响应曲线。 % ch3example1B.m clear; f_p=2400; f_s=5000; R_p=3; R_s=25; % 设计要求指标 [n, fn]=ellipord(f_p,f_s,R_p,R_s,'s'); % 计算阶数和截止频率 Wn=2*pi*fn; % 转换为角频率 [b,a]=ellip(n,R_p,R_s,Wn,'s'); % 计算H(s) f=0:100:10000; % 计算频率点和频率范围 s=j*2*pi*f; % s=jw=j*2*pi*f H_s=polyval(b,s)./polyval(a,s); % 计算相应频率点处H(s)的值 figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel('频率Hz');ylabel('幅度dB'); subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabel('频率Hz');ylabel('相角rad'); figure(2); freqs(b,a); % 也可用指令freqs直接画出H(s)的频率响应曲线。 % ch3example2A.m f_N=8000; % 采样率 f_p=2100; f_s=2500; R_p=3; R_s=25; % 设计要求指标 Ws=f_s/(f_N/2); Wp=f_p/(f_N/2); % 计算归一化频率 [n, Wn]=buttord(Wp,Ws,R_p,R_s); % 计算阶数和截止频率 [b,a]=butter(n, Wn); % 计算H(z) figure(1); freqz(b,a, 1000, 8000) % 作出H(z)的幅频相频图, freqz(b,a, 计算点数, 采样率)

基于Matlab和VR技术的移动机器人建模及仿真

文章编号:100422261(2004)0120039204 基于Matlab 和VR 技术的移动机器人建模及仿真Ξ 葛为民1,2,曹作良2,彭商贤1 (1.天津大学机械工程学院,天津300072;2.天津理工学院机械工程学院,天津300191) 摘 要:利用Matlab 建立移动机器人的动力学模型,在虚拟现实(VR )环境下,实时仿真移动机器人路径跟踪的运动特性,为基于Internet 的机器人遥操作试验搭建了仿真平台.实验结果表明,虚拟模型准确地模拟了真实移动机器人的动力学特征;通过对模型的参数修改,为实现对真实机器人的最优控制和设计提供了可信的参考方案.关键词:Matlab ;虚拟现实;移动机器人;遥操作中图分类号:TP242.2 文献标识码:A Dynamic modeling and simulation of mobile robot based on matlab and VR technology GE Wei 2min 1,2,C AO Zuo 2liang 2,PE NG Shang 2xian 1 (1.School of Mechanical Eng.,T ianjin University ,T ianjin 300072,China ;2.School of Mechanical Eng.,T ianjin Institute of T echnology ,T ianjin 300191,China ) Abstract :This paper proposes an approach that develops a dynam ic m odel of a m obile robot taking advantage of the M atlab.M eantime ,in a developed virtual reality environment ,the built m odel simulates the m otion of path tracking and obstacle av oidance.Furtherm ore ,it provides a platformfor experiments of m obile robot teleoperation.The experi 2mental results approve that ,the virtual m odel represents the dynam ic properties of real robot accurately and ,w ith the change of parameters of the virtual m odel ,it helps to find out the optim ization methods of controlling and designing the m obile robot indeed. K eyw ords :M atlab ;virtual reality ;m obile robot ;teleoperation 在当今工业现代化的高速发展时期,特别是自动化设备在各个领域的广泛应用,移动机器人(AG V )的应用越来越显示出它的重要性和优越性.AG V 的重要特征是它的可移动性,对这种可移动性的控制是AG V 研制的核心问题.课题组研制的T UT -1型AG V 采用3种传感器(磁导航传感器、CC D 摄像机、超声波传感器)跟踪磁条来对AG V 进行引导和避障,经过这3种传感器的信息融合,测算出AG V 的位置和运动方向作为反馈与给定的运动状态进行比较,来调整AG V 下一步的运动[1]. 在天津市自然科学基金的资助下,课题组利用T UT -1这个平台开展基于Internet 的AG V 遥操作系 统的研究.为模拟AG V 的运动特性,利用Matlab 进行 AG V 的动力学建模.同时,在虚拟现实环境下,利用Matlab 模型仿真AG V 的路径跟踪,研究和探索AG V 最优的控制和配置方案. 1 实验和建模过程 如图1所示,T UT -1移动机器人在室内进行导航和避障的实验[2].AG V 通过磁导航传感器和CC D 摄像机跟踪磁条引导前进,当AG V 接近墙壁时,通过超声传感器引导.AG V 将实时采集到的磁条位置信息作为反馈,与给定的磁条标准位置信息进行比较来调整 Ξ收稿日期:2003212225 基金项目:天津市自然科学基金资助项目(023615011) 第一作者:葛为民(1968— ),男,讲师,博士研究生  第20卷第1期2004年3月天 津 理 工 学 院 学 报 JOURNA L OF TIAN JIN INSTITUTE OF TECHN OLOG Y V ol.20N o.1 Mar.2004

倒立摆系统地建模及Matlab仿真

倒立摆系统的建模及Matlab仿真 1.系统的物理模型 考虑如图(1)所示的倒立摆系统。图中,倒立摆安装在一个小车上。这里仅考虑倒立摆在图面内运动的二维问题。 图(1)倒立摆系统 假定倒立摆系统的参数如下。 摆杆的质量:m=0.1g 摆杆的长度:l=1m小车的质量: M=1kg重力加速度:g=9.8m/2s 摆杆的质量在摆杆的中心。 设计一个控制系统,使得当给定任意初始条件(由干扰引起)时,最大超调量 ≤10%,调节时间ts ≤4s ,通过小车的水平运动使倒立摆保持在垂直位置。

2.系统的数学模型 2.1建立倒置摆的运动方程并将其线性化。 为简化问题,在数学模型中首先假设:1)摆杆为刚体;2)忽略摆杆与支点之间的摩擦;3)忽略小车与接触面间的摩擦。 设小车瞬时位置为z,摆心瞬时位置为(θsin l z +),在u 作用下,小车及摆均产生加速远动,根据牛顿第二定律,在水平直线远动方向的惯性力应与u 平衡,于是有 u l z dt d m dt z d M =++)sin (22 22θ 即: u ml ml z m M =-++θθθθsin cos )(2&&&&& ① 绕摆轴转动的惯性力矩与重力矩平衡,因而有 θθθsin cos )sin (22mgl l l z dt d m =??? ????+ 即: θθθθθθθsin cos sin cos cos 22g l l z =-+&&&&& ② 以上两个方程都是非线性方程,为求得解析解,需作线性化处理。由于控制的目的是保持倒立摆直 立,在试驾合适的外力条件下,假定θ很小,接近于零时合理的,则1cos ,sin ≈≈θθθ,且可忽略θθ 2&项。于是有 u ml z m M =++θ&&&& )( ③ θθg l z =+&&&& ④ 联立求解可得 u Ml Ml m M u M M mg z 1)(1 -+=+- =θθθ&&&&

相关文档
最新文档