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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
0 A; F7 f: [9 W& ~% k$ ~/ e2 Lfunction z=nd(a,v,n,r,rv)
A=a;
# M: p* ~( C8 X( P8 DV=v;
2 W- I/ f4 k& j. d, z2 Z6 _  _5 aN=n;
R=r;
7 H: j; m6 f: H% G; v1 FRV=rv;
for k=1:N
; R8 J' g- r  U    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
, z* V' l% p; x* T$ t5 z    z(1,k)=k;
; [# C7 c6 C: C- \6 D$ |/ i& ]1 L    z(2,k)=A(k);
end;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
) g& c; I" M/ h, S系统显示：
/ O$ }2 D; |4 ~! ^# e' ^' E??? Attempted to access A(2); index out of bounds because numel(A)=1.
$ k& j% f3 [( I8 X* }, _/ iError in ==> nongdu at 8
- `' O0 d, F# W" c% B  s1 D, p    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
4 G2 w$ L4 G" I) g( A# b
1 Z  p1 z- o7 P% ~1.请问为什么啊？
- N" {- H$ E1 A' i; Q4 `1 _2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
7 K2 Q7 H2 H- _9 [4 s看书上直接可以写成A（k），到底应该如何呢？
# P  [& K' y( T! F- j3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
3 [- T, b& P1 CA=a;
# B7 [  D* F6 K9 ]' MV=v;
( a# J7 g& L4 T8 ?3 HN=n;
: T) h0 d- E6 b, W% C! QR=r;
RV=rv;
for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
end
" z7 x! {2 Q- I/ Q9 o8 x4 j
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
0 h' z- U; L3 |( OV=v;
N=n;
R=r;
RV=rv;
  G0 V+ x8 [5 l/ E6 ?for k=1:N
    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)，这属于非法赋值
所以应该修改成：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
+ M2 [% d) v4 G" a1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
, W, v/ H- C) \4 u所以应该修 ...

, D! E" z8 Y7 K& b  A7 c谢谢谢  我终于明白了