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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
4 E( N. o9 D7 zfunction z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
R=r;
; e. B5 X/ R, }9 FRV=rv;
for k=1:N
3 D2 r0 Y: R2 G; d7 {- @    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
2 k  A& P" B: q& J& l! \% S    z(1,k)=k;
8 j! ^' x& L% ]- p9 j    z(2,k)=A(k);
end;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
2 l  |- q5 k0 T( L! W系统显示：
* I* {* L' \2 ~: q  j" y! V2 J??? Attempted to access A(2); index out of bounds because numel(A)=1.
( H# H# ?. U; }1 `; t% u9 pError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

+ A, j* ~1 Z0 B/ x& j- J0 q1.请问为什么啊？
( |' P, @& W! q* N2 y4 V2 l2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
7 b# V7 ~9 R+ L& x$ \3能不能有些教程或者实际的程序啊，当例子。

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
- j; U+ n5 K/ a3 A7 s1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
" L' c. E9 ~4 I所以应该修 ...

谢谢谢  我终于明白了

厚积薄发

1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
+ m" M9 R6 {3 r1 z8 a所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
. N! x* h8 F) e* }/ @; `2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
; l, D, X$ w5 q6 H+ k0 a0 V# w) H3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
4 U) u; k; C* _N=n;
R=r;
: l7 j7 T- V, \8 R- U! P1 cRV=rv;
0 K' Q$ R, b: C$ N9 k" |for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
" R, Q8 |- U& k: D1 |: @: w    z(1,k)=k;
# `( r& x" G2 i' k; x* }    z(2,k)=A(k);
$ b- G- _5 R/ L5 b( w5 B" Pend

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
; H  h$ b" F! H; N; AV=v;
N=n;
1 a( c+ o( F: B( ]# NR=r;
RV=rv;
. ]2 I8 t+ P5 a. N+ qfor k=1:N
0 R  H, \8 Q3 k8 Y" ^    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
0 g' y1 G) `2 c6 A' p  z$ u5 L    z(1,k)=k;
    z(2,k)=A(k);
end

+ B! j5 V5 Y# v4 g: q不知道这样可不可以。