Which numerical integral equation is correct?
1 次查看(过去 30 天)
显示 更早的评论
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:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/225811/image.png)
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.
First is:
X = repmat(0,size(BS))
for j = 2:400
X(j,1)=10.*(BS(j-1,1)+BS(j,1)).*0.5+X(j-1,1);
end
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!
0 个评论
回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!