Maybe I didn't understand what you are asking. Now it seems like you are only asking for a time-dependent model of a resonant system.
One simple approach is to model the string as a having a single resonant frequency (given by the length, tension, mass density, etc.) If the string is fixed at each end, the resonant mode is a pure sinusoid (or e^i(kx))
Now, when you pluck the string, the triangle shape (if that's what you want to start with) is a sum of sines and cosines- for a simple triangle function, the expansion is found here for a time series:
http://en.wikipedia.org/wiki/Triangle_wave
Your expansion is in terms of length (kx/L), though, not time (wt). Using an off-centered triangle function changes the coefficients, but the different amplitudes are fairly straightforward to calculate, and many programs will do it for you.
Now you need a dispersion relation- a way to relate the spatial wavelength to a temporal frequency. A simple one will do- for a stretched string it's here:
http://en.wikipedia.org/wiki/Dispersion_relation#Waves_on_a_string
That's the input. Now you multiply the Fourier series with the frequency response of the string. This is a basic method of linear systems analysis. Because there is only one resonant frequency, the frequency response of the string uses a (complex) Lorentzian function. You must put in the width of the resonant peak by hand- this is the same thing as saying you must put in, by hand, the dissipation.
To summarize so far, you start off with a string shape given as a Fourier series A_k * e^i(kx/L), convert it to A_w e^i(wt) by the dispersion relation, and multiply by the lineshape to give something like A_w e^i((2w-w_0) + ig(w-w_0))t, where w_0 is the resonance frequency and g the damping coefficient. So it's clear that off-resonant modes decay as A_w *e^(-g(w-w_0)t). I may have left out some factors of '2' and 'pi'.
All my reference books are in my office- Seigmann has a particularly clear derivation in 'Lasers"; I can check the specifics tomorrow.
Does this help?