C实现模糊控制算法

上传人:无*** 文档编号:100096483 上传时间:2022-06-02 格式:DOCX 页数:6 大小:22.41KB
收藏 版权申诉 举报 下载
C实现模糊控制算法_第1页
第1页 / 共6页
C实现模糊控制算法_第2页
第2页 / 共6页
C实现模糊控制算法_第3页
第3页 / 共6页
资源描述:

《C实现模糊控制算法》由会员分享,可在线阅读,更多相关《C实现模糊控制算法(6页珍藏版)》请在装配图网上搜索。

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

2、MAX100int 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,/-2-4,-3,-2,-1, 0, 1, 2,/-1-4,-3,-1,

3、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;/* 偏差 , 偏差微分以及输出值的隶属度 */intPn,Dn,Un4;long temp1,temp2;/* 隶属度的确定*/*根据PD的指定语言值获得有效隶属度*/if(P-PFF3 & PPFF3)if(P=-PFF2)Pn=-2;PF0=FMAX

4、*(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-P)/(PFF2-PFF1);else if(P=PFF3)Pn=3; PF0=FMAX*(float

5、)(PFF3-P)/(PFF3-PFF2);else if(P=PFF3)Pn=3;PF0=0;PF1=FMAX-PF0;if(D-DFF3 & DDFF3)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-DFF0);else if(D=DFF1)Dn=1;DF0=FMAX*(float)(DFF1-D)/(DFF1-DFF0);el

6、se 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-1+3;Un1=rulePn+3Dn-1+3;Un2=rulePn-1+3Dn+3;Un3=rulePn+3Dn+3;

7、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)UF2=0;elseUF1=0;if(Un1=Un3)if(UF1UF3)UF3=0;elseUF1=0;if(Un2=Un3)if(UF2UF3)UF3=0;elseUF2=0;/* 重心法

8、反模糊*/*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()int a=0,e,ec;/*int nowpoin

9、t,p1,p2=1;FILE *in,*out;in=fopen(in.txt,r);out=fopen(out.txt,w);*/while(!feof(in)while(1)/fscanf(in,%d,&nowpoint);/p1=nowpoint;/e=0-nowpoint;/ec= p1-p2;printf( 请输入 e:);scanf(%d,&e);printf( 请输入 ec:);scanf(%d,&ec);a=Fuzzy(e,ec);/fprintf(out,%d,a);/printf(%d:,p1);,e,ec);printf(e: %dec: %dprintf(a: %dn,a);/p2=p1;/fclose(in);/fclose(out);

展开阅读全文
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

copyright@ 2023-2025  zhuangpeitu.com 装配图网版权所有   联系电话:18123376007

备案号:ICP2024067431-1 川公网安备51140202000466号


本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!