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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
R=r;
RV=rv;
) n" o7 I! z# G; M9 E+ X3 ]4 {for k=1:N
' z- D) q! R" X    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
& H' ~" G1 E+ S    z(2,k)=A(k);
' L6 Z. a+ ~4 U5 n1 P- Xend;
: F. \. e3 w( Z然后z=nd(0.2,200*10^4,360,0.02,10^4);
" V$ A5 S2 ^( ^- w系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
7 {: \) a$ H7 X, G& [7 f9 qError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
" \1 O7 a/ d: b+ ^
0 F6 L6 K2 s* c1.请问为什么啊？
- o0 @; x* l9 \2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
/ M. F  s! d* T6 C看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
+ {. M. M4 M' l3 sV=v;
N=n;
R=r;
4 R  y' T' }4 }/ A- iRV=rv;
+ [! p9 a& }% Z8 O: f" U2 q% t7 ifor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
: R9 J1 v9 C. D, Q    z(2,k)=A(k);
! L8 I  S" l5 X9 d9 s6 jend
/ _3 e+ [  B9 _
! C, g: n, |5 g1 T( r2 k* ~不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
N=n;
4 ^0 S5 b0 ?6 S: ER=r;
7 E4 K- ^- H" dRV=rv;
  T) w) j/ }- G  W3 o, b! wfor k=1:N
, n$ ]# z; a4 l& s* g- b5 R& \    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
8 y0 c% X/ _3 S" v( c    z(1,k)=k;
    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)，这属于非法赋值
7 E/ |0 O& A/ V5 U# d' g9 x0 N# |所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
; K# \- [% c" Y6 H( z2 O3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

厚积薄发 发表于 2009-12-25 12:08
  H7 \5 j/ ]0 E7 w0 H1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
2 F# s! g$ \% u% Q% U; a, J9 Y. K所以应该修 ...

1 U/ Z5 X+ n. o' |, p谢谢谢  我终于明白了