View source for Equal tuning accuracy functions

{{Technical}}

There are a number of ways to quantify the tuning accuracy of equal temperaments.

== Mu function ==
The [https://www.desmos.com/calculator/wanguxom9z mu function] is defined as:

<math> \mu( x, \sigma) = \sum_{k=1}^{\infty}f \left( x, k \right) </math>

for a given edo ''x'' and weight σ > 1, where

<math>f(x, k) = \frac{|(2x \log_{2}(k) \mod 2) - 1|}{k^{\sigma}}.</math>

The function essentially sums up the relative error on all integer harmonics ''k'', weighted by the inverse power to sigma of ''k'' in order to converge to a finite value.

It is derived as follows.

For each integer harmonic ''k'', the relative error on that integer in the continuum of equal tunings follows a zigzag line where 1 is an equal division of ''k'', and 0 is an odd equal division of 2''k'' (which has the largest possible error on ''k''). Such a zigzag line takes the form of:

<math>|(2x\mod 2) - 1|</math>

for {{nowrap|''k'' {{=}} 2}}, if integer values of ''x'' are EDOs.

Equal divisions of any integer ''k'' can be found by multiplying 2''x'' by <math>\log_{2}(k).</math>

As such, finding our final function is simply a matter of summing up

<math>|(2x \log_{2}(k) \mod 2) - 1|</math>

for all integers ''k''. To make the sum finite at all values, we weight each term by 1/''k''<sup>σ</sup>, producing our final formula for ''f'', and thus for μ.

== Zeta function ==
There is a derivation of the Riemann zeta function as a variant of the previously established mu function, which is occasionally used to determine optimal octave stretches for edos, and which allows the analytic continuation of the function to allow values of sigma less than 1. 

[https://www.desmos.com/calculator/deafikrhvg] <math> \sum_{k=1}^{\infty}\frac{\cos\left(\log_{2}\left(k\right)2\pi x\right)}{k^{\sigma}} </math>

Let's clean up the function by removing the scale factors on ''x''. This just scales the function's inputs from EDO to the equal divisions of a different interval (e<sup>2π</sup>/1, to be specific), and these can be added back later to go back to EDO.

[https://www.desmos.com/calculator/26ypbwbglg] <math> \sum_{k=1}^{\infty}\frac{\cos (\ln\left(k\right)x)}{k^{\sigma}} </math>

By the complex exponential theorem, which relates trigonometric functions and the exponential function, we know that

<math> e^{ix}=\cos\left(x\right)+i\sin\left(x\right) </math>

where ''i'' is the square root of −1, and ''e'' is the natural exponential constant.

''i'' is the imaginary unit, which is on a line perpendicular to the real number line. A complex (two-dimensional) number may be written as {{nowrap|''a'' + ''bi''}}.

With this knowledge, cos(''x'') can be rewritten as Re(''e<sup>ix</sup>'')—but since we're only doing multiplication and addition and this is the only place complex numbers appear, we can just ignore the Re() and add it back later.

[https://www.desmos.com/calculator/e7wn17tzjf] <math> \mu_{c}\left(\sigma, x\right)=\sum_{k=1}^{\infty}\frac{e^{i\left(\ln\left(k\right)x\right)}}{k^{\sigma}} </math>

Since {{nowrap|''e''<sup>ln(''n'')''x''</sup> {{=}} ''n''<sup>''x''</sup>}}, the exponentials and logarithms cancel each other out (i.e. {{nowrap|''e''<sup>ln(''n'')</sup> {{=}} ''n''}}), so:

[https://www.desmos.com/calculator/f4ojwn0an4] <math> \sum_{k=1}^{\infty}\frac{k^{ix}}{k^{\sigma}} </math>

[https://www.desmos.com/calculator/6388kalfmq] <math> \sum_{k=1}^{\infty}k^{ix}k^{-\sigma} </math>

[https://www.desmos.com/calculator/l3q2dtd6xn] <math> \sum_{k=1}^{\infty}k^{-\sigma+ix} </math>

{{nowrap|−σ + ''ix''}} is just a complex number, which we may write as −''s'':

[https://www.desmos.com/calculator/esbdlxdoui <math> \mu_{d}\left(s\right)=\sum_{k=1}^{\infty}k^{-s} </math>]

where, for {{nowrap|''s'' {{=}} σ − ''ix''}}, {{nowrap|Re(μ<sub>''d''</sub>(''s'')) {{=}} μ<sub>''c''</sub>(σ, ''x'')}}, our badness function.

μ<sub>''d''</sub> is a mathematical function called the Riemann zeta function, so {{nowrap|μ<sub>''d''(s)</sub> {{=}} ζ(''s'')}}, and re-adding the Re() function gives Re(ζ(''s'')) with {{nowrap|''s'' {{=}} σ − ''ix''}}; ''x'' is the equal division and σ is the weight.