Indem du das für n nacheinander 0, 1, 2 und 3 einsetzt.
Ist aber vermutlich nicht das, was du letztlich haben möchtest...

h[0] =    0
h[1] =    1
h[2] =  -3/2
h[3] =   3/4

Gruß
Chopin

Den Frequenz- und Phasengang als Bodeplot könnte ich noch anbieten,
berechnet mit Scilab:

 
z=poly(0,'z'); //Defines z to be a polynomial variable
hz = (z-1)/(z^2+z*0.5);
sys=syslin('d',hz); //Creates sys as a discrete-time ('d') syslin model.
fmin=0.001; //Min freq in Hz
fmax=0.499; //Max freq in Hz
scf(1);clf;
bode(sys,fmin,fmax); //Plots frequency response in Bode diagram
scf(2);
[tg,fr]=group(256,sys); //Plots group delay
plot2d("ln",(fr+.001),tg);

ob das Ergebnis stimmt, weiss ich nicht, nach dem Beispiel hier zum
Befehl "flts"

http://help.scilab.org/docs/5.3.3/en_US/flts.html

habe ich das noch an den Code oben angehängt:

scf(3);
u=1:10
y=flts(u,sys); //Plots pulse response
plot2d(y)

Das scheint eher die Step-response zu sein, hab die Beispiele nochmal
angeschaut, das sollte eher hinkommen:

imprep=flts(eye(1,20),tf2ss(sys));   //Impulse response
clf();
plot(imprep,'b')

jetzt stimmen auch die Werte mit Fredericks Zahlen überein.