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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
' Z) h9 {2 n1 IA=a;
7 v+ B* k9 G& ]V=v;
N=n;
R=r;
4 o. U& k, P2 m: vRV=rv;
for k=1:N
, V  u6 x2 ?! P    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
7 P5 n7 \$ n1 I    z(2,k)=A(k);
end;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
??? 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);
5 a* s- S$ ]: W: S. v& X3 [1 p- Y
1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
, ]# t- E& T3 |$ s$ ^* j! n看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
$ b* B" T! S- e. w( X1 e5 HA=a;
$ W2 l% b) Y5 o% MV=v;
N=n;
  Y8 r/ D4 k1 C2 A! p" rR=r;
RV=rv;
for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
5 q* k/ k( K) k2 H8 c+ ^/ m    z(1,k)=k;
    z(2,k)=A(k);
8 @$ `$ C: p/ C; `end
& ^; Y3 t( b5 {% V. @$ ?
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
8 e# R8 t0 Z0 W+ m" g" e$ p! @4 Nfunction z=nd(a,v,n,r,rv)
: I9 |# g" ?" |' t+ @: S7 hA(1)=a;
9 x- _- P' C/ M& U. rV=v;
8 C( f* M) k3 }+ |* i- qN=n;
; w# M1 X5 b+ Q2 [' z8 X, cR=r;
RV=rv;
+ W0 J) t# y4 t" p: B2 [for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    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）
2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
0 T( l- q0 k7 C7 z) _2 P1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
所以应该修 ...

谢谢谢  我终于明白了