Skip to main content

plotting - Plot of continuous function cuts off


Why the following plot cuts off at integer points ?


Plot[Cos[Pi IntegerPart[x]]
  Min[Abs[FractionalPart[x]], Abs[FractionalPart[x] - 1]], {x, 0, 10}]

enter image description here


Update


It gets connected when you add Exclusions -> None as an option, but It is not clear why this function has exclusive points, since it is continuous. Also it gets connected if you increase the PlotPoints (thanks to the comments). Also why at x=9 the cut is more severe than say at x=1 or x=2 ?



Can someone explain what is the issue here, and how these two options solve(or maybe hide) the issue ?



Answer



This behavior has been evolving ever since Mathematica was created. Something different happens in V11 than in V10. (I'm unable to go back further.) I'll describe V10 first, since its behavior conforms to the OP.


IN V10, the gaps are because the functions FractionalPart and IntegerPart are discontinuous, which makes the gaps symmetric. In plotting, Mathematica does not check limits at their discontinuities to see if the expression being plotted happens to be continuous. Rather, it assumes the discontinuities in the component functions propagate to discontinuities in the plot and puts a little gap in the plot at each one.


I assume this choice (not to check limits) was made for the sake of speed. The discontinuities are identified by a time-constrained symbolic analysis of the function.


The size of the gap is determined by the sampling. With higher setting for PlotPoints, the smaller the gap. The sampling can be observed using Mesh -> All. The sampling is a result of an asymmetric recursive subdivision of "active" subintervals, which depends principally on the angles between successive line segments used to approximate the graph. As a result, the gaps are not identical.


Plot[Cos[Pi IntegerPart[x]] Min[Abs[FractionalPart[x]], 
   Abs[FractionalPart[x] - 1]], {x, 0, 10}, Mesh -> All, 
 MeshStyle -> Red]


Mathematica graphics


The gaps are not visible at ordinary size with lots of plot points:


Plot[Cos[Pi IntegerPart[x]] Min[Abs[FractionalPart[x]], 
   Abs[FractionalPart[x] - 1]], {x, 0, 10}, PlotPoints -> 1000]

Mathematica graphics


Another approach, mentioned by @kglr in a comment, is to turn off the processing with Exclusions -> None.


Plot[Cos[Pi IntegerPart[x]] Min[Abs[FractionalPart[x]], 
   Abs[FractionalPart[x] - 1]], {x, 0, 10}, Exclusions -> None]


IN V11, Plot adds sampling points close to the discontinuity. It still does not connect the lines, but if the function is continuous, the gap should be imperceptible.


plot = Plot[
  Cos[Pi IntegerPart[x]] Min[Abs[FractionalPart[x]], Abs[FractionalPart[x] - 1]],
  {x, 0, 10}, Mesh -> All, MeshStyle -> Red]
Count[plot, _Line, Infinity]

Mathematica graphics


(*  20  <-- number of continuous lines *)

One continuous line is created with Exclusions -> None:



Count[
 Plot[Cos[Pi IntegerPart[x]] Min[Abs[FractionalPart[x]], 
    Abs[FractionalPart[x] - 1]], {x, 0, 10}, Exclusions -> None],
 _Line, Infinity]
(*  1  <-- number of continuous lines *)

Comments

Post a Comment

Popular posts from this blog

plotting - How to draw lines between specified dots on ListPlot?

I would like to create a plot where I have unconnected dots and some connected. So far, I have figured out how to draw the dots. My code is the following: ListPlot[{{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 4}, {2, 5}, {3, 6}, {4, 7}, {1, 7}, {2, 8}, {3, 9}, {4, 10}, {1, 10}, {2, 11}, {3, 12}, {4,13}, {2.5, 7}}, Ticks -> {{1, 2, 3, 4}, None}, AxesStyle -> Thin, TicksStyle -> Directive[Black, Bold, 12], Mesh -> Full] I have thought using ListLinePlot command, but I don't know how to specify to the command to draw only selected lines between the dots. Do have any suggestions/hints on how to do that? Thank you. Answer One possibility would be to use Epilog with Line : ListPlot[ {{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 4}, {2, 5}, {3, 6}, {4, 7}, {1, 7}, {2, 8}, {3, 9}, {4, 10}, {1, 10}, {2, 11}, {3, 12}, {4, 13}, {2.5, 7}}, Ticks -> {{1, 2, 3, 4}, None}, AxesStyle -> Thin, TicksStyle -> Directive[Black, Bold, 12], Mesh -> Full, Epilog -> { Line[ ...
Post a Comment
Read more

equation solving - Invert and fit implicitly defined curve

I need to fit an implicitly defined curve. I thought I could get some data out of Solve , and then using FindFit . Therefore, I would like to find the relation the parametric curve defined by $F(x,y)=0$: Solve[-(1/2) + 1/2 (0.41202 BesselK[0, 0.1 Sqrt[x^2 + y^2]] + (0.101483 x BesselK[1, 0.1 Sqrt[x^2 + y^2]])/Sqrt[x^2 + y^2]) == 0, y] But I can't get an output: Solve was unable to solve the system with inexact coefficients or the system obtained by direct rationalization of inexact numbers present in the system. Since many of the methods used by Solve require exact input, providing Solve with an exact version of the system may help. >> Edit: In particular, I would like to fit the data coming from the curve with the expression of another curve, and not with a function $f(x)$. In particular, since this clearly looks like a cardioid , I would like it to fit to something like it. What other strategies could I try?
Post a Comment
Read more

dynamic - How can I make a clickable ArrayPlot that returns input?

I would like to create a dynamic ArrayPlot so that the rectangles, when clicked, provide the input. Can I use ArrayPlot for this? Or is there something else I should have to use? Answer ArrayPlot is much more than just a simple array like Grid : it represents a ranged 2D dataset, and its visualization can be finetuned by options like DataReversed and DataRange . These features make it quite complicated to reproduce the same layout and order with Grid . Here I offer AnnotatedArrayPlot which comes in handy when your dataset is more than just a flat 2D array. The dynamic interface allows highlighting individual cells and possibly interacting with them. AnnotatedArrayPlot works the same way as ArrayPlot and accepts the same options plus Enabled , HighlightCoordinates , HighlightStyle and HighlightElementFunction . data = {{Missing["HasSomeMoreData"], GrayLevel[ 1], {RGBColor[0, 1, 1], RGBColor[0, 0, 1], GrayLevel[1]}, RGBColor[0, 1, 0]}, {GrayLevel[0], GrayLevel...
Post a Comment
Read more