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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
" V" k, c0 V. F) z8 Yfunction z=nd(a,v,n,r,rv)
A=a;
V=v;
( q" ?% z7 ]( L, e, `& W( q5 kN=n;
R=r;
1 ^: B" K" @! m9 ~4 }- s# E% HRV=rv;
for k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
4 J, J) r8 ^9 q1 v! r    z(1,k)=k;
    z(2,k)=A(k);
$ K6 J8 G7 w0 g6 n% iend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
; Z7 C0 z. P& Z- U3 ?/ r! ^??? Attempted to access A(2); index out of bounds because numel(A)=1.
/ P/ b1 W* f! z" j! a( Z3 y, WError in ==> nongdu at 8
/ v% m, s* \" V. Z    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

3 P% Q& w& _& K: i1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
4 e0 \4 V% b+ M- w- l3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
- X! S; u/ A% lV=v;
# a" Q* d7 f, CN=n;
, U6 V, W/ _) o" C2 i9 b4 R+ eR=r;
: {$ [+ A8 }, k: V' Y1 _RV=rv;
for k=1:N
! D. N1 @/ U. k4 q# d0 q    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
+ `3 \' [" U: x% u    z(1,k)=k;
; m7 Z3 L" F2 y/ M. ^    z(2,k)=A(k);
end
7 A' S: m3 |0 p4 J1 R; q
5 Q/ V8 j* \/ ~" }+ W不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
: S+ Y9 L: c0 Y6 |+ k  Afunction z=nd(a,v,n,r,rv)
2 `* }& ]% {2 u6 l" UA(1)=a;
5 O+ H; f, x; n# v, M2 RV=v;
N=n;
, \- x+ }* W0 R2 [% ?R=r;
: l/ }" |8 p2 bRV=rv;
& {' Z# Y& O# ?1 F( R" nfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
0 B( X  k0 F' h  y; `% R    z(1,k)=k;
  \" O- Z9 a% H* r    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）
' T9 ^- T' ]/ Q- a2.不论是写成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)，这属于非法赋值
+ I. m7 l; h5 `& o& W0 B所以应该修 ...

5 {9 z5 G$ Q% C2 P谢谢谢  我终于明白了