二元二阶微分方程组求数值解的matlab程序并给出图像
来源:学生作业帮 编辑:搜搜考试网作业帮 分类:综合作业 时间:2024/07/28 04:01:39
二元二阶微分方程组求数值解的matlab程序并给出图像
![二元二阶微分方程组求数值解的matlab程序并给出图像](/uploads/image/z/15587477-53-7.jpg?t=%E4%BA%8C%E5%85%83%E4%BA%8C%E9%98%B6%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%E6%B1%82%E6%95%B0%E5%80%BC%E8%A7%A3%E7%9A%84matlab%E7%A8%8B%E5%BA%8F%E5%B9%B6%E7%BB%99%E5%87%BA%E5%9B%BE%E5%83%8F)
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-θ 图像
![](http://img.wesiedu.com/upload/b/c3/bc36597d55cae25278097a17672bfc8b.jpg)
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-θ 图像
![](http://img.wesiedu.com/upload/b/c3/bc36597d55cae25278097a17672bfc8b.jpg)