- 注册时间
- 2013-10-24
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 8856
- 在线时间
- 小时
|
发表于 2018-1-28 22:22:27
|
显示全部楼层
本帖最后由 kastin 于 2018-1-28 22:55 编辑
完全图可以做出来,包括边界颜色。其中红色表示三个交点,绿色表示两个交点,蓝色表示1个交点,白色表示没有交点,品红虚线表示间断线。
- Clear[f, g1, g2, x0, x1, s]
- f[x_ /; x > 0 && x != 1, a_ /; a > 0 && a != 1] :=
- Log[a, Log[a, x]] - x;
- (*使用超越函数表示极值点*)
- x0[a_ /; 0 < a <= E^(-E)] := Exp@ProductLog[-1, 1/Log[a]];
- x1[a_ /; 0 < a] := Exp@ProductLog[1/Log[a]];
- (* (*或 数值求根方式表示极值点*)
- g1[0]=0;
- g1[x_]:=1/(x Log[x] Log[#])-1&;
- g2[x_]:=-(x+Log[x])/(x Log[x])^2/Log[#]&;
- x0[a_/;0<a<=E^(-E)]:=x/.First@FindRoot[g1[x][a],{x,1/E-a}];
- x1[a_/;0<a<=E^(-E)||a>1]:=x/.First@FindRoot[g1[x][a],{x,1/E+a}]; *)
- (*三个交点*)
- s[1] = RegionPlot[
- f[x1[a], a] < Log[a, Log[a, b]] < f[x0[a], a] &&
- 0 < a < E^(-E), {a, 0, E^(-E)}, {b, 0, 1}, PlotStyle -> Red,
- BoundaryStyle -> {Thick, Green}];
- (*一个交点*)
- s[2] = RegionPlot[
- Log[a, Log[a, b]] > f[x0[a], a], {a, 0, E^(-E)}, {b, 0, E^(-E)},
- PlotStyle -> Blue, BoundaryStyle -> None];
- s[3] = RegionPlot[
- 0 < a < E^(-E) && b > 1, {a, 0, E^(-E)}, {b, 1 + 0.00001, 5},
- PlotStyle -> Blue, BoundaryStyle -> None];
- s[4] = RegionPlot[
- Log[a, Log[a, b]] < f[x1[a], a], {a, 0, 1}, {b, 0, E^(-E)},
- PlotStyle -> Blue, BoundaryStyle -> None];
- s[5] = RegionPlot[a > E^(-E) && 0 < b < E^(-E), {a, E^(-E), 5}, {b, 0, E^(-E)},
- PlotStyle -> Blue, BoundaryStyle -> None];
- s[6] = RegionPlot[(0 < a < 5 && E^(-E) < b < 1) || (E^(-E) < a < 1 &&
- 0 < b < 5), {a, 0, 5}, {b, 0, 5}, PlotStyle -> Blue,
- BoundaryStyle -> None];
- (*两个交点*)
- s[7] = RegionPlot[
- Log[a, Log[a, b]] < f[x1[a], a], {a, 1, 5}, {b, 1, 5},
- PlotStyle -> Green, BoundaryStyle -> {Thick, Blue}];
- (*间断线*)
- s[8] = ListLinePlot[{{{1, 0}, {1, 5}}, {{0, 1}, {5, 1}}},
- PlotStyle -> {{Magenta, Thick, Dashed}}];
- Show[Table[s[k], {k, 1, 8}], PlotRange -> All, AspectRatio -> 1,
- FrameLabel -> {a, b}]
复制代码
注意左下角较小范围放大图是
|
评分
-
参与人数 1 | 威望 +12 |
金币 +12 |
贡献 +12 |
经验 +12 |
鲜花 +12 |
收起
理由
|
wayne
| + 12 |
+ 12 |
+ 12 |
+ 12 |
+ 12 |
很给力! |
查看全部评分
|