sin.md

   1 # Sine
   2
   3 Sine, abbreviated *sin*, is a [trigonometric](trigonometry.md) [function](function.md) that simply said models a smooth oscillation, it is one of the most important and basic functions in geometry, [mathematics](math.md) and [physics](physics.md), and of course in [programming](programming.md). Along with [cosine](cos.md), [tangent](tan.md) and [cotangent](cot.md) it belongs to a group of functions that can be defined by ratios of sides of a right [triangle](triangle.md) depending on one of the angles in it (hence *trigonometric* -- "triangle measuring"). If some measurement looks like sine function, we say it is *harmonic*. This is very common in nature and technology, e.g. a weight on a spring goes up and down by this function, [alternating current](ac.md) voltage has the sine shape (because it is generated by a circular motion), the amount of light falling onto a surface depends on the light angle via sine function etc.
   4
   5 The function is most commonly defined using a right triangle as follows. Consider the following triangle:
   6
   7 ```
   8           /|
   9          / |
  10         /  |
  11       c/   |
  12       /    |a
  13      /     |
  14     /     _|
  15    /A____|_|
  16        b
  17 ```
  18
  19 *Sin(A)*, where *A* is the angle between side *b* and *c*, is the ratio *a* / *c*. The function can be defined in many other ways, for example it is the curve we get when tracking only one direction (e.g. horizontal) of a point moving alongside circle (well, in the end it's actually related to the triangle definition) -- imagine something making a circular motion while casting a shadow onto a nearby wall (in a direction that's in the plane of rotation), the shadow will go up and down by sine of time. It can also be defined as a solution to some [differential equations](differential_equation.md) etc.
  20
  21 The graph of the sine function is following:
  22
  23 ```
  24                   ^ sin(x)
  25                   |
  26                 1_|_
  27                   |     .--'''--.
  28       -1/2 pi     | _.''         ''._        3/2 pi
  29 .________|________.'________|________'|________|________.' --> x
  30  '._     |     _.'|0        |         |'._     |     _.'|
  31     ''--___--''  _|_      1/2 pi      pi  ''--___--''  2 pi
  32                -1 |
  33 ```
  34
  35 Basically we'll observe a sine function in places where something oscillates naturally, gradually, i.e. where something goes there and back repeatedly but in a gradual, smooth manner, slowing down as it's approaching one end, coming to stop for a while, then reversing direction, slowly accelerating the other way and so on. Very rare in nature something comes there and back with constant speed, coming to stop and reversing directions instantly.
  36
  37 **Why the fuck are there these [pi](pi.md) values on the x line???** Nubs often can't comprehend this. These pi values are values in **[radians](radian.md)**, units of measuring angles where *2 pi* is the full angle (360 degrees). In fact sine is sometimes shown with [degrees](degree.md) instead of radians (so imagine 90 degrees on the line where there is 1/2 pi etc.), but mathematicians prefer radians. **But why are there angles in the first place???** Why doesn't it go e.g. from 0 to 1 like all other nice functions? Well, it's because of the relation to geometry, remember the fucking triangle above... also if you define sine with a circle it all repeats after *2 pi*. Just draw some picture if you don't get it.
  38
  39 Some additional facts and properties regarding the sine functions are:
  40
  41 - The domain are all [real numbers](real_number.md), the [codomain](codomain.md) are real numbers in interval <-1,1> (including both bounds).
  42 - It is an [odd function](odd_function.md) (*-sin(x) = sin(-x)*).
  43 - It is periodic, with a period of 2 [pi](pi.md).
  44 - Sine is just shifted [cosine](cos.md), i.e. *sin(x) = cos(x - 1/2 pi)*
  45 - Its inverse function is [arcus sine](asin.md), abbreviated *asin*, also written as *sin^-1* -- this function tells you what argument you need to give to sin to get a specific result number. It's actually an inverse of only part of the sine function because the whole sine function can't be inverted, it isn't [bijective](bijection.md).
  46 - [Derivative](derivative.md) of *sin(x)* is *cos(x)*, the [integral](integral.md) of *sin(x) dx* is *-cos(x)*.
  47 - By adding many differently shifted and scaled sine functions we can create basically any other function, see e.g. [cosine transform](cosine_transform.md).
  48 - Sine and [cosine](cos.md) functions are used to draw [circles](circle.md). If you plot points with *x* coordinate equal to *sin(t)* and *y* coordinate equal to *cos(t)* for *t* going from 0 to *2 * pi*, you'll get a unit circle.
  49 - *sin(x)^2 + cos(x)^2 = 1*
  50
  51 Some values of the sine function are:
  52
  53 | x (rad)  | x (deg)  |  sin(x)            |
  54 |----------|----------|--------------------|
  55 | 0        | 0        | 0                  |
  56 | pi / 12  | 15       | ~0.259             |
  57 | pi / 6   | 30       | 0.5                |
  58 | pi / 4   | 45       | sqrt(2)/2 ~= 0.707 |
  59 | pi / 3   | 60       | sqrt(3)/2 ~= 0.866 |
  60 | pi / 2   | 90       | 1                  |
  61 | 2 pi     | 360      | 0                  |
  62
  63
  64
  65 ## Programming
  66
  67 In programming languages the sine function is generally available in some math library, for example in [C](c.md) the function `sin` is in `math.h`. Spare yourself bugs, **always check if your sin function expects [radians](radian.md) or degrees!**
  68
  69 **Want to make your own sine function for whatever reason (performance, curiosity, ...)?** Then firstly consider what you expect from it. If you want a small, fast and perhaps integer only `sin` function (the one we'd prefer in [LRS](lrs.md)) that doesn't need extreme accuracy, consider using a **[look up table](lut.md)**. You simply precompute the values of the sine function into a static table in memory and the function just retrieves them when called -- this is super fast. Note that you can save a lot of space by **only storing sine values between 0 and 1/2 pi**, the remaining parts of the function are just different transformations of this part. You can further save space and/or make the function work with [floats](float.md) by further [interpolating](interpolation.md) (even just linearly) between the stored values, for example if `sin(3.45)` is called and you only have values stored for `sin(3.4)` and `sin(3.5)`, you simply average them.
  70
  71 Lot of times, e.g. in many calculators where speed isn't really critical, sine is computed using [Taylor series](taylor_series.md) -- a sum of infinitely many terms of which if we take the first *N*, we get an [approximation](approximation.md) of the function (the more terms we add, the more precise we get). For sine the series is
  72
  73 *sin(x) = x - x^3 / 3! + x^5 / 5! - x^7 / 7! + ...*
  74
  75 Adding just the first 3 terms (*x - x^3 / 6 + x^5 / 120*) already gives a very accurate approximation in range <-pi/2,pi/2> (error < 0.5 %). Here is a [C](c.md) function that uses this to compute an 8bit sine (the magic numbers are made so as to incorporate pi while using power of two divisors, also note the use of many operations that will make the function relatively slow):
  76
  77 ```
  78 // x = 255 means full angle, returns 0 to 255
  79 unsigned char sin8(unsigned char x)
  80 {
  81   int a = x;
  82   char flip = 0;
  83
  84   if (a > 127)
  85   {
  86     a -= 128;
  87     flip = 1;
  88   }
  89
  90   if (a > 63)
  91     a = 128 - a;
  92
  93   int result = (411999 * a) - (a * a * a * 41);
  94
  95   a /= 4;
  96
  97   a = a * a * a * a * a;
  98
  99   result = (a + result) / 131072;
 100   return flip ? (127 - result) : (127 + result);
 101 }
 102 ```
 103
 104 If you just need a super fast and very rough sine-like value, there exists an **ugly engineering approximation** of sine that can be useful sometimes, it says that
 105
 106 *sin(x) = x, for small x*
 107
 108 Indeed, sine looks similar to a mere line near 0, but you can see it quickly diverges.
 109
 110 Very rough and fast approximations e.g. for primitive music synthesis can be done with the traditional very basic [square](square_function.md) or [triangle](triangle_function.md) functions. The following is a simple 8bit linear approximation that's more accurate than square or triangle (approximates sine with a linear function in each octant):
 111
 112 ```
 113 unsigned char sinA(unsigned char x)
 114 {
 115   unsigned char quadrant = x / 64;
 116
 117   x %= 64;
 118
 119   if (quadrant % 2 == 1)
 120     x = 63 - x;
 121
 122   x = x < 32 ? (2 * x + x) : (64 + x);
 123
 124   return quadrant <= 1 ? (128 + x) : (127 - x);
 125 }
 126 ```
 127
 128 Similar approximation can be made with a quadratic curve, the following is a modification of the above function that does this (notice that now we need at least 16 bits for the computation so the data type changed to int): { I quickly made this just now, maybe it can be improved. ~drummyfish }
 129
 130 ```
 131 int sinA(int x)
 132 {
 133   unsigned char quadrant = x / 64;
 134
 135   x %= 64;
 136
 137   if (quadrant % 2 == 1)
 138     x = 63 - x;
 139
 140   x -= 63;
 141   x = (x * x) / 32;
 142
 143   return quadrant <= 1 ? (255 - x) : x;
 144 }
 145 ```
 146
 147 Sine can also be surprisingly accurately approximated with the [smoothstep](smoothstep.md) function, which is just a polynomial *3 * x^2 - 2 * x^3*.
 148
 149 TODO: code for that
 150
 151 Furthermore there exist other nice approximations, such as the extremely accurate **Bhaskara I's approximation** (angle in radians): *sin(x) ~= (16 * x * (pi - x)) / (5 * pi^2 - 4 * x * (pi - x))*. (This formula is actually more elegant for cosine, so it may be even better to consider using that.) Here is a [C](c.md) [fixed point](fixed_point.md) implementation:
 152
 153 ```
 154 #define UNIT 1024
 155 #define PI ((int) (UNIT * 3.14159265))
 156
 157 /* Integer sine using Bhaskara's approx. Returns a number
 158 in <-UNIT, UNIT> interval. Argument is in radians * UNIT. */
 159
 160 int sinInt(int x)
 161 {
 162   int sign = 1;
 163
 164   if (x < 0) // odd function
 165   {
 166     x *= -1;
 167     sign = -1;
 168   }
 169
 170   x %= 2 * PI;
 171
 172   if (x > PI)
 173   {
 174     x -= PI;
 175     sign *= -1;
 176   }
 177
 178   x *= PI - x;
 179
 180   return sign * (16 * x) / ((5 * PI * PI - 4 * x) / UNIT);
 181 }
 182 ```