【板球仿真】模糊控制板球系统仿真【含Matlab源码 2982期】

Matlab领域

于 2024-01-15 15:15:19 发布

阅读量818

点赞数 18

分类专栏： Matlab物理应用（初级版）文章标签： matlab

本文链接：https://blog.csdn.net/weixin_63266434/article/details/135602355

版权

Matlab物理应用（初级版）专栏收录该内容

272 篇文章 8 订阅

订阅专栏

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【板球仿真】基于matlab模糊控制板球系统仿真【含Matlab源码 2982期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分源代码

clear all;
close all;

a=newfis(‘fuzzf’);

f1=0.5;
a=addvar(a,‘input’,‘e’,[-288f1,288f1]); % Parameter e
a=addmf(a,‘input’,1,‘NB’,‘trimf’,[-385f1,-288f1,-192f1]);
a=addmf(a,‘input’,1,‘NM’,‘trimf’,[-288f1,-192f1,-96f1]);
a=addmf(a,‘input’,1,‘NS’,‘trimf’,[-192f1,-96f1,0f1]);
a=addmf(a,‘input’,1,‘ZO’,‘trimf’,[-96f1,0,96f1]);
a=addmf(a,‘input’,1,‘PS’,‘trimf’,[0f1,96f1,192f1]);
a=addmf(a,‘input’,1,‘PM’,‘trimf’,[96f1,192f1,288f1]);
a=addmf(a,‘input’,1,‘PB’,‘trimf’,[192f1,288f1,385f1]);

f2=0.5;
a=addvar(a,‘input’,‘ec’,[-80f2,80f2]); % Parameter ec
a=addmf(a,‘input’,2,‘NB’,‘trimf’,[-107f2,-80f2,-53f2]);
a=addmf(a,‘input’,2,‘NM’,‘trimf’,[-80f2,-53f2,-27f2]);
a=addmf(a,‘input’,2,‘NS’,‘trimf’,[-53f2,-27f2,-0f2]);
a=addmf(a,‘input’,2,‘ZO’,‘trimf’,[-27f2,0f2,27f2]);
a=addmf(a,‘input’,2,‘PS’,‘trimf’,[0f2,27f2,53f2]);
a=addmf(a,‘input’,2,‘PM’,‘trimf’,[27f2,53f2,80f2]);
a=addmf(a,‘input’,2,‘PB’,‘trimf’,[53f2,80f2,107*f2]);

%f3=1.4;
f3=0.5;
a=addvar(a,‘output’,‘u’,[-30f3,30f3]); % Parameter u
a=addmf(a,‘output’,1,‘NB’,‘trimf’,[-40f3,-30f3,-20f3]);
a=addmf(a,‘output’,1,‘NM’,‘trimf’,[-30f3,-20f3,-10f3]);
a=addmf(a,‘output’,1,‘NS’,‘trimf’,[-20f3,-10f3,0f3]);
a=addmf(a,‘output’,1,‘ZO’,‘trimf’,[-10f3,0,10f3]);
a=addmf(a,‘output’,1,‘PS’,‘trimf’,[0f3,10f3,20f3]);
a=addmf(a,‘output’,1,‘PM’,‘trimf’,[10f3,20f3,30f3]);
a=addmf(a,‘output’,1,‘PB’,‘trimf’,[20f3,30f3,40f3]);

%Each rule is a PD rule: error=rin-yout (nagative feedback)

rulelist=[1 1 1 1 1; %Edit rule base
1 2 1 1 1;
1 3 1 1 1;
1 4 1 1 1;
1 5 2 1 1;
1 6 2 1 1;
1 7 3 1 1;

      2 1 1 1 1;
      2 2 1 1 1;
      2 3 1 1 1;
      2 4 2 1 1;
      2 5 2 1 1;
      2 6 3 1 1;
      2 7 3 1 1;
      
      3 1 2 1 1;
      3 2 2 1 1;
      3 3 3 1 1;
      3 4 3 1 1;
      3 5 4 1 1;
      3 6 5 1 1;
      3 7 5 1 1;
      
      4 1 1 1 1;
      4 2 2 1 1;
      4 3 3 1 1;
      4 4 4 1 1;
      4 5 5 1 1;
      4 6 5 1 1;
      4 7 7 1 1;
      
      5 1 3 1 1;
      5 2 3 1 1;
      5 3 4 1 1;
      5 4 5 1 1;
      5 5 5 1 1;
      5 6 6 1 1;
      5 7 6 1 1;
      
      6 1 5 1 1;
      6 2 5 1 1;
      6 3 6 1 1;
      6 4 6 1 1;
      6 5 7 1 1;
      6 6 7 1 1;
      6 7 7 1 1;
   
      7 1 5 1 1;
      7 2 6 1 1;
      7 3 6 1 1;
      7 4 7 1 1;
      7 5 7 1 1;
      7 6 7 1 1;
      7 7 7 1 1];

a=addrule(a,rulelist);
%showrule(a) % Show fuzzy rule base

a1=setfis(a,‘DefuzzMethod’,‘centroid’); % Defuzzy
writefis(a1,‘ljk’); % save to fuzzy file “ljk.fis” which can be
% simulated with fuzzy tool
a2=readfis(‘ljk’);
%plotfis(a2);

%fuzzy ljk.fis %Demo fuzzy control simulation
%ruleview(a2);

disp(‘ģ :e=[-3,+3],ec=[-3,+3]’);
Ulist=zeros(7,7);

for i=1:7
for j=1:7
e(i)=-4+i;
ec(j)=-4+j;
Ulist(i,j)=evalfis([e(i),ec(j)],a2);
end
end