| 查看: 917 | 回复: 5 | |||
[交流]
求助各位大神，参数拟合结果不对，金币全送。谢谢
|
 求助各位大神，用下面的程序做的参数拟合，拟合结果和图形都不对，希望大神们帮帮忙，金币全送。谢谢function fit_nonl12 clear all;clc format long data=[ 2561.0000 10.0000 2560.0000 15.0000 2560.0000 25.0000 2558.0000 40.0000 2557.0000 60.0000 2556.0000 100.0000 2555.0000 160.0000 2554.0000 250.0000 2552.0000 400.0000 2551.0000 630.0000 2550.0000 1000.0000 ]; xspan=data(:,2); Texp=data(:,1); k0=[2.0 30.0 6.0 0.5]; lb=[0.1 0.1 0.1 0.0]; ub=[10.0 100.0 20.0 1.0]; % OPTIONS=optimset('MaxFunEvals',1000); [k,resnorm,residual,exitflag,output,lambda,jacobian]=... lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,xspan,Texp); ci=nlparci(k,residual,jacobian); %residual; fprintf('\n\n拟合结果;\n') fprintf('\n\t参数 E1 = %.16f',k(1)) fprintf('\n\t参数 E2 = %.16f',k(2)) fprintf('\n\t参数 E3 = %.16f',k(3)) fprintf('\n\t参数 E4 = %.16f',k(4)) y=kineticsEqs(xspan,k); R2=1-sum((Texp-y).^2)./sum((Texp-mean(y)).^2); fprintf('\n\tR^2 = %.16f',R2); figure plot(xspan,kineticsEqs(xspan,k),'b',xspan,Texp,'or'),legend('计算值','实验值','Location','Best') % function f =ObjFunc(k,xspan,Texp) f=kineticsEqs(xspan,k)-Texp; function x1 =kineticsEqs(x,k) E1=k(1);E2=k(2);E3=k(3);E4=k(4); x1=E1+(((2.*pi.*E3).^(-E4)+x.^(E4).*cos(pi.*E4./2)).*E2.*(2.*pi.*E3).^(-E4))./(((2.*pi.*E3).^(-E4)+x.^(E4).*cos(pi.*E4./2)).^2+(x.*(E4).*sin(pi.*E4./2)).^2); |
回复此楼» 本帖已获得的红花（最新10朵）
why小海» 猜你喜欢
🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队
已经有0人回复
-
第一性原理计算方向2026级博士申请 PRB*1，四级484
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有161人回复
-
求助VISSIM破解版软件
已经有0人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有4人回复
» 抢金币啦！回帖就可以得到:
限广州，征女友
+2/92
-
华中科技大学2026级申请考核制博士生1名-新型共轭MOF/COF合成及其超电应用方向
+1/79
-
招贤纳博（已结束）
+1/75
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收各类博士生
+1/43
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
+2/28
-
中山大学院士团队王来源教授课题组招聘博士后
+2/18
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
+1/11
-
南京大学能源与资源学院蔡亮课题组诚招2026年申请-考核制博士生2-3名
+1/10
-
现代材料与先进制造团队硕士生、博士生以及博士毕业生招收公告
+1/9
-
东华大学 唐正 课题组诚招2026年博士研究生-有机半导体材料与器件等
+1/9
-
某外资仪器厂家急招技术支持工程师
+1/8
-
湘潭大学2026年招收 力学 博士研究生2名（“申请-考核制”、硕博连读）
+1/7
-
海南大学生物医学工程学院光免疫诊疗团队诚招神经生物学、光学、分子生物学博士
+1/5
-
长江大学武汉校区诚招工程热物理、油气、电气等新能源博士-2025
+1/5
-
长江学者团队招聘高校教师7名(地点杭州、有事业编)+博后5名
+1/5
-
南开大学齐迹课题组诚聘博士后
+1/4
-
长江大学武汉校区诚招新能源博士-2025
+1/4
-
中国科学技术大学 精准智能化学重点实验室 武建昌课题组招聘博士，博士后
+1/3
-
上海交通大学 张峻课题组招收2026年申请考核博士生1 名
+1/1
-
兰州大学物理学院韩卫华教授课题组招收 2026年博士研究生 (物理、电子以及核能源方向)
+1/1
2楼2019-03-26 13:50:09
3楼2019-03-26 19:51:16
4楼2019-03-26 20:19:44
5楼2019-03-27 20:41:08
6楼2019-03-28 11:22:35