神马作文网二元二阶微分方程组求数值解的matlab程序并给出图像
来源：学生作业帮 编辑：神马作文网作业帮 分类：综合作业 时间：2024/11/12 17:35:54
二元二阶微分方程组求数值解的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-θ 图像
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-θ 图像