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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
A=a;
" P8 D( p. e. B1 |V=v;
1 i' |1 f4 r6 i+ F) }6 [" xN=n;
R=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;
1 w4 x; T: P3 N' b2 Y然后z=nd(0.2,200*10^4,360,0.02,10^4);
& [3 x' A0 C# F系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
Error in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
/ n9 ?) N, r- h2 a
1.请问为什么啊？
6 [4 ~- T/ u1 U$ p, D) A+ f2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
3 y% H( V1 V: s8 W# ]/ E; p- ~看书上直接可以写成A（k），到底应该如何呢？
' Q( b& A' a, X2 j- d3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
5 N3 C7 X3 r1 F0 j( F" A- k' K' aA=a;
V=v;
' r# p, f( J6 D3 g% Y+ O5 TN=n;
R=r;
RV=rv;
6 u4 Y/ m0 [7 R" g5 |' Jfor k=1:N
# j, A0 g' c# p2 H& i3 Q    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
& W) g9 Q. O/ N( [* H  n2 h/ w" z: K    z(1,k)=k;
0 `1 V+ I# P; ~6 i    z(2,k)=A(k);
end
9 l0 k* z. v: G  m% b( V6 j4 {
! A) ^; d/ G  V7 e3 o' e不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
4 k/ P6 C! t4 G5 }+ t) O7 }V=v;
N=n;
1 @( P5 L$ o+ t# r4 f1 MR=r;
4 m% ~5 `7 M1 M0 C9 sRV=rv;
# \/ _% ^& ~1 p: O8 \; G+ F( {for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
" w8 s; m' e: c$ N+ k    z(1,k)=k;
    z(2,k)=A(k);
5 q' b& }! ~5 _- m( X3 X: jend

dreamofyou

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

厚积薄发

1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
1 Y4 |0 l7 p6 r所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
- N4 _# O. H8 q. ?5 i3 Q* Q, K; l3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

) C5 e$ W. {7 k! |谢谢谢  我终于明白了