In this article, a novel wavelet collocation method based on Fibonacci wavelets is proposed to solve the dual-phase-lag (DPL) bioheat transfer model in multilayer skin tissues during hyperthermia treatment. Firstly, the Fibonacci polynomials and the corresponding wavelets along with their fundamental properties are briefly studied. Secondly, the operational matrices of integration for the Fibonacci wavelets are built by following the celebrated approach of Chen and Haiso. Thirdly, the proposed method is utilized to reduce the underlying DPL model into a system of algebraic equations, which has been solved using the Newton iteration method. Towards the culmination, the effect of different parameters including the tissue-wall temperature, time-lag due to heat flux, time-lag due to temperature gradient, blood perfusion, metabolic heat generation, heat loss due to diffusion of water, and boundary conditions of various kinds on multilayer skin tissues during hyperthermia treatment are briefly presented and all the outcomes are portrayed graphically.