Volterra-Modell in Matlab

Hallo Mathematik- und Programmierfreunde,
ich wusste nicht so genau, wo ich meine Frage stellen soll, könnte auch 
in die DSP-Abteilung passen, aber weils um Programmieren geht (mit 
Matlab halt) frage ich mal hier.

Also, ich habe eine nichtlineare Übertragungskette. Davon möchte ich ein 
mathematisches Modell erstellen. Da es sich um ein diskretes System 
handelt, kann man gut mit digitalen Filtern arbeiten, und ich möchte das 
mal versuchen, als Volterra-Modell zu modellieren. Falls euch das nichts 
sagt:

im Prinzip arbeitet man beim Volterra-Modell genau gleich, wie bei einem 
FIR-Filter, aber die Impulsantwort ist mehrdimensional (angenommen der 
Eingang sei x[k]):

- es gibt ein konstantes Glied c
- es gibt ein lineares Glied: a[0]*x[k] + a[1]*x[k-1] + a[2]*x[k-2] + 
... das kennt man allgemein als FIR-Filter.
- es gibt ein quadratisches Glied: a[0][0]*x[k]^2 + a[0][1]*x[k]*x[k-1] 
+ a[0][2]*x[k]*x[k-2] + ... + a[1][1]*x[k-1]^2 + a[1][2]*x[k-1]*x[k-2] + 
...

und so weiter. Ihr seht, der Aufwand, um diese ganzen Summen 
auszuschreiben,  wächst exorbitant schnell! Ich möchte das mit Matlab 
automatisieren, aber mir Fehlt im Moment noch ein Ansatz, wie ich das 
irgendwie allgemeingültig formulieren kann, sodass ich einfach 
verschiedene Filter testen kann. Wisst ihr wie ich meine? gibt es etwa 
schon fertige Volterra-Modelle für Matlab?


Gruss Ansgar.