小波神经网络程序.docVIP

下载本文档

22
0
约8.46千字
约 10页
2017-12-19 发布于河南
举报
版权申诉

小波神经网络程序.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

小波神经网络程序

主要思想：用小波变换代替普通神经网络的激励函数小波的参数b（平移因子）和a（伸缩因子）以及小波前的系数W在网络中受到训练以逼近目标函数（信号）。附件是我简单地画了一下流程图%BP小波神经网络快速学习算法clearclcSamIn=-10:0.01:10;SamOut=2*sin(3*SamIn).*exp(-0.1*SamIn.*SamIn);% figure;plot(SamIn,SamOut);% pause% Bp小波神经网络参数初始化W=rand(15,1);B=10*rands(15,1);% 对伸缩因子A进行初始化temp=rand(15,1);i=1;while(i=15)? ? if temp(i,1)=0.5? ?? ???A(i,1)=2*temp(i,1);? ? else? ?? ???A(i,1)=8*temp(i,1)-3;? ? end? ? i=i+1;end% BP算法训练小波神经网络N=2001;g=1.01;l=0.0008;t=SamIn;study_rate=0.08;count=0;Err_NetworkOut_SamOuts=[];? ?% 全局变量NetworkOutSet=[];? ?? ?? ?? ?% 局部变量? ? for i=1:N? ?? ???NetworkOut=0;? ?? ???% 局部变量? ?? ???for j=1:15? ?? ?? ?? ?tao=(t(i)-B(j,1))/A(j,1);? ?? ?? ?? ?h=cos(1.75*tao).*exp(-tao^2/2);? ?? ?? ?? ?NetworkOut=NetworkOut+h*W(j,1);? ?? ???end? ?? ???NetworkOutSet=[NetworkOutSet NetworkOut];? ? end? ? % 计算误差函数和相对误差? ? count=count+1;? ? Err_NetworkOut_SamOut=(1/2)*sumsqr(NetworkOutSet-SamOut);? ?? ?% 误差函数? ? former=Err_NetworkOut_SamOut;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? % 记录前一次的误差? ? Err_NetworkOut_SamOuts=[Err_NetworkOut_SamOuts Err_NetworkOut_SamOut]; % 误差函数集合? ? Err_SamOut=(1/2)*sumsqr(SamOut);? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?% 函数平方和? ? ErrRate=Err_NetworkOut_SamOut/Err_SamOut;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?% 相对误差? ? % 计算增加量? ? E_diff_W=0; % 局部变量? ? E_diff_B=0; % 局部变量? ? E_diff_A=0; % 局部变量? ? E_diff_W1=[]; % 局部变量? ? E_diff_B1=[]; % 局部变量? ? E_diff_A1=[]; % 局部变量? ? for i=1:15? ?? ???for j=1:N? ?? ?? ?? ?tao=(t(j)-B(i,1))/A(i,1);? ?? ?? ?? ?h=cos(1.75*tao).*exp(-tao^2/2);? ?? ?? ?? ?f=1.75*sin(1.75*tao)+tao*cos(1.75*tao);? ?? ?? ?? ?s1=SamOut(j)-NetworkOutSet(j);? ?? ?? ?? ?s2=s1*h;? ?? ?? ?? ?s3=s1*(W(i,1)/A(i,1))*exp(-tao^2/2).*f;? ?? ?? ?? ?E_diff_W=E_diff_W+s2;? ?? ?? ?? ?E_diff_B=E_diff_B+s3;? ?? ?? ?? ?E_diff_A=E_diff_A+s3*tao;? ?? ???end? ?? ???E_diff_W=-E_diff_W;? ?? ???E_diff_B=-E_diff_B;? ?? ???E_diff_A=-E_diff_A;? ?? ???E_diff_W1=[E_diff_W1,E_diff_W];? ?? ???E_diff_B1=[E_diff_B1,E_diff_B