二维模糊系统程序及matlab仿真运行结果

《二维模糊系统程序及matlab仿真运行结果》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、题目：设计一个在=1，l]x[-l,1]上的模糊系统，使其以精度£二0.1—致地逼近函数=sin(x1^)+cos(x2^)+sin(xI^)cos(x2-zr),并进行Matlab仿真。原理：由于OOdg=sup

2、m:os(x,^)+m:os(x1^)cos(x2^)

3、=27TxeU=sup

4、-^in(x2^)--7zsin(x,^)sin(x2-7r)

5、=2tixeU由式磬v磬〜可知，当取1^=0.0075,h2=0.0075时，有0.0075x2tt+0.0075x2^<0.1满足精度要求。由于L=2,此时模糊集的个数为N=-+1=267.67取N=268,h即七和％2分

6、别在"=[-1，1]上定义268个具有三角隶属函数的模糊集所设计的模糊系统为268268z/W/2=1268268该模糊系统由268x268=71824条规则來逼近函数。程序：%FuzzyApproachingclearallclearallT=0.1;xl=-l:T:l;x2=-l:T:l;L=2;h二0.0075;N=268;fori=4:4:Nforj=4:4:Nel(i)=-l+h*(i-4);e2⑴=-l+h*(j-4);gx(i,j)=sin(el(i)A3)+co$(e2(j)A3)+$in(el(i))*co$(e2(j));endenddf=zeros(L/

7、T+1,L/T+1);cf=zeros(L/T+l,L/T+1);form=4:4:Nifm==4u1=trimf(x+L/(N-4)1);elseifm==Nul=trimf(xl,[l-L/(N-4),1,1]);elseul=trimf(x1，[e1(m-4)，el(m)，e1(m+4)]);endforn=4:4:Nifn==4u2=trimf(x2,[-1，-1，•1+L/(N-4)]);elseifn==Nu2=trimf(x2,[l-L/(N-4),1,1]);elseu2=trimf(x2,Le2(n-4),e2(n),e2(n+4)J);endfigure。)

8、;subplot(211);holdon;plot(xl，ul);xlabel('xl’);ylabel('Membershipfunction');subplot(212);holdon;plot(x2，u2);xlabel(’x2.);ylabel(’Membershipfunction');fori=l:l:L/T+lforj=l:l:L/T+ld=df(ij)+ul(i)*u2(j)；df(i,j)=d;c=cf(i，j)+gx(m,n)*u1(i)*u2(j);cf(i，j)=c;endendendendfori=l:l:L/T+lforj=l:l:L/T+lf(i，

9、j)=cf(i，j)/df(i，j);y(i,j)=sin(el(i)A3)+cos(e2(j)A3)+sin(el(i))*cos(e2(j));endendfigure(2);subplot(311);surf(xl,x2,f)；titleff(x).);subplot(312);surf(xl,x2,y);title(’g(x)’)；subplot(313);surf(xl,x2,f-y);title('Approachingerror');程序运行结果:080.6犯-0.446似H.D.8•0.4似00.20.40.60.8080.60.20■°2-0.40.8•1.1

10、00.20.40.60.8物a麵酿