Sepic电路课程设计说明书资料

1、?怎孝自羯N大亨SouthChinaUniversityofTechnology课程设计说明书课程名称：电力电子课程设计设计题目：Sepic电路的建模与仿真专业：电气工程及其自动化班级：2009级电气（4）班学号：200930213291姓名：姓培正指导教师：郭红霞华南理工大学电力学院2013年1月课程设计任务书1. .题目Sepic电路建模、仿真2. 任务建立Sepic电路的方程，编写算法程序，进行仿真，对仿真结果进行分析，合理选取电路中的各元件参数。3. 要求课程设计说明书采用A4纸打印，装订成本；内容包括建立方程、编写程序、仿真结果分析、生成曲线、电路参数分析、选定。V1=2040VV2

2、=26VI0=01AF=50kHz指导教师评语：指导教师：2013年月目录1 Sepic电路分析11.1 Sepic电路简介11.2 原理分析11.3 电力运行状态分析22 Sepic电路各元件的参数选择72.1 Sepic电路参数初值72.2 电路各元件的参数确定73 控制策略的设定114 Matlab编程仿真124.1 根据状态方程编写Matlab子程序124.2 求解算法的基本思路134.3 Matlab求解Sepic电路主程序155 通过分析仿真结果合理选取电路参数L1，L2，C1，C2185.1 参数L1的确定185.2 参数L2的确定205.3 参数C1的确定215.4 参数C2的

3、确定225.5 采用校核后的参数仿真246采用Matlab分析Sepic斩波电路的性能246.1 计算电感L2的电流IL2出现断续的次数246.2 纹波系数的计算256.3 电压调整率256.4 负载调整率266.5 电路的扰动分析277参考文献301Sepic电路分析1.1 Sepic电路简介Sepic斩波电路是开关电源六种基本DC/DC变换拓扑之一，是一种允许输出电压大于、小于或者等于输入电压的DC/DC斩波电路。其输出电压由主控开关（三极管或MOS管）的占空比控制。SEPIC变换器是一种四阶非线性系统，因具有可升降压、同极性输出、输入电流脉动小、输出易于扩展等特点，而广泛应用于升降压型直

4、流变换电路和功率因数校正电路。这种电路最大的好处是输入输出同极性。尤其适合于电池供电的应用场合，允许电池电压高于或者小于所需要的输入电压。比如一块锂电池的电压为3V4.2V,如果负载需要3.3V,那么Sepic电路可以实现这种转换。另外一个好处是输入输出的隔离，通过主回路上的电容C1实现。同时具备完全关断功能，当开关管关闭时，输出电压为0V。1.2 原理分析Sepic斩波电路的原理图如图1所示。由可控开关Q、储能电感L1、L2二极管D、储能电容C1、滤波电容C2、负载电阻R和控制电路等组成。图1、Sepic斩波电路的原理图Sepic斩波电路的基本工作原理是：当开关管Q受控制电路的脉冲信号触发而

5、导通时，V1L1Q回路C1QL2回路同时导通，L1和L2储能。V处于断态时，V1L1C1D负载（C2和R）回路及L2D负载回路同时导通，此阶段V1和L1既向负载供电，同时也向C1充电，C1储存的能量在Q处于通态时向L2转移。Sepic斩波电路的输入输出关系由下式给出:(1)OUT图2、SPEIC的开关波形（Vqi : Qi漏源电压）1.3 电力运行状态分析对于理想情况下的电路分析，储能电感Li、L2足够大，其时间常数远大于开关的周期，流过储能电感的电流Il可近似认为是线性的。电容Ci、C2也足够大，能够维持两端电压恒定。止匕外，开关管Q及二极管都具有理想的开关特性。分析电路图可以得到：1.3.

6、1 Q开通时电路运行分析MOSFET开通时的等效电路如图2所示：图3、Q开通时的等效电路图Q开通时，输入电源Vi对Li充电，储能电容Ci对L2充电，电容C2维持着负载R的两端电压。此时有VliViVL2VC1iC1I2(2)iC21.3.2 Q关断时电路运行分析MOSFET关断时的等效电路如图2所示:图4、Q关断时的等效电路图Q关断后，充在电感Li上的电荷对电容Ci放电，充在电感L2上的电荷通过二极管D对负载放电。此时有ViVL1VciV20VL2V2iC1iC2IlIl(3)稳态时，一个周期T内电感L两端电压Ul对时间的积分为零，即1.3.3 输入直流电压Vi和输出直流电压V2的关系T(4)

