24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 xingshen 的 10 个金币

xingshen

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 428
帖子: 8
在线: 8分钟
虫号: 774518
注册: 2009-05-19
性别: GG
专业: 控制理论与控制工程

[交流] 【求助】编了个模糊神经网络的程序，但是运行基本没效果，高手给看看哪出错了

clear
%clc
tic,
x=[0 0 0 0 0 0.125 0.125 0.125 0.125 0.25...
0.25 0.25 0.25 0.25 0.25 0.375 0.375 0.375...
0.375 0.5 0.5 0.5 0.5 0.5 0.5 0.625 0.625 0.625...
0.625 0.75 0.75 0.75 0.75 0.75 0.875 0.875 0.875 1.0 1.0;
0.02 0.05 0.09 0.11 0.26 0.08 0.11 0.19 0.49 0.17 0.19...
0.22 0.26 0.36 0.52 0.26 0.3 0.4 0.61 0.39 0.41 0.44 0.49...
0.57 0.71 0.51 0.59 0.65 0.85 0.66 0.7 0.75 0.83 0.95 0.84 0.92 0.99 0.97 0.99];
y=[0.643 0.50 0.357 0.214 0.071 0.786 0.571 0.286 0.071...
0.857 0.714 0.571 0.429 0.286 0.143 1.0 0.714 0.429 0.214...
1.0 0.857 0.714 0.571 0.429 0.286 1.0 0.714 0.571 0.357 1.0...
0.857 0.714 0.571 0.429 0.929 0.714 0.571 0.929 0.786];
[p1,p2]=size(x);
%利用模糊c-均值法对样本聚类确定模糊规则个数
k=7;
%初始化隶属度函数的参数m,b及输出层初始权值
for i=1:p1;
for j=1:k;
   m(i,j)=rand(1);
   b(i,j)=rand(1);
end
end
for  j=1:k;
   w(j)=rand(1);
end
%推理计算输出值
for q=1:p2;
%选用高斯函数作为隶属度
for i=1:p1;
for j=1:k;
   u(i,j)=gaussmf(x(i,q),[m(i,j),b(i,j)]);
end
end
for i=1:k;
v(i)=1;
j=1;
while  j<=p1;
      v(i)=v(i)*u(j,i);
      j=j+1;
end
end
out1(q)=w*v';
e(q)=y(q)-out1(q);
end
out=out1;
E=(1/2)*sumsqr(y-out);
EE=E;
lr=0.3;
%误差反传后的参数修正过程
r=1;
p=1;
s=1000;
while p<=s & EE>1e-4
%隶属度参数m,b输出层权值参数w的修正过程
wc=zeros(1,k);
for i=1:k;
wc(i)=lr*e(r)*v(i);
end
mc=zeros(p1,k);
for i=1:p1;
   for j=1:k;
      mc(i,j)=2*lr*e(r)*w(j)*(v(j)./u(i,j))*exp(-((x(i,r)-m(i,j)).^2)./(b(i,j).^2))*(x(i,r)-m(i,j))./(b(i,j).^2);
   end
end
bc=zeros(p1,k);
for i=1:p1;
   for j=1:k;
      bc(i,j)=2*lr*e(r)*w(j)*(v(j)./u(i,j))*exp(-((x(i,r)-m(i,j)).^2)/(b(i,j).^2))*((x(i,r)-m(i,j)).^2)/(b(i,j).^3);
   end
end
%参数修正m b w
m=m+mc;
b=b+bc;
w=w+wc;
%利用修正过的参数重新计算输出
for q=1:p2;
for i=1:p1;
   for j=1:k;
      u(i,j)=gaussmf(x(i,q),[m(i,j),b(i,j)]);
   end
end
for i=1:k;
   v(i)=1;
   j=1;
   while  j<=p1;
      v(i)=v(i)*u(j,i);
      j=j+1;
   end
end
out1(q)=w*v';
end
  out=out1;
p=p+1;
EE=(1/2)*sumsqr(y-out);
E(p)=EE;
r=r+1;
if  r>p2
   r=1;
end
e(r)=(y(r)-out(r));
end
m,b,w,E,e
epoch=1:size(E,2);
figure
plot(epoch,E,'r');
axis([0 1.5*s min(E) max(E)]);
set(gca,'xtick',0:s/10:1.5*s);
title('误差变化曲线');xlabel('步数');ylabel('误差');
figure
plot(1:39,y,'-',1:39,out,'*')
toc

回复此楼

» 猜你喜欢

求推荐期刊，重谢已经有4人回复
无聊看看时间戳打发时间已经有5人回复
无聊看看filecode 已经有7人回复
有机能跨方向申博吗已经有3人回复
求助！！已经有4人回复
kh550接枝Sio2失败求助已经有3人回复
职称论文投稿已经有12人回复
反应求助已经有6人回复
基于自然哲学类比的风化壳型稀土矿已经有16人回复
看《给阿ma的情书》有感已经有10人回复

1楼 2009-11-29 16:51:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 xingshen 的主题更新

返回列表