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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
A=a;
V=v;
N=n;
# L- t7 e$ g! F3 z4 ZR=r;
$ R1 Y/ w1 I% N' ]6 w% H& LRV=rv;
for k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
$ C9 J# n. y+ h  Q6 l3 Q    z(1,k)=k;
    z(2,k)=A(k);
end;
/ V8 e5 f' j6 e! g; V/ `& W1 L然后z=nd(0.2,200*10^4,360,0.02,10^4);
系统显示：
??? 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);
2 F# O! l6 {  B, J# I
1.请问为什么啊？
8 o# ~+ D' d; _$ j" X% _0 A8 @2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
+ Y' C( u8 j8 S( N2 t看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
7 E/ H4 o3 }+ C0 Y5 i" L9 o7 `A=a;
& J. ~% o% P, }/ D" tV=v;
3 _+ k8 K6 o0 g* S  U0 h9 @* H- mN=n;
R=r;
5 l6 _* R% u1 W, x' U$ }RV=rv;
, R* r& \/ L- |2 Zfor k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
, Z4 o& R" z5 o( m    z(1,k)=k;
2 x8 E8 S. ]9 I0 W    z(2,k)=A(k);
end
6 `  |% A' L! R* h7 Y& e* h
不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
function z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
% J' G6 W5 G4 @% H  q  \' s$ ]" r# eN=n;
R=r;
RV=rv;
for k=1:N
+ J" H' x) x2 p    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
    z(2,k)=A(k);
% |! J- e3 I' l4 \1 w/ Aend

dreamofyou

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

厚积薄发

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

dreamofyou

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

谢谢谢  我终于明白了