MATLAB绘制威布尔分布曲线

MATLAB 绘制威布尔分布曲线

威布尔分布概率密度函数：

1(/)(,,)()a a x m a x f x m a e m m

--=

威布尔分布概率分布函数： ()()1a

mx F x e -=-

其中m>0,是尺度参数也叫比例参数，a>0是形状参数。

X 是随机变量，是未知参数，表示时间延滞。

图1：设定尺度参数m 值为1，取五个形状参数a ，自变量x

代码如下：

m=[1 1 1 1 1,2];

a=[0.5 1 1.5 2.5 5,5];

x=linspace(0,5);

linecolor=['r','b','g','k','y'];

for n=1:5

y1=m(n)*a(n)*((m(n)*x).^(a(n)-1)).*(exp(-(m(n)*x).^a(n)));

y=1-exp(-(m(n)*x).^a(n));

subplot(1,2,2)

title('图1：概率分布函数');

plot(x,y);

hold on;

subplot(1,2,1)

type=linecolor(n);

title('图1：概率密度函数');

plot(x,y1,type);

hold on;

legend('m=1,a=0.5','m=1,a=1','m=1,a=1.5','m=1,a=2.5','m=1,a=5'); end

图2：设定形状参数a值为2，取五个尺度参数m，自变量x

代码如下：

m=[0.5 0.75 1 1.5 1.75,2];

a=[2 2 2 2 2.5];

x=linspace(0,5);

linecolor=['r','y','b','g','k'];

for n=1:5

y1=m(n)*a(n)*((m(n)*x).^(a(n)-1)).*(exp(-(m(n)*x).^a(n)));

y=1-exp(-(m(n)*x).^a(n));

subplot(1,2,2)

title('图2：概率分布函数');

plot(x,y);

hold on;

subplot(1,2,1)

type=linecolor(n);

title('图2：概率密度函数');

plot(x,y1,type);

hold on;

legend('m=0.5,a=2','m=0.75,a=2','m=1,a=2','m=1.5,a=2','m=1.75,a=2'); end

图3：设定尺度参数m值为1，自变量为x，a的三维概率分布图代码如下：

m=1;

[x,a]=meshgrid(0:0.05:4,0:0.05:5);

fx=m.*a.*(m.*x).^(a-1).*(exp(-(m.*x).^a));

Fx=1-exp(-(m.*x).^a);

subplot(1,2,1)

mesh(x,a,fx);

title('图3：m=1,a,x三维概率密度分布');

subplot(1,2,2)

mesh(x,a,Fx);

title('图3：m=1,a,x三维概率分布图');

图4：设定形状参数a值为2，自变量为x，m的三维概率分布图代码如下：

a=2;

[x,m]=meshgrid(0:0.05:5,0:0.05:2);

fx=m.*a.*(m.*x).^(a-1).*(exp(-(m.*x).^a)); Fx=1-exp(-(m.*x).^a);

subplot(1,2,1)

mesh(x,m,fx);

title('图4：a=2,m,三维概率密度分布'); subplot(1,2,2)

mesh(x,m,Fx);

title('图4：a=2,m,x三维概率分布图');