Wir untersuchen numerisch und teilweise analytisch die kompakte Quantenelektrodynamik auf dem Gitter mit Wilson-Fermionen. Dabei konzentrieren wir uns auf zwei wesentliche Teilprobleme der Theorie: der Einfluss von Eichfeld-Moden mit verschwindendem Impuls in der Coulomb-Phase und die Effizienz von verschiedenen Monte-Carlo-Algorithmen unter Berücksichtigung dynamischer Fermionen. Wir zeigen, dass der Einfluss der Null-Impuls-Moden auf die eichabhängigen Gitter-Observablen wie Photon- und Fermion-Korrelatoren nahe der kritischen chiralen Grenzlinie innerhalb der Coulomb Phase zu einem Verhalten führt, das vom naiv erwarteten gitterstörungstheoretischen Verhalten abweicht. Diese Moden sind auch für die Abschirmung des kritischen Verhaltens der eichinvarianten Fermion-Observablen nahe der chiralen Grenzlinie verantwortlich. Eine Entfernung dieser Null-Impuls-Moden aus den Eichfeld-Konfigurationen führt innerhalb der Coulomb-Phase zum störungstheoretisch erwarteten Verhalten der eichabhängigen Observablen. Die kritischen Eigenschaften der eichinvarianten Fermion-Observablen in der Coulomb-Phase werden nach dem Beseitigen der Null-Impuls-Moden sichtbar. Der kritische Hopping-Parameter, den man aus den invarianten Fermion-Observablen erhält, stimmt gut mit demjenigen überein, der aus den eichabhängigen Observablen extrahiert werden kann. Wir führen den zweistufigen Multiboson-Algorithmus für numerische Untersuchungen im U(1)-Gittermodell mit einer geraden Anzahl von dynamischen Fermion-Flavour-Freiheitsgraden ein. Wir diskutieren die geeignete Wahl der technischen Parameter sowohl für den zweistufigen Multiboson-Algorithmus als auch für den hybriden Monte-Carlo-Algorithmus. Wir geben theoretische Abschätzungen für die Effizienz dieser Simulationsmethoden. Wir zeigen numerisch und theoretisch, daß der zweistufige Multiboson-Algorithmus eine gute Alternative darstellt und zumindestens mit der hybriden Monte-Carlo-Methode konkurrieren kann. Wir argumentieren, daß eine weitere Verbesserung der Effizienz des zweistufigen Multiboson-Algorithmus durch eine Vergrößerung der Zahl lokaler Update-Schleifen und auch durch die Reduktion der Ordnungen der ersten und zweiten Polynome zu Lasten des sogenannten 'Reweighting' erzielt werden kann. / We investigate numerically and in part analytically the compact lattice quantum electrodynamics with Wilson fermions. We studied the following particular tasks of the theory: the problem of the zero-momentum gauge field modes in the Coulomb phase and the performance of different Monte Carlo algorithms in the presence of dynamical fermions. We show that the influence of the zero-momentum modes on the gauge dependent lattice observables like photon and fermion correlators within the Coulomb phase leads to a behaviour of these observables different from standard perturbation theory. These modes are responsible also for the screening of the critical behaviour of the gauge invariant fermion values near the chiral limit line. Within the Coulomb phase the elimination of these zero-momentum modes from gauge configurations leads to the perturbatively expected behaviour of gauge dependent observables. The critical properties of gauge invariant fermion observables upon removing the zero-momentum modes are restored. The critical hopping-parameter obtained from the invariant fermion observables coincides with that extracted from gauge dependent values. We implement the two-step multiboson algorithm for numerical investigations in the U(1) lattice model with even dynamical Wilson fermion flavours. We discuss the scheme of an appropriate choice of technical parameters for both two-step multiboson and hybrid Monte Carlo algorithms. We give the theoretical estimates of the performance of such simulation methods. We show both numerically and theoretically that the two-step multiboson algorithm is a good alternative and at least competitive with the hybrid Monte Carlo method. We argue that an improvement of efficiency of the two-step multiboson algorithm can be achieved by increasing the number of local update sweeps and also by decreasing the orders of first and second polynomials corrected for by the reweighting step.