发帖

【优惠升级】华秋PCB首单最高立减100元，SMT免费贴片！！！

[问答]

请大家帮忙看一个模糊PID控制的M程序

8133

问答对人有帮助，内容完整，我也想知道答案 0 最近在写一个模糊PID控制的matlab m程序（原程序如下，）但是我运行的时候出现了如下这种错误： ??? Error using ==> close at 82 Specified window does not exist. Error in ==> Fuzzypidm at 1 close all； ??? Attempted to access rules(1,6); index out of bounds because size(rules)=[49,5]. Error in ==> setfis at 103 fis.rule(i).weight=rules(i, numInputs+numOutputs+1); Error in ==> addrule at 44 fis=setfis(fis,'ruleList',newRuleList); Error in ==> Fuzzypidm at 97 a=addrule(a,rulelist); 复制代码本人新手，不知道哪里有问题，请地里各位前辈帮忙看看程序。不甚感激！ clear all;close all； a=newfis('fuzzpidm'); a=addvar(a,'input','e',[-3,3]); a=addmf(a,'input',1,'NB','trimf',[-3,-3,-2]); a=addmf(a,'input',1,'NM','trimf',[-3,-2,-1]); a=addmf(a,'input',1,'NS','trimf',[-2,-1,0]); a=addmf(a,'input',1,'Z','trimf',[-1,0,1]); a=addmf(a,'input',1,'PS','trimf',[0,1,2]); a=addmf(a,'input',1,'PM','trimf',[1,2,3]); a=addmf(a,'input',1,'PB','trimf',[2,3,3]); a=addvar(a,'input','ec',[-3,3]); a=addmf(a,'input',2,'NB','trimf',[-3,-3,-2]); a=addmf(a,'input',2,'NM','trimf',[-3,-2,-1]); a=addmf(a,'input',2,'NS','trimf',[-2,-1,0]); a=addmf(a,'input',2,'Z','trimf',[-1,0,1]); a=addmf(a,'input',2,'PS','trimf',[0,1,2]); a=addmf(a,'input',2,'PM','trimf',[1,2,3]); a=addmf(a,'input',2,'PB','trimf',[2,3,3]); a=addvar(a,'output','kp',[0,3]); a=addmf(a,'output',1,'Z','trimf',[0,0,1]); a=addmf(a,'output',1,'PS','trimf',[0,1,2]); a=addmf(a,'output',1,'PM','trimf',[1,2,3]); a=addmf(a,'output',1,'PB','trimf',[2,3,3]); a=addvar(a,'output','ki',[0,3]); a=addmf(a,'output',2,'Z','trimf',[0,0,1]); a=addmf(a,'output',2,'PS','trimf',[0,1,2]); a=addmf(a,'output',2,'PM','trimf',[1,2,3]); a=addmf(a,'output',2,'PB','trimf',[2,3,3]); a=addvar(a,'output','kd',[0,3]); a=addmf(a,'output',3,'Z','trimf',[0,0,1]); a=addmf(a,'output',3,'PS','trimf',[0,1,2]); a=addmf(a,'output',3,'PM','trimf',[1,2,3]); a=addmf(a,'output',3,'PB','trimf',[2,3,3]); rulelist=[1 1 6 4 5; 1 2 5 5 6; 1 3 6 6 7; 1 4 6 7 5; 1 5 6 6 5; 1 6 5 5 6; 1 7 6 4 7; 2 1 7 4 5; 2 2 6 5 5; 2 3 7 7 6; 2 4 7 7 6; 2 5 7 7 4; 2 6 6 5 5; 2 7 5 4 5; 3 1 7 4 5; 3 2 6 4 4; 3 3 7 7 4; 3 4 7 7 5; 3 5 7 7 6; 3 6 6 4 4; 3 7 5 4 5; 4 1 7 4 7; 4 2 6 4 7; 4 3 7 7 5; 4 4 4 7 5; 4 5 7 7 5; 4 6 6 4 5; 4 7 5 4 4; 5 1 7 4 4; 5 2 6 4 5; 5 3 7 7 6; 5 4 7 7 4; 5 5 7 7 5; 5 6 6 4 5; 5 7 5 4 5; 6 1 7 4 5; 6 2 6 5 5; 6 3 7 7 5; 6 4 7 7 5; 6 5 7 7 5; 6 6 6 5 6; 6 7 5 4 5; 7 1 6 4 7; 7 2 5 5 7; 7 3 6 6 5; 7 4 6 7 5; 7 5 6 6 7; 7 6 5 5 7; 7 7 6 4 5]; a=addrule(a,rulelist); a1=setfis(a,'DefuzzMethod','mom');%Defuzzy writefis(a1,'fuzzpidm'); a=readfis('fuzzpidm'); %PID controller ts=10; sys=tf(2,[144,1],'inputdelay',30);% dsys=c2d(sys,ts,'zoh'); [num,den]=tfdata(dsys,'v'); u_1=0.0;u_2=0.0;u_3=0.0;u_4=0.0; y_1=0.0; x=[0,0,0]'; error_1=0; e_1=0.0; ec_1=0.0; kp0=2.88; kd0=0.048; ki0=43.2; 复制代码 %看论坛上加zero限定值。？也不知道这里加的对不对，但是如果不加的话，for循环里面的以下变量就会出现warning：说变量的值会随着循环而变化，考虑以下preallocating 什么的... time=zero(500); rin=zero(500); kp=zero(500); ki=zero(500); kd=zero(500); u=zero(500); yout=zero(500); error=zero(500); for k=1:1:500; time(k)=kts; rin(k)=1; %using fuzzy inference to tunning PID k_pid=evalfis([e_1,ec_1],a); kp(k)=k_pid(1); ki(k)=k_pid(2); kd(k)=k_pid(3); u(k)=kp(k)x(1)+kd(k)x(2)+10.5ki(k)x(3); yout(k)=-den(2)y_1+num(2)u_4; error(k)=rin(k)-yout(k); u_4=u_3; u_3=u_2; u_2=u_1; u_1=u(k); y_1=yout(k); x(1)=error(k); x(2)=error(k)-error_1; x(3)=x(3)+error(k); e_1=x(1); ec_1=x(2); error_2=error_1; error_1=errpr(k); end showrule(a) figure(1);plot(time,rin,'b',time,yout,'r'); xlabel('time(s)');ylabel('rin,yout'); figure(2);plot(time,error,'r'); xlabel('time(s)');ylabel('error'); figure(3);plot(time,u,'r'); xlabel('time(s)');ylabel('u'); figure(4);plot(time,kp,'r'); xlabel('time(s)');ylabel('kp'); figure(5);plot(time,ki,'r'); xlabel('time(s)');ylabel('ki'); figure(6);plot(time,kd,'r'); xlabel('time(s)');ylabel('kd'); figure(7);plotmf(a,'input',1); figure(8);plotmf(a,'input',2); figure(9);plotmf(a,'output',1); figure(10);plotmf(a,'output',2); figure(11);plotmf(a,'output',3); plotfis(a); fuzzy fuzzpidm.fis 复制代码* 1
2012-4-12 16:00:28　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × babe媛该类别下有 9 个回答。邀请回答大学生电子该类别下有 8 个回答。邀请回答迷失的傻该类别下有 7 个回答。邀请回答 CHENYIXIAN 该类别下有 7 个回答。邀请回答 upc_wh 该类别下有 7 个回答。邀请回答电子工程师该类别下有 6 个回答。邀请回答卟李该类别下有 6 个回答。邀请回答 xiejianghao 该类别下有 6 个回答。邀请回答晶晶若水该类别下有 5 个回答。邀请回答冰霜晨星该类别下有 5 个回答。邀请回答 aiolia8672 该类别下有 5 个回答。邀请回答 xxpaulxx 该类别下有 4 个回答。邀请回答 lilylu 该类别下有 4 个回答。邀请回答 dplion 该类别下有 4 个回答。邀请回答 gengyanan2031 该类别下有 4 个回答。邀请回答 xichyu 该类别下有 4 个回答。邀请回答 mawenbin123 该类别下有 4 个回答。邀请回答 chenbingjy 该类别下有 4 个回答。邀请回答归来的红雀该类别下有 4 个回答。邀请回答陈崇任该类别下有 4 个回答。邀请回答举报神之小风相关推荐 • 做了一个PID控制的程序，有些问题。求高手帮忙分析 3079 • 求大神帮忙，关于PID控制或模糊控制的单片机仿真及程序 4858 • 请问怎么用LabVIEW温度模糊PID控制程序，实现对温度信号的采集与模糊pid控制，最后输出一个电压量？ 4048 • labview中的模糊控制和pid结合的问题 2632 • 关于labview使用PID控制步进电机的问题，大家帮忙看看，谢谢！ 2109 • 请问如何去设计模糊PID自整定控制器？ 2469 • 请大家帮忙看看这个问题！ 2388 • labview模糊控制pid 1612 • 如何去设计一种基于STM32控制器的模糊PID温度控制系统呢 2519 • 请大家帮忙辨认一个工业相机 3519 4个回答

