假设我先定义会用到的符号 w y z
syms w y z
要求解的非线性微分方程组如下:
w'=w+y+z
y'=w*y*z
z'=w*y-z
我先用前面定义的符号组成一个3x1的D矩阵
D1是D矩阵的第一个元素，以此类推
D1=w+y+z
D2=w*y*z
D3=w*y-z
故D=[D1;D2;D3]
按照ode45的语法，我需要先把前面的变量转换成x(1),x(2),x(3) (我卡在这边)
像是w=x(1),y=x(2),z=x(3),
然后定义[email protected](t,x)[D1;D2;D3]
接着解微分方程组
[t,xa]=ode45(f,[0 100],[0 0.5 2])
有2个问题想请问版上大大，
一个是如何将我目前有的变量转换成x(1),x(2),x(3)
另一个问题是我前面描述的方法可以求解这个微分方程组吗？(跟ode45没有很熟)
针对第一个问题，我有试着在定义变量那边多定义一个x变量。
但是因为它是1*1的符号，所以我没有办法定义y=x(2),就卡住了...
麻烦大家帮忙了，谢谢