Some years ago we had reasonable success modelling lossy and dispersive transmission lines. The base was modeling the two-port with some controlled sources plus some equivalent lumped networks. One of them was used to model the desired Zo(s)=sqrt((R+Ls)/(G+Cs)). The other, the propagation function F(s)=exp(-theta(s)), with theta(s)=l*sqrt((R+Ls)(G+Cs)) was modeled with an ideal transmission line together with a shaping network.
the Taylor series at s=0. In our approach, we took a bilinear transformation from s to z plane
s=so*(1-z)/(1+z) with so = sqrt(RG/LC)
z=0, and fixing the first series term at z=1 and z=-1 (lastly this corresponds to fitting to f=inf and f=0 and the intermediate frequency so).
This performed much better than the previous techniques without having to do any kind of optimization, i.e. the procedure was explicit.
You may google the title "An explicit method for modeling lossy and dispersive transmission lines"
Pere
Pere