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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
A=a;
V=v;
! h  U' K2 D/ I. U5 r3 M, xN=n;
. e& z& J% ]7 Q  h3 V7 kR=r;
RV=rv;
for k=1:N
, i! q/ \9 Y/ b9 i- C+ N    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    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.
. y& @, p  Y/ S3 OError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
/ m  q6 j0 y0 l" r, W- ?7 b6 O
! u, z5 n/ \2 m/ H- z7 ^" h( a1.请问为什么啊？
5 Y; ?, C5 ]# V' s4 @( g% W2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
! z3 u6 Z. @4 U  M$ {0 d! E看书上直接可以写成A（k），到底应该如何呢？
& {6 ~' }! H+ T- O9 B' Q& Y  P3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
: |7 a9 w/ Q. V2 ]A=a;
. H( \$ T. z! F) c7 O& D, W7 UV=v;
N=n;
8 r+ J2 {7 P# uR=r;
RV=rv;
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

; @4 R. b1 V- J/ [6 ?+ D6 u1 X* X不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
: x% z4 ]+ u& W5 bfunction z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
9 D1 @" ?  b! d" a1 X# W0 D/ Q/ x8 jN=n;
: s% F3 ~2 k3 H$ X5 BR=r;
RV=rv;
/ c2 u5 y2 z+ [for k=1:N
! L/ Z$ O6 \  s9 C- K" B    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）
  R1 O' I0 K; D! x3 D& j2.不论是写成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)，这属于非法赋值
所以应该修 ...

谢谢谢  我终于明白了