二元二阶微分方程组求数值解的matlab程序并给出图像

如题所述

function dz=yueh(t,z)                   %保存为 yueh.m 文件

dz=[z(2);

    10-z(1)-0.001*z(3);

    z(4);

    -0.0001*z(1)-0.001*z(3)];

-------------------------------------------------------

运行一下：

[t,z]=ode45('yueh',[0 1],[0 0 0 0]); 

% 感觉初始值不完整啊，这里假设x(0)=θ(0)=x'(0)=θ'(0)=0

plot(z(:,1),z(:,3))          %    这里是 x-θ 图像