请高手帮忙！！关于BP实现

akyi · 发表于 2006-4-11 10:34

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

为什么我编的程序运算速度刚刚开始还可以但是到了后面的速度就比较慢了？ 谢谢帮忙！！！ clear % clc tic PP =[1 0 0 1 ;0 1 0 1 ]; TT = [1 1 0 0 ]; HiddenUnitNum = 2; [InDim,SamNum] = size(PP); [OutDim,SamNum] = size(TT); SamIn=PP; SamOut=TT;
MaxEpochs=30000; lr=0.5; E0=0.001; TotalEpochs = 0;
for n = 0:29 %rand('seed',1); %rand('state',sum(100*clock)); W1=0.6*rand(HiddenUnitNum,InDim)-0.3; %初始化权值[-0.3,0.3] B1=0.6*rand(HiddenUnitNum,1)-0.3; W2=0.6*rand(OutDim,HiddenUnitNum)-0.3; B2=0.6*rand(OutDim,1)-0.3; W1Ex=[W1 B1]; W2Ex=[W2 B2]; w1=0; % w2=0; alpha=0.7; % I_H权值的动量因子 beta=0.7; %H_O权值的动量因子
SamInEx=[SamIn' ones(SamNum,1)]'; % 扩展输入矩阵 ErrHistory=[];
for i=1:MaxEpochs HiddenOut=logsig(W1Ex*SamInEx); %隐层输出 HiddenOutEx=[HiddenOut' ones(SamNum,1)]'; NetworkOut=logsig(W2Ex*HiddenOutEx); % 输出层输出 Error=SamOut-NetworkOut; % 误差 SSE=(1/(2*SamNum))*sumsqr(Error); % sumsqr(): Sum squared elements of a matrix ErrHistory=[ErrHistory SSE]; if SSE<E0,break,end %Delta2=Error; %% Delta2 = Error.*NetworkOut.*(1-NetworkOut); %%隐层神经元的误差信号 Delta1=W2'*Delta2.*HiddenOut.*(1-HiddenOut); %% dW2Ex=Delta2*HiddenOutEx'; % 调整量 dW1Ex=Delta1*SamInEx'; % %%%%%%%%%%%%%标准BP算法%%%%%%%%%%% % W1Ex=W1Ex+lr*dW1Ex; % W2Ex=W2Ex+lr*dW2Ex; %%%%%%%%%%%%%%%%%%%%加入动量项的BP算法%%%%%%%%%%% a=W1Ex; b=W2Ex; W1Ex=W1Ex+lr*dW1Ex+alpha*w1; W2Ex=W2Ex+lr*dW2Ex+beta*w2; w1=W1Ex-a; w2=W2Ex-b; %%%%%%%%%%%%%%%%%%%%%%%%%%%%% % W2=W2Ex(:,1:HiddenUnitNum); %%W1=W1Ex(:,1:InDim); end TotalEpochs = TotalEpochs + i;%TotalEpochs=TotalEpochs+1; end
% 输出结果 AverageEpochs = TotalEpochs/30 W1=W1Ex(:,1:InDim) B1=W1Ex(:,InDim+1) W2=W2Ex(:,1:HiddenUnitNum) B2=W2Ex(:,1+HiddenUnitNum)
[xx,Num]=size(ErrHistory); plot(1:Num,ErrHistory,'k--'); title('训练误差') toc

happy · 发表于 2006-4-11 21:26

什么叫做刚开始还可以，到后来就慢了，请说的具体点

账号		自动登录	找回密码
密码			我要加入

请高手帮忙！！关于BP实现

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（akyi）请高手帮忙！！关于BP实现