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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
A=a;
" ?7 T! c8 T% N9 @( ]V=v;
N=n;
R=r;
0 d6 O" Z& a8 D& \RV=rv;
for k=1:N
4 [& F+ G+ s& b4 Z    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
end;
* K: y' s3 B! X; }( ~( @7 @2 N1 g然后z=nd(0.2,200*10^4,360,0.02,10^4);
* b9 F# L( D7 g' _系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
3 q' k; f6 T7 Q+ \; Q' E  ^Error in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

0 q( W8 M( u9 O( \. }) J. M' R# d1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
# k- z* t. f+ k4 S2 n8 y3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
- ^8 A' G9 s+ C7 ?0 UA=a;
) V0 h& X6 R9 e! F0 QV=v;
N=n;
" j# {+ D: N; C  C  y7 R  \R=r;
) c' ~& R5 L4 D0 }) z7 w4 u, wRV=rv;
for k=1:N
1 R5 t5 X* G9 F6 r* V4 M# s2 {    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
end

; @# B/ g8 e- ]! ?! |9 @  K不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
# W$ H  j) {: I8 j5 Qfunction z=nd(a,v,n,r,rv)
" C- _: _# l0 H7 lA(1)=a;
V=v;
N=n;
R=r;
2 b. P8 Y7 f0 A2 ^8 xRV=rv;
for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
1 M. o+ o% @' [% e    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）
, k8 p4 i# n6 v/ Y4 v4 S3 e+ {2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
  n6 G' o5 @$ Z8 y1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
所以应该修 ...

' H& `0 B" A; B: n2 D1 \谢谢谢  我终于明白了