请教：matlab微分方程数值解法问题

东南小忽悠

式中：p(t) 是随时间任意变化的函数。。如 记录的地震波函数 ，应该如何求解u，u'，u" 随 时间 t 的变化曲线（有没有现成类似的matlab程序）

如：对于固定时间间隔，p(t)的值如下
1.154088e-002
6.341486e-003
2.165699e-003
-3.850290e-003
-7.415056e-003
-4.054916e-003
-1.494242e-003
1.891596e-003
4.060150e-003
6.296754e-003
8.556200e-003
7.682497e-003
1.081660e-003
1.192538e-003
3.817455e-003
3.654230e-003
5.246026e-003
6.217283e-003
3.837917e-003
6.381460e-004
-7.182354e-003
-1.569858e-002
-2.531741e-002
-3.075949e-002
-2.499476e-002
-1.976920e-002
-1.907252e-002

hsfy919

回复 1 # 东南小忽悠 的帖子

这是基本的振动力学问题，请参阅《振动力学》相关书籍，matlab没有提供直接求解函数，需要理解求解原理后，自己编程。不过应该是不难实现的

youxianyuan

差分。离散后解，
要不自己编 一个ODE45就好了，
这个步长你就可以控制了