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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
! [8 I# a- L2 {, u+ \( d# M9 gfunction z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
R=r;
RV=rv;
3 g8 W0 v) `$ N/ S: xfor k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
5 c/ B# x7 Q- P; B- m: ^8 Z( N) s( z    z(2,k)=A(k);
% A9 Q8 w0 A! lend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
) A$ q& H$ K; V1 t# R系统显示：
# z3 T+ e0 P1 A! a9 W??? Attempted to access A(2); index out of bounds because numel(A)=1.
4 y1 [. g9 f! o8 A* }Error in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

1.请问为什么啊？
7 D; e% Q1 l7 j* c6 y7 I2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
! R# G6 S$ q  D+ y; C; a看书上直接可以写成A（k），到底应该如何呢？
  E7 z2 ]& D' ?3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
% B8 W& C( k& p6 @R=r;
RV=rv;
6 ?% ?- A# }9 l; @8 l& L. P  wfor k=1:N
$ N9 R, H" P' J. q( l7 M; K    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
7 \9 C2 q! R5 C$ l: Y- h  G% |    z(2,k)=A(k);
) |/ n- X, |: D. W+ ?5 v: G$ m+ Xend

  \: x3 _/ _6 J: N4 M/ N不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
& y/ w. H/ @  u8 a0 {! Y4 |function z=nd(a,v,n,r,rv)
A(1)=a;
5 n5 K& A+ Q% J" Z) kV=v;
9 n$ Q+ W6 S" |: g) `, e' W; KN=n;
9 n7 O, ?3 N$ }, e0 lR=r;
, `1 [: y& f: w8 KRV=rv;
/ C* s3 S7 r2 \7 t8 M' Cfor k=1:N
" Y7 l' m& f. M3 f* |0 M    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
end

dreamofyou

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

厚积薄发

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

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
( L3 b1 y# t2 O0 W1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
, ~  w: f" {* T% o所以应该修 ...

7 T8 b6 C8 i, \! C8 {0 U谢谢谢  我终于明白了