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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
1 }! x7 z! p) |- U& a+ v( v- Z# C. xfunction z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
" @8 T/ v8 r1 A2 U/ X( s. ]R=r;
3 ]; F7 i( ~. ]8 ?RV=rv;
( e4 m7 q  p5 b5 n. }for k=1:N
8 s, T- M6 k1 C' u0 X    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
# F" o( v, n; |: m. W    z(1,k)=k;
: ~3 E2 [: ]2 L    z(2,k)=A(k);
" w8 n6 Q; a3 Cend;
$ W$ O! j  k2 W* s$ r7 T然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
0 S0 K# g% Y9 X??? Attempted to access A(2); index out of bounds because numel(A)=1.
+ R1 I' [/ ^0 O- v7 oError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
+ E, L  o# f. T( n- L' R
  ]) A: @. @: |# F% ~* L: X1.请问为什么啊？
! B# K) x7 h+ R9 r4 u, l( S# f2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
8 l( A8 z, `$ N& NR=r;
RV=rv;
for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
: \0 o+ K9 G- V6 _    z(1,k)=k;
    z(2,k)=A(k);
end
, T/ L" g+ ~1 m0 b" J( s2 l
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
) Q( W0 s1 ~* }! J9 aN=n;
R=r;
' G2 {5 F1 b: o9 b6 iRV=rv;
6 Y5 d9 W/ _6 \- g& O( l( l7 ifor k=1:N
# @2 ?: P' G8 I% m) e4 |* }. W    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
1 U4 T+ F' f+ a0 Q    z(1,k)=k;
4 ^& U6 s) n! H    z(2,k)=A(k);
% v) _4 G7 x! i2 @* y* G0 Y2 G0 [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）
: B! x% @* ]4 I3 i9 G9 w2 C0 r2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
4 [# l6 p  k0 D3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

谢谢谢  我终于明白了