求这个二元微分方程组的MATLAB程序

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/28 00:03:48

求这个二元微分方程组的MATLAB程序
求这个二元微分方程组的MATLAB程序

求这个二元微分方程组的MATLAB程序
假设x(1)=x1,x(2)=x1',x(3)=x2;把高阶导数换算成低阶倒数;得到3 个方程联立的方程组
1:x(1)'=x(2) 2:x(2)'=-2*x(2)-2*x(3) 3:x(3)'=-x(2)-2x(1)+x(3).
求解程序如下:
clear all
clc
f=@(t,x)[x(2);-2*x(2)-2*x(3);-x(2)-2*x(1)+x(3)];
x0=[1;-1;0];
[t,x]=ode45(f,[0,5],x0);
plot(t,x)
figure;
plot(x(:,1),x(:,3))
可以的到其时间曲线和相平面曲线.