声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2668|回复: 3

[经典算法] [求助]请问这个公式怎么拟合阿!谢谢!

[复制链接]
发表于 2006-4-5 17:56 | 显示全部楼层 |阅读模式

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

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

x
请教下高手们:这个公式怎么把回归的系数求出来阿!同时保证所求参数的精度!听说可以用Marquardt来算,但是我不会。请高手们指导一下。
1.jpg
回复
分享到:

使用道具 举报

 楼主| 发表于 2006-4-6 11:07 | 显示全部楼层
我有T,u,ε和σ的数据,现在想通过这些数据把a1—a6给求出来。请大侠们指导一下阿。
发表于 2006-4-6 17:32 | 显示全部楼层
一般的神经网络教材中都会有该法(麦夸脱法)的详细介绍,具体如何应用俺就不知道了!!!
发表于 2006-4-12 07:09 | 显示全部楼层

回复:(jellycc)[求助]请问这个公式怎么拟合阿!谢谢...

一个基于Levenberg-Marquardt算法的神经网络源程序

in=[ua';ub';ia';ib'];
out=[W'/200];
VV.P=[ua1';ub1';ia1';ib1'];
VV.T=[W1'/200];

[R,Q]=size(in);
S1=9;
S2=8;
S3=1;
net=newff(minmax(in),[S1 S2 S3],{'tansig' 'tansig' 'purelin'});%Levenberg-Marquardt法
% net=newff(minmax(x),[S1 S2 S3],{'tansig' 'tansig' 'purelin'},'trainlm');%Levenberg-Marquardt法
%训练参数的设计
net=init(net);
net.trainFcn= 'trainbr';
net.performFcn='sse';%性能函数
net.trainParam.goal=0.000001;%训练目标
net.trainParam.show=50;%每50步显示一次结果
net.trainParam.epochs=5000;
net.trainParam.max_fail=500;
net.trainParam.goal=0.00001 ;
net.trainParam.mu=0.055 ;
net.trainParam.mu_dec=0.1 ;
net.trainParam.mu_inc=1 ;
net.trainParam.mu_max=1e+10;
net.trainParam.mem_reduc=1 ;
net.trainParam.min_grad=1e-10;

[net,tr]=train(net,in,out);

save TrainedNN net;
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-13 14:14 , Processed in 0.159230 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表