基于OLS 的RBF 网设计算法

gghhjj

1. SamNum = 100; % 训练样本数
   
2. TestSamNum = 101; % 测试样本数
   
3. SP = 0.6; % 隐节点扩展常数
   
4. ErrorLimit = 0.9; % 目标误差
   
5. % 根据目标函数获得样本输入输出
   
6. rand('state',sum(100*clock))
   
7. NoiseVar = 0.1;
   
8. Noise = NoiseVar*randn(1,SamNum);
   
9. SamIn = 8*rand(1,SamNum)-4;
   
10. SamOutNoNoise = 1.1*(1-SamIn+2*SamIn.^2).*exp(-SamIn.^2/2);
    
11. SamOut = SamOutNoNoise + Noise;
    
12. TestSamIn = -4:0.08:4;
    
13. TestSamOut = 1.1*(1-TestSamIn+2*TestSamIn.^2).*exp(-TestSamIn.^2/2);
    
14. figure
    
15. hold on
    
16. grid
    
17. plot(SamIn,SamOut,'k+')
    
18. plot(TestSamIn,TestSamOut,'k--')
    
19. xlabel('Input x');
    
20. ylabel('Output y');
    
21. [InDim,MaxUnitNum] = size(SamIn); % 样本输入维数和最大允许隐节点数
    
22. % 计算隐节点输出阵
    
23. Distance = dist(SamIn',SamIn);
    
24. HiddenUnitOut = radbas(Distance/SP);
    
25. PosSelected = [];
    
26. VectorsSelected = [];
    
27. HiddenUnitOutSelected = [];
    
28. ErrHistory = []; % 用于记录每次增加隐节点后的训练误差
    
29. VectorsSelectFrom = HiddenUnitOut;
    
30. dd = sum((SamOut.*SamOut)')';
    
31. for k = 1 : MaxUnitNum
    
32. % 计算各隐节点输出矢量与目标输出矢量的夹角平方值
    
33. PP = sum(VectorsSelectFrom.*VectorsSelectFrom)';
    
34. Denominator = dd * PP';
    
35. [xxx,SelectedNum] = size(PosSelected);
    
36. if SelectedNum>0,
    
37. [lin,xxx] = size(Denominator);
    
38. Denominator(:,PosSelected) = ones(lin,1);
    
39. end
    
40. Angle = ((SamOut*VectorsSelectFrom) .^ 2) ./ Denominator;
    
41. % 选择具有最大投影的矢量，得到相应的数据中心
    
42. [value,pos] = max(Angle);
    
43. PosSelected = [PosSelected pos];
    
44. % 计算RBF 网训练误差
    
45. HiddenUnitOutSelected = [HiddenUnitOutSelected; HiddenUnitOut(pos,];
    
46. HiddenUnitOutEx = [HiddenUnitOutSelected; ones(1,SamNum)];
    
47. W2Ex = SamOut*pinv(HiddenUnitOutEx); % 用广义逆求广义输出权值
    
48. W2 = W2Ex(:,1:k); % 得到输出权值
    
49. B2 = W2Ex(:,k+1); % 得到偏移
    
50. NNOut = W2*HiddenUnitOutSelected+B2; % 计算RBF 网输出
    
51. SSE = sumsqr(SamOut-NNOut)
    
52. % 记录每次增加隐节点后的训练误差
    
53. ErrHistory = [ErrHistory SSE];
    
54. if SSE < ErrorLimit, break, end
    
55. % 作Gram-Schmidt 正交化
    
56. NewVector = VectorsSelectFrom(:,pos);
    
57. ProjectionLen = NewVector' * VectorsSelectFrom / (NewVector'*NewVector);
    
58. VectorsSelectFrom = VectorsSelectFrom - NewVector * ProjectionLen;
    
59. end
    
60. UnitCenters = SamIn(PosSelected);%%%%%%%%%%%
    
61. % 测试
    
62. TestDistance = dist(UnitCenters',TestSamIn);%%%%%%%%
    
63. TestHiddenUnitOut = radbas(TestDistance/SP);
    
64. TestNNOut = W2*TestHiddenUnitOut+B2;
    
65. plot(TestSamIn,TestNNOut,'k-')
    
66. k
    
67. UnitCenters
    
68. W2
    
69. B2

复制代码

nicaili

HiddenUnitOutEx = [HiddenUnitOutSelected; ones(1,SamNum)];
Error in ==> jiyuolsderbf at 46
HiddenUnitOutEx = [HiddenUnitOutSelected; ones(1,SamNum)];
麻烦楼主看看！

heiyanzi123

楼上的朋友，你的程序有点小错，HiddenUnitOutSelected = [HiddenUnitOutSelected; HiddenUnitOut(pos,];在这句话里HiddenUnitOut(pos,]是不是写的不全？麻烦你看看，是不是忘了写什么东西？谢谢

jiarunda

加上 :) 就可以了吧！

ying_327

:handshake