貌似第一种写法的解是正确的，第二种写法的解是错误的。

henri

原方程为：

在用ode45解的时候，第一种写法
2 A: s% i# u! d  Jfunction dy=t1(t,y)
dy=zeros(2,1);
y(2)=dy(1);   
& c0 v7 K/ P  I6 K+ `5 a7 P$ ~$ zdy(2)=1;


当我把第三行等式交换位置后，即第三行改为dy(1)=y(2)
( |' E$ ]. \3 C  z
最后两种写法得到的结果完全不一样，貌似第一种写法的解是正确的，第二种写法的解是错误的。请教大神，感谢!

Ferrya

+ P2 A' C3 D# g( I0 h
先去搞清楚“赋值”的含义，等号在matlab里面是赋值，不是相等，不能左右互换
a = 1; b = 2; c = 3;
a = b;
disp([a b c])

+ S# \' L9 f6 V9 r# U: v% D$ F' S
结果是     2     2     3

Ferrya

a = 1; b = 2; c = 3;
8 C/ R* O+ @! L' B. H3 qb = a;
disp([a b c])
  B! J) r8 o$ M* [
结果是     1     1     3
. K% a2 {3 d5 L& W

Ferrya

正确的应当是
5 i7 G8 }3 v+ J/ T0 a2 B  Zfunction dy=t1(t,y)
dy=zeros(2,1);
dy(1)=y(2);
% M, C; ^- v! m! P0 N( Mdy(2)=1;