Abstract
In this research, we consider the inverse problem for the wave equation under an unknown initial condition. A generalized solution to the direct problem was formulated, its correctness was established, and the stability assessment was obtained. The inverse problem was reduced to an optimization problem, where the objective function was minimized using gradient methods, including the accelerated Nesterov algorithm. The conjugate problem was constructed, and the functional gradient was computed, while the existence of the Frechet derivative was proved. For the first time, the quaternion Fourier transform (QFT) was applied to the numerical solution of a direct problem, making it possible to analyze multidimensional wave processes more efficiently. A computational experiment was carried out, which demonstrated that if there is insufficient additional information, the restoration of the initial condition is incomplete. The introduction of the second boundary condition makes it possible to significantly improve the accuracy and stability of the solution. The results confirm the importance of an integrated approach and the availability of sufficient a priori information when solving inverse problems.