Erzeugung eines „unsymmetrischen“ Dreiecks aus Sägezahn ...

Prof. Dr. R. Kessler,C:\Si05\dgln\matlab\unsymdreieck.doc, S. 2/2
% Datei rdrei1.m Rundatei für Simulink-Datei drei1.mdl
% clear; fs=1; sw=0.2;tmax=5;dt=0.01;bild=2; rdrei1;
format compact;
if bild ==1, drei1; pause; end; % Simulink-Modell auf Bildschirm
sim('drei1');
figure(bild);clf;
plot(t,saeg,':m', t,drei,'--m', t,dreid,'m'); grid on;
xlabel('-> sec, Dreieck(-1..+1):glatt, Dreieck(0..+1)gestrichelt, Sägezahn: gepunktet');
axis([0,max(t),-1.1, 1.1]);
S1=['bild ',num2str(bild)]; S2=[', fs=',num2str(fs)]; S3=[', sw=',num2str(sw)];
S4=[', dt=',num2str(dt)];
tit=[S1,S2,S3,S4]; title(tit);
ylabel('Simulink:(rot), Matlab:schwarz');
disp('das war Simulink(rot), Taste! ');pause;
% jetzt mit Matlab:
Np=floor(tmax/dt); tp=zeros(Np,1); ysp=tp; dp=tp; ddp=tp;
k=0;t=0;
while t < tmax;
ys= t*fs - floor(t*fs);
usw=ys-sw;
d=abs(usw) * ( ( sign(usw) 0)/(1- sw) );
dd=(d-0.5)*2;
k=k+1; tp(k)=t; dp(k)=d; ysp(k)=ys; ddp(k)=dd;
t=t+dt;
end;
hold on;
plot(tp,ysp,':k', tp,dp,'--k', tp,ddp,'k');