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

回复此楼

» 猜你喜欢

酰胺脱乙酰基已经有12人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全已经有3人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全已经有3人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有4人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有10人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有5人回复
有院领导为了换新车，用横向课题经费买了俩车已经有5人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有6人回复
面上项目申报已经有3人回复
博士延得我，科研能力直往上蹿已经有7人回复

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

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

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

返回列表