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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
$ o" S$ x& I& G1 z0 N$ Ffunction z=nd(a,v,n,r,rv)
! [" p5 K& h+ d/ m5 o6 T; fA=a;
V=v;
N=n;
R=r;
RV=rv;
7 e4 J+ d3 p; f& Kfor k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
+ v( ]. u3 e1 D6 u% E4 `/ E    z(1,k)=k;
    z(2,k)=A(k);
5 V% x- c! V. k) f% Y$ \end;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
. h% b* N( ~6 K系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
4 j7 |3 r' V9 P  D/ EError in ==> nongdu at 8
" p% P1 u0 A" {# B+ t. @  ~+ P5 S    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

+ C1 b1 @, a; a8 Q6 t6 c4 B: v5 x1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
9 A( ~9 P  S& s6 l0 D看书上直接可以写成A（k），到底应该如何呢？
9 j" D6 C0 ]: q4 s3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
% s, L6 H' ^: U" o7 RA=a;
' |& q8 W* d  Y( p" S) x6 aV=v;
3 a* N2 c, q* Y' K2 _, l. }N=n;
R=r;
3 ^" v0 ?" ?1 e9 PRV=rv;
( V( [8 s5 j+ Pfor 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

9 @) L; a" ?, R7 M不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
8 t' q( C- E. {( X4 H( UA(1)=a;
# s+ \8 K7 {3 i4 w( yV=v;
N=n;
R=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);
" _; f+ k- S& p& n3 e9 E5 Fend

dreamofyou

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

厚积薄发

1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
7 P! n% M' M  G5 q+ l' M所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
1 K8 Y$ v( }' L3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

谢谢谢  我终于明白了