Momentum and heat friction forces between fast ions (isotropic to the lowest order in the fast-ion rest frame) and thermal bulk plasma species (near Maxwellian) have been derived analytically for the first time using the linearized Fokker-Planck collision operator. The effects of flow-type distortions of the test particle distribution and flow-type restoring effects from field particles are both retained self-consistently. Here, the momentum and heat friction forces are defined as the momentum moments of the Fokker-Planck Coulomb collision operator weighted by the Laguerre polynomials L0(3/2)(wα/ν2Tα) and L1(3/2)(wα2/νTα2), respectively, i.e. the momentum friction force is F0α/β ≡ ∫ mαwαL0(3/2)Cα β[fα,fβ]d3wα , while the heat friction force is F1αβ ≡ ∫ mαwαL1(3/2)Cαβ[fα,fβ]d3wα. Here wα identical to ν - Vα represents the random velocity measured in the moving reference frame of the test particles, and the superscript α/β denotes test particles of species or colliding with field particles of species β. The momentum and heat friction forces have been explicitly calculated for electron-fast-ion collisions (Fke/f), fast-ion-electron collisions(Fkf/e), thermal-ion-fast-ion collisions (Fki/f) and fast-ion-thermal-ion collisions (Fkf/i), wherein the subscripts k correspond to the indices of Laguerre polynomials (k=0, 1). The collisional moment matrix is found to be asymmetric in the presence of fast ions-because the approximate self-adjointness of the Coulomb collision operator is not valid for this situation