lmsimu.m

nnToolKit 神经网络工具包是基于 MATLAB 神经网络工具箱自行开发的一组神经网络算法函数库

%此为一仿真程序
%首先读入权域值参数

function retdouble = LmSimu(ModelNo,NetPara,SimulatePara,InputFun,OutputFun,DataDir)
NNTWARN OFF

%%%% 输入参数赋值开始 %%%%%%%%%%%%%%%%%%%%%%%
% 这部分代码主要是方便用户调试用
ModelNo='1';
NetPara(1)=7;
NetPara(2)=1;
NetPara(3)=6;

SimulatePara(1)=0.495723;
SimulatePara(2)=0.286143;
SimulatePara(3)=0.849300;
SimulatePara(4)=0.839320;
SimulatePara(5)=0.932467;
SimulatePara(6)=0.428714;
SimulatePara(7)=0.100800;

DataDir='.';
InputFun = 'tansig';
OutputFun = 'purelin';
%%%% 输入参数赋值结束 %%%%%%%%%%%%%%%%%%%%%%%%
retdouble = -1;
olddir=pwd;
% 数据所在目录
cd(DataDir);

% 输入层到中间层的传递函数
if (length(InputFun)==0)
    InputFun = 'tansig';
end
% 中间层到输出层的传递函数  
if (length(OutputFun)==0)
    OutputFun = 'purelin'
end    

%读入训练好的神经网络的权域值参数
frw1=fopen(sprintf('w%s%s',ModelNo,'1.dat'),'r');
frb1=fopen(sprintf('b%s%s',ModelNo,'1.dat'),'r');
frw2=fopen(sprintf('w%s%s',ModelNo,'2.dat'),'r');
frb2=fopen(sprintf('b%s%s',ModelNo,'2.dat'),'r');

InputDim=NetPara(1);      %输入层节点数
OutputDim=NetPara(2);     %输出层节点数
MidDim=NetPara(3);        %中间层节点数

[w1,count]=fscanf(frw1,'%f',[MidDim,InputDim]);
[b1,count]=fscanf(frb1,'%f',[MidDim,1]);
[w2,count]=fscanf(frw2,'%f',[OutputDim,MidDim]);
[b2,count]=fscanf(frb2,'%f',[OutputDim,1]);

fclose(frw1);
fclose(frb1);
fclose(frw2);
fclose(frb2);

for i=1:InputDim
    p1(i)=SimulatePara(i);
end

%网络仿真
a=simuff(p1',w1,b1,InputFun,w2,b2,OutputFun);

%将仿真结果写入文件
fwresult=fopen(sprintf('result%s%s',ModelNo,'.dat'),'w');
fprintf(fwresult,'%-15.6f\n',a);
fclose(fwresult);
retdouble = a;
logsig(1);
purelin(1);
tansig(1);
cd(olddir);