ly36.m

来自「本人非线性控制理论的大作业,绝对经典,完全可以运行,全力打造,真心奉献,希望会对」· M 代码 · 共 45 行

function dx1=ly36(t,x)
dx1=zeros(1,1);
%x1=x(1)
%x2=x(2)
dx1(1)=[5-x(1)]^5;
%dx1(2)=4*x(1)^2*x(2)+x(2)*(x(1)^2+x(2)^2-2);

%[t,x1]=ode45(@liyan,[0 100],[3 4])
%figure
%hold on
%plot(t,x1)

%[t,x1]=ode45(@liyan,[0 100],[3.01 4.01])
%figure
%hold on
%grid
%plot(t,x1)

%x1'=x2
%x2'=6*sin(t)-0.1*x2-x1^5

%一个二阶微分方程：
%y''+y'+y=sin(t)

%初始条件为y(0)=5,y'(0)=6。

%过程：
%先降阶为一阶微分方程组
%y'=z
%z'=-z-y+sin(t)

%编制如下函数m文件
%function dy=liyan(t,x)
%dy=zeros(2,1);
%y=x(1)
%z=x(2)
%dy(1)=x(2);
%dy(2)=sin(t)-x(2)-x(1);

%*******************
%然后用ode45解方程
%[t,y]=ode45(@liyan,[0 20],[5 6])
%plot(t,y)就可以画出y和y'的图像