A combustion speed, for example, is estimated or evaluated, with a required accuracy, more simply than the conventional art, while reducing man-hours to produce a heat generation rate waveform of an internal combustion engine. An increase rate of a heat generation rate relative to a change in a crank angle in a heat generation rate increasing period (e.g., a first-half combustion period a ) in which the heat generation rate increases after ignition of an air-fuel mixture is defined as a heat generation rate gradient b / a that is one of characteristic values of the heat generation rate waveform. The heat generation rate gradient is estimated based on a fuel density (e.g., fuel density Á fuel@dQpeak at heat generation rate maximum time) at a predetermined time set in advance in the heat generation rate increasing period so as to produce the heat generation rate waveform using the estimated heat generation rate gradient.
展开▼