% MATLAB, Herbert Weidner
r = -0.7; %Reflexionsfaktor an beiden Enden des Wellenleiters
steps = 1/80; zeit = (0.: steps : 3.);
wave = sin( 2. * pi * zeit);
%huell=zeros(1,length(zeit)); %nur bei init !
hinlaufend = sin (2. * pi * zeit - 2. * pi * phi);
ruecklaufend = -r*sin (2. * pi * zeit + 2. * pi * phi);
sum = hinlaufend + ruecklaufend;
plot(zeit,hinlaufend, zeit,ruecklaufend, zeit,sum, zeit,huell)
ylim([-2 2])
active.BackgroundColor = [0 0 0];
axis tight
set(gcf,'nextplot','replacechildren','visible','off')
f = getframe; k=1;
[im,map] = rgb2ind(f.cdata,256,'nodither');
for phi=steps:steps:1
hinlaufend = sin (2. * pi * zeit - 2. * pi * phi);
ruecklaufend = -r*sin (2. * pi * zeit + 2. * pi * phi);
%vertauschtes VZ wegen Schreibrichtung links->rechts
sum = hinlaufend + ruecklaufend;
%huell=max(huell,sum); %nur bei init !
plot(zeit,hinlaufend, zeit,ruecklaufend, zeit,sum, zeit,huell)
ylim([-2 2]); f = getframe;
im(:,:,1,k) = rgb2ind(f.cdata,map,'nodither'); k=k+1;
end