模糊PID控制的C程序

1、精品好资料学习推荐/e:-3,3 ec:-3,3 kp:-0.3,0.3#include#define NB 0#define NM 1#define NS 2#define ZO 3#define PS 4#define PM 5#define PB 6/*/float uf(float x,float a,float b,float c);float cuf(float x,float a,float b,float c);float ufl(float x,float a,float b);float cufl(float x,float a,float b);float ufr(flo

2、at x,float a,float b);float cufr(float x,float a,float b);float ufr(float x,float a,float b);float cufr(float x,float a,float b);float fand(float a,float b);float forr(float a,float b);float FuzzyKp(float e,float ec);/主程序void main()float a,b,Kp; printf(The E is:); scanf(%f,&a); printf(The Ec is:); s

3、canf(%f,&b); Kp=FuzzyKp(a,b); printf(The parameter Kp is: %fnn,Kp);float FuzzyKp(float e,float ec)float es7;float ecs7;/*模糊推理规则的可信度通过取小点运算得到*/float form77;/(表X7)（隶属度表）int i,j;/*求得的最大值赋给formab*/int a=0,b=0;float lsd;int p;float detkp;/*/int kp77=PB,PB,PM,PS,PS,ZO,ZO, PB,PB,PM,PS,PS,ZO,NS, PM,PM,PM,PS

4、,ZO,NS,NS, PM,PM,PS,ZO,NS,NM,NM, PS,PS,ZO,NS,NS,NM,NM, PS,ZO,NS,NM,NM,NM,NB, ZO,ZO,NM,NM,NM,NB,NB;esNB=ufl(e,-3,-1);esNM=uf(e,-3,-2,0);esNS=uf(e,-3,-1,1);esZO=uf(e,-2,0,2);esPS=uf(e,-1,1,3);esPM=uf(e,0,2,3);esPB=ufr(e,1,3);ecsNB=ufl(ec,-3,-1);/ececsNM=uf(ec,-3,-2,0);ecsNS=uf(ec,-3,-1,1);ecsZO=uf(ec,

5、-2,0,2);ecsPS=uf(ec,-1,1,3);ecsPM=uf(ec,0,2,3);ecsPB=ufr(ec,1,3);for(i=0;i7;i+)float w,h,r;for(j=0;j7;j+) h=esi;r=ecsj;w=fand(h,r);formij=w; for(i=0;i7;i+)for(j=0;j7;j+) if(formabformij) a=i; b=j; lsd=formab;/es和ecs最值中的最大值p=kpab;if(p=NB)detkp=cufl(lsd,-0.3,-0.1);elseif(p=NM)detkp=cuf(lsd,-0.3,0.2,0)

6、;elseif(p=NS)detkp=cuf(lsd,-0.3,-0.1,0.1);elseif(p=ZO)detkp=cuf(lsd,-0.2,0,0.2);elseif(p=PS)detkp=cuf(lsd,-0.1,0.1,0.3);elseif(p=PM)detkp=cuf(lsd,0,0.2,0.3);elseif(p=PB)detkp=cufr(lsd,0.1,0.3);return detkp;/*以下为附属子函数*/*求隶属度（三角形）模糊化处理*/float uf(float x,float a,float b,float c)if(x=a)return (0);elseif

7、(ax)&(x=b)return( (x-a)/(b-a);elseif(bx)&(xc)return (0);/*三角形反模糊化处理(最大隶属度法)*/float cuf(float x,float a,float b,float c) float y,z;z=(b-a)*x+a;y=c-(c-b)*x;return (y+z)/2;/*梯形(左)求隶属度模糊化*/float ufl(float x,float a,float b)if(x=a) return 1;elseif(ax)&(xb)return 0;/*梯形反模糊化*/float cufl(float x,float a,float b)return b-(b-a)*x;/*梯形(右)求隶属度模糊化*/float ufr(float x,float a,float b)if(x=a)return 0;if(ax)&(x=b)return 1;/*梯形反模糊化*/float cufr(float x,float a,float b)return (b-a)*x +