The aortic blood flow is described by a set of nonlinear hyperbolic partial differential equations that account for mass and momentum conservation, and nonlinear models for the mechanical properties of the artery. Identification is used for determining the wave speed, arterial taper, and cross section: these parameters reflect the elastic characteristics of the aorta wall and control the pulsatile response. The differential equations were numerically integrated by the Lax-Wendroff scheme of Abarbanel and Goldberg [ J. Comput. Phys. 10:1–21 (1972)] that avoids nonlinear oscillations. The Gauss-Newton technique was used for the parameter identification. By reference to reported elocity and pressure input-output pairs, a parameter vector is found such that the distance in the L 2 norm between the predicted outputs and the measured functions is minimal. Calculations of the velocity and pressure waves show excellent compatibility of the model with reported experimental data: starting from arbitrary parameter estimates, which yield grossly distorted waveforms, the error is typically reduced to 7–8%. Introduction of viscoelastic behaviour for the arterial wall in the form of a Volterra integral for the cross section does not lead to significant improvement. Numerical examples are presented which prove the convergence, accuracy, and stability of the algorithm. Emphasis is placed on the computational feasibility of the proposed system identification.