C实现模糊控制算法

1、由于项目需要，需要模糊控制算法，之前此类知识为 0，经过半个多月的研究， 终于有的小进展。 开始想从强大的互联网上搜点 c 代码来研究下， 结果搜遍所有 搜索引擎都搜不到，以下本人从修改的模糊控制代码，经过自己修改后可在 vc6.0 ，运行！输入 e 表示输出误差， ec 表示误差变化率，经过测试具有很好的 控制效果，对于非线性系统和数学模型难以建立的系统来说有更好的控制效果！ 现将其公开供大家学习研究！#include #includemath.h #define PMAX 100#define PMIN -100#define DMAX 100语言值的满幅值 */#define DMIN

2、-100 #define FMAX 100 /* int PFF4=0,12,24,48;/* 输入量 D 语言值特征点 */ int DFF4=0,16,32,64;/* 输出量 U 语言值特征点 */ int UFF7=0,15,30,45,60,75,90;/* 采用了调整因子的规则表 ,大误差时偏重误差 , 小误差时偏重误差变化 */*a0=0.3,a1=0.55,a2=0.74,a3=0.89 */int rule77=/ 误差变化率 -3,-2,-1, 0, 1, 2, 3 / 误差-6,-6,-6,-5,-5,-5,-4, / -3-5,-4,-4,-3,-2,-2,-1, /

3、-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) /* 模糊运算引擎 */int U; /*偏差, 偏差微分以及输出值的精确量 */unsigned int PF2,DF2,UF4; /* 偏差,偏差微分以及输出值的隶属度 */int Pn,Dn,Un4;long temp1,temp2;/* 隶属度的确定 */*根据PD的指定语言值获得

4、有效隶属度*/if(P-PFF3 & PPFF3)if(P=-PFF2)Pn=-2;PF0=FMAX*(float)(-PFF2-P)/(PFF3-PFF2);else if(P=-PFF1)Pn=-1;PF0=FMAX*(float)(-PFF1-P)/(PFF2-PFF1);else if(P=PFF0)Pn=0; PF0=FMAX*(float)(-PFF0-P)/(PFF1-PFF0);else if(P=PFF1)Pn=1; PF0=FMAX*(float)(PFF1-P)/(PFF1-PFF0);else if(P=PFF2)Pn=2; PF0=FMAX*(float)(PFF2-

5、P)/(PFF2-PFF1);else if(P=PFF3)Pn=3; PF0=FMAX*(float)(PFF3-P)/(PFF3-PFF2);else if(P=PFF3)Pn=3; PF0=0;PF1=FMAX-PF0;if(D-DFF3 & DDFF3) if(D=-DFF2)Dn=-2;DF0=FMAX*(float)(-DFF2-D)/(DFF3-DFF2); else if(D=-DFF1)Dn=-1;DF0=FMAX*(float)(-DFF1-D)/(DFF2-DFF1);else if(D=DFF0)Dn=0;DF0=FMAX*(float)(-DFF0-D)/(DFF1-

6、DFF0);else if(D=DFF1)Dn=1;DF0=FMAX*(float)(DFF1-D)/(DFF1-DFF0);else if(D=DFF2)Dn=2; DF0=FMAX*(float)(DFF2-D)/(DFF2-DFF1);else if(D=DFF3)Dn=3; DF0=FMAX*(float)(DFF3-D)/(DFF3-DFF2);else if(D=DFF3)Dn=3;DF0=0;DF1=FMAX-DF0;/* 使用误差范围优化后的规则表 rule77*/*输出值使用13个隶属函数,中心值由UFF7指定*/* 一般都是四个规则有效 */Un0=rulePn-1+3Dn

7、-1+3;Un1=rulePn+3Dn-1+3;Un2=rulePn-1+3Dn+3;Un3=rulePn+3Dn+3;if(PF0=DF0)UF0=PF0;elseUF0=DF0;if(PF1=DF0)UF1=PF1;elseUF1=DF0;if(PF0=DF1)UF2=PF0;elseUF2=DF1;if(PF1UF1)UF1=0;elseUF0=0;if(Un0=Un2)if(UF0UF2)UF2=0;elseUF0=0;if(Un0=Un3)if(UF0UF3)UF3=0;elseUF0=0;if(Un1=Un2)if(UF1UF2)UF2=0;elseUF1=0; if(Un1=Un

8、3)if(UF1UF3)UF3=0;elseUF1=0; if(Un2=Un3)if(UF2UF3)UF3=0;elseUF2=0;/* 重心法反模糊 */*Un 原值为输出隶属函数标号，转换为隶属函数值 */ if(Un0=0)Un0=UFFUn0;elseUn0=-UFF-Un0;if(Un1=0)Un1=UFFUn1;elseUn1=-UFF-Un1;if(Un2=0)Un2=UFFUn2;elseUn2=-UFF-Un2;if(Un3=0)Un3=UFFUn3;elseUn3=-UFF-Un3;temp1=UF0*Un0+UF1*Un1+UF2*Un2+UF3*Un3; temp2=UF0+UF1+UF2+UF3;U=temp1/temp2;return U; void main。 宀inf aHoeeG inf nowpoinLP-P2巴八 F_LE Jrvouc inHfopen (=in.bd=3八 oulrfopen(o=u!-* e) scanf(=%d=QOe)八 prinff (二*鸯 ec