| 查看: 924 | 回复: 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小海» 猜你喜欢
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有109人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有22人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
» 抢金币啦！回帖就可以得到:
鲁东大学硕士研究生招生
+1/127
-
寻找承接企业/科研机构/高校：启明人才计划（上海及长三角方向）
+1/83
-
中国医学科学院放射医学研究所、中国医学科学院天津医学健康研究院彭鑫课题组招生招聘
+1/79
-
以色列理工-生物质塑料等催化转化及流体力学方向---全奖博士研究生和科研助理
+2/52
-
2026申博自荐
+1/45
-
论文即将投出，求发表，求祝福，传递好运
+1/37
-
上海交通大学-化学化工学院-邱惠斌教授课题组招聘博士后
+1/32
-
2026年博士招生--北京理工大学交叉学科（航空，力学，能动，计算机等方向）
+1/27
-
中南大学冶金与环境学院陈伟老师招收环境科学与工程2026年博士生1人
+1/26
-
鄢勇课题组2026年拟招收项目聘用人员1名，方向：1. 具身智能；2. 智能感知；3. 忆阻器
+1/22
-
今年的函评专家费收到了吗
+1/16
-
招聘2026年入学博士生
+1/7
-
国家“双一流”建设高校-南京林业大学-国家级青年人才团队招聘
+1/7
-
招若干有机合成人员 (中山大学)
+1/6
-
北京理工大学（珠海）医疗机器人实验室现招收2026年博士研究生
+1/6
-
新加坡南洋理工大学招聘博士后（固态电池，固体离子导体，转化正极方向）
+1/6
-
招收26年秋季入学博士生（北科大高精尖学院 力学超材料/机器学习/增材制造相关方向）
+1/4
-
北理工柔性电子国家杰青团队招【博士后】【科研助理（读博意向）】
+1/2
-
海南大学化学院 招聘 材料与电化学方向——研究助理，博士（2026年入学）
+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