czfee.m

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

function z=Czfee(x,d)         %计算总费用
%输入:x 一组标定值组合,d:一组容差组合
%输出:该组合下的损失费用
dh=0.001;
for i=1:7        %分别求y对x1...x7的偏导数
   for j=1:7
      x1(:,j)=x(j)*ones(2,1);
   end
   x1(:,i)=linspace(x(i),x(i)+dh,2)';
   F(:,i)=FUN(x1(:,1),x1(:,2),x1(:,3),x1(:,4),x1(:,5),x1(:,6),x1(:,7));
   g(i)=diff(F(:,i),1)/dh;%求数值导数。
end

mu=F(1,1);%???
sigma=sqrt((g.^2)*((x.*d).^2)');%????
%normcdf(n,mu,sigma)是求y小于n的概率,其中y服从正态分布
y1=normcdf(1.4,mu,sigma)-normcdf(1.2,mu,sigma)+normcdf(1.8,mu,sigma)-normcdf(1.6,mu,sigma);%零件质量损失的概率
y2=1-normcdf(1.8,mu,sigma)+normcdf(1.2,mu,sigma);%零件成废品的概率

e=ones(1,7);%e(i)=1表示等级A;e(i)=2表示等级B;e(i)=3表示等级C
for i=1:7
if d(i)==0.1/3 e(i)=1;
elseif d(i)==0.05/3 e(i)=2;
    else e(i)=3;
    end
end
f=[0,25,0;20,50,0;20,50,200;50,100,500;50,0,0;10,25,100;0,25,100];
for i=1:7
    p(i)=f(i,e(i));
end
z=1000*(1000*y1+9000*y2)+1000*(p*ones(7,1));% 'p*ones(7,1)'  has the same meaning with  'sum(p)'