Matlab code for approximating difference in heat released of a disel engine given by exprimental data using the sum of two wiebe functions which parameters are deterimend by gradinet descent.
The approximated fuction is given by the following formula:
%Short version
data = wiebe(dQ)
%Long version
[data, Qtot] = wiebe(dQ, b, deltat)
- dQ - row vector that represents experimantal data of the difference in heat released of a disel engine as data points
- b - maximum heat released - default: 6.908
- deltat - time step used for sampling of experimental data - default: 0.1
- data - row vector containting 6 elements: a1, T1, r1, a2, T2, 3
- Qtot - total heat relesed
Plot of sample input (dQ)