MATLAB學(xué)習(xí)筆記之模糊算法3

1.4.2 C語(yǔ)言代碼

#include 
#include 
#include 
#define  N  7                      //定義量化論域模糊子集的個(gè)數(shù)
//模糊語(yǔ)言定義
#define  NB  -3
#define  NM  -2
#define  NS  -1
#define  ZO  0
#define  PS  1
#define  PM  2
#define  PB  3
struct
{
  //規(guī)則表
  int rule[ N ][ N ] ;                //模糊規(guī)則表
  //隸屬度函數(shù)
  char *f_type_e ;                  //e的隸屬度函數(shù)類型
  char *f_type_de ;                  //de的隸屬度函數(shù)類型
  char *f_type_u ;                  //u的隸屬度函數(shù)類型
  float paras_e[ 3*N ] ;                //e的隸屬度函數(shù)的參數(shù)
  float paras_de[ 3*N ] ;              //de的隸屬度函數(shù)的參數(shù)
  float paras_u[ 3*N ] ;                //u的隸屬度函數(shù)的參數(shù)
  //模糊運(yùn)算參數(shù)
  float Ke ;                    //Ke=n/emax,量化論域?yàn)閇-3,-2,-1,0,1,2,3]
  float Kde ;                    //Ke=n/demax,量化論域?yàn)閇-3,-2,-1,0,1,2,3]
  float Ku ;                    //Ke=umax/n,量￠化論域?yàn)閇-3,-2,-1,0,1,2,3]
  float umax ;                    //輸出的上限
  float emax ;                    //誤差基本論域上限
  float demax ;                    //誤差變化率基本論域的上限
  //控制目標(biāo)參數(shù)
  float target ;                  //控制目標(biāo)
  float actual ;                  //實(shí)際值
  float e ;                    //誤差
  float e_pre ;                    //上一次的誤差
  float de ;                    //誤差的變化率
}_fuzzy_;
/*******************************************************
Name      :trimf
Function    :三角隸屬度函數(shù)
Parameter  :
      x:自變量
      (a,b,c):定義域
Return    :結(jié)果
*******************************************************/
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    :正態(tài)隸屬度函數(shù)
Parameter  :
      x:自變量
      ave:均值
      sigma:方差
Return    :結(jié)果
*******************************************************/
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    :梯形隸屬度函數(shù)
Parameter  :
      x:自變量
      (a,b,c,d):定義域
Return    :結(jié)果
*******************************************************/
float trapmf( float x, float a, float b, float c, float d )
{
  float u ;
  if( ( x>=a )&&( xelse if( ( x>=b )&&( x

**
最終代碼運(yùn)行結(jié)果如下圖所示。

參考于CDSN

**