Matlab程序调用时出现错误，求高手指导一下（关于ode23）

yuexiaoming

调用程序如下（ m文件）：

1. function dy = fun1(t,y,flag,u)
   
2. if y(1)>=1
   
3.     a=0.1*cos(u*t)-y(1)+1;
   
4. elseif y(1)<1&y(1)>-1
   
5.     a=0.1*cos(u*t);
   
6. elseif y(1)<=-1
   
7.     a=0.1*cos(u*t)-y(1)-1;
   
8. end
   
9. dy=zeros(2,1);
   
10. dy=[y(2);a];

复制代码

主程序：

1. u=[0.5, 2.0];
   
2. for j=1:2
   
3. [t,y] = ode23(@fun1,[0 5] , [2 0],[],u(j));
   
4. figure
   
5. plot(y(:,1),y(:,2));
   
6. end

复制代码

yuexiaoming

自己顶一下

牛小贱

正确的程序应该是这样的：
调用程序如下（ m文件）：

1. function dy = fun1(t,y,flag,u)
   
2. u=[0.5,2.0];
   
3. for j=1:2
   
4.     if y(1)>=1
   
5.     a=0.1*cos(u(j)*t)-y(1)+1;
   
6. elseif y(1)<1&y(1)>-1
   
7.     a=0.1*cos(u(j)*t);
   
8. elseif y(1)<=-1
   
9.     a=0.1*cos(u(j)*t)-y(1)-1;
   
10.     end
    
11. end
    
12. dy=zeros(2,1);
    
13. dy=[y(2);a];

复制代码

主程序：

1. [t,y] = ode23(@fun1,[0 5] ,[2 0]);
   
2. figure
   
3. plot(y(:,1),y(:,2));

复制代码

另外，补充一点：ode23和ode45在 MATLAB中的用法

ode23:使用二阶龙格-库塔法求微分方程，调用格式为：
[t，y]=ode23(‘fname’,tspan,y0)
ode45:使用四阶龙格-库塔法求微分方程，调用格式为：
[t，y]=ode45(‘fname’,tspan,y0)

其中fname为由M函数定义的线性或者非线性微分方程的句柄函数名。tspan形式为[t0,tf]，表示求解区间。y0是初始状态列向量。t和y分别给出时间向量和相应的状态向量。

希望对你有所帮助！！！