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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
* I; F8 b  F! ^function z=nd(a,v,n,r,rv)
" W/ b: |0 K; }/ R; BA=a;
. M! @" o) t/ f0 m- z- ZV=v;
& b$ \5 C8 k& e* \N=n;
R=r;
- B1 v  F( h8 z" }% Z; IRV=rv;
for k=1:N
, a! O! r2 K( P1 c3 k    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
3 c" e& d3 Q' s% w5 {( o    z(1,k)=k;
    z(2,k)=A(k);
$ l7 L6 l2 [$ Y% Xend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
) o: P% b9 o7 G# D# p系统显示：
??? Attempted to access A(2); index out of bounds because numel(A)=1.
Error in ==> nongdu at 8
2 j" B/ X1 \4 [" [; i' M1 Z0 y    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
8 W0 P) X9 ?1 C8 ]% U- v
) x" l- G) H( c9 X% q1.请问为什么啊？
: N  v5 R/ m& N: {2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
# ^# p& R4 m$ h. U看书上直接可以写成A（k），到底应该如何呢？
( X5 P1 D& M# D3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
) y9 L' T0 k( b& L' y3 o8 WA=a;
! ~4 F' z2 g5 FV=v;
+ a3 B6 M% a/ ^N=n;
R=r;
6 a0 o2 q; |, Q6 o$ R" g/ t7 bRV=rv;
; n$ _, ?! u$ [* f5 y4 tfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
8 ]6 J# l7 L% U! K1 Aend
- y5 I, b# g2 M+ \/ b2 R
/ P% D5 }4 P' {7 A  S4 M不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
% L4 N3 _" B; P6 s2 i7 nfunction z=nd(a,v,n,r,rv)
; b* D) y$ ^' B3 t* b. T( lA(1)=a;
V=v;
N=n;
% h  U& R4 v! j  XR=r;
RV=rv;
( V: n: \+ g# Nfor k=1:N
$ E* _0 a7 E4 S; i* z$ |    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
+ J" x" q5 M! B& L    z(2,k)=A(k);
& {; x% A2 L' i4 T6 ^6 eend

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）
2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
! i% W5 P# j# R, F' F3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

, h4 p1 ^, ~+ I% q; |, F谢谢谢  我终于明白了