: syms t : FT =-495/8/(-15*exp(-1/2*t^4)*exp(-17/32*t^4) + : 16*exp(-1/2*t^4)^2)*t^3*exp(-1/2*t^4)*exp(-17/32*t^4) + : 64/(-15*exp(-1/2*t^4)*exp(-17/32*t^4) + : 16*exp(-1/2*t^4)^2)*t^3*exp(-1/2*t^4)^2; : (上面部分有點過長...其實是同一條式子...>O<) 問題可能是出在這 (-15*exp(-1/2*t^4)*exp(-17/32*t^4) +16*exp(-1/2*t^4)^2) 畫圖出來看 testFT= (-15*exp(-1/2*t^4)*exp(-17/32*t^4) + 16*exp(-1/2*t^4)^2) ; ezplot(FT,[0,10]) 可以看出t在2之後都趨近於0 (-15*exp(-1/2*t^4)*exp(-17/32*t^4) +16*exp(-1/2*t^4)^2) 你在分母的這道式子在t約>5.1之後 其趨近於0的值太小 超過matlab能夠接受的範圍 因此分母都會有一為0的項 所以算出來的值會怪怪的吧 : ezplot(FT,[0,10]) : double(subs(FT,t,0)) : double(subs(FT,t,5)) : double(subs(FT,t,10)) : double(int(FT,t,0,5)) : double(int(FT,t,0,10)) : double(int(FT,t,5,10)) : 跑完的結果會有一張圖...但圖只畫t=0~5 5之後到10完全沒有圖形... : 而下面的6個double指令...結果如下 : 0 , 500 , NAN , 622.2274 , 9.9972e+003 , 9.3750e+003 : 我的問題是: : 我做圖範圍是t=0~10 但t從5之後就沒有圖形... : 我算這函數當t=10沒有值... : 不過我做積分 範圍t從 0~10 卻有結果 9.9972e+003 : 範圍t從 5~10 也有結果 9.3750e+003 : 照理說積分積得出值代表有面積...應該要有圖形曲線... : 可是在5~10的範圍內卻無法做圖...且函數值為 NAN... : 這樣真的有點奇怪...>O< : 因此希望有好心的大大可以幫忙解答...^^" : 真的由衷地感謝.... --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.103.220 ※ 編輯: suzu 來自: 140.112.103.220 (08/09 23:21) ※ 編輯: suzu 來自: 140.112.103.220 (08/09 23:21)