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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
  a: x5 S! M. D6 f9 x5 t, ufunction z=nd(a,v,n,r,rv)
A=a;
9 o6 G. F' g3 ]8 l) TV=v;
N=n;
0 I" B& ^! J& D$ }$ [- d6 XR=r;
RV=rv;
" @9 h: B! I2 X+ K+ bfor k=1:N
( Z" z2 g; |5 e6 T' w- j$ h# K    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
/ a. f; `& J9 c% n" e    z(2,k)=A(k);
end;
- _9 I+ V7 C# c3 t1 @4 |+ a* n然后z=nd(0.2,200*10^4,360,0.02,10^4);
1 l3 S+ Q2 i  M) H0 [& Z系统显示：
: ?# n5 @- N" H. v- [& O, [??? Attempted to access A(2); index out of bounds because numel(A)=1.
2 n! J( F& K0 m3 q. M% l. a' I7 CError in ==> nongdu at 8
8 q2 p1 x/ @; r3 W8 \$ _    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
, i! j1 j) `2 d! x4 M9 G( N
; s6 j' N3 G% I* G( K! N1.请问为什么啊？
; T  f3 A& _) q2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
( N. Z6 U: k* k3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
! H' |2 ^" Z# D: ~2 I' F+ bA=a;
% g# a  f8 r- ~' v: S3 K3 C  {V=v;
  ^. e# g$ M0 g9 J4 I& O1 T, vN=n;
R=r;
" H+ q) z  P" kRV=rv;
4 }2 y+ A- E9 k3 l  Afor k=1:N
4 F2 Q1 o4 \- O# b5 H* d* a1 k    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
2 p; S! r( \8 ~    z(1,k)=k;
" Y" t+ h4 N" |' z. x1 V1 U. C    z(2,k)=A(k);
end

4 T, G/ J$ x  r; j- p3 Q不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
+ E/ H- N  `2 o) T+ efunction z=nd(a,v,n,r,rv)
6 F2 a0 B( h, S, ]' Z0 d* mA(1)=a;
9 p7 l& v# o. C+ W( e: a' sV=v;
2 T6 d' w: @  d, f' O. _( sN=n;
# ~+ r2 v2 W4 w, ]% _2 OR=r;
RV=rv;
8 [8 O: N* A* t4 v, N- h( Vfor k=1:N
9 }; a$ D7 o3 T) K: ?    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
. o$ y) H, _, G+ I$ M    z(1,k)=k;
1 [; _4 u% ]/ E: {9 g+ z7 b    z(2,k)=A(k);
end

dreamofyou

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

厚积薄发

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

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
5 {+ z! I. D5 X- i7 j1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
5 K0 W  \( a1 K1 T0 d5 Y所以应该修 ...

  _$ @- F8 J$ }! ]4 V谢谢谢  我终于明白了