The graphical comparison should use pdf normalization of the histogram.
Original code, but use use an actual Weibull distribution for the data
rng default
%A=rnd(6000,1);
A = wblrnd(4,3,6000,1);
[Overal_parameters]= wblfit(A)
figure()
histogram(A,'Normalization','probability')
hold on
plot(sort(A),wblpdf(sort(A),Overal_parameters(1),Overal_parameters(2)))
legend('Observed Samples','Estimated Distribution')
xlabel('Mean wind speed (m/s)')
ylabel('Annual probabilty density')
hold off
Use pdf normalization to compare to the estimated pdf
figure()
histogram(A,'Normalization','pdf')
hold on
plot(sort(A),wblpdf(sort(A),Overal_parameters(1),Overal_parameters(2)))
legend('Observed Samples','Estimated Distribution')
xlabel('Mean wind speed (m/s)')
ylabel('Annual probabilty density')
hold off
I'm not sure what the goal is from here, but at least now the graphs comparable.
