• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

如何在matlab中将模糊控制系统中的规则转化为数值表

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

如何在matlab中将模糊控制系统中的规则转化为数值表
第一步首先要建立一个模糊控制器的规则
主要有两种方式
1 直接list清单的方式构建
2 在rule生成器中填写
本文给出了第一种方式的代码如下

clear all;
close all;
load list %这个地方是装载一个规则表,具体代码最后会附带!
% rulelist=[1 1 1 1;
% 2 2 1 1;
% 3 3 1 1;
% 4 4 1 1;
% 5 5 1 1];
a=newfis(‘fuzzpid’);%创建模糊推理系统结构
a = addvar(a,‘input’,‘e’,0.07580[-1,1]); %Parameter eNH,NM,NL,0,PL,PM,PH
a = addmf(a,‘input’,1,‘NH’,‘trimf’,[-6,-6,-3]);
a = addmf(a,‘input’,1,‘NM’,‘trimf’,[-6,-3,-2]);
a =addmf(a,‘input’,1,‘NL’,‘trimf’,[-3,-1,0]);
a =addmf(a,‘input’,1,‘ZO’,‘trimf’,[-1,0,1]);
a =addmf(a,‘input’,1,‘PL’,‘trimf’,[0,1,3]);
a = addmf(a,‘input’,1,‘PM’,‘trimf’,[2,3,6]);
a = addmf(a,‘input’,1,‘PH’,‘trimf’,[3,6,6]);

a=addvar(a,‘input’,‘ec’,0.610[-1,1]); %Parameter ec
a = addmf(a,‘input’,2,‘NH’,‘trimf’,[-6,-6,-3]);
a = addmf(a,‘input’,2,‘NM’,‘trimf’,[-6,-3,-2]);
a =addmf(a,‘input’,2,‘NL’,‘trimf’,[-3,-1,0]);
a =addmf(a,‘input’,2,‘ZO’,‘trimf’,[-1,0,1]);
a =addmf(a,‘input’,2,‘PL’,‘trimf’,[0,1,3]);
a = addmf(a,‘input’,2,‘PM’,‘trimf’,[2,3,6]);
a = addmf(a,‘input’,2,‘PH’,‘trimf’,[3,6,6]);

a = addvar(a,‘output’,’ kp’,[-3,3]); %Parameter kp
a = addmf(a,‘output’,1,‘NH’,‘trimf’,[-3,-3,-2]);
a = addmf(a,‘output’,1,‘NM’,‘trimf’,[-3,-2,-1]);
a =addmf(a,‘output’,1,‘NL’,‘trimf’,[-2,-1,0]);
a =addmf(a,‘output’,1,‘ZO’,‘trimf’,[-1,0,1]);
a =addmf(a,‘output’,1,‘PL’,‘trimf’,[0,1,2]);
a = addmf(a,‘output’,1,‘PM’,‘trimf’,[1,2,3]);
a = addmf(a,‘output’,1,‘PH’,‘trimf’,[2,3,3]);

a = addvar(a,‘output’,‘ki’,[-3,3]); %Parameter ki
a = addmf(a,‘output’,2,‘NH’,‘trimf’,[-3,-3,-2]);
a = addmf(a,‘output’,2,‘NM’,‘trimf’,[-3,-2,-1]);
a =addmf(a,‘output’,2,‘NL’,‘trimf’,[-2,-1,0]);
a =addmf(a,‘output’,2,‘ZO’,‘trimf’,[-1,0,1]);
a =addmf(a,‘output’,2,‘PL’,‘trimf’,[0,1,2]);
a = addmf(a,‘output’,2,‘PM’,‘trimf’,[1,2,3]);
a = addmf(a,‘output’,2,‘PH’,‘trimf’,[2,3,3]);

