function [koeff] = firbp3()
fa = 48000;
fu = 500;     
fo = 5000; 
N = 255;
f = fir1(N, [fu/(fa/2), fo/(fa/2)]);
l = 1;
for k = 1:N+1
    koeff(1,l)   = f(1,k);
    koeff(1,l+1) = f(1,k);
    l = l + 2;
end
dlmwrite('c:\Users\Jürgen\Studium\5.Semester\Signalprozessor\APL\Bandpass\koeffizienten.dat',koeff,'delimiter','\n','precision','%.15f')
