[E_Tobi]

Ja, irgendwie so...lokal, direkt auf der Membran gilt quasi erst mal q(t)=v(t)*A - unter mehreren Annahmen:
- die durchstrahlte Fläche A ist erst mal zeitinvariant - die ganze Membran schwingt als ganzes
- die Schallschnelle v(t) ist auf der ganzen Membran homogen, und beschrieben genau auf der Membranoberfläche

Das bedeutet schon in sich dass das nur so lange geht wie die Welle nicht auf die Membran passt - andernfalls gelten beide Annahmen nicht mehr.

So, dann hörts aber auf 
Für die Impedanz brauchts Z=Schalldruck/Schallfluss. Schalldruck ist F/A. A ist wieder bekannt, aber wie hoch ist F? Wohl so groß wie anteilig die gesamte beschleunigte Masse aktio=reaktio spielen kann, aber nur solange es in Phase ist.

Ersatzmäßig müsste man "eine unendlich lange L-C-Kette" an die Membran anschließen, mit Luft-Federn, die auf Luft-Massen wirken, die auf Luft-Federn wirken, etc...
Analog zur Strahlungsimpedanz beim funken, quasi. E-Feld, H-Feld, E-Feld, etc...
Ich habe das Gefühl die Lösung ist relativ einfach, und liegt auch direkt da, aber ich sehe sie nicht 

Naja, in der Zwischenzeit habe ich über eine zusätzliche Runde "DGL-von-hinten-aufrollen" den Bassreflex-Masse-Feder-Schwinger reingeholt, aber bisher noch ohne Nachgiebigkeit der Luft im Port:

   

Klappt gut, die doppelte Impedanzspitze mit Phasenverlauf entsteht, die Luftgeschwindigkeit im Port passt, etc