fuzzy均值2.cpp

模糊规划

// fuzzy Simulation
// Written by Microsoft Visual C++
// Copyright by UTLab @ Tsinghua University
// http://orsc.edu.cn/UTLab
//E[(a8+a9+a4+a5+a6)*(a1+a7+a8+a9)*(a2+a3)]=?

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "UTLab.h"

#define N 3000

int main(void)
{
  double t, a, a1, a2, a3, a4, a5, a6, a7, a8, a9, x[N+1],y[N+1],z[N+1], mu[N+1];
  int    i;
  FILE   *fs;   
  fs   =   fopen("a.txt","w");   
  
 
  for(i=1; i<=N; i++)  {  
	  
	  a4=myu(0.01,0.03);
	  mu[i]=triangle(a4,0.01,0.02,0.03);
	 
	  
	  a5=myu(0.03,0.06);
	  t=triangle(a5,0.03,0.05,0.06);
	  
	  if(mu[i]>t) mu[i]=t;

	  a6=myu(0.04,0.06);
	  t=triangle(a6,0.04,0.05,0.06);
	 
	  if(mu[i]>t) mu[i]=t;
	  
	  a8=myu(0.01,0.03);
	  t=triangle(a8,0.01,0.02,0.03);
	  
	  	  if(mu[i]>t) mu[i]=t;
	  
	  a9=myu(0.03,0.06);
	  t=triangle(a9,0.03,0.05,0.06);
	 
	  if(mu[i]>t) mu[i]=t;
	  
	  a1=myu(0.01,0.03);
	  t=triangle(a1,0.01,0.02,0.03);
	 
	  if(mu[i]>t) mu[i]=t;

	  a7=myu(0.01,0.04);
	  t=triangle(a7,0.01,0.03,0.04);
	 
	  if(mu[i]>t) mu[i]=t;


	  a2=myu(0.02,0.05);
	  t=triangle(a2,0.02,0.04,0.05);
	  
	  if(mu[i]>t) mu[i]=t;
	  
	  a3=myu(0.02,0.05);
	  t=triangle(a3,0.02,0.04,0.05);
	 
	  if(mu[i]>t) mu[i]=t;

     
	x[i]=(1-(1-a8)*(1-a9)*(1-a4)*(1-a5)*(1-a6))*(1-(1-a1)*(1-a7)*(1-a8)*(1-a9))*(1-(1-a2)*(1-a3));
     fprintf(fs,"%f\n",mu[i]); 
  }

  a=fuzzymean(x, mu, N, N);

  printf("Expected value is %f\n",a);

  fclose(fs); 
  return 1;

}