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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
" a% K: W& o7 U0 q1 pfunction z=nd(a,v,n,r,rv)
2 z0 M; I2 j. {2 I& {. lA=a;
: x' k* j8 d9 T+ X4 oV=v;
% c$ Q/ J9 y" M. F; O3 b6 jN=n;
R=r;
RV=rv;
( O/ I( R* x/ U1 T2 f& pfor k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
1 \7 L; p1 ]; ^6 p& Vend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
0 v8 q3 R8 C! J8 _, }系统显示：
??? 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);
) c% W- B: E/ d  o) U4 U4 ?
1.请问为什么啊？
' r1 R& \  @, [9 a7 b1 |$ ~2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
5 |1 T- B1 w% {& J. @3 rA=a;
V=v;
. g4 k" s% S/ e6 B3 eN=n;
R=r;
" K$ e% I" \4 C( E& LRV=rv;
for k=1:N
: {/ j6 I& g. C! A$ ~, _7 f7 Y    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
end
4 o% P( r) T9 j
; {" @6 K) e- v- k, n+ A# U- Q不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
: y9 k9 v; z, Q4 b, oV=v;
N=n;
9 o" g3 e6 O4 iR=r;
RV=rv;
for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
) L7 x- A: m8 n% y* ?, D! {3 T    z(1,k)=k;
    z(2,k)=A(k);
6 C4 S5 Z6 S$ Kend

dreamofyou

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

厚积薄发

1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
8 Z. h. G3 I& H3 m* e/ N+ X( R$ L所以应该修改成：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)，这属于非法赋值
所以应该修 ...

  i2 j9 t2 h8 g4 E0 A8 r谢谢谢  我终于明白了