1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
| function Bpmodeldistinguish( X,Y )
%*****************************************************
%利用matalb实现基于BP神经网络模式识别实现代码
%*****************************************************
% X是三个场景的平均特征向量所组成的矩阵（就是fea_average）
% Y是待匹配的特征向量（就是b1或b2）
%**************************************************************************
%**************** 先对X和Y进行标准化处理 *******************
%**************************************************************************
X=[1e9*X(1,:);1e3*X(2,:);X(3,:)/100];
X=X';
X0=ones(3,16);
X0(:,[1:3])=X;
X=X0;
Y=[1e9*Y(1);1e3*Y(2);Y(3)/100];
Y=Y';
Y0=ones(1,16);
Y0(1:3)=Y;
Y=Y0;
%**************************************************************************
%*******************主函数main*********************************************
%*************************START********************************************
%**************************************************************************
Y1=[1 0 0]; %输出模式
Y2=[0 1 0];
Y3=[0 0 1];
Yo=[Y1;Y2;Y3]; %输出的三种模式
n=16; %输入层神经元个数
p=8; %中间层神经元个数
q=3; %输出神经元个数
k=3 ;%训练模式个数
a1=0.2; b1=0.2; %学习系数，
%rou=0.5;%动量系数，
emax=0.01; cntmax=100;%最大误差，训练次数
[w,v,theta,r,t,mse]=bptrain(n,p,q,X,Yo,k,emax,cntmax,a1,b1);%调用函数bptrain训练网络
c=bptest(p,q,n,w,v,theta,r,Y);
c(c>0.5)=1;
c(c<0.2)=0;
disp('测试数据的识别结果:')
c
if (min(c==Y1)==1)
fprintf('该数据和场景1匹配\n');
elseif (min(c==Y2)==1)
fprintf('该数据和场景2匹配\n');
else
fprintf('该数据和场景3匹配\n');
end
%**************************************************************************
%**************************END*********************************************
%**************************************************************************
function [w,v,theta,r,t,mse]=bptrain(n,p,q,X,Yo,k,emax,cntmax,a1,b1)
%n表示输入神经元个数,p表示中间层神经元个数，q表示输出神经元个数,
%X表示输入训练模式，Yo表示标准输出，k表示训练模式的个数
%emax表示最大误差，cntmax表示最大训练次数，a1，b1表示学习系数，rou表示动量系数
%w、theta表示训练结束后输入层与中间层连接权系数和阈值，
%v、r表示训练结束后中间层与输出层连接权和阈值
%t表示训练时间
tic
w=rands(n,p);%输入层与隐含层连接权
v=rands(p,q); %隐含层与输出层连接权
theta=rands(1,p);%中间层的阈值
r=rands(1,q);%输出层的阈值
cnt=1;
mse=zeros(1,cntmax);%全局误差为零
er=0;
while ((er>emax)|(cnt<=cntmax))
E=zeros(1,q);
%循环识别模式
for cp=1:k
X0=X(cp,:);
Y0=Yo(cp,:);
%计算中间层的输入Y(j)
Y=X0*w;
%计算中间层的输出b
Y=Y-theta; %中间层阈值
for j=1:p
b(j)=1/(1+exp(-Y(j)));%中间层输出f(sj)
end
%计算输出层输出c
Y=b*v;
Y=Y-r; % 输出层阈值
for t=1:q
c(t)=1/(1+exp(-Y(t))); %输出层输出
end
%计算输出层校正误差d
for t=1:q
d(t)=(Y0(t)-c(t))*c(t)*(1-c(t));
end
%计算中间层校正误差e
xy=d*v';
for t=1:p
e(t)=xy(t)*b(t)*(1-b(t));
end
%计算下一次的中间层和输出层之间新的连接权v(i,j),阈值t2(j)
for t=1:q
for j=1:p
v(j,t)=v(j,t)+a1*d(t)*b(j);
end
r(t)=r(t)+a1*d(t);
end
%计算下一次的输入层和中间层之间新的连接权w(i,j),阈值t1(j)
for j=1:p
for i=1:n
w(i,j)=w(i,j)+b1*e(j)*X0(i);
end
theta(j)=theta(j)+b1*e(j);
end
for t=1:q
E(cp)=(Y0(t)-c(t))*(Y0(t)-c(t))+E(cp);%求当前学习模式的全局误差
end
E(cp)=E(cp)*0.5;
%输入下一模式
end
er=sum(E);%计算全局误差
mse(cnt)=er;
cnt=cnt+1;%更新学习次数
end
t=toc;
function c=bptest(p,q,n,w,v,theta,r,X)
%p,q,n表示输入层、中间层和输出层神经元个数
%w，v表示训练好的神经网络输入层到中间层、中间层到输出层权值。
%X为输入测试的模式
%c表示模式X送入神经网络的识别结果。
%计算中间层的输入Y(j)
Y=X*w;
%计算中间层的输出b
Y=Y-theta; %中间层阈值
for j=1:p
b(j)=1/(1+exp(-Y(j)));%中间层输出f(sj)
end
%计算输出层输出c
Y=b*v;
Y=Y-r; % 输出层阈值
thr1=0.01;thr2=0.5;
for t=1:q
c(t)=1/(1+exp(-Y(t))); %输出层输出
end
|
