用matlab怎样求分段二阶线性微分方程,并画出相图

heyong2002

请问各位高手,用matlab怎样求分段 二阶线性微分方程,并画出相图

分段积分求相图

heyong2002

公式dx(dx/dt)/dt + f(x)=0.1cos0.5t.其中f(x)分段线性

heyong2002

x>=1时.f(x)=x-1
-1<x<1.f(x)=0
x<=-1.f(x)=x+1

heyong2002

请问各位高手,用matlab怎样求分段 二阶线性微分方程,并画出相图
分段积分求相图公式
dx(dx/dt)/dt + f(x)=0.1cos0.5t.其中f(x)分段线性
x>=1时.f(x)=x-1
-1<x<1.f(x)=0
x<=-1.f(x)=x+1

happy

ode系列就能求解
画相图之用输出求解结果dx/dt 和 x，然后plot

heyong2002

但是每段x的初值不一样．x的取值也不同呀．
你的意思是不是对每个方程用一次ode

happy

这个函数不是可以直接表示出来吗？

heyong2002

高手,能不能把你的程序给我侃看看,参考

happy

1. function dy = fun1(t,y)
   
2. dy=zeros(2,1);
   
3. dy(1)=y(2);
   
4. dyy(y(1)>=1)=0.1*cos(0.5*t)-y(1)+1;
   
5. dyy(y(1)<1&y(1)>-1)=0.1*cos(0.5*t);
   
6. dyy(y(1)<=-1)=0.1*cos(0.5*t)-y(1)-1;
   
7. dy(2)=dyy;
   
8. 
   
9. [T,Y] = ode45(@fun1,[0 500],[2 0]);
   
10. plot(Y(:,1),Y(:,2))

复制代码

大概是这样，不一定对，没有仔细检查，你自己好好检查一下吧