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

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

TAPI编程

开发平台：

Matlab

11.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)=(-a(j,i)./delta(j));
  66. end
  67. end
  68. P=1:1:6;
  69. T=t;
  70. %figure;subplot(2,2,1);plot(P,t);title('待逼近的函数样本点');
  71. %axis([1,6,-1,1])
  72.   p=a;
  73.  r=radbas(p);
  74.  err_goal=0.99;
  75.  sc=1;
  76.  net=newrb(p,T',err_goal,sc);
  77. %NEWRB, neurons = 0, SSE = 0.999213
  78. Y=sim(net,p);
  79. %axis([1,6,-0.4,0.4]);
  80. figure;
  81. plot(P,T,'r');
  82. hold on;
  83. plot(P,Y,':*');
  84. title('RBF网络拟合曲线图11');
  85. legend('化验值','估计值');
  86. ylabel('淀粉利用率（％）');
  87. xlabel('样本个数');
  88. axis([1,6,-1,1]);