The viscoelastic properties of tissues are often assessed by analyzing the phase velocity dispersion curve computed using a 2-D Fourier transform onthe single fixed duration acoustic radiation force excitation pulse (FD-ARFEP)-induced shear wave with broadband frequencies. This study proposes a multi-frequency ARFEP (MF-ARFEP) using a single imaging transducer to generate shear waves at target frequencies of 100:100:1000 Hz to reduce noise in the phase velocity images. The MF-ARFEP is generated by summing sinusoids with target frequencies and then sampled to collect tracking pulses in-between the sampled MF-ARFEP. The MF-ARFEP is validated by imaging 6 and 70 kPa inclusions in a 18 kPa background, exvivo canine liver, and in vivo 4T1 breast cancer mouse tumor with comparison to the FD-ARFEP. The median phase velocity was similar between MF-ARFEP versus FD-ARFEP. However, the CNR was on an average 1.3 times higher in MF-ARFEP versus FD-ARFEP -derived images of inclusions and the coefficient of variation was on an average 1.1 and 2.2 times lower in MF-ARFEP versus FD-ARFEP-derived images of tumor and liver. These results demonstrate the feasibility of the MF-ARFEP excitation pulse to generate phase velocity images at 100:100:1000 Hz frequencies with reduced noise compared to the FD-ARFEP. Ongoing studies entail the translation of this new viscoelasticity imaging method for the characterization of breast cancer tumors in patients.