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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
# o% D/ A( u* p7 Z1 C  Efunction z=nd(a,v,n,r,rv)
( }# T* T* `$ XA=a;
, S7 n& l) E( h, L+ i+ I7 sV=v;
N=n;
7 k" m' Y" Z0 QR=r;
RV=rv;
for k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
end;
7 Z1 ~9 s) p& l+ G7 Q然后z=nd(0.2,200*10^4,360,0.02,10^4);
3 u$ U) f. U2 ]3 x系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
2 i0 D# k3 p4 w1 G" G) ^% D& `, FError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

1.请问为什么啊？
. b+ v" a( A5 c: q2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
! C  g' ~  P" K2 ]  r3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
R=r;
% H; q; z5 v. Y# HRV=rv;
9 O9 ^, V2 t  a! M4 e5 J7 ^) wfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
  ]$ h+ w; a0 O    z(1,k)=k;
    z(2,k)=A(k);
# l, P$ e# G8 A4 V0 R* @$ G! Zend
4 }" B+ Y: i* J6 K
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
6 P6 j6 b( j; A, P* ?3 Nfunction z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
# `  V' j# R8 B+ FN=n;
R=r;
9 \1 I0 S/ A' {RV=rv;
for k=1:N
6 d: K' J. m$ F0 y, `  P    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
) f! h8 b7 D& [" T' C3 N    z(1,k)=k;
, x9 Z. G- c6 n/ i, v1 A9 f" G# i    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)，这属于非法赋值
所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
7 s3 y# W. P$ ]+ F; v. J; X2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
7 |! _! D0 f( R" M4 [5 |3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

' c: q+ M4 d1 y$ f谢谢谢  我终于明白了