
% Da Filter implementieren mit Fdatool und code generieren.
%*************************************************************
Fs = 60000;
N    = 20;   Fc   = 20000;      
flag = 'noscale';  % Sampling Flag
win = hamming(N+1);
% Calculate the coefficients using the FIR1 function.
b  = fir1(N, Fc/(Fs/2), 'low', win, flag);
%**************************************************************
%**************************************************************


% Anregungssignal. Es handel sich hiermit um eine MLS der Ordnung 3.
seq = [1 1 -1 -1 1 -1 1];


% Filter mit dem Signal anregen und in variable input speichern
% Also seq mit b falten
input = conv(seq, b);


% Ausgangsignal nach der Kreuzkorrelation von dem Anregungssignal seq mit
% dem gefilterten Anregungssignal input. Werte in Variable ouput speichern
output = xcoor(seq, input);


% Variablen darstellen. Zu merken ist dass ich das Filter in output
% erkennen kann. Vergleichen sie selbst
figure; plot(b); grid on;
figure; plot(input); grid on;
figure; plot(output); grid on;

