各位版友好,
目前在做的事情是, 畫出不等式所圍出的空間並填滿之.
做法是先用 ContourPlot 畫出不等式的邊界範圍 (空間沒有被填滿, 只是畫上方程式),
接著使用 RegionPlot 找出不等式填滿的空間 (沒有邊界範圍, 但有被填滿的空間);
我想使用 Map 讓程式自動跑出一個有邊界又填滿的圖,
但一直卡住, 因此想請問是不是有更好的方法呢? 以下附上我正在進行的工作,
非常感謝各位的耐心.
ContourPlot[{(1 - \[Phi])*(ta*tm)^(\[Alpha]) + \[Phi]*(1 - \[Mu])*
ta - 1 == 0.1, (((1 - \[Phi])*(ta*tm)^(\[Alpha]) + \[Phi]*(1 - \[Mu])*ta -
1)/(2 ((1 - \[Phi])*(ta*tm)^(\[Alpha]) + \[Phi]*(1 - \[Mu])*ta) - 1)) == 0.1,
Ya (1 - T) - Yb (T (tm)^(2 (\[Sigma] - 1)) - 1) == 0,
(1 - \[Phi])*(ta*tm)^(\[Alpha]) + \[Phi]*(1 - \[Mu])*ta - 1 == 0},
{ta, 1, 1.2}, {tm, 1, 1.5}, Axes -> True,
AxesLabel -> {"ta", "tm"}, Frame -> False,
AxesStyle -> Arrowheads[{0, 0.02}],
BoundaryStyle -> {Red, Thickness[0.01]}]
RegionPlot[(1 - \[Phi])*(ta*tm)^(\[Alpha]) + \[Phi]*(1 - \[Mu])*ta - 1 > 0.1 &&
((1 - \[Phi])*(ta*tm)^(\[Alpha]) + \[Phi]*(1 - \[Mu])*ta - 1)/
(2 ((1 - \[Phi])*(ta*tm)^(\[Alpha]) + \[Phi]*(1 - \[Mu])*ta) - 1) < 0.1 &&
Ya (1 - T) > Yb (T (tm)^(2 (\[Sigma] - 1)) - 1), {ta, 1, 1.2}, {tm, 1, 1.5}]
__________________________________________________________________________
簡單來說, 就是要填滿三條不等式所圍出來的空間.
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.109.122.222
※ 編輯: bonus1123 來自: 140.109.122.222 (10/14 21:58)