各位版友好,
目前在做的事情是, 画出不等式所围出的空间并填满之.
做法是先用 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}]
__________________________________________________________________________
简单来说, 就是要填满三条不等式所围出来的空间.