Ok, credo di aver trovato una soluzione quanto meno funzionante, la posto di seguito nel caso a qualcuno servisse in futuro (oppure nel caso qualcuno voglia sgridarmi per l'oscenità scritta ).
Non potendo calcolarmi la risposta di TF nel dominio del tempo in un dato istante e sollecitata da un determinato ingresso attraverso lsim(), ho deciso di intraprendere la strada più complicata ed antitrasformare TF nella funzione
f(t, u).
Successivamente, la valuto nell'istante
t e nell'ingresso
u che mi vengono passati, e gli e termino la funzione ritornando tale valore.
Il punto meno piacevole è stato scoprire che il comando ilaplace(), che serve per effettuare l'antitrasformata di Laplace, non funziona con l'oggetto TF in quanto facente parte della classe
tf mentre invece ilaplace() richiede un oggetto di tipo
syms.
Ho dovuto quindi lavorare un pò di tfdata() e poly2sym() per ottenere quanto desiderato.
Ecco il codice
[num, den] = tfdata(TF);
syms s;
TF_syms = poly2sym(cell2mat(num), s) / poly2sym(cell2mat(den), s);
iTF = ilaplace(TF_syms)
% trucchetto
iTF = @(t) subs(iTF, 't', t)
% ora è possibile calcorale iTF(t) con t valore scalare a piacere
return iTF(time);