computeFloquetHamiltonian

computeFloquetHamiltonian(H, T, nt)

Compute the (effective) Floquet Hamiltonian from a time-periodic Hamiltonian.

Parameters:

• H (function_handle) – Matrix-valued Hamiltonian function handle \(H(t)\)

• T (double) – Modulation period \(T\) in [s]

• nt (double, optional) – Number of time samples over one period (uniform grid)

Returns:

Effective Floquet Hamiltonian \(H_F\) in [Hz]

Return type:

double matrix

Notes:

The evolution operator over one period is approximated by a product \(U(T) \approx \prod_k \exp\{-i\,2\pi\, H(t_k)\, \Delta t\}\), and the effective Hamiltonian is extracted as \(H_F = \frac{i}{2\pi T}\,\log U(T)\).