Pr=1300; x0=3.0;xr=20;x1=25; %datos de la curva de potencia k=2.3849; c=6.0208; %parametros weibul calculados anteriormente potencia=xlsread('sotavento_curva potencia.xlsx','B2:B27'); x=0:1:25; pot=potencia(x>=x0 & x<=xr); hold on x=x0:1:xr; plot(x,pot,'ro','markersize',3,'markerfacecolor','r') title('Ajuste de la curva de potencia de un generador'); axis([0 15 -10 1550]) xlabel('velocidad') ylabel('potencia') grid on %ajuste p=polyfit(x,pot',3); %ajuste a un polinomio cúbico yp=polyval(p,x); plot(x,yp,'k') hold off
