%-----------------------------------------------------------------------
% UC3825 в режиме управления по току (slua053)

f=logspace(0,6,1E3);	% диапазон исследуемых частот от 1 Гц до 1 МГц 	 

w=2*pi*f; 		% w — круговая частота
s=j*w;

Vin=48;		% входное напряжение
Vo=5;		% выходное напряжение
fs=1.5e6;	% частота коммутации
Io=10;		% ток нагрузки
N=0.2;		% коэфициент трансформации
Lo=0.74e-6;	% индуктивность выходного дросселя 
RLo=0.01;	% ESR выходного дросселя
Co=3e-6;	% емкость выходного конденсатора
RCo=10e-3;	% ESR выходного конденсатора
Vd=0.1;		% падение напряжения на выходном выпрямителе

Ro=Vo/Io;

% импеданс выходного конденсатора
ZCo=RCo+1./(s*Co); 

% импеданс дросселя
ZLo=RLo+s*Lo;

% комплексная нагрузка
ZRo=Ro.*ZCo./(Ro+ZCo); 		

% передаточная характеристика LC-фильтра
KLC=ZRo./(ZRo+ZLo); 

% входное сопротивление LC-фильтра
ZLCin= ZRo+ZLo;

% передаточная функция силовой части
Gvd=Vin*N*KLC;

%-----------------------------------------------------------------------
% Передаточные функции в режиме управления по пиковому току
Rs=1.5/4;			% сопротивление датчика тока
Rf=1e3;				% резистор фильтра
Ra=6.8e3;			% резистор компенсации наклона
Va=1.8;				% размах пилы генератора
Von=Vin*N-Vd-Vo;		% напряжение на дросселе в интервале импульса
Voff=Vo+Vd;				% напряжение на дросселе в интервале паузы
Sn=Von/Lo*N*Rs*Ra/(Rf+Ra);		% прямой «наклон»
Sf=Voff/Lo*N*Rs*Ra/(Rf+Ra);		% обратный «наклон»
Se=Va*fs*Rf/(Rf+Ra);

% Коэффициент передачи модулятора
Fm=fs/(Sn+Se);

% Коэффициент передачи датчика тока
Hi=Rs*Ra/(Rf+Ra);

wn=pi*fs;
Qz=-2/pi;

He=1+s/(wn*Qz)+s.^2/(wn^2);

% передаточная характеристика ток-коэффициент заполнения
Gid=Vin*N./ZLCin;

% Передаточная характеристика петли управления по току
Ti=Fm*Hi*Gid.*He;

%-----------------------------------------------------------------------
% Передаточная характеристика выход-управление
Gvc=Fm.*Gvd./(1+Ti);

%-----------------------------------------------------------------------
% инвертирующий усилитель ошибки

R9=3.3E3;
C12=560E-12;
R10=20E3;	

AV0=95;		% усиление ОУ на постоянном токе, дБ
Fu=5.5E6;	% частота единичного усиления ОУ

gm=1E-3;
Rp1=(1/gm)*10^(AV0/20);
Cp1=gm/(2*pi*Fu);

% параллельное соединение Rp1, Cp1
Zp1=Rp1./(1+s*Rp1*Cp1);

% передаточная функция ОУ
AV=gm.*Zp1;

% последовательное соединение Cz и Rz
Zc=R9+1./(s*C12); 

%передаточная функция ОУ в инвертирующей схеме
KEA=-(Zc/R10)./(1+(1+Zc/R10)./AV);

%-----------------------------------------------------------------------

% Передаточная характеристика разомкнутой петли управления
TA=KEA.*Gvc;

% модуль передаточной характеристики в дБ
magTAdB=20*log10(abs(TA)); 

% фаза передаточной характеристики в гр.
phTAdeg= angle(TA)*180/pi;		

hold on;
figure(1);
semilogx(f,magTAdB,'linewidth',3);
ylim([-20, 80]);
grid on;
hold on;
figure(2);
semilogx(f,phTAdeg,'linewidth',3);
ylim([-180, 180]);
grid on;
%-----------------------------------------------------------------------