| 查看: 1122 | 回复: 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小海» 猜你喜欢
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有118人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
» 抢金币啦！回帖就可以得到:
-
山东征女友，坐标济南
+1/169
-
温州医科大学沈建良课题组招聘2026级博士生
+1/96
-
[急] 安徽大学热电材料方向2026年（9月份入学）招收博士生1名
+1/73
-
深圳大学申请-考核制博士生招生（第二批）：介电储能及先进封装聚合物电介质方向
+1/71
-
欢迎报考中科院苏州纳米所的氢能、燃料电池方向的博士生！（长期有效）
+1/46
-
【博士后招聘｜新加坡南洋理工大学 NTU】
+1/40
-
2026年江西师范大学药学院陈芬儿院士课题组招收智能药学博士生
+1/37
-
上海交通大学燃料电池研究所招聘博士后
+1/32
-
广州医科大学附属二院联合药学院招收科研/临床型博后（35-60）绩效优渥，有留院机会
+1/31
-
西安交通大学李平课题组诚招2026年材料科学与工程专业推免研究生2名
+1/29
-
坐标南京
+1/29
-
快报名！北京信息科技大学仪器科学与光电工程学院招收博士研究生（2026）
+1/27
-
中国科大-合肥国家实验室冷原子量子网络团队招聘启事
+2/26
-
南昌大学物理与材料学院舒龙龙教授团队招聘有机合成科研助理
+1/23
-
【北京/江苏】中国科学院下属国企 | 人工智能与平台软件工程师招聘
+1/11
-
北京大学李鲜婵研究员/国家纳米科学中心朱墨桃研究员招聘联合培养博士后
+1/9
-
中科院过程工程研究所 诚招科研助理及博士后
+1/8
-
注册小木虫周年大礼包无法领取
+1/6
-
台州学院智能制造（航空工程）学院招聘飞行控制、结构强度方向博士2-3人，教学科研岗
+1/4
-
海南大学化学院—功能分子器件团队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