三参数weibull分布matlab程序

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

clear all;clc
% By assuming the sample data is under the weibull 3 distribution, to get the 3 parameters,
% a is postition parameter, b is scale parameter, c is shape parameter
sample=[31.8 34.4 32.6 32.7 32.4 32.8 31.3 32.2 34.2 31.2 29.7 30.2 30.9 33.6 31.2 32.2 30.6 32.1 33.9 33.5]

S=sort(sample);num=length(sample);
pwm(1)=mean(S);

for j=1:num-1
b=0;
for i=1:num
if num-i >= j
a=prod(1:num-i)/prod(1:j)/prod(1:num-i-j);
else
a=0;
end
b=b+a*S(i);
end
b=b/num;
b=b*prod(1:j)*prod(1:num-j-1)/prod(1:num-1);
pwm(j+1)=b;
end

a=4*(pwm(4)*pwm(1)-pwm(2)*pwm(2))/(4*pwm(4)+pwm(1)-4*pwm(2))
b=(pwm(1)-a)/gamma(log((pwm(1)-2*pwm(2))/(pwm(2)-2*pwm(4)))/log(2))
c=log(2)/log( ( pwm(1)-2*pwm(2) )/(2*(pwm(2)-2*pwm(4))) )