逆变器的电感电路双环PI控制，最后matlab求出相位裕度103度,

只看该作者 · 2008-1-21 16:04

%%%%%%%%%%% 求电流环的Kip %%%%%%%%%%%%%%%%%%%%%%%%
PI = 3.1415926;
Vout = 45;
Iout = 0.7;
Alpha = (3.225/2) / Iout; % 输出 = Alpha * 输入
Beta = (3.225/2) / Vout;
Kpwm = 48;
R = 70;
C = 10.8 * 10^(-6);
L = 1.2 * 10^(-3);
fn = 1 / (2 * PI * ((L * C) ^ 0.5));
%fc = fn / 16; %穿越频率
fc = (1/10) *20000;
s = 2 * PI * fc * i;
%求得电流环的Kip
Gio = Alpha * Kpwm * (R*C*s + 1) /(R*C*L*s^2 + L*s + R);
Gio = abs(Gio);
Kip = 1/Gio
%电流环开环单位阶越响应
SysGio = tf([Gio*Alpha*Kip*Kpwm* R*C Alpha*Kip*Kpwm], [R*L*C L R]);
step(SysGio)
%margin(SysGio);
% 求电流环闭环传递函数
SysGI = tf([Kip*Kpwm*R*C Kip*Kpwm], [R*C*L L R]);
SysHI = tf([Alpha], [1]);
SysGic = feedback(SysGI, SysHI, -1)%电流闭环函数
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Kip =
0.0708
电流闭环函数为:
Transfer function:
0.00257 s + 3.399
-----------------------------------
9.072e-007 s^2 + 0.007119 s + 77.83

只看该作者 · 2008-1-21 21:12

%电压环Kvi, Kvp参数求解

Fvcut = 250; %电压环的穿越频率

zeta = 2*PI* Fvcut*i;

Fvzero = 300; %电压环PI的零点

GvNum = (zeta + 2*PI*Fvzero) *R * Beta * (0.00257 *zeta + 3.399);

GvDen = zeta * (R*C*zeta +1) * (9.072e-007 *zeta^2 + 0.007119 *zeta + 77.83);

Gv = abs(GvNum/GvDen);

Kvp = 1/ Gv

Kvi = 2 *PI* Fvzero *Kvp

%系统的开环伯德图

Gvpi = tf([Kvp Kvi], [1 0]);

Gr = tf([Beta *R], [R*C 1]);

%Gbeta = tf([Beta], [1]);

Gvo = series(Gvpi, series(SysGic, Gr))

margin(Gvo)

逆变器的电感电路双环PI控制，最后matlab求出相位裕度103度,

相关下载

相关帖子

求电流环的Kip

电压环Kvi, Kvp参数求解

无冕之王奖章

沉静之湖泊

七世轮回

技术高手奖章