如何使用Matlab求解复杂的微分方程组

2024-11-01 11:35:57
推荐回答(1个)
回答1:

这说明此微分方程没有显示解,应用数值方法求解,例如ode45等函数

就一个例子(选自matlab帮助文档)

解入下微分方程组

1,建立方程组函数

 

function dy = rigid(t,y)
dy = zeros(3,1);    % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);

2,求解并绘图

[T,Y] = ode45(@rigid,[0 12],[0 1 1]);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')

结果