Message Boards

GROUPS:

Staff Picks Mathematics Calculus Discrete Mathematics Graphics and Visualization Wolfram Language Computer-Based Maths

Shenghui Yang

Understanding a CRUX problem: visualizing using LogLinearPlot

Shenghui Yang, WOLFRAM

Posted 10 months ago

2190 Views

5 Replies

4 Total Likes

Follow this post

Problem Given b is greater than 1, find the value of the following integral Out[]= Visualization Michael Penn discussed this eccentric integration from CRUX on this channel recently. Please use the link to watch the video. The main idea is to show that, given b > 1, the integrand is zero when x is no less than 1. For x is strictly between 0 and 1, the integrand is a piecewise function with different constant value on infinite number of partitions such as (1/ k b , 1/ k+1 b ) for k in integer domain. Eventually the final result comes down to summing up a geometric series. I feel LogLinearPlot helps to visualize the problem in a straightforward manner. For convenience, I can let b = , the base of nature logarithm and take a look at the case with finite terms. For instance, I pick up 16 terms In[]:= nTerms=16; In[]:= range=Range[nTerms] Out[]= {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16} The key to proper visualization is to set proper ticks to represent the switching positions. In Ticks option, I can add symbolic logarithm base to the x axis In[]:= xTicks=With[{val=^(-range),label=Power[b,HoldForm[-#]]&/@range},Transpose[{val,label}]]; In[]:= yTicks=range; If I plot the integrand with LogLinearPlot, the interesting pattern of ascending stair (right to left) definitely gives away the hidden secret of the problem. In[]:= plot=LogLinearPlot[Floor[Log[Floor[Ceiling[x]/x]]],{x,^-nTerms,1},PlotStyle{Black,Thickness[.007]},ExclusionsNone,Ticks{xTicks,yTicks},TicksStyleDirective["Label",11]] Out[]= In Michael’s video he mentioned there are two ways to do the integral and I use two different color stripe systems to show them. The first one is just sum up the area of the vertical color stripes. Notice that if any graphics object is added to LogLinearPlot, the x coordinate must be the log b value of the input. Say, I want to add a point at ( -3 b ,2), the actual coordinate used in Epilog-like option is (-3,2) In[]:= nthVerticalRect[n_]:=Rectangle[{-n-1,0},{-n,n-0.05}] In[]:= Show[plot,PlotStyleRed,Epilog{Opacity[0.6],{Hue[#/nTerms],nthVerticalRect[#]}&/@Most[range]}] Out[]= Because the area of each stripe is width times height, the infinite integral is just the sum of all vertical stripes. Warning: the color stripes are NOT of constant width, if they are plotted on linear scale for x axis, the widths are squeezed at exponential speed as it goes from 1 to 0 In[]:= res1=Sumn* 1 b^n - 1 b^(n+1) ,{n,1,Infinity},Assumptionsb>1 Out[]= 1 -1+b The second way is to flip the color stripe and I can sum of the area of horizontal ones. In[]:= nthHorizontalRect[n_]:=Rectangle[{-nTerms,n-1},{-n-0.05,n}] In[]:= Show[plot,Epilog{Opacity[0.6],{Hue[#/nTerms],nthHorizontalRect[#]}&/@Most[range]}] Out[]= This way I have a simpler expression for area sum because of the constant height of each color stripe and each width is a term of geometric series, thus the integral is In[]:= res2=Sum 1 b^n ,{n,1,Infinity},Assumptionsb>1 Out[]= 1 -1+b Verify two results are the same In[]:= res1res2 Out[]= True Summary ◼ LogLinearPlot is useful to visualize geometric series ◼ Ticks in LogLinearPlot can be set along with customizable labels ◼ Using graphics object in log plot needs coordinate conversion

POSTED BY: Shenghui Yang

Answer

5 Replies

Sort By:

Mariusz Iwaniuk

Mariusz Iwaniuk, engineer,hobbyist.

Posted 10 months ago

Great post,but: Problem 1. << MaTex` ($Failed) Problem 2: Integrate[Floor[Log[b, Floor[Ceiling[x]/x]]], {x, 0, Infinity}, Assumptions -> b > 1] (Integrate[Floor[Log[Floor[Ceiling[x]/x]]/Log[b]], {x, 0, \[Infinity]}, Assumptions -> b > 0])(Give me input) I executed this code on Mathematica 12.2.0.

Great post,but:

Problem 1.

<< MaTex`
(*$Failed*)

Problem 2:

Integrate[Floor[Log[b, Floor[Ceiling[x]/x]]], {x, 0, Infinity}, Assumptions -> b > 1]
(*Integrate[Floor[Log[Floor[Ceiling[x]/x]]/Log[b]], {x, 0, \[Infinity]},
  Assumptions -> b > 0]*)(*Give me input*)

I executed this code on Mathematica 12.2.0.

POSTED BY: Mariusz Iwaniuk

Answer

Shenghui Yang

Shenghui Yang, WOLFRAM

Posted 10 months ago

Please follow the link to install MaTeX and required packages. In general you do not need that to run the code because I only use it to display large formula and no computation depends on it. This is a very tricky integral with very specific strategy not included in automation. Typically breaking down domain for integration to (b^(k+1),b^k) is quite unusual and probably not used anywhere else

POSTED BY: Shenghui Yang

Answer

Mariusz Iwaniuk

Mariusz Iwaniuk, engineer,hobbyist.

Posted 10 months ago

Ok thanks.

POSTED BY: Mariusz Iwaniuk

Answer

Rohit Namjoshi

Posted 10 months ago

Mariusz, For the first problem `MaTeXInstall`. Then Needs["MaTeX`"]

POSTED BY: Rohit Namjoshi

Answer

Moderation Team

Moderation Team, WOLFRAM

Posted 10 months ago

-- you have earned *Featured Contributor Badge* Your exceptional post has been selected for our editorial column *Staff Picks* http://wolfr.am/StaffPicks and Your Profile is now distinguished by a *Featured Contributor Badge* and is displayed on the Featured Contributor Board. Thank you!