看板 Mathematica 關於我們 聯絡資訊
不好意思又要麻煩各位先進,一陣子沒用 MATHEMATICA 真的生疏到一個不像話境界... 目前在從事的工作是這樣的,將方程式代不同的參數值然後畫圖出來; 我的問題是橫坐標的設定問題,因為我的橫坐標有效的參數範圍是計算出來的。 簡單來說,整個流程大概是這樣: 參數代入--> 數個方程式求值--> 畫出目標函數,橫軸參數有效範圍由兩個方程式決定。 現在碰到的問題是,我如果直接用方程式當成橫軸的上下界, MATHEMATICA 會說:Plot::plln: ".....省略" is not a machine-size real number 這是要我把橫軸上下界改成固定數的意思吧....可是我如果一次只設定一個參數的話, 橫軸的上下界一樣打上方程式,這個圖是可以畫的。 另外的問題是因為我代入五組參數值,所以我希望看到五個圖;但如果我把橫坐標區間 改成固定的值的話,我會得到一個上面有五條線的圖... 不好意思麻煩各位,請問能不能給點提示呢? 謝謝! 以下是我目前打出來的東西。 Clear["Global`*"] f[ta_] := ta data = {f[#]} & /@ Range[1.7, 2.2, 0.1] {\[Alpha], \[Rho], \[Phi], \[Mu], a, ta, tm} = {0.7, 0.5, 0.5, 0.05, 0.4, data, 2.2} h = a (1 - \[Phi]) (ta*tm)^(\[Alpha]) + (1 - a)*\[Phi]*(1 - \[Mu])* ta - 1 h/(2 h + 1) i = \[Alpha] (1 - \[Rho])/(1 - \[Alpha] (1 - \[Rho])) \[Beta] = 1 - \[Alpha] \[Sigma] = 1/(1 - \[Rho]) T = ((tm)^(\[Alpha] + \[Sigma] - 1))/((ta)^(\[Beta] + \[Sigma] - 1)) Plot [{0.5 (1/k - 1/h)*((T - 1) ((ta - 1)*i + ta (1 - \[Phi] + \[Phi]*\[Mu])) - ((1 - \[Mu])*\[Phi]*ta - 1) ((tm)^(2 (\[Sigma] - 1)) - T)), ((tm)^(2 (\[Sigma] - 1)) - T - (T - 1) (ta + i (ta + 1)))}, {k, h/(2 h + 1), h}, Axes -> True, AxesLabel -> {"k", "f(k), g(k)"}] 再次感謝各位願意閱讀此文章。 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.109.122.132 ※ 文章網址: http://www.ptt.cc/bbs/Mathematica/M.1406706131.A.9A8.html ※ 編輯: bonus1123 (140.109.122.132), 07/30/2014 16:13:52
bonus1123:感謝戴老師解答,如果有需要的板有可以寄信給我喔! 07/30 18:32