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

回复此楼

» 猜你喜欢

有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复

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

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

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

返回列表