ラプラス変換

>>101
U--r--+---Y
　　　　 |
　　　 　c
　　　　 |
　　　GND

１次のLPF（c,r）の伝達関数を、状態変数で表した後にSCILABで求めよ。

与えられたLPFは
U=r*i + 1/c*∫idt
Y=1/c*∫idt
である。
x=∫idt （x'=i)なる状態変数を導入すれば、これらの式は
U=rx'+1/c*x
Y=1/c*x
となり、変形して

x'=-1/cr*x + 1/r*U
Y=1/c*x + 0*U
となる。（状態変数）
---------------------------------
これをSCILABのss2tfにかけると

-->r=1E6 ; c=1e-6;
-->A=[-1/(c*r)] ; B=[1/r] ; C=[1/c]; D=[0];
-->ss2tf(syslin("c",A,B,C,D))
ans =
　 1
-----
1 + s
となる。（終わり）
---------------------------------
ついでに逆変換。
-->sys=syslin("c",1/(1+%s))
sys =

1
-----
1 + s

-->{A,B,C,D}=abcd(tf2ss(sys))
D =

0.
C =

1.
B =

1.
A =

- 1.