MATLAB绘制威布尔分布曲线

MATLAB绘制威布尔分布曲线威布尔分布概率密度函数：1(/)(,,)()aaxmaxfxmaemm威布尔分布概率分布函数：()1amxFxe其中m0,是尺度参数也叫比例参数，a0是形状参数。X是随机变量，是未知参数，表示时间延滞。图1：设定尺度参数m值为1，取五个形状参数a，自变量x代码如下：m=[11111,2];a=[0.511.52.55,5];x=linspace(0,5);linecolor=['r','b','g','k','y'];forn=1:5y1=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);holdon;subplot(1,2,1)type=linecolor(n);title('图1：概率密度函数');plot(x,y1,type);holdon;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.50.7511.51.75,2];a=[22222.5];x=linspace(0,5);linecolor=['r','y','b','g','k'];forn=1:5y1=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);holdon;subplot(1,2,1)type=linecolor(n);title('图2：概率密度函数');plot(x,y1,type);holdon;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三维概率分布图');