a = addvar(a,‘output’,‘kd’,[-3,3]); %Parameter kd
a = addmf(a,‘output’,3,‘NH’,‘trimf’,[-3,-3,-2]);
a = addmf(a,‘output’,3,‘NM’,‘trimf’,[-3,-2,-1]);
a =addmf(a,‘output’,3,‘NL’,‘trimf’,[-2,-1,0]);
a =addmf(a,‘output’,3,‘ZO’,‘trimf’,[-1,0,1]);
a =addmf(a,‘output’,3,‘PL’,‘trimf’,[0,1,2]);
a = addmf(a,‘output’,3,‘PM’,‘trimf’,[1,2,3]);
a = addmf(a,‘output’,3,‘PH’,‘trimf’,[2,3,3]);
rulelist=list;
a=addrule(a,rulelist);
a=setfis(a,‘DefuzzMethod’,‘centroid’);
writefis(a,‘fuzzpid’);

a=readfis(‘fuzzpid’);
figure(1);
plotmf(a,‘input’,1);
figure(2);
plotmf(a,‘input’,2);
figure(3);
plotmf(a,‘output’,1);
figure(4);
plotmf(a,‘output’,2);
figure(5);
plotmf(a,‘output’,3);
figure(6);
plotfis(a);

fuzzy fuzzpid;
showrule(a)
%%附带的7*7的表
7 7 1 7 1 1 1
7 6 1 7 5 1 1
7 5 2 6 7 1 1
7 4 2 6 7 1 1
7 3 3 5 7 1 1
7 2 4 4 6 1 1
7 1 4 4 3 1 1
6 7 1 7 3 1 1
6 6 1 7 5 1 1
6 5 2 6 7 1 1
6 4 3 5 6 1 1
6 3 3 5 6 1 1
6 2 4 4 5 1 1
6 1 5 4 4 1 1
5 7 2 7 4 1 1
5 6 2 6 5 1 1
5 5 2 5 6 1 1
5 4 3 5 6 1 1
5 3 4 4 5 1 1
5 2 5 3 5 1 1
5 1 5 3 4 1 1
4 7 2 6 4 1 1
4 6 2 6 5 1 1
4 5 3 5 5 1 1
4 4 4 4 5 1 1
4 3 5 3 5 1 1
4 2 6 2 5 1 1
4 1 6 2 4 1 1
3 7 3 6 4 1 1
3 6 3 5 4 1 1
3 5 4 4 4 1 1
3 4 5 3 4 1 1
3 3 5 3 4 1 1
3 2 6 2 4 1 1
3 1 7 1 4 1 1
2 7 3 4 1 1 1
2 6 4 4 5 1 1
2 5 5 3 3 1 1
2 4 6 3 3 1 1
2 3 6 2 3 1 1
2 2 6 1 3 1 1
2 1 7 1 1 1 1
1 7 4 4 3 1 1
1 6 4 4 5 1 1
1 5 6 3 7 1 1
1 4 6 2 7 1 1
1 3 6 2 7 1 1
1 2 7 1 6 1 1
1 1 7 1 3 1 1
第二步 建立simulink


具体的工具就在常用模块就可以找到,在搜索地方打出fuzzy可以得到模糊模块


注意!!!需要对使用readfis函数读取一下之前建立的模糊系统

第三步 在simulink 中通过 Tools -> System Test -> Launch System Test 进入系统测试界面。

进去之后在菜单栏选择 Insert -> Test Element -> Simulink ,选择第二步中建立的模型 fuzzy_model
选择之间建的simulink模型的路径

之后首先点test vector,选择new新建二个,命名input1和input2,
设置两个输入变量,因为我们是遍历选择需要将所有的模糊区间选到
我这里面是-6到6每次变化1

随后进入测试变量,选择new新建一个outlet

对properties画圈的地方设置输入变量和输出变量,最后点运行,等待一到两分钟就行了。

最后一步是调取参数表的数据
在command 界面输入:
data_biao=stresults.ResultsDataSet.output;
temp=zeros(169,1);
for i=1:169
temp(i)=data_biao{i}(6);
end
table_data=reshape(temp,13,13);
这样就可以得到你想要的参数表了
我的例子的参数表如下


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
matlab画一个局部放大的图中图发布时间:2022-07-18
下一篇:
matlab基础学习——常用数学函数以及数据类型之间的转换发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap