How to realize Riemann Siegel Theta Function with MATLAB?
5 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
yifei wang
am 1 Jun. 2020
Kommentiert: 兰花 西
am 24 Aug. 2020
How to realize Riemann Siegel Theta Function with MATLAB?
0 Kommentare
Akzeptierte Antwort
John D'Errico
am 1 Jun. 2020
Bearbeitet: John D'Errico
am 2 Jun. 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 Kommentare
兰花 西
am 24 Aug. 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 !
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Startup and Shutdown finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!