cuk.rar_Cuk Converter_biggestdx9_cuk_电力电子

clear all;clc %清除空间，清除控制窗口内容 V1=48*0.9+48*rand(1)*0.2;T=1/50000; d=0.556; %定义参数 ton=T*d; toff=T*(1-d); %Q1关断时间 y10=[0;0;0;0]; %IL1,IL2和UC1,UC2的初值为0 n=5;m=2000; %设置循环次数 t01=zeros(n,1); t02=zeros(m,1); %设定初始条件 tt=[];yy=[]; for i=1:n t01(i)=(i-1)*m*T; %计算第i个循环之前经过了多长时间 for j=1:m t02(j)=(j-1)*T; %计算第j个循环之前经过了多少时间 [t,y1]=ode45('cuk1',linspace(0,ton,6),y10); %计算导通时间内的电压、电流 tt=[tt;t+t01(i)+t02(j)];yy=[yy;y1]; %将计算结果记录在tt，yy两个矩阵中 y20=y1(end,:); %设置末值为下次计算的初值 [t,y2]=ode45('cuk2',linspace(0,toff,6),y20); %计算关断时间内的电压、电流 y10=y2(end,:); %设置末值为下一个循环的初值 tt=[tt;t+t01(i)+t02(j)+ton]; %将计算结果记录在tt，yy两个矩阵中 yy=[yy;y2]; %设置末值为下一个循环的初值 end end figure; %第一张图画电感电流iL1的波形 plot(tt,yy(:,1) ,'b') xlabel('t/s');ylabel('iL1/A'); figure; %第二张图画电感电流iL2的波形 plot(tt,yy(:,2) ,'b') xlabel('t/s');ylabel('iL2/A'); figure; %第三张图画输出电压uC1的波形 plot(tt,yy(:,3),'b') xlabel('t/s');ylabel('vC1/V'); figure %第四张图画输出电压uC2的波形 plot(tt,yy(:,4),'b') xlabel('t/s');ylabel('vC2/V');