32 lines
1.7 KiB
Matlab
Executable File
32 lines
1.7 KiB
Matlab
Executable File
% ##################################################################################
|
|
% ## 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 #####
|