解ODE 繪圖收斂的問題
請問各位；我解一個ODEAy''+By=C*V
V是一個外力作用項假設3*3的話
C=
我用for寫出 當t<10 則 V=sin(3*t)
其他時間V＝0
這是畫出來的圖理論上t超過10之後會開始收斂且震盪時域
我做出來是這樣子沒錯！
不過我改寫
if t<5
V=sin(3*t)
else
V=0
end這時候就五秒後沒有收斂出現,反而會跟外力持續作用的圖類似(從頭震盪到尾無收斂出現)
若我是給20 3040這種的都會在超過該時間開始收斂並震盪到時間末端
給5 15 6 16等等的就又沒收斂= ="
請問是我IF那邊有寫錯嗎？謝謝！！ 可否给出代表这两种(震荡收敛及持续无收敛)图形!?
还有建议简繁转换下(http://www.j4.com.tw/big-gb/), 不然看懂的人或许较少!
[ 本帖最后由 ChaChing 于 2010-3-27 22:42 编辑 ] 恩恩多谢提醒,我附上图形,
分别为t<50 ,t<10 ,t<5,t<0.5
谢谢您！
=========
顺便附上我的程式档
副程式:testx.m
function dx=testx(t,x)
globalM K P n
if t<10
V=sin(3*t)
else
V=0
end
dx=
主程式：
clear all
globalM K P n
n=3
m1=1 ;
m2=1 ;
m3=1;
k1=2;
k2=2;
k3=1;
k4=2;
M=;
K=
P=;
tinterv=0:0.1:50
yinit=
=ode45('testx',tinterv,yinit)
plot(t,x(:,1))
[ 本帖最后由 inoran 于 2010-3-28 17:59 编辑 ] 一时未看出问题, 初步猜测是否为IC影响!?
较有空时, 再研究看看! 或待高人路过 恩！我也觉得奇怪,也就是给t<10, 20, 30, 40sec
图形就会在该秒之後有收敛,若给t<5,15,25,35
或者不是十的倍数时 ,图就会怪怪的
也谢谢楼上这几次热心的回答！ 处理了许多帖了, 有点累了! 该休息准备明天上班, 明晚再看看!
页:
[1]