Rappresentazione con Scilab
Inviato: 7 apr 2015, 22:14
Un saluto a tutti 
Mi è richiesto di rappresentare con Scilab la trasformata inversa della funzione:
![\begin{aligned}
F(s) & = \left ( {4.8356\cdot10^{10} + 10638.30s \over s^3 + 4549709.865s^2 + 2.901353965\cdot10^{10}s} \right ) \\
& = \underbrace{{1.6667 \over s}}_{\text{Forced response}} + \underbrace{{1.0996\cdot 10^{-6} \over s + 4543323.893} - {1.6667 \over s + 6385.97}}_{\text{Natural response}}\\
\mathcal{L}^{-1}[F(s)] & =F(t)\\
& = 1.6667 + 1.0996\cdot 10^{-6}\text{e}^{-4543323.893 t } - 1.6667\text e^{-6385.97 t}\\
F(100 \mu) &= 0.78654247
\end{aligned}](/forum/latexrender/pictures/975e51254e6da58e022f21c9fbcd065d.png "\begin{aligned}
F(s) & = \left ( {4.8356\cdot10^{10} + 10638.30s \over s^3 + 4549709.865s^2 + 2.901353965\cdot10^{10}s} \right ) \\
& = \underbrace{{1.6667 \over s}}_{\text{Forced response}} + \underbrace{{1.0996\cdot 10^{-6} \over s + 4543323.893} - {1.6667 \over s + 6385.97}}_{\text{Natural response}}\\
\mathcal{L}^{-1}[F(s)] & =F(t)\\
& = 1.6667 + 1.0996\cdot 10^{-6}\text{e}^{-4543323.893 t } - 1.6667\text e^{-6385.97 t}\\
F(100 \mu) &= 0.78654247
\end{aligned}")
Per un tempo da 1 ns a 5ms con passi 100 ns.
Ho provato con il codice seguente:
E l'immagine ottenuta è questa:
Forse sarà ora che mi compri un paio di occhiali nuovi, ma secondo me MATLAB dà un risultato più coerente con i calcoli fatti:
Dove sbaglio con Scilab ?
Ringrazio in anticipo chi avrá voglia di aiutarmi.
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 1135 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 1135 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.
