求助：差分方程在matlab的实现问题

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
  c9 w- b, \/ a; T6 O/ C" L# GA=a;
5 J3 j6 f  S5 A- l+ NV=v;
% V" e. T, Z4 q3 D' Y. \N=n;
, Z- @8 {+ J: z% K8 IR=r;
RV=rv;
( |  ?$ j3 R, _+ Sfor k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
/ ^7 S% m* A6 n6 j2 A    z(1,k)=k;
5 h# M. p3 e7 o* W1 G8 a    z(2,k)=A(k);
# I4 g% O. s% f7 c. z3 b' yend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
* J% B: \4 z% \0 u0 {& r9 Z系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
) Q" H& H, |8 L" `" aError in ==> nongdu at 8
- x* ~& m" j/ k    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

: X7 k" _5 S" f$ m1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
; |9 ]5 b% N; i: W看书上直接可以写成A（k），到底应该如何呢？
9 _7 W. m: ]. `3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
( a3 x; A# d1 a; h% @A=a;
, b+ z5 m4 Z% n: B# qV=v;
3 m% W/ ^7 @% WN=n;
R=r;
RV=rv;
for k=1:N
- n. v: ^' _0 G# `. C    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
3 L6 g( j5 Y: X& o: i    z(1,k)=k;
    z(2,k)=A(k);
end
! [! g' ~# B! d0 L
; h) c9 M3 r, i% l  U0 [, G不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
- D" x* v  a* H% G8 {$ r, OV=v;
# l# }! X4 z$ y& MN=n;
# T1 @" M& Z7 ]$ s8 ?9 Q3 NR=r;
; m. T, }* h6 n9 Z" d' W; k2 [RV=rv;
% B- x+ {! r  d& e, Dfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
* z9 ]9 p, j7 p, _) P    z(1,k)=k;
    z(2,k)=A(k);
- J3 ~( I* H6 n) y' Y% c+ i6 M$ nend

dreamofyou

唉，谢谢，我再试试，matlab里边用差分方**的很诡异啊。

厚积薄发

1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
/ V  c' b) Q# |0 r5 h. Q9 @; \" O1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
所以应该修 ...

谢谢谢  我终于明白了