我们在前面的多摆模型中,利用多摆的微分方程模型,求解出了多摆每时每刻的位置随时间的变化。当然那是一个高度复杂的非线性模型,难以上手分析。...这篇文章,我们首先利用一个二阶的线性模型进行求解,并引入微分方程定性分析中常用的工具——相图。...首先考虑下面这个经典的二阶阻尼振动方程:
将它整理为线性系统,如dx=Ax形式的样子:
矩阵A是一个二阶矩阵。我们取k=0.925,c=0.3。...(负实,无虚部)
% A=[0,1;
% -0.96,-2];
%只有阻尼(一个为0一个为负)
% A=[0,1;
% 0,-2];
%只有弹簧(实部都为0)
% A=[0,1;
%...4,4])
box on
xlabel('x');ylabel('dx');
function [F,Output]=Fdydx(x,y,Input)
%形式为Y'=F(x,Y)的方程,参见数值分析求解常系数微分方程相关知识