文档介绍：课程设计说明书
课程名称：电力电子课程设计
设计题目：Cuk电路的建模与仿真
专 业： 电气工程及其自动化
班 级： 06电气（4）班
学 号： 2006301方程的M函数为：
function y=cuk2(t,x)
%参数设定
L1=15e-5;
L2=15e-5;
C1=47e-6;
C2=47e-6;
v1=48*+rand(1)*48*;
r=60;
%微分方程M函数主体
y=[(v1-x(3))/L1;
-x(4)/L2;
x(1)/C1;
(x(2)-x(4)/r)/C2];
%输入x(1)为iL1,输入x(2)为iL2;输入x(3)为uC2;输入x(4)为uC1
（3）求解微分方程的主程序如下：
% 参数设定
f=50000; %频率为50kHz
T=1/f; %周期为T
d=; %占空比为d
ton=T*d; %Q1导通时间
toff=T*(1-d); %Q1关断时间
y10=[0,0,0,0]; %iL1,iL2和uC1,uC2的初值为0
n=4000; %循环次数为4000次
% 求解微分方程
clear; %清空内存
f=50000; %频率
T=1/f; %周期
n=3;m=2000 %设置循环次数
t01=zeros(m,1);
t02=zeros(n,1);
y10=[0,0]; %设定初始条件
d=; %设置初始占空比
u1=36+24*rand(1) %输入电压，36-60的随机数
tt=[],yy=[]
for j=1:n
ton=T*d %计算导通时间
toff=(1-d)*T %计算关断时间
t02(j)=(j-1)*m*T %计算第j个循环之前经过了多长时间
for i=1:m
t01(i)=(i-1)*T; %计算第i个循环之前经过了多少时间
[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(1); %第一张图画电感电流iL1的波形
hold on;
plot(tt,yy(:,1) ,'b')
title('电感电流iL1波形 ');
xlabel('时间t(s)');
ylabel('电流iL1(A)');
figure(2); %第二张图画输出电压uC2的波形
hold on;
plot(tt,yy(:,4),'b')
title('输出电压uC1 波形');
xlabel('时间t(s)');
ylabel('电压uC1(V)');
figure(3); %第三张图画电感电流iL2的波形
hold on;
plot(tt,yy(:,2) ,'b')
title('电感电流iL2波形 ');
xlabel('时间t(s)');
ylabel('电流iL2(A)');
figure(4); %第四张图画输出电压uC1的波形
hold on;
plot(tt,yy(:,3),'b')
title('输出电压uC2 波形');
xlabel('时间t(s)');
ylabel('电压uC2(V)');
(1)电流的波形
（2）电流的波形
（3