Sepic电路建模、仿真设计

1、课程设计任务书1 题目sepic电路建模、仿真2 任务建立sepic电路的方程，编写算法程序，进行仿真，对仿真结果进行分析，合理选取电路中的各元件参数。3 要求课程设计说明书采用a4纸打印，装订成本；内容包括建立方程、编写程序、仿真结果分析、生成曲线、电路参数分析、选定。 v1=2040v v2=26vi0=0 1af=50khz目录1 sepic电路分析11.1 sepic电路简介11.2 原理分析11.3 电力运行状态分析22 sepic电路各元件的参数选择72.1 sepic电路参数初值72.2 电路各元件的参数确定73 控制策略的设定114 matlab编程仿真124.1根据状态方程编

2、写matlab子程序124.2 求解算法的基本思路134.3 matlab求解sepic电路主程序155 通过分析仿真结果合理选取电路参数l1，l2，c1，c2185.1参数l1的确定185.2参数l2的确定205.3参数c1的确定215.4参数c2的确定225.5 采用校核后的参数仿真246 采用matlab分析sepic斩波电路的性能246.1 计算电感l2的电流il2出现断续的次数246.2 纹波系数的计算256.3 电压调整率256.4 负载调整率266.5 电路的扰动分析277 参考文献301 sepic电路分析1.1 sepic电路简介 sepic斩波电路是开关电源六种基本dc/d

3、c变换拓扑之一，是一种允许输出电压大于、小于或者等于输入电压的dc/dc斩波电路。其输出电压由主控开关(三极管或mos管)的占空比控制。sepic变换器是一种四阶非线性系统, 因具有可升降压、同极性输出、输入电流脉动小、输出易于扩展等特点, 而广泛应用于升降压型直流变换电路和功率因数校正电路。这种电路最大的好处是输入输出同极性。尤其适合于电池供电的应用场合，允许电池电压高于或者小于所需要的输入电压。比如一块锂电池的电压为3v 4.2v，如果负载需要3.3v，那么sepic电路可以实现这种转换。另外一个好处是输入输出的隔离，通过主回路上的电容c1实现。同时具备完全关断功能，当开关管关闭时，输出电

4、压为0v。1.2 原理分析sepic斩波电路的原理图如图1所示。由可控开关q、储能电感l1、l2二极管d、储能电容c1、滤波电容c2、负载电阻r和控制电路等组成。图1、sepic斩波电路的原理图sepic斩波电路的基本工作原理是：当开关管q受控制电路的脉冲信号触发而导通时，v1l1q回路c1ql2回路同时导通，l1和l2储能。v处于断态时，v1l1c1d负载（c2和r）回路及l2d负载回路同时导通，此阶段v1和l1既向负载供电，同时也向c1充电，c1储存的能量在q处于通态时向l2转移。sepic斩波电路的输入输出关系由下式给出： （1） 图2、speic的开关波形（vq1q1漏源电压）1.3

5、电力运行状态分析对于理想情况下的电路分析，储能电感l1、l2足够大，其时间常数远大于开关的周期，流过储能电感的电流il可近似认为是线性的。电容c1、c2也足够大，能够维持两端电压恒定。此外，开关管q及二极管都具有理想的开关特性。分析电路图可以得到：1.3.1 q开通时电路运行分析mosfet开通时的等效电路如图2所示：图3、q开通时的等效电路图q开通时，输入电源v1对l1充电，储能电容c1对l2充电，电容c2维持着负载r的两端电压。此时有 （2）1.3.2 q关断时电路运行分析 mosfet关断时的等效电路如图2所示：图4、q关断时的等效电路图q关断后，充在电感l1上的电荷对电容c1放电，充在

6、电感l2上的电荷通过二极管d对负载放电。此时有 （3）1.3.3 输入直流电压v1和输出直流电压v2的关系稳态时，一个周期t内电感l两端电压ul对时间的积分为零，即 （4） 当q处于通态时，电感l1、l2两端的电压分别为、，当q处于关断时，电感l1、l2两端的电压分别为、。将数据代入式4得： （5）求解得： （6）稳态时，电容c的电流在一个周期t内的平均值应为零，也就是其对时间的积分为零，即 （7）当q处于通态时，流过电容c1、c2的电流分别为、，当q处于关断时，流过电容c1、c2的电流分别为、。将数据代入式7得： （8）求解得： （9）由式6知，所以可通过控制占空比的大小来控制输出电压v2的

