Contents
% Raquel Fernandez Gomez
EJERCICIO 2
Apartado d)
[NUMEROS,TEXTO,RESTO]=xlsread('sotaventogaliciaanual.xlsx'); k0=mean(NUMEROS); c0=std(NUMEROS)^2; BINS=0:25; y=histc(NUMEROS,BINS); frecuencias=y/sum(y); a0=[c0 k0]; x=0:25; x=x'; f=@(a,x) (a(1)/a(2))*((x/a(2)).^(a(1)-1)).*exp(-(x/a(2)).^a(1)); %distribución de weibull af=nlinfit(x,frecuencias,f,a0); %función de potencia [NUMEROS,TEXTO,RESTO]=xlsread('sotavento_curva potencia.xlsx'); velocidad=NUMEROS(:,1); pot=NUMEROS(:,2); P=polyfit(velocidad,pot,3); g=@(x) f(af,x).*polyval(P,x); Potencia_media=quad(g,0,10)
Error using ^ Inputs must be a scalar and a square matrix. To compute elementwise POWER, use POWER (.^) instead. Error in apartd (line 8) c0=std(NUMEROS)^2;