Question

MATLAB中ode45方法求解微分方程组

在MATLAB中将图片中的方程采用ode45来解，现在等式右边的式子不采用三角函数（采用三角函数的已会），而是采用对应的离散点数据输入，请问如何编程实现？如有不明白请追问，谢谢！

Answer 1

使用ode45解微分方程组，可以用下列格式来求解。

[t,x]=ode45(odefun,tspan,x0)

式中：odefun——给定微分方程组的自定义函数；tspan——自变量t的范围，即【t0,tf】;x0——x的初值向量，即x1（0），x1'（0），x2（0），x2'（0），x3（0），x3'（0）

由于没有给出自变量t的范围和x的初值，所以不能得到其x1，x2，x3的值。

追问

ODE45这个方法我会，上面的例子给出了初始条件和时间后，我也知道怎么编程序，结果也能算对，但上面图片的等式右边是一个与时间相关的具体函数表达式嘛，我现在的问题是如果右边是关于时间的离散点数据作为输入，这样的话，如何实现？？感谢！

追答

你的意思是指，2.0*sin(3.75*t)、-2.0*cos（2.2*t）和1.0*sin(2.8*t)为各为一组离散点数据。然后再求相应每组的x1，x2，x3值。如是这样，可以通过循环语句来实现。

Answer 2

我想问一下图片中的方程组，右边若采用三角函数该怎么编程序求解，能不能指导一下，谢谢