IIR Filter auf Matlab realisieren.

[filterMonster (Gast)]

Ich schaffe es nicht einen IIR Filter zu ploten. Ein funktionierenden 
Frequenzganz ist möglich siehe code

h1 = zeros(1,8);
h1(1,1:3) = [0.5 1 0.5];  

N = 1024;
H1 = fft(h1,N);
f = linspace(-1/2,1/2,N);
roots(h1)
figure(1)
plot(f,fftshift(abs(H1)),'g');
hold on;
grid on;
xlabel('f/f_a'); ylabel('H(f)');
title('Frequenzgang');

figure(2)
plot(f,fftshift(20*log10(abs(H1))),'g');
hold on;
grid on;
xlabel('f/f_a'); 
ylabel('H(f)');
title('Frequenzgang dB');

Ein FIR Filter funktioniert wunderbar. Teste den Code. Ich kann aber 
keinen IIR FIlter programmieren. ich habe es mit tf probiert davon kann 
man aber keine fft machen. wie kann ich den einen gebrochen Rationale 
Funktionen. Einen BRuch darstellen. Dann könnte ich den IIR Filter 
realisieren... Ich weiß gerade nicht weiter.

Eigentlich will ich nur diese stelle auf einen Gebrochen rationale 
funktion ändern siehe untern

h1 = zeros(1,8);
h1(1,1:3) = [0.5 1 0.5];

Das ich an dieser Stelle des Codes einen bruch realisieren kann

$$\frac{ b0 *z^{ -1 } + b1*z^{ -2 }}{1 + a1*z^{ -1 } + a2*z^{ -2 } }$$

[Detlef _A (Gast)]

Monster,

so geht das.

Cheers
Detlef

clear

fb=[1 -1 1];
fa=[1 -1.5 1];
N = 1024;
H1=freqz(fb,fa,N);
f=1:N;
figure(1)
plot(f,(abs(H1)),'g');
xlabel('f/f_a'); ylabel('H(f)');
title('Frequenzgang');

figure(2)
plot(f,(20*log10(abs(H1))),'g');
xlabel('f/f_a');
ylabel('H(f)');
title('Frequenzgang dB');

return