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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
- o: ]- k$ `3 r, Gfunction z=nd(a,v,n,r,rv)
- v9 c# t2 ~+ Y/ ]6 eA=a;
V=v;
N=n;
8 p' O" N' n# [% x6 Q3 P' K* iR=r;
RV=rv;
for 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);
end;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
" ?7 W+ q: ?  \; v+ q, u0 \??? Attempted to access A(2); index out of bounds because numel(A)=1.
( X+ L1 ^% @$ D! z& s& `5 g& cError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
# w8 S  V# Q$ Q7 Q5 s看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
7 n2 h6 @) D1 i3 W1 p9 VA=a;
6 {# W' v$ r! QV=v;
N=n;
3 n1 |. L! K9 p" }) PR=r;
RV=rv;
' z  H$ A! G% q& xfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
  g7 s8 ^0 v8 _2 l' m    z(1,k)=k;
    z(2,k)=A(k);
2 C+ U8 e7 g0 n6 E, Gend
8 V1 @* Z7 |; [+ o8 j4 ]  i+ T
8 q$ F6 _, P6 I& W- i9 `% b不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
, [9 ~# K4 A/ iN=n;
R=r;
7 }- y0 c+ C" f1 t8 d) Q0 \, e6 eRV=rv;
6 q/ w! E0 y( ~: Z2 F! U' _; pfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
+ F7 M! i  @) ^1 S; ]! [3 _2 A    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）
' a9 T+ P' ]5 R# m1 P4 \2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
+ O& H2 b+ h: n  U- R" o1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
所以应该修 ...

2 R% i0 N# H+ W: F+ O1 N谢谢谢  我终于明白了