【请教】ODE中如何设置条件语句？

tangxi

我用ODE解运动方程：Mx''+Cx'+Kx=f
下面是引用的程序

function tangxioder
clear;clc
x0=zeros(6,1);
[t,x]=ode45(@tangxi,[0,100],x0);
plot(t,x(:,1))

function dx=tangxi(t,x)
f=[sin(t);0;0];
M=eye(3);
C=eye(3)*0.1;
K=eye(3)-0.5*diag(ones(2,1),1)-0.5*diag(ones(2,1),-1);
dx=[x(4:6);inv(M)*(f-C*x(4:6)-K*x(1:3))];

我的问题是：我把f改成这样行不行？
if t<=3,
     f1=t*100;
elseif t>3&t<=6,
     f1=(6-t)*100;
else
     f1=0;
end
f=[f1;0;0];
如果行的话，该如何安插？我试过放在原来的位置取代原来的 f 好像行不通

想了一整天了，烦请各位高手指教，先谢过了！

[ 本帖最后由 tangxi 于 2008-11-11 09:54 编辑 ]

sigma665

原来的f是3*1的

tangxi

不好意思，我改了一下

ch_j1985

将t分为三段求三次，不知可否？

tangxi

我尝试了您说的，语言编排上老出错
而且分三段理论上行的通吗？

无水1324

这里面是对时间进行分段，上面的做法是不合理的
需要调用ｅｖｅｎｔ这个函数，你查一下ｏｄｅ４５
就知道了

tangxi

感谢院长，我按您说的试试