Maple计算积分胜出Mathematica一局（案例）

最新推荐文章于 2022-08-17 11:09:09 发布

stereohomology

最新推荐文章于 2022-08-17 11:09:09 发布

阅读量3.6k

点赞数

分类专栏：学习学习文章标签： mathematica maple 微积分分段函数

本文链接：https://blog.csdn.net/stereohomology/article/details/52371915

版权

学习学习专栏收录该内容

670 篇文章 14 订阅

订阅专栏

——虽然Maple的界面用起来比较不称心，我从来不敢否认它的强大。这不，碰上一个简单定积分计算的例子：Maple可以计算出来，mathematica不支持。

问题

一个包含 向下取整 函数的定积分：

\int 10 x ⌊ 1 x ⌋ d x = ?

$\int_{0}^1x\left\lfloor \dfrac{1}x \right\rfloor{\rm d}x=?$

解答

一些研究

分段函数piecewise求导数、积分原函数都不能用光滑连续函数类似的办法。通常解决较难问题的一般方法是 把待求解的陌生或难的问题，通过一些变换，转化成已经求解过的问题的形式、或者相对熟悉的容易解决的问题的形式

变换和转化有不同的类型，取决于思路、经验、视野。

我想到的是，原始问题中的定积分，实际就是求面积。

被积函数绘制出来时分段的如下形式；绘出积分对应的面积到蓝色部分：

<<MaTeX`


Plot[x Floor@(1/x)//Evaluate,
{x,0.01,1},
PlotStyle->Directive[Black,Thin],
Filling->Axis,FillingStyle->Blue,
AxesOrigin->{0,0}, AspectRatio->1,ImageSize->800,WorkingPrecision->30,
PlotPoints->10000,
FrameStyle->BlackFrame,Frame->True,

FrameTicks->{ {Table[{x,MaTeX[x,"DisplayStyle"->True,Magnification->3/2]},{x, Range[0,5]/5}],None},
{Table[{x,MaTeX[x,"DisplayStyle"->True,Magnification->3/2]},{x, 1/Range[1,7]}],None}}

]

这里写图片描述

边长为1的正方形面积减去一系列小直角三角形的面积之和即可。这些小三角形的上直角边长度实际对应于下面的值：

{1 - 1 2, 1 2 - 1 3, 1 3 - 1 4, \dots, 1 k - 1 k + 1, \dots}

$\left\{1-\dfrac{1}{2},\dfrac{1}{2}-\dfrac{1}{3},\dfrac{1}{3}-\dfrac{1}{4},\cdots,\dfrac{1}{k}-\dfrac{1}{k+1},\cdots\right\}$

侧直角边长为1减去下角点纵坐标。所有的下角点都在直线 $y=1-x$ 上。从而容易得到蓝色部分的总面积，也就是定积分的数值：

S BLUE = 1 - \sum k = 1 + \infty 1 2 k ( k + 1 ) 2 = π 2 12

${\Large S}_{\text{BLUE}}=1-\sum _{k=1}^{+\infty } \frac{1}{2 k (k+1)^2}=\dfrac{\pi^2}{12}$

这个级数很奇妙。利用这个结果： $\displaystyle\sum\limits_{k=1}^{+\infty}\dfrac{1}{k^2}=\dfrac{\pi^2}{6}$ ，可以得到：

\sum k = 1 + \infty 1 k ( k + 1 ) 2 = \sum k = 1 + \infty ((1 k - 1 k + 1) - 1 ( k + 1 ) 2)

$\sum\limits_{k=1}^{+\infty}\dfrac{1}{k(k+1)^2}=\sum\limits_{k=1}^{+\infty}\left(\left(\dfrac{1}{k}-\dfrac{1}{k+1}\right)-\dfrac{1}{(k+1)^2}\right)$

= \sum k = 1 + \infty (1 k - 1 k + 1) - \sum k = 1 + \infty 1 ( k + 1 ) 2 = 1 - (π 2 6 - 1) = 2 - π 2 6

$=\sum\limits_{k=1}^{+\infty}\left(\dfrac{1}{k}-\dfrac{1}{k+1}\right)-\sum\limits_{k=1}^{+\infty}\dfrac{1}{(k+1)^2}=1-\left(\dfrac{\pi^2}{6}-1\right)=2-\dfrac{\pi^2}{6}$
最终结果从而可以知道。

把图从新画一下，

Plot[{1, x Floor@(1/x) // Evaluate}, {x, 0.1^5, 1}, 
 PlotStyle -> {Red, Directive[None, Thin]}, 
 Filling -> {1 -> {{2}, {Red}}}, AxesOrigin -> {0, 0}, 
 AspectRatio -> 1, ImageSize -> 800, WorkingPrecision -> 30, 
 PlotPoints -> 10000, FrameStyle -> BlackFrame, Frame -> True, 
 FrameTicks -> {{Table[{x, 
      MaTeX[x, "DisplayStyle" -> True, Magnification -> 3/2]}, {x, 
      Range[0, 5]/5}], 
    None}, {Table[{x, 
      MaTeX[x, "DisplayStyle" -> True, Magnification -> 3/2]}, {x, 
      1/Range[1, 7]}], None}}]