% ##################################################################################
% ##  Loesung:  Entwurf eines Differenzieres (Remez-Verfahren)                    ##
% ##################################################################################

% #####  Teilaufgabe a:  Differenzierer entwerfen  #####
% ### 1:
f = [0 1];   m = [0 1];   b1 = remez(22,f,m);
f = [0 0.9]; m = [0 1];   b2 = remez(21,f,m);

figure; [H1 W] = freqz(b1,1,512); plot(W/pi,abs(H1) ); grid; xlabel('Omega/pi');
ylabel('|H(exp(j*Omega))|'); title('Filterordnung: 22, Typ I');
figure; [H2 W] = freqz(b2,1,512); plot(W/pi,abs(H2) ); grid; xlabel('Omega/pi');
ylabel('|H(exp(j*Omega))|'); title('Filterordnung: 21, Typ III');
figure; stem([0:22], b1); axis([0 22 -0.3 0.5]);       grid; xlabel('k');
ylabel('Impulsantwort');     title('Filterordnung: 22, Typ I');
figure; stem([0:21], b2); axis([0 21 -0.3 0.5]);       grid; xlabel('k');
ylabel('Impulsantwort');     title('Filterordnung: 21, Typ III');

% #####  Teilaufgabe b:  Uebertragungsfunktion vergleichen  #####
f = [0 1];   m = [0 1];   b2 = remez(21,f,m,'d');
f = [0 0.9]; m = [0 1];   b1 = remez(22,f,m,'d');

figure; [H1 W] = freqz(b1,1,512); plot(W/pi,abs(H1) ); grid; xlabel('Omega/pi');
ylabel('|H_D(exp(j*Omega))|'); title('Filterordnung: 22, Typ II');
figure; [H2 W] = freqz(b2,1,512); plot(W/pi,abs(H2));  grid; xlabel('Omega/pi');
ylabel('|H_D(exp(j*Omega))|'); title('Filterordnung: 21, Typ IV');
figure; stem([0:22], b1); axis([0 22 -0.4 0.5]);       grid; xlabel('k');
ylabel('Impulsantwort');       title('Filterordnung: 22, Typ II');
figure;stem([0:21], b2); axis([0 21 -0.4 0.5]);        grid; xlabel('k');
ylabel('Impulsantwort');       title('Filterordnung: 21, Typ IV');
% #####  EOF  #####