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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
" B1 x8 N" u+ Q5 `/ nA=a;
6 p8 U  S+ r8 ]2 FV=v;
N=n;
R=r;
RV=rv;
* b/ L/ U! t7 t5 Bfor k=1:N
; \$ d0 l* ]$ o    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
/ K' v9 z" a- f, \( Z    z(1,k)=k;
$ x; M; Q; n5 P$ M, ]4 Y    z(2,k)=A(k);
* |2 N; C$ a# K# A  R2 pend;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
; }1 W% b" S' b??? Attempted to access A(2); index out of bounds because numel(A)=1.
: ~$ |# B+ w! f% O4 F1 ^+ K. Q0 vError in ==> nongdu at 8
7 G) w  Q% r  x    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);

1.请问为什么啊？
4 l6 k% e" |* O3 \; y8 \* E; }2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
- l) A# |+ ~8 m& D* q看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
V=v;
+ [. s% T: Y3 B' ^9 XN=n;
R=r;
2 _, w! ]" z% ~6 L+ _  QRV=rv;
- M+ s& m1 ~7 ]& s5 Q1 vfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
! W5 s2 Y+ N6 S' t9 X1 e: R- U    z(1,k)=k;
* u& k/ X# z$ n" s    z(2,k)=A(k);
2 @$ A/ W# s0 Tend

不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
" k* {& ?. S8 \" w1 LA(1)=a;
( X0 K8 d- t- q' B2 U3 AV=v;
N=n;
4 F  z* T- Z: T0 U  Z- O. |R=r;
# |2 Y" _3 K$ O  M$ `" O- VRV=rv;
for 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);
end

dreamofyou

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

厚积薄发

1.出现error:A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);因为你的等式右边出现了A(k)，这属于非法赋值
1 j5 F! S0 W* D% Y, M9 O所以应该修改成：A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2）
/ L5 i( K; D) y! Q4 U1 }- }- C2.不论是写成A(k)，还是A(k，1)，只要达到能表示向量的目的就行了，因为差分方程研究的是一个向量，而非一个矩阵
3.例子：这个就比较多了，根据差分方程的不同，主要改变的是上面的函数中的条件式子，但要始终注意赋值的合法性。

dreamofyou

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

$ N+ \) \! s2 l+ f: `' t' S: i% J& F. o谢谢谢  我终于明白了