开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > TAPI编程 > 查看源码
建模1.m
资源名称：rbf 网络设计实例.rar　[点击查看]
上传用户：yetwld
上传日期：2010-01-26
资源大小：82k
文件大小：3k
源码类别：

TAPI编程

开发平台：

Matlab

建模1.m：源码内容
  1. data=[-0.22914 -0.059577 -0.88853 -0.25353 -0.28422 0.096958;
  2. -0.48264 0.63645 0.025257 -0.36327 0.47082 -0.087841;
  3. -0.4228 0.28586 0.34107 0.048845 -0.76611 0.18407;
  4. -0.38717 -0.46772 0.14487 -0.23132 -0.089108 -0.74088;
  5. -0.40248 -0.5382 0.19561 -0.22917 0.24525 0.63041;
  6. -0.47317 -0.035252 -0.18518 0.83387 0.20625 -0.051838];
  7. x=data(:,1:5);t=data(:,6);
  8. c1=x(1:3,:);  %初值
  9. %第一次分类
  10. for i=1:1:6
  11.      for j=1:1:3
  12.     d1(i,j)=(x(i,:)-c1(j,:))*((x(i,:)-c1(j,:))');
  13. end
  14. end
  15. % x1   d1=[0    1.9658  1.9924;     
  16.  % x2    1.9658 0      1.9261;
  17.  % x3    1.9924 1.9261   0
  18.  % x4    1.298 1.5735  1.1445
  19.  % x5    1.7154 1.4841  1.8008
  20.  % x6    1.9779 1.9987  1.9443]     
  21. %类1（x1)
  22. %类2（x2,x5)
  23. %类3（x3,x4,x6)
  24. %中心c1=[-0.22914  -0.059577  -0.88853  -0.25353  -0.28422
  25. %       -0.48264 0.63645    0.025257    -0.36327    0.47082
  26. %       -0.4228     0.28586    0.34107    0.048845    -0.76611]
  28. %第二次聚类
  29. %新的中心c
  30. c2=[x(1,:);0.5*(x(2,:)+x(5,:));(x(3,:)+x(4,:)+x(6,:))./3];
  31. %c2=[ -0.22914   -0.059577  -0.88853  -0.25353  -0.28422
  32. %     -0.44256   0.049124  0.11043     -0.29622   0.35803
  33. %    -0.42772   -0.072373   0.10025   0.21713 -0.21632]
  34. %求d2
  35. for i=1:1:6
  36.      for j=1:1:3
  37.     d2(i,j)=(x(i,:)-c2(j,:))*((x(i,:)-c2(j,:))');
  38. end
  39. end
  40. %x1      d2=[0     1.4696   1.2434
  41. %x2          1.9658 0.37103   1.3201
  42. %x3          1.9924 1.4924   0.51694
  43. %x4          1.298 0.47553   0.37723
  44. %x5          1.7154 0.37103   0.63896
  45. %x6          1.9779 1.3956   0.64385]
  46. %聚类结果
  47. %类1(x1)
  48. %类2(x2,x5）
  49. %类3(x3,x4,x6)
  50. %
  51. %可知，c不再变化，故，分类结束
  52. %
  53. %最后的结果：c＝[ -0.22914   -0.059577  -0.88853  -0.25353  -0.28422
  54. %                -0.44256   0.049124  0.11043     -0.29622   0.35803
  55. %                -0.42772   -0.072373   0.10025   0.21713 -0.21632]
  56. %  求RBF基函数的宽度delta
  57. dd2=d2';
  58. dsum=sum(dd2);
  59. delta=[dsum(1) 0.5*(dsum(2)+dsum(5)) (dsum(2)+dsum(4)+dsum(6))./3];
  60. %delta=[2.713 3.1912 3.275]
  61. %据上述可知，隐含层数（采用高斯核函数）为3，输出层为线性输出
  62. for j=1:1:3
  63. for i=1:1:6
  64. a(j,i)=((x(i,:)-c2(j,:)))*((x(i,:)-c2(j,:))');
  65.     a(j,i)=exp(-a(j,i)./delta(j));
  66. end
  67. end
  68. %
  69. lr=0.001;max_epoch=500;err_goal=0.99;
  70. w=rands(3,1);
  71. y=a'*w;
  72. E=t-y;
  73. SSE=sumsqr(E);
  74. for epoch=1:1:max_epoch
  75.     if SSE<err_goal,epoch=epoch-1;break,end
  76.     w=w+lr*a*E;
  77.     y=a'*w;
  78.     E=t-y;
  79.     SSE=sumsqr(E);
  80. end
  81. y=a'*w;
  82. %标准化处理
  83. %ymax=max(y);
  84. %ymin=min(y);
  85. %i=ones(6,1);
  86. %y=0.9*i-0.8*(i*ymax-y)./(i*(ymax-ymin));
  87. %
  88. figure;
  89. plot(t,'r');
  90. hold on;
  91. plot(y,':*');
  92. title('RBF网络拟合曲线图');
  93. legend('化验值','估计值');
  94. ylabel('淀粉利用率（％）');
  95. xlabel('样本个数');
  96. axis([1,6,-1,1]);