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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
2 ^7 R* b1 F4 _( Hfunction z=nd(a,v,n,r,rv)
A=a;
' t7 U# ^% F  f/ O! b3 A* O! nV=v;
" V8 n7 A+ v- r. S: IN=n;
8 Q$ p0 n" [6 C0 D/ NR=r;
! V" c/ J- g7 Q9 U( BRV=rv;
! ]# R( s' I+ ifor k=1:N
8 ~) r8 g: _7 }0 ?* a5 x8 ]" r( {8 z    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
. A0 c4 w9 {- r3 ]5 T9 ~' r    z(1,k)=k;
9 o0 P4 K. b. M3 }    z(2,k)=A(k);
end;
: u, k! _& Q  p然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
( I) p& m. n# i5 ^' J??? Attempted to access A(2); index out of bounds because numel(A)=1.
, g$ X! @) U4 `1 i5 S3 H' |Error in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

1.请问为什么啊？
7 v6 i" ?( n6 t* O! R8 G2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
4 O2 _6 n' n+ b3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
5 {/ q+ t7 y* T9 XV=v;
N=n;
8 J# Y  G& e' A* O$ d1 sR=r;
RV=rv;
7 t+ v6 y+ p# G" @; Rfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
8 I( O8 |) ?% u0 w" P. ~    z(1,k)=k;
    z(2,k)=A(k);
end
0 g  r/ @- H' N( P+ u# O
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
# |& P- ]8 L9 F! nfunction z=nd(a,v,n,r,rv)
; N1 w4 F. c% D- D6 i! nA(1)=a;
V=v;
N=n;
: ^) j7 D! n7 SR=r;
RV=rv;
+ i/ F6 E+ M: B( h" }2 Afor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
' E8 n+ e# D! d) R9 ?( ~. K    z(1,k)=k;
    z(2,k)=A(k);
0 n$ d( V; P+ u5 D' \8 {end

dreamofyou

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

厚积薄发

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

dreamofyou

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

" M$ _2 H1 x1 R6 m谢谢谢  我终于明白了