解ODE 繪圖收斂的問題

inoran

請問各位；
我解一個ODE  Ay''+By=C*V
V是一個外力作用項  假設3*3的話
C=[1,0,0]
我用for寫出 當t<10 則 V=sin(3*t)
              其他時間V＝0
這是畫出來的圖理論上t超過10之後會開始收斂且震盪  時域[0 50]
我做出來是這樣子沒錯！
不過我改寫
if t<5
   V=sin(3*t)
else
V=0
end這時候就五秒後沒有收斂出現,反而會跟外力持續作用的圖類似(從頭震盪到尾無收斂出現)
若我是給20 30  40這種的都會在超過該時間開始收斂並震盪到時間末端
給5 15 6 16等等的就又沒收斂= ="
請問是我IF那邊有寫錯嗎？  謝謝！！

ChaChing

可否给出代表这两种(震荡收敛及持续无收敛)图形!?

还有建议简繁转换下(http://www.j4.com.tw/big-gb/), 不然看懂的人或许较少!

[ 本帖最后由 ChaChing 于 2010-3-27 22:42 编辑 ]

inoran

恩恩  多谢提醒,我附上图形,

分别为t<50 ,t<10 ,t<5,t<0.5

谢谢您！
=========

顺便附上我的程式档

副程式:testx.m
function dx=testx(t,x)
global  M K P n
if t<10
   V=sin(3*t)
else
   V=0                  
end
dx=[x(n+1:end,1);-inv(M)*K*x(1:n,1)+inv(M)*P*V]

主程式：
clear all
global  M K P n
n=3
m1=1 ;
m2=1 ;
m3=1;
k1=2;
k2=2;
k3=1;
k4=2;
M=[m1,0,0;0,m2,0;0,0,m3];
K=[k1+k2,-k2,0;-k2,k2+k3,-k3;0,-k3,k3+k4]
P=[1;0;0];
tinterv=0:0.1:50
yinit=[0 0 0 0 0 0]
[t,x]=ode45('testx',tinterv,yinit)
plot(t,x(:,1))

[ 本帖最后由 inoran 于 2010-3-28 17:59 编辑 ]

ChaChing

一时未看出问题, 初步猜测是否为IC影响!?
较有空时, 再研究看看! 或待高人路过

inoran

恩！我也觉得奇怪,也就是给t<10, 20, 30, 40sec
图形就会在该秒之後有收敛,若给t<5,15,25,35

或者不是十的倍数时 ,图就会怪怪的
也谢谢楼上这几次热心的回答！

ChaChing

处理了许多帖了, 有点累了! 该休息准备明天上班, 明晚再看看!