The Riemann zeta function and tuning: Difference between revisions

Line 38:

A link to the graph of zeta can be found at [https://samuelj.li/complex-function-plotter/#abs(zeta(i*2*pi*real(z)%2Fln(2)%2Bimag(z))) Zeta in Samuelj Plotter]. (In the top left menu, make sure that "Enable Checkerboard" is unticked and "Invert Gradient" and "Continuous Gradient" are ticked.) The function has been reoriented to place EDO size along the horizontal axis and weight along the vertical axis, and also scaled by <math>\tfrac{2\pi}{\ln(2)}</math> to ensure that the real number line aligns with edos. One can see that with higher weights, the function approaches a cyclic function with a period of 1; this corresponds to the prime 2 dominating more and more extremely as other harmonics are weighted less with higher weights. You can see this easier by raising the entire expression to an absurdly high power, such as 100. Note, however, that this visualization is inaccurate beyond a couple hundred: around 146.5, 324.5 and 473.5, and in many cases after, there appear to be zeroes that are not on the critical line; this is an artifact of the way the function is approximated and is the ultimate reason why the Riemann hypothesis remains unsolved. These actually correspond to zeroes that are very close together but on the critical line.

You may also view the graph of zeta along the critical line on Desmos: [https://www.desmos.com/calculator/dstp7wnidf Zeta in Desmos]. This makes it easier to see peaks, but only works for ~~sigma~~=1/2.

You may also view the graph of zeta along the critical line on Desmos: [https://www.desmos.com/calculator/dstp7wnidf Zeta in Desmos]. This makes it easier to see peaks, but only works for {{nowrap|σ {{=}} {{sfrac|1|2}}}}.

=== Plots ===

Line 45:

Using the [http://functions.wolfram.com/webMathematica/FunctionPlotting.jsp?name=RiemannSiegelZ online plotter] we can plot Z in the regions corresponding to scale divisions, using the conversion factor {{nowrap|''t'' {{=}} {{sfrac|2π|ln(2)}}''x''}}, for ''x'' a number near or at an edo number. Hence, for instance, to plot 12 plot around 108.777, to plot 31 plot around 281.006, and so forth. An alternative plotter is the applet [http://web.viu.ca/pughg/RiemannZeta/RiemannZetaLong.html here].

If you have access to {{w|Mathematica}}, which has Z, zeta and theta as a part of its suite of initially defined functions, you can do even better. Below is a Mathematica-generated plot of Z({{frac|2π''x''|ln(2)}}) in the region around 12edo:

If you have access to {{w|Mathematica}}, which has Z, zeta and theta as a part of its suite of initially defined functions, you can do even better. Below is a Mathematica-generated plot of Z({{frac|2π''x''|ln(2)}}) in the region around [[12edo]]:

[[File:plot12.png|alt=plot12.png|plot12.png]]

Line 51:

The peak around 12 is both higher and wider than the local maximums above 11 and 13, indicating its superiority as an approximation of JI. Note also that the peak occurs at a point slightly larger than 12; this indicates the octave is slightly compressed in the zeta tuning for 12. The size of a step in octaves is 1/''x'', and hence the size of the octave in the zeta peak value tuning for ''N''edo is ''N''/''x''; if ''x'' is slightly larger than ''N'' as here with {{nowrap|''N'' {{=}} 12}}, the size of the zeta tuned octave will be slightly less than a pure octave. Similarly, when the peak occurs with ''x'' less than ''N'', we have stretched octaves.

For larger edos, the width of the peak narrows, but for strong edos the height more than compensates, measured in terms of the area under the peak (the absolute value of the integral of Z between two zeros.) Note how ~~270~~ completely dominates its neighbors:

For larger edos, the width of the peak narrows, but for strong edos the height more than compensates, measured in terms of the area under the peak (the absolute value of the integral of Z between two zeros.) Note how [[270edo]] completely dominates its neighbors:

[[File:plot270.png|alt=plot270.png|plot270.png]]

@@ Line 38: / Line 38: @@
 A link to the graph of zeta can be found at [https://samuelj.li/complex-function-plotter/#abs(zeta(i*2*pi*real(z)%2Fln(2)%2Bimag(z))) Zeta in Samuelj Plotter]. (In the top left menu, make sure that "Enable Checkerboard" is unticked and "Invert Gradient" and "Continuous Gradient" are ticked.) The function has been reoriented to place EDO size along the horizontal axis and weight along the vertical axis, and also scaled by <math>\tfrac{2\pi}{\ln(2)}</math> to ensure that the real number line aligns with edos. One can see that with higher weights, the function approaches a cyclic function with a period of 1; this corresponds to the prime 2 dominating more and more extremely as other harmonics are weighted less with higher weights. You can see this easier by raising the entire expression to an absurdly high power, such as 100. Note, however, that this visualization is inaccurate beyond a couple hundred: around 146.5, 324.5 and 473.5, and in many cases after, there appear to be zeroes that are not on the critical line; this is an artifact of the way the function is approximated and is the ultimate reason why the Riemann hypothesis remains unsolved. These actually correspond to zeroes that are very close together but on the critical line.
-You may also view the graph of zeta along the critical line on Desmos: [https://www.desmos.com/calculator/dstp7wnidf Zeta in Desmos]. This makes it easier to see peaks, but only works for sigma=1/2.
+You may also view the graph of zeta along the critical line on Desmos: [https://www.desmos.com/calculator/dstp7wnidf Zeta in Desmos]. This makes it easier to see peaks, but only works for {{nowrap|σ {{=}} {{sfrac|1|2}}}}.
 === Plots ===
@@ Line 45: / Line 45: @@
 Using the [http://functions.wolfram.com/webMathematica/FunctionPlotting.jsp?name=RiemannSiegelZ online plotter] we can plot Z in the regions corresponding to scale divisions, using the conversion factor {{nowrap|''t'' {{=}} {{sfrac|2π|ln(2)}}''x''}}, for ''x'' a number near or at an edo number. Hence, for instance, to plot 12 plot around 108.777, to plot 31 plot around 281.006, and so forth. An alternative plotter is the applet [http://web.viu.ca/pughg/RiemannZeta/RiemannZetaLong.html here].
-If you have access to {{w|Mathematica}}, which has Z, zeta and theta as a part of its suite of initially defined functions, you can do even better. Below is a Mathematica-generated plot of Z({{frac|2π''x''|ln(2)}}) in the region around 12edo:
+If you have access to {{w|Mathematica}}, which has Z, zeta and theta as a part of its suite of initially defined functions, you can do even better. Below is a Mathematica-generated plot of Z({{frac|2π''x''|ln(2)}}) in the region around [[12edo]]:
 [[File:plot12.png|alt=plot12.png|plot12.png]]
@@ Line 51: / Line 51: @@
 The peak around 12 is both higher and wider than the local maximums above 11 and 13, indicating its superiority as an approximation of JI. Note also that the peak occurs at a point slightly larger than 12; this indicates the octave is slightly compressed in the zeta tuning for 12. The size of a step in octaves is 1/''x'', and hence the size of the octave in the zeta peak value tuning for ''N''edo is ''N''/''x''; if ''x'' is slightly larger than ''N'' as here with {{nowrap|''N'' {{=}} 12}}, the size of the zeta tuned octave will be slightly less than a pure octave. Similarly, when the peak occurs with ''x'' less than ''N'', we have stretched octaves.
-For larger edos, the width of the peak narrows, but for strong edos the height more than compensates, measured in terms of the area under the peak (the absolute value of the integral of Z between two zeros.) Note how 270 completely dominates its neighbors:
+For larger edos, the width of the peak narrows, but for strong edos the height more than compensates, measured in terms of the area under the peak (the absolute value of the integral of Z between two zeros.) Note how [[270edo]] completely dominates its neighbors:
 [[File:plot270.png|alt=plot270.png|plot270.png]]