1

我们正在尝试使用 matlab 的 ode23 求解器对 DC/DC 降压转换器进行建模。当我们尝试运行我们的代码时，会出现以下错误：

??? Error using ==> odearguments at 91
The last entry in tspan must be different
from the first entry.

Error in ==> ode23 at 171
[neq, tspan, ntspan, next, t0, tfinal,
tdir, y0, f0, odeArgs, odeFcn, ...

Error in ==> buck2 at 13
      [t,x] = ode23(@event1,[t0 tf],x0);

当我们取出修改初始条件数组的代码时，代码运行没有错误，但不会产生预期的结果。

这是我们的代码：

function buck2
close all
clear all
clc

t0 = 0; 
tf = 1;
x0 = [0 0];  % Initial conditions

for i = 1:10
    if (1 <= i <= 4),
      [t,x] = ode23(@event1,[t0 tf],x0);
      nt = length(t);          
    else
      [t,x] = ode23(@event2,[t0 tf],x0);
      nt = length(t);   
    end   
    t0 = t(nt);
    x0 = x(nt);
end

plot(t,x)

function xdot = event1(t,x)
L = (12.12e-6);
C = (19.5e-6);
RC = (2.5*19.5e-6);
xdot = [(24/L) - (x(2)/L); (x(1)/C) - (x(2)/RC)];

function xdot = event2(t,x)
L = (12.12e-6);
C = (19.5e-6);
RC = (2.5*19.5e-6);
xdot = [-(x(2)/L); (x(1)/C) - (x(2)/RC)];
4

1 回答 1

0

这是问题所在：

在循环中，第一次迭代：

你打电话给以下

[t,x] = ode23(@event1,[t0 tf],x0);

t0 = 0；和 tf = 1;

然后进一步向下循环：

t0 = t(nt); %so t0 = 1;

下一个 for 循环迭代：

[t,x] = ode23(@event1,[t0 tf],x0);

换句话说：

[t,x] = ode23(@event1, [1 1] ,x0);

解决方案：

修改t0 = t(nt);或更新 tftf = t0 +1;

更新：

另外，您应该更正以下内容x0 = x(nt,:);

于 2012-04-09T03:29:55.063 回答