Relation between Fourier phase and Analytic Signal Phase
Discussing the relationship between the Fourier transform and the Hilbert transform, let us consider an example function. It can be any square integrable function; for readability of the text, we select a function with both transforms expressible via simple analytical functions.
Let a signal s be a sinc function of time, s(t) = sinc(t). The Fourier transform of the sinc function is a rect function
$$
\mathscr{F}\{s(t)\}(ω)={{1\over{\sqrt{2π}}}·rect\left({ω\over{2π}}\right)}
$$
which is a function of frequency: OP's FT_s is a function of frequency, or, for the discrete case, an array with indices representing frequency bins.
The Hilbert transform of the sinc function
$$
\mathscr{H}\{s(t)\}(t)={{1-cos(t)}\over{t}}
$$
is the function of time, because the Hilbert transform is not a transform of the kind that maps time domain functions to frequency domain functions -- as the Fourier transform does -- but &quot;a specific linear operator that takes a function, u(t) of a real variable and produces another function of a real variable H(u)(t)&quot; (citation from the Wikipedia article on the Hilbert transform).
Comparing these formulas, one sees that the question &quot;if there is a relation between the phase of the Fourier Transform and the phase of the Analytic signal?&quot; is rather pointless.
FT_s, as defined by OP, is an array of values against frequency bins, and atan(Imag(FT_s),Real(FT_s)) is also an array indexed by the same bin enumerables, corresponding to a continuous-time frequency. Therefore, atan(Imag(FT_s),Real(FT_s)) is an array of phases of the signal frequency components.
A_s is a function of time, and atan(Imag(A_s),Real(A_s)) is also the function of time. It is the instantaneous phase of an analytic signal, and it is not the phase of the Fourier transform frequency component, not in any way.
The OP-defined two continuous-variable &quot;phase functions&quot; and corresponding discrete &quot;phase arrays&quot;, phaseFourier and phaseAnalytic, are defined on different domains: phaseFourier is defined on a frequency domain, phaseAnalyticis defined on a time domain. Still, you may want to see the graphs of these function for a sinc(t) signal.
For s(t) = sinc(t), the Fourier transform (OP's FT_s) values is a real constant up to a frequency of 1/(2π) and a zero for higher frequencies. Therefore, the array of phases of frequency components (OP's phaseFourier) is zero up to 1/(2π) and undefined in the region where the Fourier transform amplitude is zero. For sinc(t), one can safely assign all phases of FT components to zero in the entire frequency domain. The FT_s graph is trivial, and shown in every textbook and course notes on signal processing.
The instantaneous phase (OP's phaseAnalytic) graph may be shown less frequently, so here you are:

Summing up, one may agree that, talking about the relationship between the Fourier transform and the Hilbert transform, more appropriate is to indicate that the Hilbert transform is a multiplier operator. The multiplier of H is
$$
σ_{\mathscr{H}}(ω)=-i·sgn(ω)
$$
(see https://en.wikipedia.org/wiki/Hilbert_transform#Relationship_with_the_Fourier_transform).