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

dreamofyou

为算一道题我写了如下程序：（保存为nongdu.m)
function z=nd(a,v,n,r,rv)
4 k7 E  ^/ A0 |A=a;
( t+ a. I1 F# c- R5 QV=v;
N=n;
/ j( l. t  |6 P: g* A5 @, cR=r;
RV=rv;
for k=1:N
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
    z(1,k)=k;
+ P" A+ ~, g4 ]6 m7 J7 V/ k    z(2,k)=A(k);
end;
然后z=nd(0.2,200*10^4,360,0.02,10^4);
! l+ D2 k7 F  P' y系统显示：
) E* [- i' f5 I# ???? Attempted to access A(2); index out of bounds because numel(A)=1.
7 `1 F2 `& `, f9 e/ \) |' UError in ==> nongdu at 8
    A(k)=A(k+1)+(R*RV+V*A(k))*RV/(V^2);
0 v* _% V$ W, H! w# d. N
9 M* Y, F1 A1 T( N# \6 I1.请问为什么啊？
, `5 L7 d% J& E2.请问差分方程在matlab里边要怎么写啊，我以前写的时候都是用一维的矩阵写的，如A（n，1）。
看书上直接可以写成A（k），到底应该如何呢？
3能不能有些教程或者实际的程序啊，当例子。

leo12ok

function z=nd(a,v,n,r,rv)
A=a;
9 p: F' g- \6 U% ~" EV=v;
N=n;
% O1 I; }, o4 Z9 |, l8 c$ UR=r;
RV=rv;
' R+ C( @0 o0 Y% \  x' @for k=1:N
    A(k+1)=A(k)-(R*RV+V*A(k))*RV/(V^2);
, a* z6 _' c8 w( x% P4 z. p    z(1,k)=k;
    z(2,k)=A(k);
end

不知道这样可不可以。

splendid1987

我试了一下，这样改一下就行了：
' f8 b( J0 s6 ^& O$ Dfunction z=nd(a,v,n,r,rv)
A(1)=a;
V=v;
3 U( K1 v8 e( `* o2 p) X5 x/ @' yN=n;
, p/ n3 [$ s' Q6 [0 JR=r;
* j# p2 E( [5 v! N0 KRV=rv;
1 S0 j* @2 d: J/ K( ffor 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);
+ `# E, a- e8 H4 M. @  x# yend

dreamofyou

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

厚积薄发

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

dreamofyou

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

% e9 e% |6 G- a+ ]7 M% ~  L% o谢谢谢  我终于明白了