关于离散微分方程的求解问题

semeng

大家看看下面的代码
就是将 速度在每个时间间隔内的积分然后+原来的位移=新的位移。
方程表示为
x (Tk+1)=| (Tk->Tk+1) v(T) dT + x (Tk)
代码如下：
        t=[1:1:10]; % Time
        v=[65 95 110 150 170 210 240 265 272]; % velocity
        x(1)=0; % To determine the positon
        for k=[1:9]
            x(k+1)=trapz(t(k:k+1),v(k:k+1))+x(k);
        end
        x
总是有错误

[ 本帖最后由 eight 于 2008-4-28 17:18 编辑 ]

sogooda

贴个错误提示来看看

semeng

??? Index exceeds matrix dimensions.
Error in ==> chapter8 at 74
            x(k+1)=trapz(t(k:k+1),v(k:k+1))+x(k);

ch_j1985

原帖由 semeng 于 2008-4-26 09:15 发表

                               
登录/注册后可看大图

大家看看下面的代码
就是将 速度在每个时间间隔内的积分然后+原来的位移=新的位移。
方程表示为
x (Tk+1)=| (Tk->Tk+1) v(T) dT + x (Tk)
代码如下：
        t=[1:1:10]; % Time
        v=[65 95 110 150 17 ...

速度和时间的维数不同
这段代码怎么样？

1. 
   
2. clc
   
3. t=1:1:9; % Time
   
4. v=[65 95 110 150 170 210 240 265 272]; % velocity
   
5. x(1)=0; % To determine the positon
   
6. for k=1:8
   
7.     x(k+1)=trapz(t(1:k),v(1:k))+x(k);
   
8. end
   
9. x
   

复制代码

semeng

Thank you, I got it

eight

原帖由 semeng 于 2008-4-26 09:15 发表

                               
登录/注册后可看大图

大家看看下面的代码
就是将 速度在每个时间间隔内的积分然后+原来的位移=新的位移。
方程表示为
x (Tk+1)=| (Tk->Tk+1) v(T) dT + x (Tk)
代码如下：
        t=[1:1:10]; % Time
        v=[65 95 110 150 17 ...

下次发帖前，先看置顶的会员守则，把问题描述清楚。另外，你的错误信息在置顶帖也有收录，你完全可据此自行解决，不要浪费我们的心血