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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
$ k; U9 D1 r5 c& ~; IA=a;
! s% n; t& \* o1 bV=v;
. a* h6 d. u' Z: j8 RN=n;
R=r;
( Q' a2 S! r8 Q- `  ZRV=rv;
2 b4 t3 _  I! O; b% `for k=1:N
' I' Y' A8 m. A) A. K6 _2 W1 q    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
4 }) k% w; M( x+ B6 s$ ]    z(1,k)=k;
, N. f' Y6 T  p8 ~7 d- g0 J9 r, O    z(2,k)=A(k);
9 q; ]5 ]$ b4 X- i' hend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
5 N( E% ]( Z6 z$ c! m2 f+ M系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
Error in ==> nongdu at 8
3 |( L4 x$ o0 O    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

5 m% C# Y" \2 _; m0 O/ }1.请问为什么啊？
7 f0 R& H4 B! S! G2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
+ j2 @' i  d$ z4 Q看书上直接可以写成A（k），到底应该如何呢？
7 k, U5 N/ O7 Q2 P  p' d3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
7 m7 A  |: u& P" z+ x0 Z2 R; M# o( iV=v;
$ x8 n- @. y& m* e/ Y7 XN=n;
7 O0 D, @$ l) ZR=r;
- d5 w; g8 S, s, i1 ^: WRV=rv;
2 y) y& d& ?+ ~' N' f  _for k=1:N
6 w! ~) y; ?) q/ _    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
0 z5 d* t( M! Q  w0 u    z(1,k)=k;
    z(2,k)=A(k);
end

& |7 d; G/ l+ q! e* [7 a; t/ p不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
& F7 j8 b: P; b( wfunction z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
  }2 r# R& `3 N. d! F  JN=n;
R=r;
9 ~9 _- {6 s- z5 KRV=rv;
5 q/ i+ Q$ h3 H  l+ a9 f% W: Vfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
- K7 Y" O2 \0 ~: o$ H9 e# k) Z    z(2,k)=A(k);
7 m  ~' S* `/ t( M5 Pend

dreamofyou

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

厚积薄发

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

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
所以应该修 ...

谢谢谢  我终于明白了