Matlab Funktion schätzen?

Hallo Forengemeinde,

ich habe ein Feder-Massesystem bei dem die Feder eine Masse hält und an 
der Decke befestigt ist. Kann sich galube jeder vorstellen.
Mein kleines M-File für ein Feder-Masseystem sieht so aus:
%Parameter
m=50; %kg
k=10; %N/m

Die Masse soll von einer Kraft Fa angeregt werden. Diese Kraft ist 
sinusförmig. Eine weitere Kraft Fe soll nach(!) der Anregung die Kraft 
Fa kompensieren.

demnach gilt:

Leider ist da was schiefgelaufen daher für die Formel einfach auf den 
Anhang schauen ;)

Diese Gl lässt sich natürlich leicht nach Fe umstellen. Ich habe dann 
mal jede Gleichung in ein Subsystem gepackt und siehe da wie zu erwarten 
ich kann am Ausgang keine Auslenkung feststellen. (Siehe beigefügtes 
Modell).

In der Realität kenne ich aber die Kraft Fa gar nicht, das heißt ich 
muss das Fe nur aus der aktuellen Auslenkung und der Funktion selber 
erstellen.



Wie könnte ich vorgehen, um das Ziel (keine Schwingung der Masse) trotz 
Unkenntniss der Anregung zu erhalten?