- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
(完整版)模糊PID 控制的C 程序
(完整版)模糊PID 控制的C 程序
//e:[-3,3] ec:[-3,3] kp:[-0。3,0.3]
#include 〈stdio.h
#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(float 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:”); scanf(”%f”,b); Kp=FuzzyKp(a,b);
printf(”The parameter Kp is: %f\n\n,Kp); float FuzzyKp(float e,float ec)
{float es[7]; float ecs[7];
/******模糊推理规则的可信度通过取小点运算得到*****/ float form[7][7];//(表X7)(隶属度表)
int i,j;
/*************求得的最大值赋给form[a][b]*********/
int a=0,b=0; float lsd;int p; float detkp;
/**************************/ int kp[7][7]={{PB,PB,PM,PS,PS,ZO,ZO},
{PB,PB,PM,PS,PS,ZO,NS},
{PM,PM,PM,PS,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}};
es[NB]=ufl(e,—3,-1); es[NM]=uf(e,—3,-2,0); es[NS]=uf(e,-3,—1,1);
es[ZO]=uf(e,-2,0,2); es[PS]=uf(e,—1,1,3); es[PM]=uf(e,0,2,3); es[PB]=ufr(e,1,3);
ecs[NB]=ufl(ec,—3,-1);//ec ecs[NM]=uf(ec,—3,—2,0); ecs[NS]=uf(ec,-3,-1,1); ecs[ZO]=uf(ec,—2,0,2); ecs[PS]=uf(ec,-1,1,3); ecs[PM]=uf(ec,0,2,3); ecs[PB]=ufr(ec,1,3);
for(i=0;i〈7;i++)
{
float w,h,r; for(j=0;j〈7;j++)
{
h=es[i]; r=ecs[j]; w=fand(h,r); form[i][j]=w;
}
}
for(i=0;i7;i++)
{
for(j=0;j〈7;j++)
{ if(form[a][b]form[i][j])
{
a=i; b=j;
}
}
}
lsd=form[a][b];//es和ecs最值中的最大值p=kp[a][b];
if(p==NB)
detkp=cufl(lsd,—0 。 3,—0.1); else if(p==NM) detkp=cuf(lsd,—0。3,0。2,0); else if(p==NS)
detkp=cuf(lsd,-0.3,-0.1,0.1); else if(p==ZO) detkp=cuf(lsd,-0。2,0,0.2); else if(p==PS)
detkp=cuf(lsd,-0.1,0。1,0。3); else i
文档评论(0)