答案对人有帮助，有参考价值 0 a1=setfis(a,'DefuzzMethod','mom');%Defuzzywritefis(a1,'fuzzpidm'); a1改成a

2012-4-12 16:01:17 评论举报路过秋天

答案对人有帮助，有参考价值 0 路过秋天发表于 2012-4-12 16:01 a1=setfis(a,'DefuzzMethod','mom');%Defuzzywritefis(a1,'fuzzpidm'); 谢谢你的回答可是改了，还是不行。。。 ??? Error using ==> parsrule at 181 Output MF index is too high Error in ==> readfis at 231 out=parsrule(out,txtRuleList,'indexed'); Error in ==> Fuzzypidm at 102 a=readfis('fuzzpidm');

2012-4-12 16:02:23 评论举报神之小风

答案对人有帮助，有参考价值 0 嗯我习惯性的把fuzzy和PID控制分成2个函数来写分别编译2个M文件给你一个FUZZY PI控制应该是好用的你对照着看一下吧程序就不给你细看了 %Fuzzy Tunning PI Controlclear all;close all;a=newfis('fuzzpid'); a=addvar(a,'input','e',[-1,1]); %Parameter ea=addmf(a,'input',1,'N','zmf',[-1,-1/3]);a=addmf(a,'input',1,'Z','trimf',[-2/3,0,2/3]);a=addmf(a,'input',1,'P','smf',[1/3,1]); a=addvar(a,'input','ec',[-1,1]); %Parameter eca=addmf(a,'input',2,'N','zmf',[-1,-1/3]);a=addmf(a,'input',2,'Z','trimf',[-2/3,0,2/3]);a=addmf(a,'input',2,'P','smf',[1/3,1]); a=addvar(a,'output','kp',1/3[-10,10]); %Parameter kpa=addmf(a,'output',1,'N','zmf',1/3[-10,-3]);a=addmf(a,'output',1,'Z','trimf',1/3[-5,0,5]);a=addmf(a,'output',1,'P','smf',1/3[3,10]); a=addvar(a,'output','ki',1/30[-3,3]); %Parameter kia=addmf(a,'output',2,'N','zmf',1/30[-3,-1]);a=addmf(a,'output',2,'Z','trimf',1/30[-2,0,2]);a=addmf(a,'output',2,'P','smf',1/30[1,3]); rulelist=[1 1 1 2 1 1;1 2 1 2 1 1;1 3 1 2 1 1; 2 1 1 3 1 1;2 2 3 3 1 1;2 3 3 3 1 1; 3 1 3 2 1 1; 3 2 3 2 1 1;3 3 3 2 1 1];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);plotfis(a); fuzzy fuzzpid;showrule(a)ruleview fuzzpid; ========================================= %Fuzzy PI Controlclose all;clear all; warning off;a=readfis('fuzzpid'); %Load fuzzpid.fis ts=0.001;sys=tf(133,[1,25,0]);dsys=c2d(sys,ts,'z');[num,den]=tfdata(dsys,'v'); u_1=0;u_2=0;y_1=0;y_2=0;e_1=0;ec_1=0;ei=0; kp0=0;ki0=0;for k=1:1:1000time(k)=kts; yd(k)=1;%Using fuzzy inference to tunning PIk_pid=evalfis([e_1,ec_1],a);kp(k)=kp0+k_pid(1);ki(k)=ki0+k_pid(2);u(k)=kp(k)e_1+ki(k)ei; y(k)=-den(2)y_1-den(3)y_2+num(2)u_1+num(3)u_2;e(k)=yd(k)-y(k);%%%%%%%%%%%%%%Return of parameters%%%%%%%%%%%%%%%u_2=u_1;u_1=u(k); y_2=y_1;y_1=y(k); ei=ei+e(k)ts; % Calculating I ec(k)=e(k)-e_1;e_1=e(k);ec_1=ec(k);endfigure(1);plot(time,yd,'r',time,y,'k:','linewidth',2);xlabel('time(s)');ylabel('yd,y');legend('ideal position','position tracking');figure(2);subplot(211);plot(time,kp,'r','linewidth',2);xlabel('time(s)');ylabel('kp');subplot(212);plot(time,ki,'r','linewidth',2);xlabel('time(s)');ylabel('ki');figure(3);plot(time,u,'r','linewidth',2);xlabel('time(s)');ylabel('Control input'); 复制代码