7、0ULdt当Q处于通态时，电感Li、L2两端的电压分别为Vi、Vci,当Q处于关断时,电感Li、L2两端的电压分别为ViVciV2、V20将数据代入式4得:(5)tonMtoffViV2Vci0tonVci骁(V2)0求解得:Vci Vi(6)稳态时，电容C的电流在一个周期T内的平均值应为零，也就是其对时间的积分为零，即T0icdt0当Q处于通态时，流过电容Ci、C2的电流分别为I2、必,当Q处于关上。将数据代入式7得:R断时，流过电容Ci、C2的电流分别为Ii、IiI2tonI2toffI10VVton(片)toff(IiI2力0RR求解得:IiV2(9)由式6知，V2V1,所以可通过控制占

8、空比的大小来控制输出电压V2的1大小。即当ton<toff时，<0.5,V2<V1,电路属于降压式；当ton=toff时，=0.5,V2=V1；当ton>toff时，>0.5,V2>V1,电路属于升压式。1.3.4 电路的状态方程由图2、3可知，等效电路与开关Q的状态有关，所以Sepic斩波电路可分为Q通态和Q断态两个状态来分析。1）当Q处于通态，系统的微分方程组如下所示也1VdtL1diL2VC1dtL2(10)dVC1iL2dtC1dV2Vc2C2C2dtRC22)当Q处于断态，系统的微分方程组如下所示diLi"dTdiL2dtdVcidtdV