7、大小。即当tontoff时，0.5， v2toff时，0.5， v2v1，电路属于升压式。1.3.4 电路的状态方程由图2、3可知，等效电路与开关q的状态有关，所以sepic斩波电路可分为q通态和q断态两个状态来分析。1）当q处于通态，系统的微分方程组如下所示 （10）2）当q处于断态，系统的微分方程组如下所示 （11）3）当q处于断态时，充在电感l1上的电荷对电容c1放电，充在电感l2上的电荷通过二极管d对负载放电，即此过程有可能会出现电感电流的断续。由于电感l1直接与电源相连，故一般来说l1的电流不会出现断续现象，下面主要讨论电感l2出现断续后，微分方程组的变化。电感l2断续后，即，此时微

8、分方程组如下所示： （12）设，将其代入式10、11，合并后如下：（13）其中，u=1表示q处于导通状态，u=0表示q处于关断状态。此外，u=0同时令，即表示q关断时电感l2的电流出现断续时的状态。2 sepic电路各元件的参数选择2.1 sepic电路参数初值 题目中给定，输入电压v1=2040v，输出电压v2=26v，负载电流i0=0 1a，开关管q的控制端的信号频率f=50khz，即周期t=2*10-5s。为了简化电路计算，更好地描述电路运行状态，现作如下假设：（1）电源电压为40v时为最差状态。（2）电路能达到满载电流1a。（3）忽略开关管的正向导通压降和二极管的正向压降。（4）忽略线

9、路电阻和电磁振荡所造成的能量损耗。2.2 电路各元件的参数确定2.2.1 负载电阻rl的确定负载电阻rl按式14确定 （14） 求得负载电阻rl=26.2.2.2 电感l1、l2的确定speic使用两个电感l1和l2，这两个电感可以绕在同一个磁芯上，因为在整个开关周期内加在它们上面的电压是一样的。使用耦合电感比起使用两个独立的电感可以节省pcb的空间并且可以降低成本。确定电感的一个好规则就是，在最小输入电压下，使得纹波电流峰峰值大约等于最大输入电流的40%。在数值相同的电感l1和l2中流动的纹波电流由下面公式算出:电感由15式求得 （15）f为开关频率，max是最小vin时的占空比。维持电感发

10、挥作用的电感峰值电流还没有饱和，可由下式计算如果l1和l2绕在同一个磁芯上，因为互感作用上式中的电感值就可用2l代替。电感值可这样计算 2.2.3 储能电容c1的确定储能电容c1的选择主要看rms电流（有效电流），可由下式得出sepic电容必须能够承受跟输出功率有关的有效电流。这种特性使sepic特别适用于流过电容的有效电流（跟电容技术有关）相对较小的小功率应用。sepic电容的电压额定值必须大于最大输入电压。c1的纹波电压的峰峰值可以这样计算 （16）取=0.4v得=28.261uf。满足需要的有效电流的电容在c1上一般不会产生太大的纹波电压，因此峰值电压通常接近输入电压。2.2.4 滤波电

11、容c2的确定在sepic中，当电源开关q1打开时，电感充电，输出电流由输出电容提供。因此输出电容会出现很大的纹波电流。选定的输出电容必须能够提供最大的有效电流。输出电容上的有效电流是图5、输出纹波电压esr、esl和大容量的输出电容直接控制输出纹波。如图4所示，假设一半纹波跟esr有关，另外一半跟容量有关，因此 (17)输出电容必须满足有效电流、esr和容量的需求。取纹波电压为2%的输出电压2.2.5输出二极体的选择选择能够承受峰值电流和反向电压的二极体。在speic中，二极体的峰值电流跟开关的峰值电流iq1peak相同。二极体必须承受的最小反向峰值电压是vrd1=vin+vout=66v （

12、18）跟升压转换器相似，二极体的平均电流跟输出电流相同。二极体的功耗等于输出电流乘以二极体的正向压降。为了提高效率建议使用肖特基二极体。2.2.6功率mosfet的选择最小阈值电压vth（min）、导通电阻rds(on)、栅漏电荷qgd和最大漏源电压vds（max）是选择mosfet的关键参数。逻辑电平或子逻辑电平阈值mosfet应该根据栅极电压使用。峰值开关电压等于vinvout。峰值开关电流由下式计算 (19)流过开关的rms电流由下式给出 (20)mosfet的散耗功率pq1大概是 (21) pq1，mosfet总的功耗包括导通损耗（上式第一项）和开关损耗（上式第二项）。igate为栅极

13、驱动电流。rds(on)值应该选最大工作结温时的值，一般在mosfet资料手册中给出。要确保导通损耗加上开关损耗不会超过封装的额定值或者超过散热设备的极限。2.2.7 编程计算所需参数在下面编程计算过程中，所需的电路参数如表1所示：表1、sepic斩波电路各元件参数值电路元件负载电阻（）电感l1（mh）电感l2（mh）电容c1(f)电容c2(f)频率f（khz）数值260.4350.43528.26143.480503 控制策略的设定由知，由于v2=26v，v1=2040v，即有=0.3940.565。v1初值为40v，即占空比的初值为0.394.由于输入不稳定，要想得到稳定的输出，需要对占空

