matlab绘制威布尔分布曲线

MATLAB绘制威布尔分布曲线 威布尔分布概率密度函数： 威布尔分布概率分布函数：

其中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: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);?????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三维概率??布图');