%apardado b
dw=xlsread('sotaventogaliciaanual.xlsx');
%interpolar si fuese necesario
if any(isnan(dw)) %si hay algún NaN
    x=1:length(dw);
    i=find(~isnan(dw));
    velocidad=interp1(x(i),dw(i),x);
end

%histo
x=0.5:1:max(dw);
horas=hist(dw,x);

%convierte a frecuencias y ajusta a la función de Weibull
frec=horas/sum(horas);
f=@(a,x) (a(1)/a(2))*((x/a(2)).^(a(1)-1)).*exp(-(x/a(2)).^a(1));
a0=[2 8];  %valor inicial parámetros
af=nlinfit(x,frec,f,a0);
hold on

%diagrama de frecuencias
bar(x,frec,'c')

%para representar la curva de ajuste
x=linspace(0,max(dw),100);
y=f(af,x);
plot(x,y,'r')

title('Ajuste a la función Weibull')
xlabel('Vel')
ylabel('Frec')
hold off

Published with MATLAB® R2019b