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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
3 T0 ^' `" O  nR=r;
RV=rv;
for k=1:N
* e+ T  K3 n) r; P9 J    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
9 S# H( j1 x1 K  [. C' F    z(1,k)=k;
8 m+ g, z; T% d* |    z(2,k)=A(k);
end;
% |) N  i* d( M6 ?4 M然后z=nd(0.2,200*10^4,360,0.02,10^4);
; A. r8 p  Q8 a" @8 o系统显示：
8 E4 Y7 t) |6 Z% t, ~  H7 V??? Attempted to access A(2); index out of bounds because numel(A)=1.
Error in ==> nongdu at 8
( ]2 M; R3 Q2 ?8 R. |    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

8 R# z4 @, O5 c! J. f& o' i1.请问为什么啊？
  O. N* j7 |! ^5 k8 c2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
* k+ a4 L8 V  N) z3能不能有些教程或者实际的程序啊，当例子。

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
9 s9 K1 r1 U4 Q0 u. g$ @; ^* R1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
5 j7 M" A/ w6 P3 c所以应该修 ...

谢谢谢  我终于明白了

厚积薄发

1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
  O" L5 k4 q. P7 |0 i! G4 g" _所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
* H* S2 B- S, P8 _3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
2 J7 A/ E& ~; I, U6 UA(1)=a;
V=v;
8 c3 i/ V# B6 g5 i/ e7 p0 X5 FN=n;
( X( m5 ?0 v1 rR=r;
( H' w9 E1 C4 Q) |5 @; i  {RV=rv;
: q0 A/ u4 [8 @( |& X1 pfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
9 Y, V' g! C9 q$ h/ A2 u    z(2,k)=A(k);
end

leo12ok

function z=nd(a,v,n,r,rv)
; t5 T8 g( D( w  b3 eA=a;
- \. t) W. O% d' i# L, d- sV=v;
  G" x* i. o0 g& uN=n;
R=r;
, a$ `' i' |. w/ q$ P6 L& dRV=rv;
for k=1:N
& H) |. B- \. b: a" m# d/ f    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
end

: C3 p! R3 t% V! ?4 a不知道这样可不可以。