mubiaof.m

#include "stdio.h" #include "math.h" #include "malloc.h" #include "string.h" #define m 6 #defin

function Z=mubiaof(x,mu,sigma)
C=[1.7 2.25 4.5 1.12];%托运单价
v=[0.0012 0.0015 0.003 0.0008];%单位质量的体积
for i=1:4
    a=normcdf(x(i+4),mu(i),sigma(i))-normcdf(0,mu(i),sigma(i));
    a=1./(sqrt(2.*pi).*sigma(i)).*a.*C(i).*x(i);
%     b0=1./(sqrt(2.*pi).*sigma(i));
%     b1=sigma(i).^2;
%     b2=exp(-mu(i).^2./(2.*b1));
%     b3=exp(-(x(i+4)-mu(i)).^2./(2.*b1));
%     b4=normcdf((x(i+4)-mu(i))./sigma(i),0,1)-normcdf(-mu(i)./sigma(i),0,1);
%     a=b0.*b1.*(b2-b3)+sigma(i).*mu(i).*b4;
    tmp(i)=(1./(sqrt(2.*pi).*sigma(i))).*(normcdf(24000,mu(i),sigma(i))-normcdf(x(i+4),mu(i),sigma(i))).*C(i).*x(i+4)+a;
end
% s=x(9:12);
% m=[x(13:16);x(17:20);x(21:24)];
% if m*v'<=9.084
%     if sum(m)<= s
%         if sum(m')<=8000 
             Z=sum(tmp);%满足约束条件才可输出
%              return;
%          end
%      end
%  end
% end
% Z=0;%不满足约束条件则置零