clc
syms L Vdc w w1 kp1 ki1 kp2 ki2 kid Lf Cf R real;
Hi=kp1+ki1/(i*w-i*w1)
Hv=kp2+ki2/(i*w-i*w1)
Tp=(Hi+i*kid)*Hv*Vdc;
Z=1/((i*w*Cf)+(i*w*Lf*R)/(i*w*Lf+R));
Y=1/Z;
Zp=(Hi*Vdc+i*w*L)/(1+Y*(Hi*Vdc+i*w*Lf)+Tp);
U=simplify(real(Zp))
V=simplify(imag(Zp))
mag=sqrt(U^2+V^2);%----公式8
mag=subs(mag,L,17.9e-3);mag=subs(mag,Vdc,300000);mag=subs(mag,kp1,0.075e-3);mag=subs(mag,ki1,0.094);
mag=subs(mag,w1,120*pi);mag=subs(mag,kp2,11.1e-3);mag=subs(mag,ki2,8.388);mag=subs(mag,kid,23e-6);
mag=subs(mag,Lf,0.43e-6);mag=subs(mag,Cf,17.7e-6);mag=subs(mag,R,10);
magdB=20*log10(mag); %幅频特性
xn=logspace(1,3,10000);
jizhun=subs(magdB,w,120*pi);
plot(xn,subs(magdB,w,xn)/jizhun)
- 1
- 2
前往页