In an attempt to simulate guitar sounds, a discrete model was developed that consisted of real strings with stiffness and internal damping, the player's finger, the resonance box, and the modelization of the sound pressure. The vibrating string equation was solved in the time domain using the finite‐difference method. The player's figure was modeled as a force density term in the string equation, which accounted for the initial plucking conditions. The resonance box was modeled as impedancelike boundary conditions at one end of the string. These end conditions were obtained under the assumption that the driving force exerted by the string at its ends excites a set of second‐order resonances in parallel, the first resonance of the plate being coupled with the air resonance of the box. Finally, the sound pressure was modeled as the contribution of a set of monopoles, where each monopole corresponds to one given resonance. The radiation efficiency of each monopole was derived from a curve fitting on a guitar sound‐pressure spectrum measured in an anechoic chamber. The waveforms obtained are in good agreement with experimental waveforms measured on real guitars. Sound examples will be presented high‐lighting the realism of the simulation.