.png)
I have no idea what you are doing.
Try this:
filename = 'solar_spectrum.xlsx';
[D,S,R] = xlsread(filename);
Phw = D(:,1);
hw = D(:,2);
E_g = hw >= 0.5; % Logical Index — Choose ‘E_g’ (Here: 0.5)
S = trapz(hw(E_g), Phw(E_g))
figure
plot(hw, Phw)
hold on
patch([hw(E_g) fliplr(hw(E_g))], [Phw(E_g) fliplr(Phw(E_g))], 'y')
hold off
grid
producing:
S =
2.0305
The plot shows the approximate area covered by the integration.
Experiment to get the result you want.
