matlab图像区域填充的原理_Matlab函数区域填充fill函数

x=0.2;p=130;H=150;L=100;s=30;k=0.5;r=0.3;u=0.7;t0=(p-L)/p;

GH=normcdf((p-L*x).*((1-x).^(-1)),H,s)

GL=normcdf((p-L*x).*((1-x).^(-1)),L,s)

GHP=normcdf(p,H,s)

GLP=normcdf(p,L,s)

GL1=normcdf(L,L,s)

GH1=normcdf(L,H,s)

t=t0:0.01:1

z1=(L-(p-(p-L)./t)).*(1-k+k*r)+u.*(p-(p-(p-L)./t)).*(r.*(1-GLP)+(1-k).*(1-r).*(1-GHP))

y1=(L-(p-(p-L)./t)).*(k*r+1-k)+u*(p-(p-(p-L)./t)).*(r*(1-GL)+(1-r)*(1-GH)-k*(1-r).*(1-GHP))+u.*(L-(p-(p-L)./t)).*(r.*(GL-GL1)+(1-r).*(GH-GH1))

plot(t,z1,'-',t,y1,':')

legend('one-stage GB mechanism','both-stage GB

mechanism')

xlabel('(p-\theta_1)/(p-c)'),ylabel('change in

profit')

title('the comparison between one- and both-stage GB

mechanism')

hold on

% 只涂从0开始到交点位置​

t_ = t0:0.01:0.7608;

z1_=(L-(p-(p-L)./t_)).*(1-k+k*r)+u.*(p-(p-(p-L)./t_)).*(r.*(1-GLP)+(1-k).*(1-r).*(1-GHP));

y1_=(L-(p-(p-L)./t_)).*(k*r+1-k)+u*(p-(p-(p-L)./t_)).*(r*(1-GL)+(1-r)*(1-GH)-k*(1-r).*(1-GHP))+u.*(L-(p-(p-L)./t_)).*(r.*(GL-GL1)+(1-r).*(GH-GH1));

fill([t_,fliplr(t_)], [y1_,fliplr(z1_)],'r')