clear;clc;
Uin=300;                                        %unit: V
Uomin=40000;                              %unit: V
Uomax=120000;                           %unit: V

P=1600;
fr=60000;
m=10;
Io=0.03;                                          %unit: mA
Uonom=50000;

n = Uin/(Uonom/m)                       %变压器变比
Mmin= n*(Uomin)/m/Uin
Mmax= n*(Uomax)/m/Uin
k=2.5;
Qmax=1/(k*Mmax)*sqrt(k+Mmax*Mmax/(Mmax*Mmax-1))
Q=0.9*Qmax

fmin=fr/(sqrt(1+k*(1-1/(Mmax*Mmax))))
tmax=1/fmin
fnmin=fmin/fr
fmax=fr/(sqrt(1+k*(1-1/Mmin)))
tmin=1/fmax
fnmax=fmax/fr

Ro=(Uonom/Io)/1000   %unit: kohm
Rac1 = (((n^2)*(Ro))/(2*(m^2)))*1000   %unit: ohm,不考虑倍压电路压降

Urpp=0.001*Uonom;
Co = (10^9)*1.3*(m*Io)/(2*3*fmin*Urpp)  %unit:  nF
C = Co/(10^9);
temp1 = Ro/(8*(m*m));
temp2 = (1/(4*m*m*fr*C))*((m^3)/24+(m^2)/8+m/3);
temp3 = (1/(8*(fr^2)*(C^2)*Ro))*((m*m/24+m/8+1/3)^2);
Rac2 = n*n*(temp1 + temp2+ temp3)

Rac = Rac1
Zo=Q*Rac

Cr=1/(2*pi*fr*Zo)
Lr=Zo/(2*pi*fr)
Lm=k*Lr

Im=Uin/(4*fmax*(Lr+Lm))

Coss=20e-12
Td=83.2e-9
Ip=(2*Coss)*Uin/Td
Io=P/Uonom;