基于模糊逻辑的二阶系统控制算法出错指教

ylhit

二阶系统的参数为num=20;den=[1.6 4.4 1]。我用的是matlab6.5 出现的错误是：evalfis使用错误。
小弟新学模糊控制，望高手指点。

[ 本帖最后由 eight 于 2008-4-10 10:21 编辑 ]

sigma665

把错误提示贴出来

ylhit

a = -1.000000, b = 0.000000, c = 1.000000, d = 0.000000
??? Illegal parameters in fisTrapezoidMf() --> c > d

Error in ==> C:\MATLAB6p5p1\toolbox\fuzzy\fuzzy\evalfismex.dll
Error in ==> C:\MATLAB6p5p1\toolbox\fuzzy\fuzzy\evalfis.m
On line 80  ==> [output,IRR,ORR,ARR] = evalfismex(input, fis, numofpoints);

Error in ==> C:\Documents and Settings\bai\桌面\flc.m
On line 62  ==>     u=ku*evalfis(in,a);

ylhit

自己重新检查了一下，上面的错误是由于建立fis时出错。以下是修改之后的程序。
ke=30;kd=0.2;ku=1;
num=20;
den=[1.6 4.4 1];
[a1,b,c,d]=tf2ss(num,den);
x=[0;0];%the init  2*1 vector
T=0.01;h=T; %step
N=250;
R=1.5*ones(1,N);

a=newfis('simple');
a=addvar(a,'input','e',[-6,6]);
a=addmf(a,'input',1,'NB','trapmf',[-6,-6,-5,-3]);
a=addmf(a,'input',1,'NS','trapmf',[-5,-3,-2,0]);
a=addmf(a,'input',1,'ZR','trimf',[-2,0,2]);
a=addmf(a,'input',1,'PS','trapmf',[0,2,3,5]);
a=addmf(a,'input',1,'PB','trapmf',[3,5,6,6]);
a=addvar(a,'input','de',[-6,6]);
a=addmf(a,'input',2,'NB','trapmf',[-6 -6 -5 -3]);
a=addmf(a,'input',2,'NS','trapmf',[-5 -3 -2 0]);
a=addmf(a,'input',2,'ZR','trimf',[-2 0 2]);
a=addmf(a,'input',2,'PS','trapmf',[0 2 3 5]);
a=addmf(a,'input',2,'PB','trapmf',[3 5 6 6]);
a=addvar(a,'output','u',[-3,3]);
a=addmf(a,'output',1,'NB','trapmf',[-3 -3 -3 -2]);
a=addmf(a,'output',1,'NS','trimf',[-2 -1 0]);
a=addmf(a,'output',1,'ZR','trimf',[-1 0 1]);
a=addmf(a,'output',1,'PS','trimf',[0 1 2]);
a=addmf(a,'output',1,'PB','trapmf',[2 3 3 3]);

rr=[5 5 4 4 3
    5 4 4 3 3
    4 4 3 3 2
    4 3 3 2 2
    3 3 2 2 1];
r1=zeros(prod(size(rr)),3);k=1;
for i=1:size(rr,1)
    for j=1:size(rr,2)
        r1(k,:)=[i,j,rr(i,j)];
        k=k+1;
    end
end
[r,s]=size(r1);
r2=ones(r,2);
rulelist=[r1,r2];
a=addrule(a,rulelist);

e=0;de=0;


for k=1:N
    e1=ke*e;
    de1=kd*de;
    if e1>=6
        e1=6;
    elseif e1<=-6
        e1=-6;
    end
    if de1>=6
        de1=6;
    elseif de1<=-6
        de1=-6;
    end
    %fuzzy controler
    in=[e1;de1]
    u=ku*evalfis(in,a);
    uu(1,k)=u;
    %under control system
    k0=a1*x+b*u;
    k1=a1*(x+h*k0/2)+b*u;
    k2=a1*(x+h*k1/2)+b*u;
    k3=a1*(x+h*k2)+b*u;
    x=x+(k0+2*k1+2*k2+k3)*h/6;
    y=c*x+d*u;
    yy(1,k)=y;
   
    e1=e;
    e=y-R(1,k);
    de=(e-e1)/T;
end

kk=[1:N]*T;
figure(1);
plot(kk,R,'k',kk,yy,'r');grid on
xlabel('time(s)');ylabel('output');

chg555

我也遇到了上述问题，但程序是书上的，建立FIS时应该没错，能解释一下具体原因吗？
a=newfis('a');
a=addvar(a,'input','wd',[0,30]);
a=addmf(a,'input',1,'len','trapmf',[0,0,10,20]);
a=addmf(a,'input',1,'re','trapmf',[10,20,30,30]);
a=addvar(a,'input','pr',[0,3]);
a=addmf(a,'input',2,'nor','trapmf',[0,1,2]);
a=addmf(a,'input',2,'high','trapmf',[1,2,3,3]);
a=addvar(a,'output','zl',[-10,10]);
a=addmf(a,'output',1,'fu','trimf',[-10,-5,0]);
a=addmf(a,'output',1,'ling','trimf',[-5,0,5]);
a=addmf(a,'output',1,'zheng','trimf',[0,5,10]);
rulelist=[1 2 3 1 1;2 2 1 1 1;0 1 2 1 2];
a=addrule(a,rulelist);
>> evalfis([20 2;10 3],a)
a = 0.000000, b = 1.000000, c = 2.000000, d = 0.000000
??? Error using ==> evalfismex
Illegal parameters in fisTrapezoidMf() --> c > d
Error in ==> evalfis at 84
[output,IRR,ORR,ARR] = evalfismex(input, fis, numofpoints);

[ 本帖最后由 ChaChing 于 2009-10-30 23:35 编辑 ]