這麼簡單的東西

y"""" = (y""")"

y""" = (y"")"

y"" = (y")"

y" = (y)"

Y = [y y" y"" y"""]

Y" = [y" y"" y""" y""""]

Y" = [ 0 1 0 0; * Y +[0 * U

0 0 1 0; 0

0 0 0 1; 0

-5 -4 -3 -3] 1]

U = wxp(-3t) + [sin(4t-pi/3)]exp(-5t)

初值 [0.2 0.5 0.5 1]

如果要求的是y

則

z = [1 0 0 0]Y

系統的a = [ 0 1 0 0; b = [0 c =[1 0 0 0] d =0

0 0 1 0; 0

0 0 0 1; 0

-5 -4 -3 -3] 1]

自己寫下程式很簡單的。

simulink檔案這裡沒法上傳。