Ejercicio 2
Contents
Apartado b)
velocidad=xlsread('sotaventogalicianual.xlsx'); x=0.5:1:max(velocidad); promedio=mean(velocidad); destd=std(velocidad); interv=hist(velocidad,x); frec=interv/sum(interv); f=@(a,x) (a(1)/a(2))*((x/a(2)).^(a(1)-1)).*exp(-(x/a(2)).^a(1)); % f=@(k,x) (k/c)*((x/c).^(k-1)).*exp(-(x/c).^k); a0=[promedio destd^2]; af=nlinfit(x,frec,f,a0) hold on bar(x,frec,'c') x=linspace(0,max(velocidad),100); y=f(af,x); plot(x,y,'r')
af =
2.3849 6.0208
