| 查看: 1044 | 回复: 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小海» 猜你喜欢
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有80人回复
-
物理学I论文润色/翻译怎么收费?
已经有81人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
» 抢金币啦！回帖就可以得到:
-
春天就该逛吃逛吃
+1/75
-
天津城建大学招收物理、化学或者生物调剂生一名，前往天津大学联合培养
+1/46
-
武汉高校接受研究生调剂，材料、化学、纺织方向
+2/44
-
广东唯一石化院校——资源与环境专硕招生
+1/43
-
球差相关知识
+1/41
-
哈尔滨工业大学（深圳）-何自开教授团队诚招化学工程与技术专业博士研究生
+1/31
-
招收荧光探针/有机合成/生物学 博士后--薪金: 30万-50万
+1/28
-
赣南师范大学智能制造与未来能源学院2026年硕士调剂信息
+1/16
-
天津医科大学基础医学院张恒课题组博士后招聘
+1/9
-
中国科学院上海硅酸盐研究所招聘科研助理（纳米材料与生物应用方向）
+1/8
-
南京林业大学化学工程学院柏惺峰课题组诚聘科研助理、水杉英才E类/F类青年教师或博后
+1/7
-
宁波东方理工大学余鹏课题组2026年上海交大联培博士生、博士后、科研助理招聘
+1/6
-
一志愿中国海洋大学095137，初试成绩312，不挑学校和地区，有学上就行！
+1/6
-
【第三轮招生】澳科大诚招2026年秋季药剂学硕士研究生（3月31日下午18:00截止）
+1/5
-
上海师范大学化学与材料科学学院任新意副研究员招收调剂学生3-4名（有机化学专业）
+1/5
-
调剂 22408 340分
+1/4
-
【博士招生】天津理工大学国家杰青王铁课题组招收2026年博士研究生
+1/4
-
深圳信息职业大学招聘联合培养硕博研究生（急）
+1/3
-
中国科大仿生界面材料科学全国重点实验室张振教授/江雷院士诚聘博士后
+1/3
-
广东工业大学环境科学与工程学院招收资源与环境专业博士研究生
+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
