1/s vs. bilineare Transformation

Hi,

In einem sehr einfachen (analogem!) Simulink Modell erstelle ich ein 
Signal mit 5 Sinusen (Frequenzen 34, 334, 106, 149, 448) gefolgt von 
einem Integrator. Die Simulation dauert eine 1-1/1000 Sekunden und 
"simout" wird mit 1/1000 sub-gesampled.

Das gleiche Signal erstelle ich in einem 1000-Elemente Vektor und 
filtere das Signal mit einem Integrator per bilinearer Transformation:

[numd, dend] = bilinear([1], [1 0], 1000);
x2 = filter(numd, dend, x);

Zeige ist mir nun das Spektrum beider Signale per fft, so ist dieses 
nicht ident:

stem(abs([fft(simout.signals.values) , fft(x_I)]))
http://snag.gy/8srWk.jpg

Hat jemand eine Idee wo der Fehler liegen koennte?
Danke,
Peter