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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
+ J/ d  y# a3 k/ a- Z/ ?+ dfunction z=nd(a,v,n,r,rv)
A=a;
# O, O, ?+ |% v, AV=v;
8 @4 }% x6 W2 ~* GN=n;
; b2 t' M4 z6 N8 [" YR=r;
RV=rv;
for k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
! `5 Y8 A6 m2 h; m* Y6 q( v4 H" |    z(2,k)=A(k);
end;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
2 ~5 k  A' E! O& d* r??? Attempted to access A(2); index out of bounds because numel(A)=1.
$ S% i. s' O7 C2 @* o9 n( GError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

, d% w; G' a1 h1.请问为什么啊？
0 a7 l! s$ f. }. Z9 O2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
4 w" F5 y7 V1 A3 j3 H看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
6 o( M9 r, [$ t2 r. Y: gA=a;
" x. H+ k6 m! [* K+ R/ p3 gV=v;
7 E, t0 a5 F* B* g8 `: xN=n;
" ~$ e6 _; E! c& ~6 X% eR=r;
RV=rv;
for k=1:N
: i4 g" t& ^; ^  E    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
2 f. E: w( A, U2 u3 V0 ~1 `$ o    z(2,k)=A(k);
, \: x" N+ H; h( Mend
# O2 e& D5 a5 M/ q0 c$ o# K
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
& s0 ?! q( u' D! e- x- mN=n;
R=r;
" x( v7 Q3 _/ ~- o3 q0 }; P; PRV=rv;
* x4 s" {# Q% f6 k! pfor k=1:N
+ z$ G4 s( t/ L    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
& e7 A4 c$ R8 c6 b9 S    z(2,k)=A(k);
! ^6 p" G" c# V/ i& @end

dreamofyou

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

厚积薄发

1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
2 @2 L7 Z9 g! K( P所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
: q' \. B% m; f( j0 a5 V) Y2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
6 q+ e, g7 K. H! C3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

0 E8 f2 R. O+ D谢谢谢  我终于明白了