请教大神，这个是为什么呢？

paprill

" o9 x8 _- d5 B* a) B" H
请教大神，这个是为什么呢？？ode45 无法执行赋值，因为左侧的大小为 2×4，右侧的大小为 2×5？非常感谢！
6 H  B4 H9 L3 {: m# c; C
; d# M0 n% D  v  q+ y  tfunction  f=boat1(t,x)
if(norm(x)>1e-5)   %限制条件，保证在x足够小的情况下可以结束计算
; Y2 t# Y/ ^. Y    f=[-2*x(1)/sqrt(x(1).^2+x(2).^2);1-2*x(2)/sqrt(x(1).^2+x(2).^2)]
% ^5 u( l+ a: m1 ]/ nelse
    f=[0,0]

$ |. x; N" v, |$ e6 P6 g' Fend


y0=[1000 0];      %微分方程的初值
[t,x]=ode45(@boat1,[0,1000],y0)      %求解微分方程的数值解
subplot(1,2,1),plot(t,x)；         %画出关于x(1)和x(2)的系统时间响应曲线
# n! d) r& y0 X- m9 t! ?plot(x(:,1),x(:,2))          %画出关于x(1)和x(2)的系统相平面曲线
axis([-1 1000,0 1000])          %限定区间
5 h. \% w9 p) J5 p. S$ X6 u. ]text(10,15,'B')         %标注B点
text(1000,0,'A')         %标注A点
* j5 A' W) h5 N, _4 C7 D  d3 x" m% A, Axlabel('X')
ylable(‘Y’)
$ J+ Z) N- M) G5 R2 X: M: T
. S' D4 P! w9 w7 I: p
无法执行赋值，因为左侧的大小为 2×4，右侧的大小为 2×5。

* H  h: H) Y' P; k; Z/ ~0 P出错 ode45 (line 488)
1 F0 U4 Q' k' [( f4 n# M# u
        yout(:,idx) = yout_new;

Demyar

你的代码有点问题。