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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
: x/ P8 o( k* y& }function z=nd(a,v,n,r,rv)
& F1 K7 a! ~% O' j, KA=a;
/ U! I. Q3 ^4 N3 Z/ U; @" pV=v;
N=n;
1 {+ [$ k( B/ k; E# O; l8 pR=r;
RV=rv;
: J5 ^+ U/ J9 r; T1 l3 Y  N1 t  Afor k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
9 y& q0 A# D+ Z, M    z(1,k)=k;
    z(2,k)=A(k);
end;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
, `: m1 n' [& b8 j$ T??? Attempted to access A(2); index out of bounds because numel(A)=1.
Error in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
. E$ h& x1 w& p" A8 T1 I8 M
1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
- c: q5 [0 b" n# h& {3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
5 K6 v0 \2 o- T& T, _# D( SA=a;
V=v;
# X3 P; }+ S3 ~! M5 \8 qN=n;
1 p. U/ ^4 g2 ?" ?6 b3 {/ IR=r;
RV=rv;
for k=1:N
1 y# b$ V+ ~* h+ j  s& X. @    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
1 {+ Q' l3 Y9 W% N) E    z(1,k)=k;
    z(2,k)=A(k);
$ {  u1 r) c3 Y$ G; Pend
2 u  B: e2 E* t* I9 V
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
, W- S( Z6 q6 i9 q6 aV=v;
N=n;
/ Y+ k2 Y" t  YR=r;
6 ]. L  T) O; ]2 B/ z& }+ p# URV=rv;
, u9 T; q& }( i( B+ @# Wfor k=1:N
1 u- a+ U7 R6 ?    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
0 J+ J+ O6 S3 O; tend

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
1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
所以应该修 ...

( ]4 {. Y$ B, b: B# j谢谢谢  我终于明白了