Which numerical integral equation is correct?
1 view (last 30 days)
I would have a question in connection with numerical integration in Matlab R2015b. I would like to integrate the ceilometer measured attenuated backscatter coefficient. The spatial range of backscatter coefficioent is 10 m. I should calculate the following equation:
where dz is the range (dz = 10) and β is the backscatter coefficient. The backscatter coefficient values are in range x time matrix. I tried to calculate the integral value in two different ways.
X = repmat(0,size(BS))
for j = 2:400
X_sum = sum(X);
The second is:
TR = cumtrapz(BS);
TR_sum = sum(TR);
The resulted X and X_sum are one order of magnitude larger than TR and TR_sum. Could someone write me which calculation is correct?
Thank you for your help in advance!