How to realize Riemann Siegel Theta Function with MATLAB?

1 visualización (últimos 30 días)
How to realize Riemann Siegel Theta Function with MATLAB?

Respuesta aceptada

John D'Errico
John D'Errico el 1 de Jun. de 2020
Editada: John D'Errico el 2 de Jun. de 2020
You start by writing the formula for it. That can be found on virtually the first line in this link.
However in MATLAB, the gamma function seems to be defined only for a real variable, unless you use the symbolic toolbox version of gamma. So this should work:
t = (1:5)';
[t,angle(double(gamma(sym(1/4 + i*t/2)))) - log(pi)/2*t]
ans =
1 -1.76754795281229
2 -2.52591091881613
3 -2.99456469601083
4 -3.29063503121648
5 -3.45962037536346
That is valid for real arguments t. I'd need to do some reading to know if this works for complex arguments, though I assume it does not, since the wiki link I gave indicates it is valid for real arguments. Hopefully you care only about real arguments, and are uninterested in an extension to the complex plane.
  4 comentarios
yifei wang
yifei wang el 2 de Jun. de 2020
I got it. You are a great help. Thank you.
兰花 西
兰花 西 el 24 de Ag. de 2020
Excuse me, when I set the parameter to a relatively large imaginary number(-1.93097838020833 + 38142819340462.3 i) in the above formula, I cannot get a valid result. Do you know how to solve it? Thank you !

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Startup and Shutdown en Help Center y File Exchange.

Productos


Versión

R2017b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by