9、c2dtVi Vc2 VciLiVC2iL1CiL1 iL2C2VC2RC2(ii)3)当Q处于断态时，充在电感Li上的电荷对电容C1放电，充在电感L2上的电荷通过二极管D对负载放电，即此过程有可能会出现电感电流的断续。由于电感L1直接与电源相连，故一般来说L1的电流不会出现断续现象，下面主要讨论电感L2出现断续后，微分方程组的变化。电感L2断续后，即日0,此时微分方程组如下所示：diLiViVc2VcidtLiL2(dVCiiLidtCid%iLiVc2dtC2RC2设Y(i)k,Y(2)L2,Y(3)Vci,丫”,将其代入式i0、ii,合并后如下：Y(1)(Y(3)Y(4)/Lu(Y(3)

10、Y(4)/LV"L?Y(2)Y(4)/L2u(Y(3)Y(4)/L2?_(13)Y(3)Y(1)/Ciu(Y(1)Y(2)/Ci卜)?Y(4)(Y(1)Y(2)/C2Y(4)/(RC2)u(Y(1)Y(2)/C2其中，u=1表示Q处于导通状态，u=0表示Q处于关断状态。此外，u=0同时令Y(2)0,即表示Q关断时电感L2的电流出现断续时的状态。2Sepic电路各元件的参数选择2.1 Sepic电路参数初值题目中给定，输入电压V1=2040V,输出电压V2=26V,负载电流Io=01A,开关管Q的控制端的信号频率F=50kHZ,即周期T=2*10-5So为了简化电路计算，更好地描述电路

11、运行状态，现作如下假设：(1)电源电压为40V时为最差状态。(2)电路能达到满载电流1A。(3)忽略开关管的正向导通压降和二极管的正向压降。(4)忽略线路电阻和电磁振荡所造成的能量损耗。2.2 电路各元件的参数确定2.2.1 负载电阻Rl的确定负载电阻Rl按式14确定Rl运(14)Io求得负载电阻Rl=26Q.2.2.2 电感Li、L2的确定SPEIC使用两个电感L1和L2,这两个电感可以绕在同一个磁芯上，因为在整个开关周期内加在它们上面的电压是一样的。使用耦合电感比起使用两个独立的电感可以节省PCB的空间并且可以降低成本。确定电感的一个好规则就是，在最小输入电压下，使得纹波电流峰峰值大约等于

12、最大输入电流的40%在数值相同的电感L1和L2中流动的纹波电流由下面公式算出:IL I in 40%I outVou 40%Vin (min)26140%A 0.52A20电感由15式求得L2 LVin (min)IL fmax200.52 5026 mH 0.435mH26 20(15)f为开关频率,omax 是最小Vin时的占空比维持电感发挥作用的电感峰值电流还没有饱和，可由下式计算:I L1( peak) I outVoutVin(min)40%、， (1+)=1.56A2I L2( peak) I out(1 + 鲁=1.2A如果L1和L2绕在同一个磁芯上，因为互感作用上式中的电感值就

13、可用2L代替。电感值可这样计算:L1L2LVin (min)22 IL f max202 0.52 50 2626 mH 0.218mHi 202.2.3 储能电容Ci的确定储能电容Ci的选择主要看RMS电流（有效电流）,可由下式得出：I C1 I outVoutVin (min)SEPIC电容必须能够承受跟输出功率有关的有效电流。这种特性使SEPIC特别适用于流过电容的有效电流（跟电容技术有关）相对较小的小功率应用。SEPIC电容的电压额定值必须大于最大输入电压。Ci的纹波电压的峰峰值可以这样计算：VC1outmaxC1f(16)取Vci=0.4V得Ci=28.261uF。满足需要的有效电流

14、的电容在Ci上一般不会产生太大的纹波电压，因此峰值电压通常接近输入电压。2.2.4 滤波电容C2的确定在SEPIC中，当电源开关Q1打开时，电感充电，输出电流由输出电容提供。因此输出电容会出现很大的纹波电流。选定的输出电容必须能够提供最大的有效电流。输出电容上的有效电流是：.VOUTA 1.140A1 Cout(RMS)1 OUTV VIN (min)图5、输出纹波电压ESRESL和大容量的输出电容直接控制输出纹波。如图4所示，假设一半纹波跟ESR有关，另外一半跟容量有关，因此ESRV-0.5I L1(peak) IL 2( peak)0.2 26 0.52.760.942C2I OUT输出电

15、容必须满足有效电流、Vripple0.5 f(17)ESR和容量的需求。取纹波电压为2%的输出电压' outVripple 0.5 fmax126 l3 F0.02 26 0.5 50 103 4643.480 F2.2.5 输出二极体的选择选择能够承受峰值电流和反向电压的二极体。在SPEIC中，二极体的峰值电流跟开关的峰值电流iQIpeak相同。二极体必须承受的最小反向峰值电压是：VRD1=Vin+Vout=66V（18）跟升压转换器相似，二极体的平均电流跟输出电流相同。二极体的功耗等于输出电流乘以二极体的正向压降。为了提高效率建议使用肖特基二极体。2.2.6 功率MOSFET的选择

16、最小阈值电压Vth（min）、导通电阻RdS（ON）、栅漏电荷QGD和最大漏源电压VDS（max）是选择mosfet的关键参数。逻辑电平或子逻辑电平阈值mosfet应该根据栅极电压使用。峰值开关电压等于Vin+Vouto峰值开关电流由下式计算：I Q1( peak)IL1(peak) IL2(peak) 2.76A(19)流过开关的Rms电流由下式给出：I Q1(RMS)1 out(VOUTVIN (min) ) VOUTVIN (min)7.733A(20)MOSFET的散耗功率Pqi大概是:2P21I Q1(RMS)RdS(ON)max(VIN (min) Vout ) I Q1(peak

17、)QgdfIgate(21)Pqi,MOSFET总的功耗包括导通损耗（上式第一项）和开关损耗（上式第二项）。Igate为栅极驱动电流。Rds（on）值应该选最大工作结温时的值，一般在mosfet资料手册中给出。要确保导通损耗加上开关损耗不会超过封装的额定值或者超过散热设备的极限。2.2.7 编程计算所需参数在下面编程计算过程中，所需的电路参数如表1所示：表1、Sepic斩波电路各元件参数值电感Li电感L2电容Ci电容C2频率F元件(Q)(mH)(mH)(N)(N)(kHZ)260.4350.43528.26143.480503控制策略的设定由V2“WiVi知，-2,由于V2=26V,V1=20

18、40V,即有toff1ViV2=0.3940.565。Vi初值为40V,即占空比的初值为0.394.由于输入不稳定，要想得到稳定的输出，需要对占空比拟定相应的控制策略。本例采用的控制策略为：在每一次循环的结尾处对占空比d（i）作一定的调整，满足下式di1dik26uoi（22）其中，k取0.00003,uoi为每次计算后的输出电压V2,26为理想输出电压。当uoi26时，k26uoi0,即对占空比进行正向的调整，占空比增大,由式uoid i-Vi可知，输出的1uo i增大，即缩小与26的差距。反之,当uoi26时，k26uoi0,即对占空比进行负向的调整，占空比减小，输出的uoi也变小，使输出

19、逼近26。由上面分析可知，只要k取值足够小，循环次数n足够大，可以使输出电压uoi稳定在26V附近。虽然k的取值越小，精度越高，但是k取值变小的同时也要求迭代次数n变大，否则迭代n次后还没有收敛结束。另一方面，迭代次数n的增大使得程序运行时间变长，本例选择k=0.00003,n=2500,程序运行一次的时间约为5s。根据对输出电压平均值进行调制的方式不同，斩波电路可有三种控制方式：1）保持开关周期T不变，调节开关导通时间ton,称为脉冲宽度调制。2）保持开关导通时间ton不变，改变开关周期T,称为频率调制。3）ton和T都可调，使占空比改变，称为混合型。根据题目要求，开关频率F固定为50kHZ

20、,故应采取第一种控制方式。4Matlab编程仿真在电力电子技术教学中,通常利用PowerPoint,Authorware,VB,Flash等软件辅助教学，这些软件虽然可以提供一些比较生动的动画和波形，但其并非仿真软件，关于电力电子器件的物理概念不突出，且应用灵活性较差。EWB、PSpiceProtel、Matlab、SABER、PLECS等专用软件体现了很好的灵活性，能够很好地满足教学要求。在这些软件中，Matlab/Simulink软件构建仿真系统较灵活，被广泛应用于电力电子技术教学中；PLECS软件能为系统级电路仿真提供一个与Simulink模型完全无缝的结合，在电力电子系统和电力驱动器的

21、模拟上可以进行简化，其仿真速度快，开关转换理想化，稳定性好。因此，Matlab和PLECS软件相结合进行电力电子系统仿真是一种十分理想的选择。下面采用Matlab对Sepic进行仿真计算并分析仿真结果。4.1根据状态方程编写Matlab子程序根据上述的电路分析可知，电路可能出现三种状态，每一种状态对应着不同的微分方程组。根据综合后的微分方程式13,可建立对应于电路的三种工作状态的子函数。综合后的微分方程所下所示Y(1)(Y(3)Y(4)/Liu(Y(3)Y(4)/LiV"Li?Y(2)Y(4)/L2u(Y(3)丫(4)人?Y(3)Y(1)/Ciu(Y(1)Y(2)/Ci?Y(4)(Y

22、(1)Y(2)/C2Y(4)/(RC2)u(Y(1)Y(2)/C21）建立子函数fun1.mQ开通时，输入电源Vi对Li充电，储能电容Ci对L2充电，电容C2维持着负载R的两端电压。此时根据式23,令u=1,可得此时的微分方程，子函数语句如下：functiondy2=fun1(t,y);globalV1RL1C1C2L2;%定义全局变量dy2=V1/L1;y(3)/L2;-y(2)/ci；-y(4)/(C2*R);%状态一的微分方程2)建立子函数fun2.mQ关断后，充在电感Li上的电荷对电容Ci放电，充在电感L2上的电荷通过二极管D对负载放电。此时根据式23,令u=0,可得此时的微分方程，子

23、函数语句如下：functiondy=fun2(t,y);globalV1RLiCiC2L2;%定义全局变量dy=(Vl-y(3)-y(4)/Ll；-y(4)/L2;y(i)/ci;(y(1)+y(2)-y(4)/R)/C2;%状态二的微分方程3)建立子函数fun3.mQ关断时电感L2对负荷放电，放电结束后电流出现断续，此时根据式23,令u=0、Y(2)0可得此时的微分方程，子函数语句如下：functiondydt=fun3(t,y)globalViRLiCiC2L2;%定义全局变量dydt=(Vi-y(3)-y(4)/Li;0；y(i)/ci;(y(i)-y(4)/R)/C2；%状态三的微分方

24、程4.2 求解算法的基本思路基于Matlab编程米用的思路如下所小算法流程图4.3 Matlab求解Sepic电路主程序在4.1中建立了三个Matlab的子程序funl.m、fun2.m、fun3.m,分别对应了电路的三种状态。其中funl.m为开关管Q导通时的微分方程模型，fun2.m为开关管Q关断时的微分方程模型，fun3.m为开关管Q关断时电感L2放电结束后的微分方程模型。根据4.2的算法流程图编写的Matlab主程序如下：主程序：clear;globalV1RL1C1C2L2dfnm;%定义全局变量m=0;%记录电感L2的电流断续的次数R=26;L1=0.435e-3;L2=0.435

25、e-3;C1=28.261e-6;C2=43.480e-6;f=50000;T=1/f;n=2500;%迭代的次数V1=20+20*rand(1);%输入电压为2040Vd=0.394*ones(n,1);%定义占空比初值为0.394的一组向量uo=zeros(n,1);%定义输出初值1的一组向量yy=0,0,0,0;%微分方程求解的初值fori=1:n;ton=d(i)*T;toff=T-ton;%ton为导通时间，toff为关断时间t,y1=ode45('funT,linspace(0,ton,6),yy);yy=y1(end,:);%将导通时间的末值作为关断时间的初值t,y2=o

26、de45('fun2',linspace(ton,T,6),yy);yy=y2(end,:);%将关断时间的末值作为下一次导通时间的初值ify2(end,2)<0%如果电流末值小于零，重新计算关断时间内的电压、电流fora=1:length(y2)%找出iL2=0的点ify2(a,2)<0b=a;break,endendm=m+1;%记录电流断续次数nnmm=size(y2);toff1=toff*(b-1.5)/(nn-1);%根据iL2=0的点计算toffltoff2=toff-toff1;%根据心=0的点toff2t1,y21=ode45('fun2&

27、#39;,linspace(0,toff1,6),yy);%计算toffl内的电压、电流y21(end,1)=0;%设置末值为下次计算的初值t2,y22=ode45('fun3',linspace(0,toff2,6),y21(end,:);%计算toff2内的电压、电流t=t1;t2+toff1;y2=y21;y22;end%输出为uo%输出电感L1电流%输出电感L2电流%输出电容C1电压%每次迭代调整占空比uo(i)=(y1(end,4)+y2(end,4)/2;I1(i)=y2(end,1);I2(i)=y2(end,2);VC1(i)=y2(end,3);k=0.000

28、03;d(i+1)=d(i)+k*(26-uo(i);ifd(i+1)>=1;d(i+1)=0.999;elseifd(i+1)<=0;d(i+1)=0.001;endenduo(i)%画图figure。);%画出输出电压V2波形plot(linspace(0,1/100000*n,n),uo);title('输出电压V2的波形图')xlabel('t(s)')ylabel('V2(V)')gridon;figure(2);plot(linspace(0,1/100000*n,n),I1);title('I1的波形图'

29、)xlabel('t(s)')ylabel('I1(A)')gridon;figure(3);plot(linspace(0,1/100000*n,n),I2);title('I2的波形图')xlabel('t(s)')ylabel('I2(A)')gridonfigure(4);plot(linspace(0,1/100000*n,n),VC1);%画出电感L1电流波形%画出电感L2电流波形%画出电容C1电压波形title('VC1的波形图')xlabel('t(s)')ylabe

30、l('VC1(V)')gridon;程序运行后，结果如下所示：其中V2输出为25.9920V。50100.0050.010.0150.020.025t(s)输出电压V2的波形图-1000.0050.010.0150.020.025t(s)ii的波形图输出电压V2电感L1电流IL1)t(s)t(s)电感L2电流IL2t(s)电容C1电压VC15通过分析仿真结果合理选取电路参数Li,L2,Ci,C2预设参数为L1=0.435mH,L2=0.435mH,C1=28.261uF,C2=43.480uF,现通过Matlab仿真结果对比分析波形，对参数进行校核。5.1 参数L1的确定1）减

31、小L1,L1取0.2mH,输出电压V2波形如下，止匕时V2不能很好地收敛,波形出现畸变。2）增大L1,L1取2mH,输出电压V2波形如下，此时V2得到很好地收敛，收敛速度较快。3）增大L1,L1取20mH,输出电压V2波形如下，此时V2得到很好地收敛，振荡次数更少，超调量也得到一定的下降，收敛速度更快。由上面的分析可知，一方面L1太小会使V2不能收敛；另一方面，增大L1可以加快V2的收敛速度。此外，当L1取值足够大后，再增大L1,此时收敛速度变化不明显。综上，可将L1取值为10mH。5.2 参数L2的确定1）减小L2,L2取值为0.1mH,此时输出电压收敛速度变快，但是超调量较大。2）增大L2

32、,当L2取为4mH,此时输出电压收敛速度变慢，但超调量明显减小3）继续增大L2,当L2取值为10mH时，超调量虽然较小，但输出电压V2不能稳定在一个较小的范围。综上，增大L2可以减小超调量，但是收敛速度变慢，且L2太大时会使得输出电压V2不能稳定在一个较小的范围。本例取L2=2mH05.3 参数Ci的确定1）减小Ci,C1取5uF,输出电压V2波形如下，此时V2收敛速度加快，但是V2最终的脉动较大。2)增大C1,C1取50uF,输出电压V2波形如下，此时V2收敛速度缓慢，只要仿真时间足够，V2能得到较好的收敛。输出电压V2的波形图500 00.0050.010.0150.020.02540O2

33、Ot(s)综上，增大C1会使V2的收敛过程变慢，但C1过小会使V2脉动变大。综合两方面因素，本例C1不改变，即取初值C1=28.261uF。5.4 参数C2的确定403510输出电压V2的波形图51)减小C2,取C2=0.5uF,此时超调量明显变小，但是波形发生畸变。0.0050.010.0150.020.025t(s)VC2,取C2=80 uF,此时波形发生畸变，V2不能很好收敛。V2）取C2=10uF,此时波形超调量相对C1=43.480uF时有所下降，波形同时也得到较好的收敛，但是收敛速度变慢。3）继续增大C2可以减小由上面仿真知，C2过大及过小都会使得波形发生畸变，减小V2超调量，但是

34、收敛速度会变慢。本例折中取C2=30uFo综上，通过Matlab仿真校核后，各电路元件参数如下所示：元件电感Li(mH)电感L2(mH)电容Ci(nD电容C2(小»10228.261305.5 采用校核后的参数仿真采用该参数进行仿真，V2的波形如下图所示输出电压V2电感L1电流电感L2电流电容C1电压6采用Matlab分析Sepic斩波电路的性能6.1 计算电感L2的电流IL2出现断续的次数在程序的条件判断语句ify2(end,2)<0内，加入m=m+1语句来记录Il2出现断续的次数，最后输出V2=25.9983,m=45。即整个收敛过程，Il2出现了45次断续。换言之，电路元

35、件按本例所选的参数下，不能保证IL2连续。6.2 纹波系数的计算纹波系数为纹波电压的幅值与直流电压幅度之比，是一种动态特性指标。用公式表示如下:件址少米心输出电压峰值-输出电压谷值纹波系数=+、-100%2直流电压设t=0.02s时，V2已收敛结束，采用语句axis(0.02,0.021,25.9,26.1)l取V2波形图如下所示,输出电压V2的波形图26.126.052625.9525.9 0.020.0210.022t(s)0.0230.024借助语句x,y=ginput(1)读出，输出电压峰值二26.0401,输出电压谷值=25.9463。则纹波系数为人冲三跖26.0401-25.9463纹波系数=100%=0.18%<1%226可知输出电压波形摆动较小，电压稳定性高6.3 电压调整率电压调整率为在一正常的固定负载下，由输入电压变化所造成其输出电压偏差率的百分比，是一种静态特性。计算公式为(24)电压调整率=Uo(max)Uo(min)100%Uo(normal)本例负载确定为R=26Q,输出电流为Io=1A。分别将输入电压设为Vimax=40V、Vlmin=20V进行仿真。1）输入Vimax=40V,此时输出V2=Uo（max=25.9985V。2）输入Vimin=20V,此时输出V2=Uo（min）=25