matlab离散点数字微分,如何解已知离散数据的微分方程？

如下函数：已知x-t，t=0：188，x是一组离散点；要求解微分方程y''=m*x'+a*x-p-(m+n)*y'-(a+b)*y ，初始条件:x=0,y=0,y'=0；希望得到y''与t的figure.请问各位大虾如何解之。以下是我写的fcn函数，令z=y'，希望通过自己编写的四阶龙哥库塔算法求解，但是得到的结果效果不理想，附件中有我的相关数据和解法。还请各位高手给出正解~

function f=fcn(t,y,z)

q=0.150*160*pi;

a=28000/(q*160*pi);

b=25/(q*160*pi);

m=0.025/q;

n=0.0345/q;

p=200/(q*160*pi);

x=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2.2e-005, 0.002272,0.008222,0.017872,0.0312,0.046,0.0608,0.0756,0.0904,0.1052,0.12,0.13525,0.15337,0.17758,0.21122,0.25739,0.31871,0.39695,0.49282,0.60638,0.73763,0.88657,1.0532,1.2375,1.4393,1.6579,1.8919,2.1394,2.3983,2.6661,2.9404,3.2191,3.5,3.7815,4.0626,4.3419,4.6186,4.8915,5.1599,5.4229,5.6802,5.9315,6.1767,6.4157,6.6482,6.8742,7.0936,7.3061,7.5117,7.7103,7.9016,8.