24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ethan1129

铁杆木虫 (著名写手)

应助: 1 (幼儿园)
金币: 8632
散金: 1453
红花: 3
帖子: 1530
在线: 244.7小时
虫号: 2344158
注册: 2013-03-13
性别: GG
专业: 人工晶体

[求助] Auto2Fit拟合问题

想根据试验数据拟合经验公式中的几个参数，目前一共是3个经验参数需要拟合。但是经验公式中有一项是个积分形式，想问一下，Auto2Fit软件还可以拟合吗？希望知道的给位，能够讲解一下。小弟接触Auto2Fit不到2周。谢谢大家了

回复此楼

» 猜你喜欢

基元I理论下三大核心空间现象精准推导与细节解析已经有0人回复
基于基元 I 统一理论的反重力理论推导已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有114人回复
基于基元I统一理论的量子力学本源推导已经有1人回复
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费已经有20人回复
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）已经有1人回复
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）已经有0人回复

1楼 2016-03-13 20:38:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lu_yu_lan

新虫 (初入文坛)

应助: 3 (幼儿园)
金币: 101
散金: 5
红花: 1
帖子: 43
在线: 26.4小时
虫号: 4935903
注册: 2016-08-22

引用回帖:

6楼: Originally posted by ethan1129 at 2016-03-13 21:21:52
刚才的代码有点小问题，请帮忙试一下这个。
Variable x,y;
Parameter a,b,c,d;
Constant h=6.626196*10^(-34);
Constant f=3*10^(10);
Constant k=1.3806488*10^(-23);
Constant m=h*f/k;
Function y=a*exp ...

用OpenLu求解：

CODE:

!!!using["IMSL","luopt","math"]; //使用命名空间

Constant(::h,f,k,m)= h=6.626196*10^(-34), f=3*10^(10), k=1.3806488*10^(-23), m=h*f/k;

jf(t)= 3.27*10^(-6)+3.25*10^(-9)*t;

f(a,b,c,d:i,s,x: h,f,k,m,xArray,max)=

{

   i=0,s=0,while{i<max,

     x=xArray(i,0),

     s=s+[a*exp(-3*d*QDAGS(@jf,203.0,x))+b*(1+2/(exp(a*m/(2*x))-1))+c*(1+3/(exp(a*m/(3*x))-1)+3/(exp(a*m/(3*x))-1)^2) - xArray(i,1)]^2,

     i++

   },

   s

};

main(::xArray,max)=

{

       xArray=matrix{   //存放实验数据

"

203        768.463

 253        767.922

 303        767.39

 353        766.649

 403        765.744

 428        765.285

 453        764.664

 478        764.044

 503        763.397

 528        762.797

 553        762.171

 603        761.024

 653        759.803

 "

       },

    len[xArray,0,&max], //用len函数取矩阵的行数

    ClearImslErr(),     //清空IMSL错误输出

    ERSET(0,0,0), //关闭IMSL所有警告

    Opt[@f],      //Opt函数全局优化

    ERSET(0,2,2), ERSET(0,1,0)   //恢复IMSL警告

};

结果：
786.5100171936438 -14.84244262943891 -0.808280883032634 -3.486072801491857 8.191536241381738e-002

莫非这是楼主想要的结果，但结果并非最优解。用OpenLu在负数范围内拟合可得7#1stopt的最优解。