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

matlab用ode45求解微分方程组~

这个方程要解出来还要初始条件吧
不妨设(x(0),y(0))=(A,B);
还有用ode45只能得到数值解，不能给出解析式表达
那么先建立一个函数文件
test_fun.m
function dx=tsst_fun(x,t)
dx(1)= c-b*x(1)-w*x(1)*x(2);
dx(2)=w*x(1)*x(2)-(b+m)*x(2);
那么调用ode45
x0=(0,0);
t0=0:0.01：10；
；%根据你所想要的求得值设定t0，间隔是任意的，与求解所用的步长无关，
[x,t]=ode45(@test_fun,t0,x0);
得到了
这样可以么？

例如：求解下列微分方程组

求解步骤：
①自定义函数 rigid。m
function dy = rigid(t,y)dy = zeros(3,1); % a column vectordy(1) = y(2) * y(3);dy(2) = -y(1) * y(3);dy(3) = -0.51 * y(1) * y(2);
②在执行窗口下，执行下列命令
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);[t,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
③求解结果

使用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值。如是这样，可以通过循环语句来实现。

---

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