1.4.2 C語言代碼

#include 
#include 
#include 
#define  N  7                      //定義量化論域模糊子集的個數
//模糊語言定義
#define  NB  -3
#define  NM  -2
#define  NS  -1
#define  ZO  0
#define  PS  1
#define  PM  2
#define  PB  3
struct
{
  //規則表
  int rule[ N ][ N ] ;                //模糊規則表
  //隸屬度函數
  char *f_type_e ;                  //e的隸屬度函數類型
  char *f_type_de ;                  //de的隸屬度函數類型
  char *f_type_u ;                  //u的隸屬度函數類型
  float paras_e[ 3*N ] ;                //e的隸屬度函數的參數
  float paras_de[ 3*N ] ;              //de的隸屬度函數的參數
  float paras_u[ 3*N ] ;                //u的隸屬度函數的參數
  //模糊運算參數
  float Ke ;                    //Ke=n/emax,量化論域為[-3,-2,-1,0,1,2,3]
  float Kde ;                    //Ke=n/demax,量化論域為[-3,-2,-1,0,1,2,3]
  float Ku ;                    //Ke=umax/n,量￠化論域為[-3,-2,-1,0,1,2,3]
  float umax ;                    //輸出的上限
  float emax ;                    //誤差基本論域上限
  float demax ;                    //誤差變化率基本論域的上限
  //控制目標參數
  float target ;                  //控制目標
  float actual ;                  //實際值
  float e ;                    //誤差
  float e_pre ;                    //上一次的誤差
  float de ;                    //誤差的變化率
}_fuzzy_;
/*******************************************************
Name      :trimf
Function    :三角隸屬度函數
Parameter  :
      x:自變量
      (a,b,c):定義域
Return    :結果
*******************************************************/
float trimf( float x, float a, float b, float c )
{
  float u ;
  if( ( x>=a )&&( x<=b ) )
    u = ( x-a )/( b-a ) ;
  else if( ( x>b )&&( x<=c ) )
    u = ( c-x )/( c-b ) ;
  else
    u = 0 ;
  return u ;
}
/*******************************************************
Name      :gaussmf
Function    :正態隸屬度函數
Parameter  :
      x:自變量
      ave:均值
      sigma:方差
Return    :結果
*******************************************************/
float gaussmf( float x, float ave, float sigma )
{
  float u ;
  if( sigma<0 )
    printf("方差不允許小于0\\n");
  else
    u = exp( -pow( ( ( x-ave )/sigma ), 2 ) ) ;
  return u ;
}
/*******************************************************
Name      :trapmf
Function    :梯形隸屬度函數
Parameter  :
      x:自變量
      (a,b,c,d):定義域
Return    :結果
*******************************************************/
float trapmf( float x, float a, float b, float c, float d )
{
  float u ;
  if( ( x>=a )&&( xelse if( ( x>=b )&&( x

**
最終代碼運行結果如下圖所示。

參考于CDSN

**