2012-4-12 16:23:33 评论举报路过秋天

答案对人有帮助，有参考价值 0 你好，a=readfis（）语句出错的问题解决了吗？如果解决了能告诉我方法吗？我也遇到了同样的问题

2017-3-10 13:45:21 评论举报郝静仪

撰写答案

你正在撰写答案

如果你是对答案或其他答案精选点评或询问，请使用“评论”功能。

声明：本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人，不代表电子发烧友网立场。文章及其配图仅供工程师学习之用，如有内容图片侵权或者其他问题，请联系本站作侵删。侵权投诉

我要提问

精选推荐

想请教一下图中是simulink的什么模块，需要这种三段斜率函数模块但没找到在哪

1182 浏览 1 评论
怎么修改matlab与comsol联合的代码啊

2110 浏览 1 评论
200圆！求助大佬给一份VSG并网和离网模式的simulink仿真

2370 浏览 0 评论
MATLAB(3)--矩阵的引用(sub2ind、ind2sub、reshape函数使用)

3193 浏览 0 评论
求问MATLAB仿真KK接收器出现的问题？

3515 浏览 0 评论

快速回复 返回顶部 返回列表

关注微信公众号

电子发烧友网

电子发烧友论坛

社区合作: 刘勇; 联系电话：15994832713; 邮箱地址：liuyong@huaqiu.com

