| 查看: 1235 | 回复: 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小海» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有171人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
» 抢金币啦！回帖就可以得到:
-
坐标无锡，诚征女友
+3/786
-
大连理工大学大连市复杂工业场景具身智能重点实验室科研助理招聘启事
+2/174
-
在职博士是不是没戏了，想读，没招
+1/168
-
强耀生物-热销成品蛋白-现货不用等~
+1/94
-
环境控制与良率强绑定：手套箱的工程化进阶路径
+1/93
-
中山大学深圳校区刘朝兴课题组拟招收2027级推免生
+1/93
-
石河子大学能源与材料学院李航副教授课题组招收硕、博士
+1/77
-
事业编制！天津师范大学诚聘副高/讲师职位
+2/36
-
同济大学环境学院肖倩课题组招收博士研究生（资格审核制）
+1/36
-
求助邢其毅第四版基础有机化学PDF版
+1/34
-
清华大学深圳国际研究生院宋桥课题组招收2027年普博生
+1/13
-
清华大学深圳国际研究生院宋桥课题组诚聘博士后
+1/12
-
🔥27届硕士招生｜中国计量科学研究院｜生物质谱测量技术与计量
+1/10
-
澳大利亚皇家墨尔本理工大学（RMIT University）招收博士生及 CSC 联合培养博士生
+1/7
-
推荐一个比较实用的生物信息学学习网站——The Omics Hub
+1/6
-
贵州大学2027有机化学申请考核制博士招生
+1/5
-
中南大学材料学院李明教授课题组招收推免硕士生
+1/4
-
27年应届畜牧学学硕，申请生物相关或医学博士，六级484分，JAFC二改小修
+1/2
-
给不想从事本专业的学弟学妹一点建议2
+1/2
-
ICBMBA 2026，生物功能材料与生物医学应用国际会议诚邀参与！
+1/2
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