RBF预测.asv
上传用户:yetwld
上传日期:2010-01-26
资源大小:82k
文件大小:2k
源码类别:

TAPI编程

开发平台:

Matlab

  1. %RBF预测模型
  2. %标准化后的测试数据集t_data
  3. t_data=[0.1 0.68668 0.67143 0.79999 0.79091 0.40004;
  4. 0.36667 0.58001 0.1 0.7 0.71818 0.20001;
  5. 0.26 0.47334 0.1 0.29997 0.2091 0.80009;
  6. 0.36667 0.9 0.9 0.29997 0.13637 0.37504;
  7. 0.26 0.84668 0.67143 0.1 0.42727 0.37504;
  8. 0.36667 0.58001 0.44286 0.49999 0.1 0.55006;
  9. 0.15333 0.47334 0.44286 0.7 0.42727 0.60006;
  10. 0.1 0.84668 0.67143 0.29997 0.5 0.1;
  11. 0.15333 0.42 0.21429 0.49999 0.5 0.55006;
  12. 0.20667 0.79335 0.21429 0.59999 0.5 0.32503;
  13. 0.1 0.42 0.21429 0.9 0.5 0.45005;
  14. 0.1 0.84668 0.32857 0.59999 0.5 0.27502;
  15. 0.20667 0.47334 0.32857 0.29997 0.13637 0.50005;
  16. 0.1 0.68668 0.67143 0.49999 0.24546 0.20001;
  17. 0.42 0.58001 0.21429 0.9 0.9 0.42504;
  18. 0.31334 0.58001 0.44286 0.49999 0.31818 0.25002;
  19. 0.15333 0.42 0.1 0.19999 0.35454 0.55006;
  20. 0.20667 0.47334 0.32857 0.29997 0.31818 0.27502;
  21. 0.15333 0.68668 0.44286 0.29997 0.31818 0.40004;
  22. 0.20667 0.20667 0.21429 0.39999 0.28183 0.52506;
  23. 0.26 0.79335 0.21429 0.49999 0.57273 0.9;
  24. 0.42 0.36667 0.1 0.59999 0.35454 0.30003;
  25. 0.47334 0.36667 0.1 0.59999 0.57273 0.35003;
  26. 0.1 0.47334 0.67143 0.7 0.42727 0.49894;
  27. 0.42 0.58001 0.67143 0.49999 0.24546 0.47505;
  28. 0.31334 0.1 0.32857 0.9 0.79091 0.8501;
  29. 0.1 0.52667 0.21429 0.9 0.5 0.50005;
  30. 0.52667 0.55867 0.21429 0.1 0.28183 0.42504;
  31. 0.9 0.58001 0.55715 0.1 0.17273 0.32503;
  32. 0.15333 0.68668 0.62572 0.29997 0.2091 0.57506];
  33. %初始化数据
  34. tt=t_data(:,6);x=t_data(:,1:5);tt=tt';
  35. %随机选取中心
  36. c=x;
  37. %定义delta平方为样本各点的协方差之和
  38. delta=cov(x');
  39. delta=sum(delta);
  40. %隐含层输出H
  41. for i=1:1:30
  42.   for j=1:1:30
  43.      R(i,j)=((x(i,:)-c(j,:)))*((x(i,:)-c(j,:))');
  44.      R(i,j)=exp(-R(i,j)./delta(j));
  45.   end
  46. end
  47. p=R;
  48. %建模
  49. %r=radbas(p);
  50. err_goal=0.005;
  51. sc=3;
  52. net=newrb(p,tt,err_goal,sc,50,5);
  53. %测试
  54. ty=sim(net,p);
  55. tE=tt-ty;
  56. tSSE=sse(tE);
  57. tMSE=mse(tE);
  58. %预测(测试)曲线
  59. figure;
  60.  plot(tt,'-+');
  61.  hold on;
  62.  plot(ty,'r:*');
  63.  legend('化验值','预测值');
  64.  title('BRF网络模型输出预测曲线');
  65.  xlabel('输入样本点');
  66.  ylabel('淀粉利用率(%)');
  67. axis([1,30,-2,2]);