Laplace transformatie numeriek bereken

[armageddon_2k1]

Okee, er is het volgende aan de hand.
Vanwege mijn interesse in snarenmodellering (gitaar) interesse in Systeem en Regeltechniek ben ik bezig met wat testprogrammaatjes maken in C# die enkele dingen numeriek voor me moeten berekenen. Zo ben ik bezig met een programmaatje waar ik simpele blokdiagrammen kan maken (a la Simulink, maar dan simpel) en zo inputs en outputs te plotten. Dit gaat aardig, alleen stuit ik nu op Laplace transformaties die erg veel gebruikt worden in deze diagrammen. Laplace kan ik met de hand makkelijk inverteren, maar ik kan me niet indenken hoe ik numeriek een Laplace output kan uitrekenen. Het liefst wil ik dat ik een input bijvoorbeeld met de joystick ingeef, en dat deze dan door een transferfuntie gehaald wordt en ik zo een output krijg.

Ik heb al op internet gezocht en onder andere de FFT (Fast Fourier Transform) Methode gevonden om het op te lossen, maar daar kan ik verder weinig over vinden.

Heeft er iemand hints, tips, links?

[mrClass]

Ik weet in ieder geval wel dat je een convolutie (En laplace kun je implementeren met een convolutie) kan uitvoeren door de fourier getransformeerde van het laplace filter te vermenigvuldigen met de fourier getransformeerde van je signaal. Na deze vermenigvuldiging kun je weer terug fourier transformeren om je laplace getransformeerde signaal te krijgen.

Dit principe heet het "convolutie theorema". Met het convolutie theorema kun je ook de inverse van een laplace transformatie doen (via het frequentie domein = fourier domein).

zo dus: invFFT(FTT(laplace_convolutie_masker) * FFT(signaal)) = convolutie_met_laplace

[armageddon_2k1]

Bedankt, ik zal er eens op zoeken.
Ik heb wel door dat al deze theorie toch redelijk aan de grens is van mijn huidige kunnen maarja, ik leer het wel van.

[caipirinha]

misschien heb je hier nog iets aan:

uit de labview 7.1 signal processing toolbox help

The discrete version of the Laplace transform is computationally expensive. An efficient strategy for the discrete Laplace transform is based on the Fractional Fast Fourier Transform (FFFT).

[armageddon_2k1]

Ja, de FFT is mijn inderdaad bekend (zie TS) maar evenzo bedankt.
Ik ga eens verder kijken op internet. Het punt is alleen dat veel van dit soort theorie behandeld wordt in boeken welke niet direct op internet te vinden zijn.

Wellicht eens een hoogleraar vragen om een duwtje in de goede richting.