matlab中simulink的编码

已知y''''+3y'''+3y''+4y'+5y=exp(-3t)+[sin(4t+pi/3)]exp(-5t)
其中y(0)=1,y'(0)=y''(0)=0.5,y'''(0)=0.2
要求，写出程序码和思路（越详细越好……）
解释一下，exp(-3t)是以e为低的指数函数，变量是-3t，后面得同理，pi是π，我怕大家看不清楚，在此声明一下，写得好的还会加赏金的~
用simulink搭建模型

这么简单的东西
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文件这里没法上传。

具体设计方法楼上已经很明白了，我只好将simulink设计图画出来，一目了然。

注意为积分器设置初值。