社区管理: elecfans短短; 微信：elecfans_666; 邮箱：users@huaqiu.com

【优惠升级】华秋PCB首单最高立减100元，SMT免费贴片！！！

返回 Matlab论坛

回复

关闭

站长推荐 /6

快速回复 返回顶部 返回列表

- 技术社区: HarmonyOS技术社区

RISC-V MCU技术社区

FPGA开发者技术社区

- OpenHarmony开源社区: OpenHarmony开源社区

- 嵌入式论坛: ARM技术论坛

STM32/STM8技术论坛

嵌入式技术论坛

单片机/MCU论坛

RISC-V技术论坛

瑞芯微Rockchip开发者社区

FPGA|CPLD|ASIC论坛

DSP论坛

- 电路图及DIY: 电路设计论坛

DIY及创意

电子元器件论坛

专家问答

- 电源技术论坛: 电源技术论坛

无线充电技术

- 综合技术与应用: 机器人论坛

USB论坛

电机控制

模拟技术

音视频技术

综合技术交流

上位机软件（C/Python/Java等）

- 无线通信论坛: WIFI技术

蓝牙技术

天线|RF射频|微波|雷达技术

- EDA设计论坛: PCB设计论坛

DigiPCBA论坛

Protel|AD|DXP论坛

PADS技术论坛

Allegro论坛

multisim论坛

proteus论坛|仿真论坛

KiCad EDA 中文论坛

DFM|可制造性设计论坛

- 测试测量论坛: LabVIEW论坛

Matlab论坛

测试测量技术

传感技术

- 招聘/交友/外包/交易/杂谈: 项目外包

供需及二手交易

工程师杂谈|交友

招聘|求职|工程师职场

- 官方社区: 发烧友官方/活动

华秋商城

华秋电路

time

recommend

hot

post

—
—
—

版
块
导
航