14、比拟定相应的控制策略。本例采用的控制策略为：在每一次循环的结尾处对占空比d（i）作一定的调整，满足下式 （22）其中，k取0.00003，为每次计算后的输出电压v2，26为理想输出电压。当，即对占空比进行正向的调整，占空比增大，由式可知，输出的增大，即缩小与26的差距。反之，当，即对占空比进行负向的调整，占空比减小，输出的也变小，使输出逼近26。由上面分析可知，只要k取值足够小，循环次数n足够大，可以使输出电压稳定在26v附近。虽然k的取值越小，精度越高，但是k取值变小的同时也要求迭代次数n变大，否则迭代n次后还没有收敛结束。另一方面，迭代次数n的增大使得程序运行时间变长，本例选择k=0.00

15、003，n=2500，程序运行一次的时间约为5s。根据对输出电压平均值进行调制的方式不同，斩波电路可有三种控制方式：1）保持开关周期t不变，调节开关导通时间ton，称为脉冲宽度调制。2）保持开关导通时间ton不变，改变开关周期t，称为频率调制。3）ton和t都可调，使占空比改变，称为混合型。根据题目要求，开关频率f固定为50khz，故应采取第一种控制方式。4 matlab编程仿真在电力电子技术教学中, 通常利用power point，authorw are,，vb,，flash 等软件辅助教学, 这些软件虽然可以提供一些比较生动的动画和波形, 但其并非仿真软件, 关于电力电子器件的物理概念不突

16、出, 且应用灵活性较差。ewb、pspice、protel、matlab、sa-ber、plecs 等专用软件体现了很好的灵活性, 能够很好地满足教学要求。在这些软件中，matlab/simulink 软件构建仿真系统较灵活，被广泛应用于电力电子技术教学中；plecs 软件能为系统级电路仿真提供一个与simulink模型完全无缝的结合，在电力电子系统和电力驱动器的模拟上可以进行简化，其仿真速度快，开关转换理想化，稳定性好。因此,，matlab和plecs 软件相结合进行电力电子系统仿真是一种十分理想的选择。下面采用matlab对sepic进行仿真计算并分析仿真结果。4.1根据状态方程编写mat

17、lab子程序根据上述的电路分析可知，电路可能出现三种状态，每一种状态对应着不同的微分方程组。根据综合后的微分方程式13，可建立对应于电路的三种工作状态的子函数。综合后的微分方程所下所示 （23）1）建立子函数fun1.mq开通时，输入电源v1对l1充电，储能电容c1对l2充电，电容c2维持着负载r的两端电压。此时根据式23，令u=1，可得此时的微分方程，子函数语句如下：function dy2=fun1(t,y);global v1 r l1 c1 c2 l2;%定义全局变量dy2=v1/l1; y(3)/l2; -y(2)/c1;-y(4)/(c2*r); %状态一的微分方程2）建立子函数f

