function y=weibuer(a,b,n)
%韋布爾分布,a,b表示參數,b不能為0.n表示數據量,一般要大于1024
%a=1時,是指數分布
%a=2時,是瑞利分布
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x=junyun(n);
for i=1:n
    if (x(i)==0)
        x(i)=0.0001;
    else
        continue;
    end
end
u=-log(x);
y=b*u.^(1/a);