clear;
RE = 1E3;  % Innenwiderstand der Quelle in Ohm
RA = 1E4;  % Lastwiderstand in Ohm
C  = 1E-5; % Glaettungskapazitaet in Farrad
UF = 0.3;  % Flussspannung der Dioden

UEmax = 3; % Amplitude der Eingangsspannung in V
f  = 100;   % Frequenz des Eingangssignals in Hz
dt = 2E-4; % Abtastintervall in s
N  = 100;  % Abtastwerte je Periode
NP = 3;    % Anzahl der simulierten Perioden
dt = 1/(f*N);

ua(1)=0;   % Anfangsspannung des Kondensators in V

t= (0:N*NP)*dt;
ue = 3* sin(2*pi*f*t);

for n=1:N*NP
    if ue(n)-ua(n)>2*UF         % Arbeitsbereich I
        ic = (ue(n)-ua(n)-2*UF)/RE - ua(n)/RA;
    elseif -ue(n)-ua(n)>2*UF    % Arbeitsbereich II
        ic = (-ue(n)-ua(n)-2*UF)/RE - ua(n)/RA;
    else                        % Arbeitsbereich III
        ic=-ua(n)/RA;
    end;
    ua(n+1)=ua(n)+dt/C*ic;
end;

plot(t, ue, 'LineWidth', 3, t, ua, 'LineWidth', 3);
set(gca,'FontSize',24); 
xlabel('t in s'); ylabel('u in V'); 
title('Aufgabe 7.2: Brueckengleichrichter mit Glaettungskondensator');