- 相關(guān)推薦
模糊控制算法c程序
由于項(xiàng)目需要，需要模糊控制算法，之前此類知識(shí)為0，經(jīng)過半個(gè)多月的研究，終于有的小進(jìn)展。開始想從強(qiáng)大的互聯(lián)網(wǎng)上搜點(diǎn)c代碼來研究下，結(jié)果搜遍所有搜索引擎都搜不到，以下本人從修改的模糊控制代碼，經(jīng)過自己修改后可在vc6.0，運(yùn)行！輸入e表示輸出誤差，ec表示誤差變化率，經(jīng)過測試具有很好的控制效果，對(duì)于非線性系統(tǒng)和數(shù)學(xué)模型難以建立的系統(tǒng)來說有更好的控制效果！現(xiàn)將其公開供大家學(xué)習(xí)研究！
#include
#include
#define PMAX 100
#define PMIN -100
#define DMAX 100
#define DMIN -100
#define FMAX 100 /*語言值的滿幅值*/
int PFF[4]={0,12,24,48};
/*輸入量D語言值特征點(diǎn)*/
int DFF[4]={0,16,32,64};
/*輸出量U語言值特征點(diǎn)*/
int UFF[7]={0,15,30,45,60,75,90};
/*采用了調(diào)整因子的規(guī)則表,大誤差時(shí)偏重誤差,小誤差時(shí)偏重誤差變化*/ /*a0=0.3,a1=0.55,a2=0.74,a3=0.89 */
int rule[7][7]={
//誤差變化率 -3,-2,-1, 0, 1, 2, 3 // 誤差
{-6,-6,-6,-5,-5,-5,-4,}, // -3
{-5,-4,-4,-3,-2,-2,-1,}, // -2
{-4,-3,-2,-1, 0, 1, 2,}, // -1
{-4,-3,-1, 0, 1, 3, 4,}, // 0
{-2,-1, 0, 1, 2, 3, 4,}, // 1
{ 1, 2, 2, 3, 4, 4, 5,}, // 2
{ 4, 5, 5, 5, 6, 6, 6}}; // 3
/**********************************************************/ int Fuzzy(int P,int D) /*模糊運(yùn)算引擎*/
{
int U; /*偏差,偏差微分以及輸出值的精確量*/
unsigned int PF[2],DF[2],UF[4]; /*偏差,偏差微分以及輸出值的隸屬度*/
int Pn,Dn,Un[4];
long temp1,temp2;
/*隸屬度的確定*/
/*根據(jù)PD的指定語言值獲得有效隸屬度*/
if(P>-PFF[3] && P
{
if(P
{
Pn=-2;
PF[0]=FMAX*((float)(-PFF[2]-P)/(PFF[3]-PFF[2])); }
else if(P
{
Pn=-1;
PF[0]=FMAX*((float)(-PFF[1]-P)/(PFF[2]-PFF[1])); }
else if(P
{
Pn=0;
PF[0]=FMAX*((float)(-PFF[0]-P)/(PFF[1]-PFF[0])); }
else if(P
{
Pn=1; PF[0]=FMAX*((float)(PFF[1]-P)/(PFF[1]-PFF[0])); }
else if(P
{
Pn=2; PF[0]=FMAX*((float)(PFF[2]-P)/(PFF[2]-PFF[1])); }
else if(P
{
Pn=3; PF[0]=FMAX*((float)(PFF[3]-P)/(PFF[3]-PFF[2])); }
}
else if(P
{
Pn=-2; PF[0]=FMAX;
}
else if(P>=PFF[3])
{
Pn=3; PF[0]=0;
}
PF[1]=FMAX-PF[0];
if(D>-DFF[3] && D
{
if(D
{
Dn=-2;DF[0]=FMAX*((float)(-DFF[2]-D)/(DFF[3]-DFF[2])); }
else if(D
{
Dn=-1;
DF[0]=FMAX*((float)(-DFF[1]-D)/(DFF[2]-DFF[1])); }
else if(D
{
Dn=0;
DF[0]=FMAX*((float)(-DFF[0]-D)/(DFF[1]-DFF[0])); }
else if(D
{
Dn=1;
DF[0]=FMAX*((float)(DFF[1]-D)/(DFF[1]-DFF[0])); }
else if(D
{
Dn=2; DF[0]=FMAX*((float)(DFF[2]-D)/(DFF[2]-DFF[1])); }
else if(D
{
Dn=3; DF[0]=FMAX*((float)(DFF[3]-D)/(DFF[3]-DFF[2])); }
}
else if(D
{
Dn=-2;
DF[0]=FMAX;
}
else if(D>=DFF[3])
{
Dn=3;
DF[0]=0;
}
DF[1]=FMAX-DF[0];
/*使用誤差范圍優(yōu)化后的規(guī)則表rule[7][7]*/
/*輸出值使用13個(gè)隸屬函數(shù),中心值由UFF[7]指定*/ /*一般都是四個(gè)規(guī)則有效*/
Un[0]=rule[Pn-1+3][Dn-1+3];
Un[1]=rule[Pn+3][Dn-1+3];
Un[2]=rule[Pn-1+3][Dn+3];
Un[3]=rule[Pn+3][Dn+3];
if(PF[0]
UF[0]=PF[0];
else
UF[0]=DF[0];
if(PF[1]
UF[1]=PF[1];
else
UF[1]=DF[0];
if(PF[0]
UF[2]=PF[0];
else
UF[2]=DF[1];
if(PF[1]
UF[3]=PF[1];
else
UF[3]=DF[1];
/*同隸屬函數(shù)輸出語言值求大*/
if(Un[0]==Un[1])
{
if(UF[0]>UF[1])
UF[1]=0;
else
UF[0]=0;
}
if(Un[0]==Un[2])
{
if(UF[0]>UF[2])
UF[2]=0;
else
UF[0]=0;
}
if(Un[0]==Un[3])
{
if(UF[0]>UF[3])
UF[3]=0;
else
UF[0]=0;
}
if(Un[1]==Un[2])
{
if(UF[1]>UF[2])
UF[2]=0;
else
UF[1]=0;
}
if(Un[1]==Un[3])
{
if(UF[1]>UF[3])
UF[3]=0;
else
UF[1]=0;
}
if(Un[2]==Un[3])
{
if(UF[2]>UF[3])
UF[3]=0;
else
UF[2]=0;
}
/*重心法反模糊*/
/*Un[]原值為輸出隸屬函數(shù)標(biāo)號(hào)，轉(zhuǎn)換為隸屬函數(shù)值*/ if(Un[0]>=0)
Un[0]=UFF[Un[0]];
else
Un[0]=-UFF[-Un[0]];
if(Un[1]>=0)
Un[1]=UFF[Un[1]];
else
Un[1]=-UFF[-Un[1]];
if(Un[2]>=0)
Un[2]=UFF[Un[2]];
else
Un[2]=-UFF[-Un[2]];
if(Un[3]>=0)
Un[3]=UFF[Un[3]];
else
Un[3]=-UFF[-Un[3]];
temp1=UF[0]*Un[0]+UF[1]*Un[1]+UF[2]*Un[2]+UF[3]*Un[3]; temp2=UF[0]+UF[1]+UF[2]+UF[3];
U=temp1/temp2;
return U;
}
void main()
{
int a=0,e,ec;
/*int nowpoint,p1,p2=1;
FILE *in,*out;
in=fopen(
out=fopen(
while(1)
{
//fscanf(in,
//e=0-nowpoint;
//ec= p1-p2;
printf(
scanf(
printf(
scanf(
a=Fuzzy(e,ec);
//fprintf(out,
//printf(
printf(
//p2=p1;
}
//fclose(in);
//fclose(out);
}
【模糊控制算法c程序】相關(guān)文章：
基于粒子群算法的污水處理模糊控制研究04-28
基于蟻群算法的結(jié)構(gòu)面產(chǎn)狀模糊C均值聚類分析04-26
基于遺傳算法的航空發(fā)動(dòng)機(jī)模糊優(yōu)化控制05-01
發(fā)動(dòng)機(jī)試驗(yàn)臺(tái)測控系統(tǒng)的模糊控制算法研究04-28
基于修正因子智能權(quán)函數(shù)的汽車ABS模糊控制算法仿真研究04-27
基于模糊理論的漸消濾波算法05-01
擾動(dòng)值模糊推理的三I算法04-29
算法與程序設(shè)計(jì)的教案04-28