| 查看: 1152 | 回复: 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小海» 猜你喜欢
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有6人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有158人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有23人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
宝宝起名
+1/482
-
华中科技大学牛冉研究员课题组诚招硕士研究生、研究助理、博士研究生以及博士后
+1/79
-
南京大学化学学院高亚飞课题组招聘金属有机化学2027级博士研究生及研究助理（可转组内
+1/75
-
南京大学高亚飞课题组招聘金属有机化学2027级博士研究生及研究助理（可转组内博士）
+1/75
-
以色列理工-生物质塑料等催化或多相反应方向2-3名---全奖博士研究生和科研助理
+2/56
-
安徽大学材料学院招聘博士后（热电材料、热功能材料、超级电容器、气敏传感方向）
+1/46
-
美国圣路易斯大学招聘生物信息学博士后
+1/43
-
《穿越持科研神辅，杀懵外星万舰》第一章 洞中醒来是异人
+1/36
-
研究所招聘科研/实验助理一名
+1/31
-
研究所招聘科研/实验助理一名
+1/31
-
深圳大学杨楚罗/黄忠衍团队招聘博士后
+1/30
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
+1/27
-
浙江师范大学膜分离技术团队申利国教授招聘博士后研究人员
+1/23
-
2氯5甲基吡啶和2氯6三氯甲基吡啶大生产工艺需求
+1/20
-
澳大利亚昆士兰大学（UQ）】全奖PhD招生：碳捕集与液流电池储能方向
+1/14
-
北航-国家级人才团队-诚聘2027博士
+1/3
-
计算机方向 文章 辅助
+1/2
-
西南大学·食品学院 招收博士后（多糖/蛋白方向）
+2/1
-
【Scoupus英文期刊】Biofunctional Materials生物功能性材料期刊长期征稿中
+1/1
-
深圳大鹏湾实验室杰青团队招收材料学科研助理
+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