18、un2.mq关断后，充在电感l1上的电荷对电容c1放电，充在电感l2上的电荷通过二极管d对负载放电。此时根据式23，令u=0，可得此时的微分方程，子函数语句如下：function dy=fun2(t,y);global v1 r l1 c1 c2 l2; %定义全局变量dy=(v1-y(3)-y(4)/l1; -y(4)/l2; y(1)/c1; (y(1)+y(2)-y(4)/r)/c2; %状态二的微分方程3)建立子函数fun3.mq关断时电感l2对负荷放电，放电结束后电流出现断续，此时根据式23，令u=0、可得此时的微分方程，子函数语句如下：function dydt=fun3 (t,y

19、)global v1 r l1 c1 c2 l2; %定义全局变量dydt=(v1-y(3)-y(4)/l1; 0; y(1)/c1;(y(1)-y(4)/r)/c2; %状态三的微分方程4.2 求解算法的基本思路 基于matlab编程采用的思路如下所示算法流程图4.3 matlab求解sepic电路主程序 在4.1中建立了三个matlab的子程序fun1.m、fun2.m、fun3.m，分别对应了电路的三种状态。其中fun1.m为开关管q导通时的微分方程模型，fun2.m为开关管q关断时的微分方程模型，fun3.m为开关管q关断时电感l2放电结束后的微分方程模型。根据4.2的算法流程图编写的

20、matlab主程序如下：主程序：clear;global v1 r l1 c1 c2 l2 d f n m; %定义全局变量m=0; %记录电感l2的电流断续的次数r=26;l1=0.435e-3;l2=0.435e-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; %微分方程求解的初值for i=1:n;ton=d(i)*

21、t;toff=t-ton; %ton为导通时间，toff为关断时间 t,y1=ode45(fun1,linspace(0,ton,6),yy); yy=y1(end,:); %将导通时间的末值作为关断时间的初值 t,y2=ode45(fun2,linspace(ton,t,6),yy); yy=y2(end,:); %将关断时间的末值作为下一次导通时间的初值 if y2(end,2)0 %如果电流末值小于零，重新计算关断时间内的电压、电流 for a=1:length(y2) %找出il20的点 if y2(a,2)=1;d(i+1)=0.999;elseif d(i+1)=0;d(i+1)=

22、0.001;endenduo(i)%画图 figure(1); plot(linspace(0,1/100000*n,n),uo); %画出输出电压v2波形 title(输出电压v2的波形图) xlabel(t(s) ylabel(v2(v) grid on;figure(2); plot(linspace(0,1/100000*n,n),i1); %画出电感l1电流波形 title( i1的波形图) xlabel(t(s) ylabel(i1(a)grid on;figure(3); plot(linspace(0,1/100000*n,n),i2); %画出电感l2电流波形 title(

23、i2的波形图) xlabel(t(s) ylabel(i2(a) grid onfigure(4); plot(linspace(0,1/100000*n,n),vc1); %画出电容c1电压波形 title( vc1的波形图) xlabel(t(s) ylabel(vc1(v) grid on;程序运行后，结果如下所示：其中v2输出为25.9920v。 输出电压v2 电感l1电流il1电感l2电流il2 电容c1电压vc15 通过分析仿真结果合理选取电路参数l1，l2，c1，c2预设参数为l1=0.435mh，l2=0.435mh，c1=28.261uf，c2=43.480uf，现通过mat

24、lab仿真结果对比分析波形，对参数进行校核。5.1参数l1的确定1）减小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取值为

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

26、会使v2的收敛过程变慢，但c1过小会使v2脉动变大。综合两方面因素，本例c1不改变，即取初值c1=28.261uf。5.4参数c2的确定1）减小c2，取c2=0.5 uf，此时超调量明显变小，但是波形发生畸变。2）取c2=10 uf，此时波形超调量相对c1=43.480uf时有所下降，波形同时也得到较好的收敛，但是收敛速度变慢。3）继续增大c2，取c2=80 uf，此时波形发生畸变，v2不能很好收敛。由上面仿真知，c2过大及过小都会使得波形发生畸变，减小c2可以减小v2超调量，但是收敛速度会变慢。本例折中取c2=30uf。综上，通过matlab仿真校核后，各电路元件参数如下所示：元件电感l1（

27、mh）电感l2（mh）电容c1(f)电容c2(f)数值10228.261305.5 采用校核后的参数仿真采用该参数进行仿真，v2的波形如下图所示 输出电压v2 电感l1电流 电感l2电流 电容c1电压6 采用matlab分析sepic斩波电路的性能6.1 计算电感l2的电流il2出现断续的次数 在程序的条件判断语句if y2(end,2)0内，加入m=m+1语句来记录il2出现断续的次数，最后输出v2 =25.9983，m =45。即整个收敛过程，il2出现了45次断续。换言之，电路元件按本例所选的参数下，不能保证il2连续。6.2 纹波系数的计算纹波系数为纹波电压的幅值与直流电压幅度之比，是

28、一种动态特性指标。用公式表示如下：设t=0.02s时，v2已收敛结束，采用语句axis(0.02,0.021,25.9,26.1)截取v2波形图如下所示， 借助语句x,y=ginput(1)读出，输出电压峰值=26.0401，输出电压谷值=25.9463。则纹波系数为可知输出电压波形摆动较小，电压稳定性高6.3 电压调整率电压调整率为在一正常的固定负载下，由输入电压变化所造成其输出电压偏差率的百分比，是一种静态特性。计算公式为 （24） 本例负载确定为r=26，输出电流为i0=1a。分别将输入电压设为v1max=40v、v1min=20v进行仿真。1）输入v1max=40v，此时输出v2=u0（max）=25.9985v。2）输入v1min=20v，此时输出v2=u0（min）=25.9938v。 将数据代入式30后，求得电压调整率为0.018%，换言之，在输入变化较大时，输出变化十分小。系统引入了反馈环节，每迭代一次都对占空比进行调整，控制策略在这里得到了很好体现。6.4 负载调整率负载调整率为在正常的固定输入电压下，由负载电流变化所造成其输出电压偏差率的百分比，也是