: 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)
推
08/10 12:12, , 1F
08/10 12:12, 1F
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 2 之 2 篇)：