Plotting a periodic function

asked 2022-08-29 02:54:24 +0200

dazedANDconfused
908 ●3 ●15 ●29

updated 2022-08-29 13:01:59 +0200

tmonteil
27333 ●31 ●203 ●514 http://wiki.sagemath.o...

I want to plot a function that is abs(x) for -1<= x <= 1 and repeats to create a sawtooth shaped function. With the code

def h(x):
    while x>1:
        x=x-2
    return abs(x)
h(4.3)

The Sage sell server tells me the value is 0.300000000000000, which is what I want as h(4.3)=h(2.3)=h(.3)=abs(.3) However, when I tried to plot h(x) with

def h(x):
    while x>1:
        x=x-2
    return abs(x)
plot(h(x),x,0,5)

The plot looks like abs(x), so that h(4.3) now appears to be 4.3. What has gone wrong? How can I properly plot the h(x) I want?

EDIT: A points plot gives me the output I am expecting, which is different than plot(h(x),x,0,5).

def h(x):
    while x>1:
        x=x-2
    return abs(x)

 points([(x, h(x)) for x in srange(0,5,.01)], pointsize=20)

edit retag flag offensive close merge delete

add a comment

2 Answers

Sort by » oldest newest most voted

answered 2022-08-29 13:01:45 +0200

tmonteil
27333 ●31 ●203 ●514 http://wiki.sagemath.o...

The way you define h makes it a Python function, so that, when you call h at the floating-point number 4.3, it is evalued by entering twice the while loop and then taking the absolute value of 4.3-2-2, that is 0.3.

However, when you call h at the symbolic expression x, the while loop is not entered since:

sage: bool(x>1)
False

hence, the abs function is called at the symbol x and what is returned is the abs(x) symbolic expression, see:

sage: abs
<built-in function abs>
sage: type(abs)
<class 'builtin_function_or_method'>

sage: abs(x)
abs(x)
sage: type(abs(x))
<class 'sage.symbolic.expression.Expression'>

Hence plot(h(x),x,0,5) is the same as plot(abs(x),x,0,5). The syntax plot(foo,x,0,5) is for plotting symbolic function, which is OK since abs(x) is a symbolic function, but it is not what you expected. If you want the plot function to evaluate the function foo at various floating-points between a and b to make the graphics, the syntax is just plot(foo,a,b), without mentionning any symbol.

Hence, you should simply call

sage: plot(h,0,5)

edit flag offensive delete link

Comments

Good explanation! It clears up some misunderstandings I had.

dazedANDconfused ( 2022-08-29 16:18:32 +0200 )edit

add a comment

answered 2022-08-29 17:19:12 +0200

slelievre
17654 ●22 ●160 ●348 http://carva.org/samue...

updated 2022-08-30 15:41:53 +0200

Instead of plot(h(x), x, -5, 5), try plot(h, (-5, 5)).

Note: the function may not be what you want for negative x.

One could also define h as a symbolic function without a while loop.

Here are two examples of doing that.

Using the floor function:

sage: h(x) = abs(x - 2 * floor((x + 1) / 2))
sage: plot(h, (-5, 5), aspect_ratio=1)
sage: plot(h(x), (-5, 5), aspect_ratio=1)

Using trigonometric functions:

sage: h(x) = arccos(cos(x*pi))/pi
sage: plot(h, (-5, 5), aspect_ratio=1)
sage: plot(h(x), (-5, 5), aspect_ratio=1)

In this case, both plot(h, (-5, 5)) and plot(h(x), (-5, 5)) give a correct plot.

edit flag offensive delete link

Comments

Yes, I was going to add another while loop which kept adding 2 for values <-1 if I needed to graph the negative portion. The h(x) you've defined, though, is different than my h(x) because in mine h(x)>=0.

dazedANDconfused ( 2022-08-30 00:08:29 +0200 )edit

Right, I had defined h(x) = x - 2 * floor((x + 1) / 2) but we want abs of that. Fixed now.

slelievre ( 2022-08-30 15:04:34 +0200 )edit

Nice, +1: I especially like the one using trig functions.

dazedANDconfused ( 2022-08-30 23:33:13 +0200 )edit

add a comment

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Plotting a periodic function

2 Answers

Comments

Comments

Your Answer

Question Tools

Stats

Related questions

Plotting a periodic function edit

2 Answers

Comments

Comments

Your Answer

Question Tools

Stats

Related questions

Plotting a periodic function