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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
% q( i7 V1 D$ D/ FR=r;
/ Y: o6 g5 j  U6 y) T2 ]RV=rv;
" i* ~1 A0 \. I4 v' Qfor k=1:N
% Y- X- D* k# X( l, @: N' D    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
% G/ m7 P2 y; S+ k. u8 a    z(1,k)=k;
    z(2,k)=A(k);
3 H( B' |# l) ^/ f# g& G# uend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
4 d6 Y6 M0 G% _* A1 y( T# zError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

* ~+ v: _* X# C: G  x% H1.请问为什么啊？
2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
  y' ~; r+ n* q, W( N) C9 ]0 [看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
R=r;
RV=rv;
3 w; ]" P, M, p6 @5 D0 kfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
" R: E* O; ]4 M8 h! n6 J    z(1,k)=k;
    z(2,k)=A(k);
# r" Z$ r1 f" G2 w7 ~* }2 C( m5 [end
) x5 T" q% X1 l+ g
% n# n3 @, F3 u# }) _% B; N不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
$ c) c) H7 a1 y3 \) dA(1)=a;
  _* G0 Q  g' Q; fV=v;
N=n;
R=r;
. s) o. A* ]! ?- NRV=rv;
for k=1:N
0 q3 r: @) O; T$ B" M    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
+ @; n5 a3 o+ X    z(1,k)=k;
    z(2,k)=A(k);
2 ?: ]' S& P( i3 J3 j. Z  w9 Eend

dreamofyou

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

厚积薄发

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

dreamofyou

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

谢谢谢  我终于明白了