Rappresentazione con Scilab

[simo85]

Un saluto a tutti

Mi è richiesto di rappresentare con Scilab la trasformata inversa della funzione:



Per un tempo da 1 ns a 5ms con passi 100 ns.

Ho provato con il codice seguente:

Codice: Seleziona tutto

s=poly(0,'s');
N = [5*(0.22*s + 1000^2)/(1.034e-4)]
D = [s^3 + 4549709.865*s^2 + 2.901353965e10*s]
H = syslin('c',N./D);
t = (1e-9:100e-9:0.005);
plot2d([t',t'],[(csim('step',t,tf2ss(H)))',0*t'])

E l'immagine ottenuta è questa:

scilab.png (11.4 KiB) Visto 1134 volte

Forse sarà ora che mi compri un paio di occhiali nuovi, ma secondo me MATLAB dà un risultato più coerente con i calcoli fatti:

matlab.png (11.59 KiB) Visto 1134 volte

Codice: Seleziona tutto

syms s;
F = (5*(1.1*s + 1000^2)/(1.034e-4))/(s^3 + 4549709.865*s^2 + 2.901353965e10*s);
t = (0:100e-9:100e-6);
ezplot(ilaplace(F), [0,0.0001]), axis square, grid on, xlabel('time (t)');

Dove sbaglio con Scilab ?

Ringrazio in anticipo chi avrá voglia di aiutarmi.

[simo85]

so che non interessa a nessuno ma nel codice di matlab qui sopra c'era un coefficiente sbagliato al denominatore.

Codice: Seleziona tutto

syms s;
F = (5*(0.22*s + 1000^2)/(1.034e-4))/(s^3 + 4549709.865*s^2 + 2.901353965e10*s);
ezplot(ilaplace(F), [0,0.0001]), axis square, grid on, xlabel('time (t)');

Ad ogni modo l'espressione in Scilab e corretta e non riesco a risolvere.

[simo85]

Per il momento ho risolto direttamente introducendo l'antitrasformata, cosi:

Codice: Seleziona tutto

t=1e-9:100e-9:100e-6;
f = 1.66667+1.09963*10^(-6) *exp(-4.54332*10^6*t)-1.66667 *exp(-6385.97*t);
plot2d(t,f,1);


Ma l'intervallo di tempo non è quello richiesto... E se lo imposto cosi:

Codice: Seleziona tutto

t = (1e-9:100e-9:0.005);

Scilab da errore sul terzo argomento di plot2d.. Mah, devo rivedermi l'help della funzione plot2d.
Ma non alle 2:30 di mattina..