Matlab实现风电场风速两参数weibull(威布尔)分布

close all clc clear all d=[3.5,3.8,4.0,4.3,4.5,4.7,4.8,5.0,5.2,5.4,5.5,5.7,6.0,6.1,6.3,6.5,17,7.3,7.7,10,12] ; %试验数据 PDF=@(v,a,b) (a./b.^a).*v.^(a-1).*exp(-(v./b).^a); start=[0.19 3]; %两个初值 （初值要求很高） phat=mle(d,'pdf',PDF,'start',start); k=phat(1); c=phat(2); v=0:0.1:20; fv=(k/c)*(v/c).^(k-1).*exp(-(v/c).^k); plot(v,fv);%风速分布 xlabel('风速v(m/s)'); ylabel('概率密度f(v)'); %weibull分布参数估计及概率表达式求解完毕