Simulation Thyristor an induktiver Last an Wechs - Home.hs ...

Prof. Dr. R. Kessler, HS-Karlsruhe, C:\ro\Si05\Andy\tephys\Bahm4\Thyri\Thyri3.doc, S. 2/3
++++++++++++++++ jetzt Matlab: +++++++++++++++++++++++++++++++++
% Datei Thyri2.m
% clear;L=1;R=1;del=90;f=1;a0=2; dt=0.002; tmax=5;bild=1;Thyri2;
format compact;
% Tephys-Datei Thyri_3.txt (September 1997, Heft R23,Seite 52
%t=t+dt
%u0=a0*sin(2*pi*f*t)
%uSt0=ja(u0)*ja(sin(2*pi*f*t-del*pi/180))*ja(180-del)*ja(del)
%uCSt=uCSt+(uSt0-uCSt)*dt*f*50
%uSt=(uSt0-uCSt)*ja(uSt0-uCSt)
%ialt=i
%FF=ja(ja(uSt)+FF*ja(ialt))
%i=FF*(i+(u0-R*i)*dt/L)
Np=floor(tmax/dt);
tp= zeros(Np,1); u0p=tp; ip=tp; FFp=tp; uStp=tp; uSt0=tp;
%Startwerte:
i=0; uCSt=0; FF=0; t=0;k=0;
while t < tmax
u0=a0*sin(2*pi*f*t); % u0 = Netzspannung
uSt0= (u0>0)*((sin(2*pi*f*t-del*pi/180)) >0) *(180>del)*(del>0) ; % uSt0=Steuersignal,
% wird 1, wenn u0 >0 ist und wenn phasenverschobenem sinus > 0 wird,
% del muss im Bereich 0< del < 180 liegen dies ist neu , in Tephys noch nicht drin
uCSt=uCSt+(uSt0-uCSt)*dt*f*50; % Das Rechteck uSTt speist einen RC-Tiefpass
uSt= (uSt0-uCSt)* ((uSt0-uCSt) > 0); % ust sind die positiven Impulse des Hochpasses RC
ialt=i;
FF=( ((uSt > 0) + FF * (ialt >0 ))) >0;
i=FF * (i+ (u0 - R*i)* dt/L);
% Plotwerte speichern:
k=k+1; tp(k)=t; FFp(k)=FF; ip(k)=i; u0p(k)=u0;
uSt0p(k)=uSt0; uCStp(k)=uCSt; uStp(k)=uSt;
t=t+dt;
end;
figure(bild); clf; ofs=2; fa='k';
plot(tp,u0p,fa, tp,uSt0p+1*ofs,fa, tp,uStp+2*ofs,fa, tp,FFp+3*ofs,fa, tp,ip+4*ofs,fa ); grid on;
text(max(t),0,' u0'); text(max(t),ofs,' uSt0'); text(max(t),2*ofs,' uSt');
text(max(t),3*ofs,' FF'); text(max(t),4*ofs,' i');
S1=['bild ',num2str(bild)]; S2=[', L=',num2str(L)]; S3=[', R=',num2str(R)];
S4=[', del=',num2str(del)]; S5=[',a0=',num2str(a0)];S6=[', f=',num2str(f)]; S7=[', dt=',num2str(dt)];
tit=[S1,S2,S3,S4,S5,S6,S7]; title(tit);