An algorithm for quantitative description of cardiac cycle dependent variation of integrated backscatter (cyclic variation) has been developed and is shown to be suitable for analysis of nonsinusoidal data typical of ultrasonic tissue characterization measurements from myocardium in vivo. The algorithm produces estimates of the magnitude of variation and of the time delay relative to the electrocardiographically recorded QRS-complex. To validate the algorithm, 246 integrated backscatter measurements were analyzed both manually and by the automated method. The magnitude and time delay estimates from the two methods correlated closely. With a separate set of data, the algorithm produced reasonable descriptions of the cyclic variation for 89 of 101 integrated backscatter measurements. Only modest computational power is required for effective implementation of this algorithm, facilitating inclusion of online automated analysis capabilities in quantitative ultrasonic tissue characterization systems.