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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
+ Y5 G. I# o( R5 x  Kfunction z=nd(a,v,n,r,rv)
- J" F  e/ e0 a$ U1 D' }A=a;
7 Z! k# [. p5 B+ X# K% oV=v;
$ M8 K) w4 m4 m0 c+ d0 q9 E# ]; ]N=n;
) s  l, O* B+ F9 y7 ~4 \; G) ~R=r;
RV=rv;
6 [: ~" A7 l+ U; C1 ffor k=1:N
0 @, F8 o( v, {5 t% H* R    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
$ m5 F' j  p: O" R+ v( A    z(2,k)=A(k);
5 q& N9 b2 H$ |. Rend;
7 z8 L+ {; C& C& Y+ o" E+ E然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
) N) ~9 W3 `8 V% g+ u0 O) W1 z  B+ g??? Attempted to access A(2); index out of bounds because numel(A)=1.
- X# e% I8 U  U5 xError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
4 s) n8 r7 {% m+ O* A$ ]3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
  Y$ b0 R; w4 D5 I8 N3 SA=a;
V=v;
3 J) u& c. F" A& tN=n;
R=r;
. Z3 q. I! f' y' [RV=rv;
3 j& b/ ^; _1 s7 }8 @/ C8 i4 O9 Cfor k=1:N
- ?2 l; t4 n  R    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
2 s* W! |' ]! W' a8 i0 a    z(1,k)=k;
    z(2,k)=A(k);
end
1 D' l+ E" n* u% B4 }
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
" e4 z& x. v% @4 r; B; w" yN=n;
R=r;
# r* _" G: B( y& P7 ?8 rRV=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);
( t6 d/ A* Q( Nend

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）
9 Q+ n$ k9 A$ i1 a( M2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
+ p  P5 A0 W2 @1 X: A6 |3 {# @  u9 D1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
9 e& m, f" P  N1 G5 G# j所以应该修 ...

5 O0 J( E+ X" d, i1 o% ^谢谢谢  我终于明白了