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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
* N, ]! e" I. M2 b% aA=a;
V=v;
$ ?6 ?' S/ E, u) j2 ~6 m" jN=n;
- l7 o, h. f3 J% _6 K5 @. i6 e$ dR=r;
" ]& M8 n3 F0 jRV=rv;
2 a+ J* `5 M8 ?( ?for k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
  h  i0 O5 u1 L. p8 a- I' V    z(1,k)=k;
# ]% _9 l( D1 B4 X/ h    z(2,k)=A(k);
/ F0 m* N3 [9 l6 V0 Yend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
# x: H- `8 Y& T! Q- U系统显示：
, z, V' B3 u7 w4 F6 O??? 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);
* i( }  \1 ]) \$ i
/ X: z/ e5 ?1 B1.请问为什么啊？
4 E" S) `0 y  g$ z) I) L2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
V=v;
! ?8 n& m6 q. q; E# @. U' uN=n;
4 a0 R5 N9 k' `  y8 v: ~& _R=r;
RV=rv;
for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
5 H. G: m- }3 Z+ X; y4 [% g# I1 e    z(1,k)=k;
    z(2,k)=A(k);
end

! u2 s6 H. f- P不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
) z! L( a. _/ _- w/ k+ b" GV=v;
( z4 }2 n% \* `8 z$ Y& ~5 BN=n;
R=r;
& D# G" L* j3 c; o  O0 ^RV=rv;
for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
1 v+ A) F% o" L+ \# C    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)，这属于非法赋值
! z3 x! \  r! _+ s: e: B# n4 x所以应该修改成：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
& e- |4 g' I7 O8 E- }* q1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
所以应该修 ...

/ I$ Y( l7 u1 Q/ a. d7 J# o/ E/ M谢谢谢  我终于明白了